Changeset 64a4950 in git

Timestamp:

Apr 16, 2014, 3:57:35 PM (11 years ago)

Author:

pavel <pavel@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

a0431dc

Parents:

3c40083

Message:

add Particle reference to HectorHit

Files:

• classes/DelphesClasses.h (modified) (1 diff)
• modules/Hector.cc (modified) (2 diffs)
• modules/TreeWriter.cc (modified) (1 diff)

classes/DelphesClasses.h

@@ -404 +404 @@
   Float_t S; // distance to the interaction point [m]
 
+  TRef Particle; // reference to generated particle
+
   static CompBase *fgCompare; //!
   const CompBase *GetCompare() const { return fgCompare; }

modules/Hector.cc

@@ -112 +112 @@
     pz = candidateMomentum.Pz();
 
-    if(candidateMomentum.Eta() <= fEtaMin || fDirection*pz <= 0.0) continue;
+    if(TMath::Abs(candidateMomentum.Eta()) <= fEtaMin || TMath::Sign(pz, Double_t(fDirection)) != pz) continue;
 
     x = 1.0E3 * candidatePosition.X();
@@ -129 +129 @@
 
     H_BeamParticle particle(candidate->Mass, candidate->Charge);
-    particle.set4Momentum(candidateMomentum);
+//    particle.set4Momentum(candidateMomentum);
+    particle.set4Momentum(candidateMomentum.Px(), candidateMomentum.Py(), 
+                          candidateMomentum.Pz(), candidateMomentum.E());
     particle.setPosition(x, y, tx, ty, z);
 

modules/TreeWriter.cc

@@ -658 +658 @@
     entry->Y = position.Y();
     entry->S = position.Z();
+
+    entry->Particle = candidate->GetCandidates()->At(0);
   }
 }