Changeset 29b722a in git for classes

Timestamp:

Dec 9, 2021, 8:42:12 PM (3 years ago)

Author:

christinaw97 <christina.wang@…>

Children:

a09b75f

Parents:

d612dec

Message:

Added DelphesCscClusterFormula

Location:

classes

Files:

• DelphesCscClusterFormula.cc (added)
• DelphesCscClusterFormula.h (added)
• DelphesHepMC2Reader.cc (modified) (4 diffs)
• DelphesHepMC3Reader.cc (modified) (3 diffs)

classes/DelphesHepMC2Reader.cc

@@ -438 +438 @@
 {
   Candidate *candidate;
+  Candidate *candidateDaughter;
   map<int, pair<int, int> >::iterator itMotherMap;
   map<int, pair<int, int> >::iterator itDaughterMap;
@@ -446 +447 @@
     candidate = static_cast<Candidate *>(allParticleOutputArray->At(i));
 
+
     if(candidate->M1 > 0)
     {
@@ -477 +479 @@
         candidate->D1 = -1;
         candidate->D2 = -1;
+        const TLorentzVector &decayPosition = candidate->Position;
+        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
       }
       else
@@ -482 +486 @@
         candidate->D1 = itDaughterMap->second.first;
         candidate->D2 = itDaughterMap->second.second;
-      }
-    }
-  }
-}
-
-//---------------------------------------------------------------------------
+        candidateDaughter = static_cast<Candidate *>(allParticleOutputArray->At(candidate->D1));
+        const TLorentzVector &decayPosition = candidateDaughter->Position;
+        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
+
+      }
+    }
+  }
+}
+
+//---------------------------------------------------------------------------

classes/DelphesHepMC3Reader.cc

@@ -458 +458 @@
   TObjArray *array;
   Candidate *candidate;
+  Candidate *candidateDaughter;
   TParticlePDG *pdgParticle;
   int pdgCode;
@@ -573 +574 @@
         candidate->D1 = -1;
         candidate->D2 = -1;
+        const TLorentzVector &decayPosition = candidate->Position;
+        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
       }
       else
@@ -578 +581 @@
         candidate->D1 = itDaughterMap->second.first;
         candidate->D2 = itDaughterMap->second.second;
-      }
-    }
-  }
-}
-
-//---------------------------------------------------------------------------
+        candidateDaughter = static_cast<Candidate *>(allParticleOutputArray->At(candidate->D1));
+        const TLorentzVector &decayPosition = candidateDaughter->Position;
+        candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
+      }
+    }
+  }
+}
+
+//---------------------------------------------------------------------------