Fork me on GitHub

Changeset f01dee6 in git for examples


Ignore:
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

Location:
examples
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • examples/Pythia8/configParticleGun.cmnd

    ra3a556c rf01dee6  
    33! 1) Settings used in the main program.
    44
    5 Main:numberOfEvents = 10000         ! number of events to generate
     5Main:numberOfEvents = 100000         ! number of events to generate
    66Main:timesAllowErrors = 3          ! how many aborts before run stops
    77Main:spareFlag1 = on                ! true means particle gun
    8 Main:spareMode1 = 11               ! 1-5 - di-quark, 21 - di-gluon, 11 - single electron, 13 - single muon, 22 - single photon
     8Main:spareMode1 = ID               ! 1-5 - di-quark, 21 - di-gluon, 11 - single electron, 13 - single muon, 22 - single photon
    99Main:spareParm1 = 10000           ! max pt
    1010
  • examples/Validation.C

    ra3a556c rf01dee6  
    1 /*
    2 This macro shows how to compute jet energy scale.
    3 root -l examples/Example4.C'("delphes_output.root", "plots.root")'
    4 
    5 The output ROOT file contains the pT(MC)/pT(Reco) distributions for
    6 various pT(Reco) and |eta| bins. The peak value of such distribution is
    7 interpreted as the jet energy correction to be applied for that
    8 given pT(Reco), |eta| bin.
    9 
    10 This can be done by modifying the "ScaleFormula" input parameter to
    11 the JetEnergyScale module in the delphes_card_XXX.tcl
    12 
    13 e.g  a smooth function:
    14 
    15   set ScaleFormula { sqrt(3.0 - 0.1*(abs(eta)))^2 / pt + 1.0) }
    16 
    17 or a binned function:
    18 
    19   set ScaleFormula {(abs(eta) > 0.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 50.0)  * (1.10) +
    20                     (abs(eta) > 0.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 100.0) * (1.05) +
    21                     (abs(eta) > 0.0 && abs(eta) <= 2.5) * (pt > 100.0)               * (1.00) +
    22                     (abs(eta) > 2.5 && abs(eta) <= 5.0) * (pt > 20.0 && pt <= 50.0)  * (1.10) +
    23                     (abs(eta) > 2.5 && abs(eta) <= 5.0) * (pt > 50.0 && pt <= 100.0) * (1.05) +
    24                     (abs(eta) > 2.5 && abs(eta) <= 5.0) * (pt > 100.0)               * (1.00)}
    25 
    26 Be aware that a binned jet energy scale can produce "steps" in the corrected
    27 jet pt distribution ...
    28 */
    29 
    301#ifdef __CLING__
    312R__LOAD_LIBRARY(libDelphes)
     
    304275            if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta > 0.0 && eta < 2.5)
    305276            {
    306                 if (eta < 1.5) {histos->at(bin).cenResolHist->Fill((bestGenMomentum.E()-recoMomentum.E())/bestGenMomentum.E());}
    307                 else if (eta < 2.5) {histos->at(bin).fwdResolHist->Fill((bestGenMomentum.E()-recoMomentum.E())/bestGenMomentum.E());}
     277                if (eta < 1.5) {histos->at(bin).cenResolHist->Fill(recoMomentum.Pt()/bestGenMomentum.Pt());}
     278                else if (eta < 2.5) {histos->at(bin).fwdResolHist->Fill(recoMomentum.Pt()/bestGenMomentum.Pt());}
    308279            }
    309280        }
     
    660631  gDirectory->cd(0);
    661632
    662 /*
    663633  ///////////
    664634  // Muons //
     
    796766  C_ph->SaveAs(phRes+".eps");
    797767
    798 */
    799768  //////////
    800769  // Jets //
    801770  //////////
    802771
    803  
    804772  // PFJets Energy Resolution
    805773  std::vector<resolPlot> plots_pfjets;
     
    856824  }
    857825
    858 /*
     826
    859827  //  gr.Write();
    860828  histos_el.first->Write();
     
    875843  //gr_eltrack.Write();
    876844  //gr_eltower.Write();
    877 */
     845
    878846
    879847  C_el1->Write();
    880848  C_el2->Write();
    881849  C_jet->Write();
    882 /*
     850
    883851  C_mu->Write();
    884852  C_ph->Write();
    885 */
     853
    886854  gr_pfjets.Write();
    887855  gr_calojets.Write();
Note: See TracChangeset for help on using the changeset viewer.