Changeset 6fbd089 in git for modules/PileUpMerger.cc

Timestamp:

Jun 2, 2016, 3:36:31 PM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

9fb48ad

Parents:

ab3bdd9

Message:

fill gen vertex

File:

• modules/PileUpMerger.cc (modified) (5 diffs)

modules/PileUpMerger.cc

@@ -116 +116 @@
   TParticlePDG *pdgParticle;
   Int_t pid, nch, nvtx = -1;
-  Float_t x, y, z, t, vx, vy;
+  Float_t x, y, z, t, vx, vy, vz, vt;
   Float_t px, py, pz, e, pt;
   Double_t dz, dphi, dt, sumpt2;
@@ -136 +136 @@
   vx = 0.0;
   vy = 0.0;
+  vz = 0.0;
+  vt = 0.0;
+  
   numberOfParticles = fInputArray->GetEntriesFast();
   nch = 0;
@@ -147 +150 @@
     t = candidate->Position.T();
     pt = candidate->Momentum.Pt();
+    
+    vz += z+dz;
+    vt += t+dt;
+    
     candidate->Position.SetZ(z + dz);
     candidate->Position.SetT(t + dt);
@@ -162 +169 @@
     vx /= numberOfParticles;
     vy /= numberOfParticles;
+    vz /= numberOfParticles;
+    vt /= numberOfParticles;
+  
   }
 
@@ -168 +178 @@
 
   vertex = factory->NewCandidate();
-  vertex->Position.SetXYZT(vx, vy, dz, dt);
+  vertex->Position.SetXYZT(vx, vy, vz, vt);
   vertex->ClusterIndex = nvtx;
   vertex->ClusterNDF = nch;
   vertex->SumPT2 = sumpt2;
-  vertex->GenSumPT2 = sumpt2;  
+  vertex->GenSumPT2 = sumpt2; 
   
   fVertexOutputArray->Add(vertex);