Changeset 50d7259 in git
- Timestamp:
- Feb 27, 2015, 11:16:07 AM (10 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 76ab6be
- Parents:
- b55dc12
- git-author:
- Michele Selvaggi <michele.selvaggi@…> (02/27/15 11:11:12)
- git-committer:
- Michele Selvaggi <michele.selvaggi@…> (02/27/15 11:16:07)
- Location:
- modules
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/TauTagging.cc
rb55dc12 r50d7259 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 -
modules/TauTagging.h
rb55dc12 r50d7259 31 31 32 32 #include "classes/DelphesModule.h" 33 #include "ExRootAnalysis/ExRootResult.h" 34 #include "ExRootAnalysis/ExRootFilter.h" 35 #include "ExRootAnalysis/ExRootClassifier.h" 33 36 34 37 #include <map> … … 76 79 }; 77 80 81 82 //------------------------------------------------------------------------------ 83 84 class TauTaggingPartonClassifier : public ExRootClassifier 85 { 86 public: 87 88 TauTaggingPartonClassifier(const TObjArray *array); 89 90 Int_t GetCategory(TObject *object); 91 92 Double_t fEtaMax, fPTMin; 93 94 const TObjArray *fParticleInputArray; 95 }; 96 97 78 98 #endif
Note:
See TracChangeset
for help on using the changeset viewer.