Changes in / [771843a:b5e2583] in git
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/PdgCodeFilter.cc
r771843a rb5e2583 76 76 fInvert = GetBool("Invert", false); 77 77 78 // no pileup79 fRequireNotPileup = GetBool("RequireNotPileup", false);80 81 78 fRequireStatus = GetBool("RequireStatus", false); 82 79 fStatus = GetInt("Status", 1); … … 130 127 if(fRequireStatus && (candidate->Status != fStatus)) continue; 131 128 if(fRequireCharge && (candidate->Charge != fCharge)) continue; 132 if(fRequireNotPileup && (candidate->IsPU >0 )) continue;133 129 134 130 pass = kTRUE; -
modules/PdgCodeFilter.h
r771843a rb5e2583 55 55 Bool_t fRequireCharge; //! 56 56 Int_t fCharge; //! 57 Bool_t fRequireNotPileup; //! 57 58 58 59 59 std::vector<Int_t> fPdgCodes; -
modules/StatusPidFilter.cc
r771843a rb5e2583 153 153 { 154 154 // PT threshold 155 155 156 fPTMin = GetDouble("PTMin", 0.5); 156 157 // keep or remove pileup particles158 fRequireNotPileup = GetBool("RequireNotPileup", false);159 157 160 158 // import input array … … 229 227 if(!pass || (candidate->Momentum.Pt() < fPTMin && !(is_b_hadron || is_b_quark || is_tau_daughter || is_W_daughter)) ) continue; 230 228 231 // not pileup particles232 if(fRequireNotPileup && (candidate->IsPU >0)) continue;233 234 229 fOutputArray->Add(candidate); 235 230 } -
modules/StatusPidFilter.h
r771843a rb5e2583 51 51 Double_t fPTMin; //! 52 52 53 Bool_t fRequireNotPileup; //!54 55 53 TIterator *fItInputArray; //! 56 54 -
readers/DelphesProMC.cpp
r771843a rb5e2583 125 125 candidate->Status = status; 126 126 127 candidate->IsPU=0;128 if (mutableParticles->barcode(i)>0) candidate->IsPU=1; // pileup particle129 130 127 candidate->M1 = mutableParticles->mother1(i); 131 128 candidate->M2 = mutableParticles->mother2(i); … … 234 231 cout << "** Reading " << argv[i] << endl; 235 232 236 // use 64 bit237 //inputFile = new ProMCBook(argv[i], "r", true);238 239 //use 32 bit zip (faster but limitted to 64k events)240 233 inputFile = new ProMCBook(argv[i], "r"); 241 234
Note:
See TracChangeset
for help on using the changeset viewer.