Fork me on GitHub

Changeset 315 in svn for trunk


Ignore:
Timestamp:
Mar 10, 2009, 11:16:44 PM (16 years ago)
Author:
Xavier Rouby
Message:

from TRootGenParticle to GenParticle in B-tag

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/JetsUtil.cc

    r310 r315  
    190190
    191191
    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)
     192void JetsUtil::RunJetBtagging(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchJet,const vector<fastjet::PseudoJet> & sorted_jets,const TSimpleArray<GenParticle>& NFCentralQ, const vector<int> &NTrackJet, const vector<float> &EHADEEM)
    193193{
    194194  TRootJet *elementJet;
    195195  TLorentzVector JET;
    196196  for (unsigned int i = 0; i < sorted_jets.size(); i++) {
    197 
    198197    JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
    199198    if(JET.Pt() > DET->PTCUT_jet)
     
    210209      }
    211210  } // for itJet : loop on all jets
    212  
    213211}
    214212
     
    218216  TLorentzVector JET;
    219217  float charge=0;
     218
    220219  for (unsigned int i = 0; i < sorted_jets.size(); i++) {
    221220    JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
     
    225224      if(
    226225         ( 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 ) &&
    228226         ( DET->NumTracks(charge,TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi()) == 1 ) &&
    229227         ( JET.Pt() > DET->PTCUT_taujet)
     
    236234      } // if tau jet
    237235    } // if JET.eta < tracker - tau_cone : Tau jet identification
    238    
    239236  } // for itJet : loop on all jets
    240  
    241  
    242 }
    243 
     237}
Note: See TracChangeset for help on using the changeset viewer.