Changeset 48f7a77 in git

Timestamp:

Apr 28, 2016, 11:41:56 AM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

632da30

Parents:

72e0320

Message:

added compute beamspot function in propagator

Location:

modules

Files:

• ParticlePropagator.cc (modified) (2 diffs)
• ParticlePropagator.h (modified) (2 diffs)

modules/ParticlePropagator.cc

@@ -66 +66 @@
 {
 }
+
+
+//------------------------------------------------------------------------------
+
+TLorentzVector ParticlePropagator::BeamSpotPosition(const TObjArray *array)
+{
+  Candidate *candidate;
+  Bool_t passed = false;
+  TLorentzVector bs;
+ 
+  fItInputArray->Reset();
+  while((candidate = static_cast<Candidate*>(fItInputArray->Next())) && !passed)
+  {
+    if(candidate->IsPU == 0) passed = true;
+    bs = candidate->Position;
+  }
+  
+  return bs;
+  
+}
+
+
 
 //------------------------------------------------------------------------------
@@ -124 +146 @@
   const Double_t c_light = 2.99792458E8;
 
+  const TLorentzVector &bs = BeamSpotPosition(fInputArray);
+
   fItInputArray->Reset();
   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))

modules/ParticlePropagator.h

@@ -35 +35 @@
 class TClonesArray;
 class TIterator;
+class TLorentzVector;
 
 class ParticlePropagator: public DelphesModule
@@ -61 +62 @@
   TObjArray *fMuonOutputArray; //!
 
+  TLorentzVector BeamSpotPosition(const TObjArray *array);
+
   ClassDef(ParticlePropagator, 1)
 };