Changes in modules/IdentificationMap.cc [341014c:95aa610] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/IdentificationMap.cc
r341014c r95aa610 17 17 */ 18 18 19 19 20 /** \class IdentificationMap 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 47 #include <algorithm> 48 #include <stdexcept> 47 49 #include <iostream> 48 50 #include <sstream> 49 #include <stdexcept>50 51 51 52 using namespace std; … … 78 79 79 80 fEfficiencyMap.clear(); 80 for(i = 0; i < size /3; ++i)81 for(i = 0; i < size/3; ++i) 81 82 { 82 83 formula = new DelphesFormula; 83 formula->Compile(param[i *3 + 2].GetString());84 pdg = param[i *3].GetInt();85 fEfficiencyMap.insert(make_pair(pdg, make_pair(param[i *3 + 1].GetInt(), formula)));84 formula->Compile(param[i*3 + 2].GetString()); 85 pdg = param[i*3].GetInt(); 86 fEfficiencyMap.insert(make_pair(pdg, make_pair(param[i*3 + 1].GetInt(), formula))); 86 87 } 87 88 … … 128 129 Double_t pt, eta, phi, e; 129 130 TMisIDMap::iterator itEfficiencyMap; 130 pair <TMisIDMap::iterator, TMisIDMap::iterator> range;131 pair <TMisIDMap::iterator, TMisIDMap::iterator> range; 131 132 DelphesFormula *formula; 132 133 Int_t pdgCodeIn, pdgCodeOut, charge; … … 135 136 136 137 fItInputArray->Reset(); 137 while((candidate = static_cast<Candidate 138 while((candidate = static_cast<Candidate*>(fItInputArray->Next()))) 138 139 { 139 140 const TLorentzVector &candidatePosition = candidate->Position; … … 143 144 pt = candidateMomentum.Pt(); 144 145 e = candidateMomentum.E(); 145 146 146 147 pdgCodeIn = candidate->PID; 147 148 charge = candidate->Charge; … … 170 171 { 171 172 // change PID of particle 172 if(pdgCodeOut != 0) candidate->PID = charge *pdgCodeOut;173 if(pdgCodeOut != 0) candidate->PID = charge*pdgCodeOut; 173 174 fOutputArray->Add(candidate); 174 175 break;
Note:
See TracChangeset
for help on using the changeset viewer.