- Timestamp:
- Jul 4, 2013, 1:44:56 PM (11 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 6e0bbe3
- Parents:
- a9c6df69
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
readers/DelphesPythia8.cpp
ra9c6df69 r474eb76 5 5 #include <signal.h> 6 6 7 #include " pythia->h"7 #include "Pythia.h" 8 8 9 9 #include "TROOT.h" … … 40 40 Int_t pid, status; 41 41 Double_t px, py, pz, e; 42 Double_t x, y, z ;42 Double_t x, y, z, t; 43 43 44 44 pdg = TDatabasePDG::Instance(); … … 50 50 pid = particle.id(); 51 51 status = particle.status(); 52 m1 = particle.mother1(); m2 = particle.mother2();53 d1 = particle.daughter1(); d2 = particle.daughter2();54 52 px = particle.px(); py = particle.py(); pz = particle.pz(); e = particle.e(); 55 53 x = particle.xProd(); y = particle.yProd(); z = particle.zProd(); t = particle.tProd(); … … 62 60 candidate->Status = status; 63 61 64 candidate->M1 = m1- 1;65 candidate->M2 = m2- 1;66 67 candidate->D1 = d1- 1;68 candidate->D2 = d2- 1;62 candidate->M1 = particle.mother1() - 1; 63 candidate->M2 = particle.mother2() - 1; 64 65 candidate->D1 = particle.daughter1() - 1; 66 candidate->D2 = particle.daughter2() - 1; 69 67 70 68 pdgParticle = pdg->GetParticle(pid); … … 74 72 candidate->Momentum.SetPxPyPzE(px, py, pz, e); 75 73 76 candidate->Position.SetXYZT(x, y, z, 0.0);74 candidate->Position.SetXYZT(x, y, z, t); 77 75 78 76 allParticleOutputArray->Add(candidate); … … 106 104 char appName[] = "DelphesPythia8"; 107 105 stringstream message; 108 FILE *inputFile = 0;109 106 TFile *outputFile = 0; 110 107 TStopwatch readStopWatch, procStopWatch; … … 115 112 DelphesFactory *factory = 0; 116 113 TObjArray *stableParticleOutputArray = 0, *allParticleOutputArray = 0, *partonOutputArray = 0; 117 Int_t i;118 114 Long64_t eventCounter, errorCounter; 119 115 Long64_t numberOfEvents, timesAllowErrors;
Note:
See TracChangeset
for help on using the changeset viewer.