Changes in modules/JetFlavorAssociation.cc [341014c:254db49] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/JetFlavorAssociation.cc
r341014c r254db49 17 17 */ 18 18 19 19 20 /** \class JetFlavorAssociation 20 21 * … … 31 32 #include "classes/DelphesFormula.h" 32 33 34 #include "ExRootAnalysis/ExRootResult.h" 35 #include "ExRootAnalysis/ExRootFilter.h" 33 36 #include "ExRootAnalysis/ExRootClassifier.h" 34 #include "ExRootAnalysis/ExRootFilter.h" 35 #include "ExRootAnalysis/ExRootResult.h" 36 37 38 #include "TMath.h" 39 #include "TString.h" 40 #include "TFormula.h" 41 #include "TRandom3.h" 42 #include "TObjArray.h" 37 43 #include "TDatabasePDG.h" 38 #include "TFormula.h"39 44 #include "TLorentzVector.h" 40 #include "TMath.h"41 #include "TObjArray.h"42 #include "TRandom3.h"43 #include "TString.h"44 45 45 46 #include <algorithm> 47 #include <stdexcept> 46 48 #include <iostream> 47 49 #include <sstream> 48 #include <stdexcept>49 50 50 51 using namespace std; … … 52 53 //------------------------------------------------------------------------------ 53 54 54 class PartonClassifier 55 class PartonClassifier: public ExRootClassifier 55 56 { 56 57 public: 58 57 59 PartonClassifier() {} 58 60 Int_t GetCategory(TObject *object); … … 65 67 Int_t PartonClassifier::GetCategory(TObject *object) 66 68 { 67 69 // select parton in the parton list 68 70 69 71 Candidate *parton = static_cast<Candidate *>(object); … … 85 87 //------------------------------------------------------------------------------ 86 88 87 class ParticleLHEFClassifier 89 class ParticleLHEFClassifier: public ExRootClassifier 88 90 { 89 91 public: 92 90 93 ParticleLHEFClassifier() {} 91 94 Int_t GetCategory(TObject *object); … … 187 190 //------------------------------------------------------------------------------ 188 191 189 void JetFlavorAssociation::Process() 190 { 192 void JetFlavorAssociation::Process(){ 191 193 192 194 Candidate *jet; … … 199 201 if(partonArray == 0) return; 200 202 201 if(fParticleLHEFInputArray) 203 if(fParticleLHEFInputArray) 202 204 { 203 205 fParticleLHEFFilter->Reset(); … … 225 227 Candidate *tempParton = 0, *tempPartonHighestPt = 0; 226 228 int pdgCode, pdgCodeMax = -1; 227 229 228 230 TIter itPartonArray(partonArray); 229 231 TIter itPartonLHEFArray(partonLHEFArray); … … 239 241 if(pdgCodeMax < pdgCode) pdgCodeMax = pdgCode; 240 242 } 241 243 242 244 if(!fParticleLHEFInputArray) continue; 243 245 244 246 itPartonLHEFArray.Reset(); 245 247 while((partonLHEF = static_cast<Candidate *>(itPartonLHEFArray.Next()))) 246 248 { 247 if(parton->Momentum.DeltaR(partonLHEF->Momentum) < 0.001 && parton->PID == partonLHEF->PID && partonLHEF->Charge == parton->Charge) 248 { 249 break; 249 if(parton->Momentum.DeltaR(partonLHEF->Momentum) < 0.001 && 250 parton->PID == partonLHEF->PID && 251 partonLHEF->Charge == parton->Charge) 252 { 253 break; 250 254 } 251 255 … … 326 330 while((partonLHEF = static_cast<Candidate *>(itPartonLHEFArray.Next()))) 327 331 { 328 if(parton->Momentum.DeltaR(partonLHEF->Momentum) < 0.01 && parton->PID == partonLHEF->PID && partonLHEF->Charge == parton->Charge) 332 if(parton->Momentum.DeltaR(partonLHEF->Momentum) < 0.01 && 333 parton->PID == partonLHEF->PID && 334 partonLHEF->Charge == parton->Charge) 329 335 { 330 336 isGoodCandidate = false;
Note:
See TracChangeset
for help on using the changeset viewer.