Changeset 9be57d9 in git for modules/TauTagging.cc
- Timestamp:
- Feb 27, 2015, 1:41:42 PM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- fa33983
- Parents:
- d5cae2b (diff), 532016b (diff), 76ab6be (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the(diff)
links above to see all the changes relative to each parent. - File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/TauTagging.cc
rd5cae2b r9be57d9 33 33 #include "classes/DelphesFactory.h" 34 34 #include "classes/DelphesFormula.h" 35 36 #include "ExRootAnalysis/ExRootResult.h"37 #include "ExRootAnalysis/ExRootFilter.h"38 #include "ExRootAnalysis/ExRootClassifier.h"39 35 40 36 #include "TMath.h" … … 53 49 using namespace std; 54 50 55 //------------------------------------------------------------------------------56 57 class TauTaggingPartonClassifier : public ExRootClassifier58 {59 public:60 61 TauTaggingPartonClassifier(const TObjArray *array);62 63 Int_t GetCategory(TObject *object);64 65 Double_t fEtaMax, fPTMin;66 67 const TObjArray *fParticleInputArray;68 };69 51 70 52 //------------------------------------------------------------------------------ … … 79 61 { 80 62 Candidate *tau = static_cast<Candidate *>(object); 81 Candidate *daughter = 0; 63 Candidate *daughter1 = 0; 64 Candidate *daughter2 = 0; 65 82 66 const TLorentzVector &momentum = tau->Momentum; 83 Int_t pdgCode, i ;67 Int_t pdgCode, i, j; 84 68 85 69 pdgCode = TMath::Abs(tau->PID); … … 98 82 for(i = tau->D1; i <= tau->D2; ++i) 99 83 { 100 daughter = static_cast<Candidate *>(fParticleInputArray->At(i)); 101 pdgCode = TMath::Abs(daughter->PID); 102 if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15 || pdgCode == 24) return -1; 84 daughter1 = static_cast<Candidate *>(fParticleInputArray->At(i)); 85 pdgCode = TMath::Abs(daughter1->PID); 86 if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15) return -1; 87 else if(pdgCode == 24) 88 { 89 if(daughter1->D1 < 0) return -1; 90 for(j = daughter1->D1; j <= daughter1->D2; ++j) 91 { 92 daughter2 = static_cast<Candidate*>(fParticleInputArray->At(j)); 93 pdgCode = TMath::Abs(daughter2->PID); 94 if(pdgCode == 11 || pdgCode == 13) return -1; 95 } 96 97 } 103 98 } 104 99
Note:
See TracChangeset
for help on using the changeset viewer.