Changeset 06a83b3 in git for examples

Timestamp:

Oct 7, 2016, 10:54:51 AM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

386f526, cb159ed

Parents:

f3be533

Message:

fixed very high energy b and tau tagging, and light mistag collection

File:

• examples/Validation.cpp (modified) (7 diffs)

examples/Validation.cpp

@@ -263 +263 @@
         //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
 
+
+
           // take the closest parton candidate
           if(TMath::Abs(pdgID) == 5)
@@ -277 +279 @@
             if( !(jet->BTag & (1 << 0)) ) continue;
           }
-
-           if(TMath::Abs(pdgID) == 1)
+        
+          if(TMath::Abs(pdgID) == 1)
           {
             Jet *jet = (Jet *)recoObj;
-            if( !(jet->BTag & (1 << 0)) ) continue;
+            if( !(jet->BTag & (1 << 0)) ) continue;
           }
 
@@ -1349 +1351 @@
   TClonesArray *branchPFJet = treeReaderJet->UseBranch("PFJet");
   TClonesArray *branchCaloJet = treeReaderJet->UseBranch("CaloJet");
-
+  TClonesArray *branchJet = treeReaderJet->UseBranch("Jet");
+  
   TClonesArray *branchParticleBJet = treeReaderBJet->UseBranch("Particle");
   TClonesArray *branchPFBJet = treeReaderBJet->UseBranch("Jet");
@@ -2557 +2560 @@
     for (k = 0; k < etaVals.size()-1; k++)
     {
-
-       h_recbjet_lmis_pt = GetEffPt<Jet>(branchPFJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
+       
+       h_recbjet_lmis_pt = GetEffPt<Jet>(branchJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
        gr_recbjet_lmis_pt[k] = TGraphErrors(h_recbjet_lmis_pt);
 
@@ -2574 +2577 @@
     for (k = 0; k < ptVals.size(); k++)
     {
-       h_recbjet_lmis_eta = GetEffEta<Jet>(branchPFJet, branchParticleJet, "Jet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
+       h_recbjet_lmis_eta = GetEffEta<Jet>(branchJet, branchParticleJet, "Jet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
        gr_recbjet_lmis_eta[k] = TGraphErrors(h_recbjet_lmis_eta);
 
@@ -2687 +2690 @@
     {
 
-       h_rectaujet_mis_pt = GetTauEffPt<Jet>(branchPFJet, branchParticleJet, "TauJet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
+       h_rectaujet_mis_pt = GetTauEffPt<Jet>(branchJet, branchParticleJet, "TauJet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
        gr_rectaujet_mis_pt[k] = TGraphErrors(h_rectaujet_mis_pt);
 
@@ -2703 +2706 @@
     for (k = 0; k < ptVals.size(); k++)
     {
-       h_rectaujet_mis_eta = GetTauEffEta<Jet>(branchPFJet, branchParticleJet, "TauJet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
+       h_rectaujet_mis_eta = GetTauEffEta<Jet>(branchJet, branchParticleJet, "TauJet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
        gr_rectaujet_mis_eta[k] = TGraphErrors(h_rectaujet_mis_eta);