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