Changeset 315 in svn

Timestamp:

Mar 10, 2009, 11:16:44 PM (16 years ago)

Author:

Xavier Rouby

Message:

from TRootGenParticle to GenParticle in B-tag

File:

• trunk/src/JetsUtil.cc (modified) (5 diffs)

trunk/src/JetsUtil.cc

@@ -190 +190 @@
 
 
-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)
+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)
 {
   TRootJet *elementJet;
   TLorentzVector JET;
   for (unsigned int i = 0; i < sorted_jets.size(); i++) {
-
     JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
     if(JET.Pt() > DET->PTCUT_jet)
@@ -210 +209 @@
       }
   } // for itJet : loop on all jets
-  
 }
 
@@ -218 +216 @@
   TLorentzVector JET;
   float charge=0;
+
   for (unsigned int i = 0; i < sorted_jets.size(); i++) {
     JET.SetPxPyPzE(sorted_jets[i].px(),sorted_jets[i].py(),sorted_jets[i].pz(),sorted_jets[i].E());
@@ -225 +224 @@
       if(
          ( Energie_tau_central/JET.E() > DET->TAU_energy_frac ) &&
-         //( DET->NumTracks(charge,TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi(),DET->TAU_track_scone) == 1 ) &&
          ( DET->NumTracks(charge,TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi()) == 1 ) &&
          ( JET.Pt() > DET->PTCUT_taujet)
@@ -236 +234 @@
       } // if tau jet
     } // if JET.eta < tracker - tau_cone : Tau jet identification
-    
   } // for itJet : loop on all jets
-  
-  
-}
-
+}