Changeset f01dee6 in git for examples/Validation.C

Timestamp:

Jun 28, 2016, 4:55:13 PM (8 years ago)

Author:

Alexandre Mertens <alexandre.mertens@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

90ccc5d

Parents:

a3a556c

Message:

Adding script for running the validation

File:

• examples/Validation.C (modified) (6 diffs)

examples/Validation.C

@@ -1 @@
-/*
-This macro shows how to compute jet energy scale.
-root -l examples/Example4.C'("delphes_output.root", "plots.root")'
-
-The output ROOT file contains the pT(MC)/pT(Reco) distributions for
-various pT(Reco) and |eta| bins. The peak value of such distribution is
-interpreted as the jet energy correction to be applied for that
-given pT(Reco), |eta| bin.
-
-This can be done by modifying the "ScaleFormula" input parameter to
-the JetEnergyScale module in the delphes_card_XXX.tcl
-
-e.g  a smooth function:
-
-  set ScaleFormula { sqrt(3.0 - 0.1*(abs(eta)))^2 / pt + 1.0) }
-
-or a binned function:
-
-  set ScaleFormula {(abs(eta) > 0.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 50.0)  * (1.10) +
-                    (abs(eta) > 0.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 100.0) * (1.05) +
-                    (abs(eta) > 0.0 && abs(eta) <= 2.5) * (pt > 100.0)               * (1.00) +
-                    (abs(eta) > 2.5 && abs(eta) <= 5.0) * (pt > 20.0 && pt <= 50.0)  * (1.10) +
-                    (abs(eta) > 2.5 && abs(eta) <= 5.0) * (pt > 50.0 && pt <= 100.0) * (1.05) +
-                    (abs(eta) > 2.5 && abs(eta) <= 5.0) * (pt > 100.0)               * (1.00)}
-
-Be aware that a binned jet energy scale can produce "steps" in the corrected
-jet pt distribution ...
-*/
-
 #ifdef __CLING__
 R__LOAD_LIBRARY(libDelphes)
@@ -304 +275 @@
             if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta > 0.0 && eta < 2.5) 
             {
-                if (eta < 1.5) {histos->at(bin).cenResolHist->Fill((bestGenMomentum.E()-recoMomentum.E())/bestGenMomentum.E());}
-                else if (eta < 2.5) {histos->at(bin).fwdResolHist->Fill((bestGenMomentum.E()-recoMomentum.E())/bestGenMomentum.E());}
+                if (eta < 1.5) {histos->at(bin).cenResolHist->Fill(recoMomentum.Pt()/bestGenMomentum.Pt());}
+                else if (eta < 2.5) {histos->at(bin).fwdResolHist->Fill(recoMomentum.Pt()/bestGenMomentum.Pt());}
             }
         }
@@ -660 +631 @@
   gDirectory->cd(0);
 
-/*
   ///////////
   // Muons //
@@ -796 +766 @@
   C_ph->SaveAs(phRes+".eps");
 
-*/
   //////////
   // Jets //
   //////////
 
-  
   // PFJets Energy Resolution
   std::vector<resolPlot> plots_pfjets;
@@ -856 +824 @@
   }
 
-/*
+
   //  gr.Write();
   histos_el.first->Write();
@@ -875 +843 @@
   //gr_eltrack.Write();
   //gr_eltower.Write();
-*/
+
 
   C_el1->Write();
   C_el2->Write();
   C_jet->Write();
-/*
+
   C_mu->Write();
   C_ph->Write();
-*/
+
   gr_pfjets.Write();
   gr_calojets.Write();