Changeset 341014c in git for modules/BTagging.cc
- Timestamp:
- Feb 12, 2019, 9:29:17 PM (6 years ago)
- Branches:
- ImprovedOutputFile, Timing, llp, master
- Children:
- 6455202
- Parents:
- 45e58be
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/BTagging.cc
r45e58be r341014c 17 17 */ 18 18 19 20 19 /** \class BTagging 21 20 * … … 34 33 #include "classes/DelphesFormula.h" 35 34 35 #include "TDatabasePDG.h" 36 #include "TFormula.h" 37 #include "TLorentzVector.h" 36 38 #include "TMath.h" 39 #include "TObjArray.h" 40 #include "TRandom3.h" 37 41 #include "TString.h" 38 #include "TFormula.h"39 #include "TRandom3.h"40 #include "TObjArray.h"41 #include "TDatabasePDG.h"42 #include "TLorentzVector.h"43 42 44 43 #include <algorithm> 45 #include <stdexcept>46 44 #include <iostream> 47 45 #include <sstream> 46 #include <stdexcept> 48 47 49 48 using namespace std; … … 66 65 void BTagging::Init() 67 66 { 68 map< Int_t, DelphesFormula *>::iterator itEfficiencyMap;67 map<Int_t, DelphesFormula *>::iterator itEfficiencyMap; 69 68 ExRootConfParam param; 70 69 DelphesFormula *formula; … … 78 77 79 78 fEfficiencyMap.clear(); 80 for(i = 0; i < size /2; ++i)79 for(i = 0; i < size / 2; ++i) 81 80 { 82 81 formula = new DelphesFormula; 83 formula->Compile(param[i *2 + 1].GetString());82 formula->Compile(param[i * 2 + 1].GetString()); 84 83 85 fEfficiencyMap[param[i *2].GetInt()] = formula;84 fEfficiencyMap[param[i * 2].GetInt()] = formula; 86 85 } 87 86 … … 106 105 void BTagging::Finish() 107 106 { 108 map< Int_t, DelphesFormula *>::iterator itEfficiencyMap;107 map<Int_t, DelphesFormula *>::iterator itEfficiencyMap; 109 108 DelphesFormula *formula; 110 109 … … 124 123 Candidate *jet; 125 124 Double_t pt, eta, phi, e; 126 map< Int_t, DelphesFormula *>::iterator itEfficiencyMap;125 map<Int_t, DelphesFormula *>::iterator itEfficiencyMap; 127 126 DelphesFormula *formula; 128 127 129 128 // loop over all input jets 130 129 fItJetInputArray->Reset(); 131 while((jet = static_cast<Candidate *>(fItJetInputArray->Next())))130 while((jet = static_cast<Candidate *>(fItJetInputArray->Next()))) 132 131 { 133 132 const TLorentzVector &jetMomentum = jet->Momentum;
Note:
See TracChangeset
for help on using the changeset viewer.