Changeset 286 in svn

Timestamp:

Mar 3, 2009, 12:44:32 AM (16 years ago)

Author:

Xavier Rouby

Message:

charge added in Taujets

Location:

trunk

Files:

• interface/JetsUtil.h (modified) (1 diff)
• src/JetsUtil.cc (modified) (2 diffs)

trunk/interface/JetsUtil.h

@@ -81 +81 @@
   vector<fastjet::PseudoJet> RunJets(const vector<fastjet::PseudoJet>&  input_particles);
   void RunJetBtagging(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchJet,const vector<fastjet::PseudoJet> & sorted_jets,const TSimpleArray<TRootGenParticle> & NFCentralQ);
-  void RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TLorentzVector> & TrackCentral);
+  //void RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TLorentzVector> & TrackCentral);
+  void RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TRootTracks> & TrackCentral);
 
  private:

trunk/src/JetsUtil.cc

@@ -157 +157 @@
 }
 
-void JetsUtil::RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TLorentzVector> & TrackCentral)
+//void JetsUtil::RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TLorentzVector> & TrackCentral)
+void JetsUtil::RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TRootTracks> & TrackCentral)
 {
   TRootTauJet *elementTauJet;
   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());
@@ -168 +170 @@
       if(
          ( Energie_tau_central/JET.E() > DET->TAU_energy_frac ) &&
-         ( DET->NumTracks(TrackCentral,DET->TAU_track_pt,JET.Eta(),JET.Phi()) == 1 ) &&
+         //( 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)
          ) {
         elementTauJet = (TRootTauJet*) branchTauJet->NewEntry();
         elementTauJet->Set(JET);
+        elementTauJet->Charge = charge;
       } // if tau jet
     } // if JET.eta < tracker - tau_cone : Tau jet identification