Changeset 310 in svn for trunk/src

Timestamp:

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

Author:

severine ovyn

Message:

add Ehad/Eem

Location:

trunk/src

Files:

• JetsUtil.cc (modified) (8 diffs)
• LHCOConverter.cc (modified) (1 diff)

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

trunk/src/LHCOConverter.cc

@@ -174 +174 @@
         particle = (TRootParticle*) branch->At(i);
         double temp =  sqrt(particle->PT*particle->PT + particle->Pz * particle->Pz);
-        double jmass = (particle->E - temp) * (particle->E + temp); // prevents some numerical sensitivity
         TLorentzVector pmu;
         pmu.SetPtEtaPhiE(particle->PT, particle->Eta, particle->Phi, particle->E);
-        jmass = pmu.M();
+        float jmass = pmu.M();
         float ntrk = 0.0;
         float btag =0;
         double ratioE = 0;
 
-        if(lhcoID == lhcoElectronID) { TRootElectron elec(*((TRootElectron*) branch->At(i))); ntrk = elec.Charge; }
-        else if (lhcoID == lhcoMuonID) { TRootMuon muon(*((TRootMuon*) branch->At(i))); ntrk = muon.Charge; }
+        if(lhcoID == lhcoPhotonID){TRootPhoton gam(*((TRootPhoton*) branch->At(i))); jmass=0; ratioE = gam.EHoverEE;}
+        if(lhcoID == lhcoElectronID) { TRootElectron elec(*((TRootElectron*) branch->At(i))); ntrk = elec.Charge; ratioE = elec.EHoverEE;}
+        else if (lhcoID == lhcoMuonID) { TRootMuon muon(*((TRootMuon*) branch->At(i))); ntrk = muon.Charge;  ratioE = muon.EHoverEE;}
         else if (lhcoID == lhcoTauJetID) { TRootTauJet taujet(*((TRootTauJet*) branch->At(i))); ntrk = taujet.Charge; }
-        else if (lhcoID == lhcoJetID) { TRootJet jet(*((TRootJet*) branch->At(i))); ntrk = jet.NTracks; btag = jet.Btag;}
+        else if (lhcoID == lhcoJetID) { TRootJet jet(*((TRootJet*) branch->At(i))); ntrk = jet.NTracks; btag = jet.Btag; ratioE = jet.EHoverEE;}
         if(lhcoID != lhcoETmisID) {
         outfile << fixed << setprecision(3)