- Timestamp:
- Mar 10, 2009, 11:16:44 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/JetsUtil.cc
r310 r315 190 190 191 191 192 void JetsUtil::RunJetBtagging(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchJet,const vector<fastjet::PseudoJet> & sorted_jets,const TSimpleArray< TRootGenParticle>& NFCentralQ, const vector<int> &NTrackJet, const vector<float> &EHADEEM)192 void JetsUtil::RunJetBtagging(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchJet,const vector<fastjet::PseudoJet> & sorted_jets,const TSimpleArray<GenParticle>& NFCentralQ, const vector<int> &NTrackJet, const vector<float> &EHADEEM) 193 193 { 194 194 TRootJet *elementJet; 195 195 TLorentzVector JET; 196 196 for (unsigned int i = 0; i < sorted_jets.size(); i++) { 197 198 197 JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E()); 199 198 if(JET.Pt() > DET->PTCUT_jet) … … 210 209 } 211 210 } // for itJet : loop on all jets 212 213 211 } 214 212 … … 218 216 TLorentzVector JET; 219 217 float charge=0; 218 220 219 for (unsigned int i = 0; i < sorted_jets.size(); i++) { 221 220 JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E()); … … 225 224 if( 226 225 ( Energie_tau_central/JET.E() > DET->TAU_energy_frac ) && 227 //( DET->NumTracks(charge,TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi(),DET->TAU_track_scone) == 1 ) &&228 226 ( DET->NumTracks(charge,TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi()) == 1 ) && 229 227 ( JET.Pt() > DET->PTCUT_taujet) … … 236 234 } // if tau jet 237 235 } // if JET.eta < tracker - tau_cone : Tau jet identification 238 239 236 } // for itJet : loop on all jets 240 241 242 } 243 237 }
Note:
See TracChangeset
for help on using the changeset viewer.