Changeset cc8716b in git for classes/DelphesHepMC2Reader.cc

Timestamp:

Mar 7, 2022, 8:51:45 PM (3 years ago)

Author:

GitHub <noreply@…>

Branches:

master

Children:

7dbc149

Parents:

4e8e72b

git-author:

Christinaw97 <cwang5@…> (03/07/22 20:51:45)

git-committer:

GitHub <noreply@…> (03/07/22 20:51:45)

Message:

Update to handle CMS endcap muon detector showers for long-lived particles (#103)

• add CscCluster modules and classes
• add DelphesCscClusterFormula
• add example/ExampleCscCluster.py

Co-authored-by: christinaw97 <christina.wang@…>

File:

• classes/DelphesHepMC2Reader.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
       {
         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
+        
+
+      }
+    }
+  }
+}
+
+//---------------------------------------------------------------------------