Changeset 310 in svn for trunk/src/JetsUtil.cc

Timestamp:

Mar 10, 2009, 11:52:16 AM (16 years ago)

Author:

severine ovyn

Message:

add Ehad/Eem

File:

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

trunk/src/JetsUtil.cc

@@ -120 +120 @@
 }
 
-vector<fastjet::PseudoJet> JetsUtil::RunJets(const vector<fastjet::PseudoJet>&  input_particles,  const vector<TRootTracks> & TrackCentral,  vector<int> &NTrackJet)
+vector<fastjet::PseudoJet> JetsUtil::RunJets(const vector<fastjet::PseudoJet>&  input_particles,  const vector<TRootTracks> & TrackCentral,  vector<int> &NTrackJet, vector<float> &EHADEEM,D_CaloTowerList list_of_active_towers)
 {
   inclusive_jets.clear();
@@ -135 +135 @@
       sorted_jets = sorted_by_pt(inclusive_jets);
 
-      //check number of tracks in jets 
+      //Bin tracks to make the link 
       float  iEtaTrack[TrackCentral.size()];
       float  iPhiTrack[TrackCentral.size()];
@@ -145 +145 @@
       for (unsigned int i = 0; i < sorted_jets.size(); i++) 
           {
+            //check number of tracks in jets 
             numTrackJet=0;
             vector<fastjet::PseudoJet> constituents = clust_seq.constituents(sorted_jets[i]);
@@ -155 +156 @@
                }
             NTrackJet.push_back(numTrackJet);
+            //now, get EHad over EEm
+            float EmVal=0,HadVal=0;
+            for (unsigned int j = 0; j < constituents.size(); j++)
+                {
+                  D_CaloTower calConsti(list_of_active_towers.getElement(constituents[j].eta(),constituents[j].phi()));
+                  EmVal  += calConsti.getEem();
+                  HadVal += calConsti.getEhad();
+                }
+            EHADEEM.push_back(HadVal/EmVal);
        }
     }
@@ -161 +171 @@
 }
 
-
-
-void JetsUtil::RunJetBtagging(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchJet,const vector<fastjet::PseudoJet> & sorted_jets,const TSimpleArray<TRootGenParticle>& NFCentralQ, const vector<int> &NTrackJet)
+vector<fastjet::PseudoJet> JetsUtil::RunJetsResol(const vector<fastjet::PseudoJet>&  input_particles)
+{
+  inclusive_jets.clear();
+  sorted_jets.clear();
+  // run the jet clustering with the above jet definition
+  if(input_particles.size()!=0)
+    {
+      fastjet::ClusterSequence clust_seq(input_particles, jet_def);
+      // extract the inclusive jets with pt > 5 GeV
+      double ptmin = 5.0;
+      inclusive_jets = clust_seq.inclusive_jets(ptmin);
+
+      // sort jets into increasing pt
+      sorted_jets = sorted_by_pt(inclusive_jets);
+    }
+  return sorted_jets;
+}
+
+
+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)
 {
   TRootJet *elementJet;
@@ -175 +202 @@
         elementJet->Set(JET);
         elementJet->NTracks = NTrackJet[i];
+        elementJet->EHoverEE = EHADEEM[i];
 
         // b-jets
@@ -185 +213 @@
 }
 
-void JetsUtil::RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TRootTracks> & TrackCentral, const vector<int> &NTrackJet)
+void JetsUtil::RunTauJets(ExRootTreeWriter *treeWriter, ExRootTreeBranch *branchTauJet,const vector<fastjet::PseudoJet> & sorted_jets,const vector<PhysicsTower> & towers, const vector<TRootTracks> & TrackCentral, const vector<int> &NTrackJet, const vector<float> &EHADEEM)
 {
   TRootTauJet *elementTauJet;
@@ -205 +233 @@
         elementTauJet->NTracks=NTrackJet[i];
         elementTauJet->Charge = charge;
+        elementTauJet->EHoverEE = EHADEEM[i];
       } // if tau jet
     } // if JET.eta < tracker - tau_cone : Tau jet identification