Changes in modules/Merger.cc [341014c:cab38f6] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Merger.cc
r341014c rcab38f6 17 17 */ 18 18 19 19 20 /** \class Merger 20 21 * … … 32 33 #include "classes/DelphesFormula.h" 33 34 35 #include "ExRootAnalysis/ExRootResult.h" 36 #include "ExRootAnalysis/ExRootFilter.h" 34 37 #include "ExRootAnalysis/ExRootClassifier.h" 35 #include "ExRootAnalysis/ExRootFilter.h"36 #include "ExRootAnalysis/ExRootResult.h"37 38 39 #include "TMath.h" 40 #include "TString.h" 41 #include "TFormula.h" 42 #include "TRandom3.h" 43 #include "TObjArray.h" 38 44 #include "TDatabasePDG.h" 39 #include "TFormula.h"40 45 #include "TLorentzVector.h" 41 #include "TMath.h"42 #include "TObjArray.h"43 #include "TRandom3.h"44 #include "TString.h"45 46 46 #include <algorithm> 47 #include <algorithm> 48 #include <stdexcept> 47 49 #include <iostream> 48 50 #include <sstream> 49 #include <stdexcept>50 51 51 52 using namespace std; … … 88 89 89 90 fMomentumOutputArray = ExportArray(GetString("MomentumOutputArray", "momentum")); 90 91 91 92 fEnergyOutputArray = ExportArray(GetString("EnergyOutputArray", "energy")); 92 93 } … … 96 97 void Merger::Finish() 97 98 { 98 vector< TIterator *>::iterator itInputList;99 vector< TIterator * >::iterator itInputList; 99 100 TIterator *iterator; 100 101 … … 112 113 Candidate *candidate; 113 114 TLorentzVector momentum; 114 Double_t sumPT, sumE; 115 vector< TIterator *>::iterator itInputList;115 Double_t sumPT, sumE; 116 vector< TIterator * >::iterator itInputList; 116 117 TIterator *iterator; 117 118 118 119 DelphesFactory *factory = GetFactory(); 119 120 120 121 momentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0); 121 122 sumPT = 0; … … 129 130 // loop over all candidates 130 131 iterator->Reset(); 131 while((candidate = static_cast<Candidate 132 while((candidate = static_cast<Candidate*>(iterator->Next()))) 132 133 { 133 134 const TLorentzVector &candidateMomentum = candidate->Momentum; … … 142 143 143 144 candidate = factory->NewCandidate(); 144 145 145 146 candidate->Position.SetXYZT(0.0, 0.0, 0.0, 0.0); 146 147 candidate->Momentum = momentum; 147 148 148 149 fMomentumOutputArray->Add(candidate); 149 150 150 151 candidate = factory->NewCandidate(); 151 152 152 153 candidate->Position.SetXYZT(0.0, 0.0, 0.0, 0.0); 153 154 candidate->Momentum.SetPtEtaPhiE(sumPT, 0.0, 0.0, sumE); 154 155 155 156 fEnergyOutputArray->Add(candidate); 156 157 }
Note:
See TracChangeset
for help on using the changeset viewer.