Fork me on GitHub

Changes in / [a7c9002:837fa70] in git


Ignore:
Files:
1 deleted
13 edited

Legend:

Unmodified
Added
Removed

  • cards/FCC/FCChh.tcl

    ra7c9002 r837fa70  
    1515
    1616set ExecutionPath {
    17 
    1817  ParticlePropagator
    1918
    … …  
    3332  Calorimeter
    3433  EFlowMerger
    35   EFlowFilter
    36  
     34
    3735  PhotonEfficiency
    3836  PhotonIsolation
    … …  
    6462
    6563  UniqueObjectFinder
    66  
    67   GenParticleFilter
    6864
    6965  TreeWriter
    … …  
    414410}
    415411
    416 ######################
    417 # EFlowFilter
    418 ######################
    419 
    420 module PdgCodeFilter EFlowFilter {
    421   set InputArray EFlowMerger/eflow
    422   set OutputArray eflow
    423  
    424   add PdgCode {11}
    425   add PdgCode {-11}
    426   add PdgCode {13}
    427   add PdgCode {-13}
    428 }
    429 
    430412
    431413###################
    … …  
    591573module Isolation PhotonIsolation {
    592574  set CandidateInputArray PhotonEfficiency/photons
    593   set IsolationInputArray EFlowFilter/eflow
     575  set IsolationInputArray EFlowMerger/eflow
    594576
    595577  set OutputArray photons
    … …  
    611593module Isolation ElectronIsolation {
    612594  set CandidateInputArray ElectronFilter/electrons
    613   set IsolationInputArray EFlowFilter/eflow
     595  set IsolationInputArray EFlowMerger/eflow
    614596
    615597  set OutputArray electrons
    … …  
    631613module Isolation MuonIsolation {
    632614  set CandidateInputArray MuonMomentumSmearing/muons
    633   set IsolationInputArray EFlowFilter/eflow
     615  set IsolationInputArray EFlowMerger/eflow
    634616
    635617  set OutputArray muons
    … …  
    655637
    656638  (pt <= 10.0)                       * (0.00) +
    657   (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 500)      * (0.001) + \
    658   (abs(eta) < 2.5)                   * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
    659   (abs(eta) < 2.5)                   * (pt > 20000.0)               * (0.000) + \
    660   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.00075) + \
    661   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
    662   (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0)               * (0.000) + \
     639  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 500) * (0.001) + \
     640  (abs(eta) < 2.5)                   * (pt > 500.0)            * (0.001)*(1.0 - pt/50000.) + \
     641  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \
     642  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0)            * (0.00075)*(1.0 - pt/50000.) + \
    663643  (abs(eta) > 4.0) * (0.00)}
    664644
    … …  
    666646
    667647  (pt <= 10.0)                       * (0.00) +
    668   (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 500)      * (0.04) + \
    669   (abs(eta) < 2.5)                   * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
    670   (abs(eta) < 2.5)                   * (pt > 20000.0)               * (0.000) + \
    671   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.03) + \
    672   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
    673   (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0)               * (0.000) + \
     648  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 500) * (0.04) + \
     649  (abs(eta) < 2.5)                   * (pt > 500.0)            * (0.04)*(1.0 - pt/50000.) + \
     650  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
     651  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0)            * (0.03)*(1.0 - pt/50000.) + \
    674652  (abs(eta) > 4.0) * (0.00)}
    675653
    676654  add EfficiencyFormula {5} {
    677655
    678   (pt <= 10.0)                                                       * (0.00) +
    679   (abs(eta) < 2.5)                    * (pt > 10.0 && pt < 500)      * (0.85) +
    680   (abs(eta) < 2.5)                    * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
    681   (abs(eta) < 2.5)                    * (pt > 20000.0)               * (0.000) +
    682   (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.64) +
    683   (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
    684   (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0)               * (0.000) +
    685   (abs(eta) >= 4.0) * (0.00)}
     656  (pt <= 10.0)                       * (0.00) +
     657  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 500) * (0.85) + \
     658  (abs(eta) < 2.5)                   * (pt > 500.0)            * (0.85)*(1.0 - pt/50000.) + \
     659  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) + \
     660  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0)            * (0.64)*(1.0 - pt/50000.) + \
     661  (abs(eta) > 4.0) * (0.00)}
    686662
    687663}
    … …  
    736712  add EfficiencyFormula {0} {
    737713
    738   (pt <= 10.0)                                                       * (0.00) +
    739   (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0)    * (0.01) + \
    740   (abs(eta) < 2.5)                   * (pt > 5000.0 && pt < 34000.0) * (0.01)  *(8./9. - pt/30000.) + \
    741   (abs(eta) < 2.5)                   * (pt > 34000.0)                * (0.000) + \
    742   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0)    * (0.0075) + \
    743   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
    744   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0)                * (0.00) + \
     714  (pt <= 10.0)                                                    * (0.00) +
     715  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0) * (0.01) + \
     716  (abs(eta) < 2.5)                   * (pt > 5000.0)              * (0.01)  *(8./9. - pt/45000.) + \
     717  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
     718  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0)              * (0.0075)*(8./9. - pt/45000.) + \
    745719  (abs(eta) > 4.0)                   * (0.00)}
    746720
    747721  add EfficiencyFormula {11} {
    748722
    749   (pt <= 10.0)                                                       * (0.00) +
    750   (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0)    * (0.005) + \
    751   (abs(eta) < 2.5)                   * (pt > 5000.0 && pt < 34000.0) * (0.005)  *(8./9. - pt/30000.) + \
    752   (abs(eta) < 2.5)                   * (pt > 34000.0)                * (0.000) + \
    753   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0)    * (0.00375) + \
    754   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
    755   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0)                * (0.00) + \
     723  (pt <= 10.0)                                                    * (0.00) +
     724  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0) * (0.005) + \
     725  (abs(eta) < 2.5)                   * (pt > 5000.0)              * (0.005)    *(8./9. - pt/45000.) + \
     726  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
     727  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0)              * (0.00375)*(8./9. - pt/45000.) + \
    756728  (abs(eta) > 4.0)                   * (0.00)}
    757729
    758730  add EfficiencyFormula {15} {
    759731
    760   (pt <= 10.0)                                                       * (0.00) +
    761   (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0)                      * (0.6)                      + \
    762   (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0)                   * (0.6) *(8./9. - pt/30000.) + \
    763   (abs(eta) < 2.5)                   * (pt > 34000.0)                * (0.000) + \
    764   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0)    * (0.45) + \
    765   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
    766   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0)                * (0.00) + \
    767   (abs(eta) > 4.0)                                                   * (0.00)}
     732  (pt <= 10.0)                                                    * (0.00) +
     733  (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0)                   * (0.6)                      + \
     734  (abs(eta) < 2.5) * (pt > 5000.0)                                * (0.6) *(8./9. - pt/45000.) + \
     735  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
     736  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0)              * (0.45)*(8./9. - pt/45000.) + \
     737  (abs(eta) > 4.0)                                                * (0.00)}
    768738
    769739}
    … …  
    783753}
    784754
    785 ######################
    786 # GenParticleFilter
    787 ######################
    788 
    789 # store only interesting particles
    790 module StatusPidFilter GenParticleFilter {
    791   set InputArray Delphes/allParticles
    792   set OutputArray filteredParticles
    793   set PTMin 0.0
    794 
    795 }
    796755
    797756##################
    … …  
    802761# add Branch InputArray BranchName BranchClass
    803762  add Branch Delphes/allParticles Particle GenParticle
    804   add Branch GenParticleFilter/filteredParticles FilteredParticle GenParticle
    805763
    806764  add Branch GenJetFinder/jets GenJet Jet

  • cards/delphes_card_ATLAS.tcl

    ra7c9002 r837fa70  
    2121  Calorimeter
    2222  EFlowMerger
    23   EFlowFilter
    24  
     23
    2524  PhotonEfficiency
    2625  PhotonIsolation
    … …  
    406405}
    407406
    408 ######################
    409 # EFlowFilter
    410 ######################
    411 
    412 module PdgCodeFilter EFlowFilter {
    413   set InputArray EFlowMerger/eflow
    414   set OutputArray eflow
    415  
    416   add PdgCode {11}
    417   add PdgCode {-11}
    418   add PdgCode {13}
    419   add PdgCode {-13}
    420 }
    421 
    422407###################
    423408# Photon efficiency
    … …  
    443428module Isolation PhotonIsolation {
    444429  set CandidateInputArray PhotonEfficiency/photons
    445   set IsolationInputArray EFlowFilter/eflow
     430  set IsolationInputArray EFlowMerger/eflow
    446431
    447432  set OutputArray photons
    … …  
    454439}
    455440
     441#################
     442# Electron filter
     443#################
     444
     445module PdgCodeFilter ElectronFilter {
     446  set InputArray HCal/eflowTracks
     447  set OutputArray electrons
     448  set Invert true
     449  add PdgCode {11}
     450  add PdgCode {-11}
     451}
    456452
    457453#####################
    … …  
    478474module Isolation ElectronIsolation {
    479475  set CandidateInputArray ElectronEfficiency/electrons
    480   set IsolationInputArray EFlowFilter/eflow
     476  set IsolationInputArray EFlowMerger/eflow
    481477
    482478  set OutputArray electrons
    … …  
    512508module Isolation MuonIsolation {
    513509  set CandidateInputArray MuonEfficiency/muons
    514   set IsolationInputArray EFlowFilter/eflow
     510  set IsolationInputArray EFlowMerger/eflow
    515511
    516512  set OutputArray muons

  • cards/delphes_card_ATLAS_PileUp.tcl

    ra7c9002 r837fa70  
    2323  EFlowMergerAllTracks
    2424  EFlowMerger
    25   EFlowFilter
    26  
     25
    2726  NeutrinoFilter
    2827  GenJetFinder
    … …  
    373372}
    374373
    375 ######################
    376 # EFlowFilter
    377 ######################
    378 
    379 module PdgCodeFilter EFlowFilter {
    380   set InputArray EFlowMergerAllTracks/eflow
    381   set OutputArray eflow
    382  
    383   add PdgCode {11}
    384   add PdgCode {-11}
    385   add PdgCode {13}
    386   add PdgCode {-13}
    387 }
    388 
    389 
    390374#############
    391375# Rho pile-up
    … …  
    524508module Isolation PhotonIsolation {
    525509  set CandidateInputArray PhotonEfficiency/photons
    526   set IsolationInputArray EFlowFilter/eflow
     510  set IsolationInputArray EFlowMergerAllTracks/eflow
    527511  set RhoInputArray Rho/rho
    528512
    … …  
    559543module Isolation ElectronIsolation {
    560544  set CandidateInputArray ElectronEfficiency/electrons
    561   set IsolationInputArray EFlowFilter/eflow
     545  set IsolationInputArray EFlowMergerAllTracks/eflow
    562546  set RhoInputArray Rho/rho
    563547
    … …  
    594578module Isolation MuonIsolation {
    595579  set CandidateInputArray MuonEfficiency/muons
    596   set IsolationInputArray EFlowFilter/eflow
     580  set IsolationInputArray EFlowMergerAllTracks/eflow
    597581  set RhoInputArray Rho/rho
    598582

  • cards/delphes_card_CMS.tcl

    ra7c9002 r837fa70  
    2121  Calorimeter
    2222  EFlowMerger
    23   EFlowFilter
    24  
     23
    2524  PhotonEfficiency
    2625  PhotonIsolation
    … …  
    318317
    319318  set EnergyMin 1.0
    320   set EnergySignificanceMin 1.0
     319  set EnergySignificanceMin 2.0
    321320
    322321  set SmearTowerCenter true
    … …  
    419418}
    420419
    421 ######################
    422 # EFlowFilter
    423 ######################
    424 
    425 module PdgCodeFilter EFlowFilter {
    426   set InputArray EFlowMerger/eflow
    427   set OutputArray eflow
    428  
    429   add PdgCode {11}
    430   add PdgCode {-11}
    431   add PdgCode {13}
    432   add PdgCode {-13}
    433 }
    434 
    435 
    436420###################
    437421# Photon efficiency
    … …  
    457441module Isolation PhotonIsolation {
    458442  set CandidateInputArray PhotonEfficiency/photons
    459   set IsolationInputArray EFlowFilter/eflow
     443  set IsolationInputArray EFlowMerger/eflow
    460444
    461445  set OutputArray photons
    … …  
    492476module Isolation ElectronIsolation {
    493477  set CandidateInputArray ElectronEfficiency/electrons
    494   set IsolationInputArray EFlowFilter/eflow
     478  set IsolationInputArray EFlowMerger/eflow
    495479
    496480  set OutputArray electrons
    … …  
    526510module Isolation MuonIsolation {
    527511  set CandidateInputArray MuonEfficiency/muons
    528   set IsolationInputArray EFlowFilter/eflow
     512  set IsolationInputArray EFlowMerger/eflow
    529513
    530514  set OutputArray muons

  • cards/delphes_card_CMS_PileUp.tcl

    ra7c9002 r837fa70  
    2626  EFlowMergerAllTracks
    2727  EFlowMerger
    28   EFlowFilter
    29  
     28
    3029  NeutrinoFilter
    3130  GenJetFinder
    … …  
    437436}
    438437
    439 ######################
    440 # EFlowFilter
    441 ######################
    442 
    443 module PdgCodeFilter EFlowFilter {
    444   set InputArray EFlowMergerAllTracks/eflow
    445   set OutputArray eflow
    446  
    447   add PdgCode {11}
    448   add PdgCode {-11}
    449   add PdgCode {13}
    450   add PdgCode {-13}
    451 }
    452438
    453439##########################
    … …  
    667653module Isolation PhotonIsolation {
    668654  set CandidateInputArray PhotonEfficiency/photons
    669   set IsolationInputArray EFlowFilter/eflow
     655  set IsolationInputArray EFlowMergerAllTracks/eflow
    670656  set RhoInputArray Rho/rho
    671657
    … …  
    702688module Isolation ElectronIsolation {
    703689  set CandidateInputArray ElectronEfficiency/electrons
    704   set IsolationInputArray EFlowFilter/eflow
     690  set IsolationInputArray EFlowMergerAllTracks/eflow
    705691  set RhoInputArray Rho/rho
    706692
    … …  
    738724module Isolation MuonIsolation {
    739725  set CandidateInputArray MuonEfficiency/muons
    740   set IsolationInputArray EFlowFilter/eflow
     726  set IsolationInputArray EFlowMergerAllTracks/eflow
    741727  set RhoInputArray Rho/rho
    742728

  • cards/delphes_card_ILD.tcl

    ra7c9002 r837fa70  
    2323  Calorimeter
    2424  EFlowMerger
    25   EFlowFilter
    26  
     25
    2726  PhotonEfficiency
    2827  PhotonIsolation
    … …  
    351350}
    352351
    353 ######################
    354 # EFlowFilter
    355 ######################
    356 
    357 module PdgCodeFilter EFlowFilter {
    358   set InputArray EFlowMerger/eflow
    359   set OutputArray eflow
    360  
    361   add PdgCode {11}
    362   add PdgCode {-11}
    363   add PdgCode {13}
    364   add PdgCode {-13}
    365 }
    366 
    367352
    368353###################
    … …  
    506491module Isolation PhotonIsolation {
    507492  set CandidateInputArray PhotonEfficiency/photons
    508   set IsolationInputArray EFlowFilter/eflow
     493  set IsolationInputArray EFlowMerger/eflow
    509494
    510495  set OutputArray photons
    … …  
    540525module Isolation ElectronIsolation {
    541526  set CandidateInputArray ElectronEfficiency/electrons
    542   set IsolationInputArray EFlowFilter/eflow
     527  set IsolationInputArray EFlowMerger/eflow
    543528
    544529  set OutputArray electrons
    … …  
    576561module Isolation MuonIsolation {
    577562  set CandidateInputArray MuonEfficiency/muons
    578   set IsolationInputArray EFlowFilter/eflow
     563  set IsolationInputArray EFlowMerger/eflow
    579564
    580565  set OutputArray muons

  • classes/DelphesHepMCReader.cc

    ra7c9002 r837fa70  
    5151using namespace std;
    5252
    53 static const int kBufferSize = 16384;
     53static const int kBufferSize  = 1024;
    5454
    5555//---------------------------------------------------------------------------

  • classes/DelphesLHEFReader.cc

    ra7c9002 r837fa70  
    4848using namespace std;
    4949
    50 static const int kBufferSize = 16384;
     50static const int kBufferSize  = 1024;
    5151
    5252//---------------------------------------------------------------------------

  • examples/Example1.py

    ra7c9002 r837fa70  
    1010
    1111ROOT.gSystem.Load("libDelphes")
    12 
    13 try:
    14   ROOT.gInterpreter.Declare('#include "classes/DelphesClasses.h"')
    15   ROOT.gInterpreter.Declare('#include "external/ExRootAnalysis/ExRootTreeReader.h"')
    16 except:
    17   pass
    1812
    1913inputFile = sys.argv[1]

  • examples/Validation.cpp

    ra7c9002 r837fa70  
    250250      eta = TMath::Abs(genMomentum.Eta());
    251251
    252 
    253252      if(eta > etamax || eta < etamin ) continue;
    254253
    255       //cout<<"b parton: "<<pt<<endl;
    256254      if (particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
    257255      {
    … …  
    281279            if( !(jet->BTag & (1 << 0)) ) continue;
    282280          }
    283 
    284             if(TMath::Abs(pdgID) == 1)
     281       
     282          if(TMath::Abs(pdgID) == 1)
    285283          {
    286284            Jet *jet = (Jet *)recoObj;
    … …  
    301299          }
    302300        }
     301
    303302        histGenPt->Fill(pt);
    304         if(deltaR < 0.3 && bestRecoMomentum.Pt()> 0.20*pt ) { histRecoPt->Fill(pt);}
     303        if(deltaR < 0.3) { histRecoPt->Fill(pt); }
    305304
    306305      }
    … …  
    13531352  TClonesArray *branchCaloJet = treeReaderJet->UseBranch("CaloJet");
    13541353  TClonesArray *branchJet = treeReaderJet->UseBranch("Jet");
    1355 
     1354 
    13561355  TClonesArray *branchParticleBJet = treeReaderBJet->UseBranch("Particle");
    13571356  TClonesArray *branchPFBJet = treeReaderBJet->UseBranch("Jet");
    … …  
    14051404
    14061405  TPaveText *pave = new TPaveText(0.0, 0.89, 0.94, 0.94,"NDC");
    1407   pave->SetTextAlign(30);
     1406  pave->SetTextAlign(kHAlignRight);
    14081407  pave->SetTextFont(132);
    14091408  pave->SetBorderSize(0);
    … …  
    14401439  const int n_ptbins = ptVals.size();
    14411440
     1441
    14421442  //////////////////////////
    14431443  // Tracking performance //
    … …  
    19061906  TGraphErrors gr_pfele_res_eta[n_ptbins];
    19071907
    1908   TGraphErrors gr_trkele_res_e[n_etabins];
    1909   TGraphErrors gr_trkele_res_eeta[n_ptbins];
    1910 
    19111908  std::vector<resolPlot> plots_pfele_res_e[n_etabins], plots_pfele_res_eta[n_ptbins];
    1912   std::vector<resolPlot> plots_trkele_res_e[n_etabins], plots_trkele_res_eeta[n_ptbins];
    19131909
    19141910  TCanvas *c_pfele_res_e[n_etabins];
    … …  
    19261922     gr_pfele_res_e[k] = EresGraph(&plots_pfele_res_e[k]);
    19271923
    1928      HistogramsCollection(&plots_trkele_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkele");
    1929      GetEres<Track>(&plots_trkele_res_e[k], branchTrackElectron, branchParticleElectron, 11, etaVals.at(k), etaVals.at(k+1), treeReaderElectron);
    1930      gr_trkele_res_e[k] = EresGraph(&plots_trkele_res_e[k]);
    1931 
    19321924     s_etaMin = Form("%.1f",etaVals.at(k));
    19331925     s_etaMax = Form("%.1f",etaVals.at(k+1));
    … …  
    19381930
    19391931     addResoGraph(mg_pfele_res_e[k], &gr_ecal_res_e[k], leg_pfele_res_e[k], markerStyles.at(0), colors.at(0), "ECAL");
    1940      addResoGraph(mg_pfele_res_e[k], &gr_trkele_res_e[k], leg_pfele_res_e[k], markerStyles.at(1), colors.at(1), "Track");
     1932     addResoGraph(mg_pfele_res_e[k], &gr_trkele_res_pt[k], leg_pfele_res_e[k], markerStyles.at(1), colors.at(1), "Track");
    19411933     addResoGraph(mg_pfele_res_e[k], &gr_pfele_res_e[k], leg_pfele_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    19421934
    … …  
    19691961     gr_pfele_res_eta[k] = EresGraphVsEta(&plots_pfele_res_eta[k]);
    19701962
    1971      HistogramsCollectionVsEta(&plots_trkele_res_eeta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "trkele", 0.0, 2.0);
    1972      GetEresVsEta<Track>(&plots_trkele_res_eeta[k], branchTrackElectron, branchParticleElectron, 11, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderElectron);
    1973      gr_trkele_res_eeta[k] = EresGraphVsEta(&plots_trkele_res_eeta[k]);
    1974 
    19751963     s_e = Form("e^{ #pm}, E = %.0f GeV",ptVals.at(k));
    19761964     if(ptVals.at(k) >= 1000.) s_e = Form("e^{ #pm}, E = %.0f TeV",ptVals.at(k)/1000.);
    … …  
    19811969
    19821970     addResoGraph(mg_pfele_res_eta[k], &gr_ecal_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(0), colors.at(0), "ECAL");
    1983      addResoGraph(mg_pfele_res_eta[k], &gr_trkele_res_eeta[k], leg_pfele_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
     1971     addResoGraph(mg_pfele_res_eta[k], &gr_trkele_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
    19841972     addResoGraph(mg_pfele_res_eta[k], &gr_pfele_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    19851973
    … …  
    20132001  TGraphErrors gr_pfpi_res_eta[n_ptbins];
    20142002
    2015   TGraphErrors gr_trkpi_res_e[n_etabins];
    2016   TGraphErrors gr_trkpi_res_eeta[n_ptbins];
    2017 
    20182003  std::vector<resolPlot> plots_pfpi_res_e[n_etabins], plots_pfpi_res_eta[n_ptbins];
    2019   std::vector<resolPlot> plots_trkpi_res_e[n_etabins], plots_trkpi_res_eeta[n_ptbins];
    20202004
    20212005  TCanvas *c_pfpi_res_e[n_etabins];
    … …  
    20332017     gr_pfpi_res_e[k] = EresGraph(&plots_pfpi_res_e[k]);
    20342018
    2035      HistogramsCollection(&plots_trkpi_res_e[k], TMath::Log10(ptMin), TMath::Log10(ptMax), "trkpi");
    2036      GetEres<Track>(&plots_trkpi_res_e[k], branchTrackPion, branchParticlePion, 211, etaVals.at(k), etaVals.at(k+1), treeReaderPion);
    2037      gr_trkpi_res_e[k] = EresGraph(&plots_trkpi_res_e[k]);
    2038 
    2039 
    20402019     s_etaMin = Form("%.1f",etaVals.at(k));
    20412020     s_etaMax = Form("%.1f",etaVals.at(k+1));
    … …  
    20462025
    20472026     addResoGraph(mg_pfpi_res_e[k], &gr_hcal_res_e[k], leg_pfpi_res_e[k], markerStyles.at(0), colors.at(0), "HCAL");
    2048      addResoGraph(mg_pfpi_res_e[k], &gr_trkpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(1), colors.at(1), "Track");
     2027     addResoGraph(mg_pfpi_res_e[k], &gr_trkpi_res_pt[k], leg_pfpi_res_e[k], markerStyles.at(1), colors.at(1), "Track");
    20492028     addResoGraph(mg_pfpi_res_e[k], &gr_pfpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    20502029
    … …  
    20772056     gr_pfpi_res_eta[k] = EresGraphVsEta(&plots_pfpi_res_eta[k]);
    20782057
    2079      HistogramsCollectionVsEta(&plots_trkpi_res_eeta[k], etaMin, etaMax, 0.5*ptVals.at(k), 2.0*ptVals.at(k), "trkpi", 0.0, 2.0);
    2080      GetEresVsEta<Track>(&plots_trkpi_res_eeta[k], branchPion, branchParticlePion, 211, 0.5*ptVals.at(k), 2.0*ptVals.at(k), treeReaderPion);
    2081      gr_trkpi_res_eeta[k] = EresGraphVsEta(&plots_trkpi_res_eeta[k]);
    2082 
    2083 
    20842058     s_e = Form("#pi^{ #pm}, E = %.0f GeV",ptVals.at(k));
    20852059     if(ptVals.at(k) >= 1000.) s_e = Form("#pi^{ #pm}, E = %.0f TeV",ptVals.at(k)/1000.);
    … …  
    20892063
    20902064     addResoGraph(mg_pfpi_res_eta[k], &gr_hcal_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(0), colors.at(0), "HCAL");
    2091      addResoGraph(mg_pfpi_res_eta[k], &gr_trkpi_res_eeta[k], leg_pfpi_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
     2065     addResoGraph(mg_pfpi_res_eta[k], &gr_trkpi_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
    20922066     addResoGraph(mg_pfpi_res_eta[k], &gr_pfpi_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    20932067
    … …  
    24442418    c_recpho_eff_eta->Print(figPath+"img_recpho_eff_eta.png","png");
    24452419
    2446 
    24472420    /////////////////////////////////////////
    24482421    // B-jets  Efficiency/ mistag rates   ///
    … …  
    25092482    c_recbjet_eff_eta->Print(figPath+"img_recbjet_eff_eta.png","png");
    25102483
     2484
     2485
    25112486    // ------ c - mistag  ------
    25122487
    … …  
    25852560    for (k = 0; k < etaVals.size()-1; k++)
    25862561    {
    2587 
     2562       
    25882563       h_recbjet_lmis_pt = GetEffPt<Jet>(branchJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
    25892564       gr_recbjet_lmis_pt[k] = TGraphErrors(h_recbjet_lmis_pt);
    … …  
    26142589
    26152590    mg_recbjet_lmis_pt->Draw("APE");
    2616     DrawAxis(mg_recbjet_lmis_pt, leg_recbjet_lmis_pt, ptMin, ptMax, 0.0, 0.5, "p_{T} [GeV]", "light - mistag rate (%)", true, false);
     2591    DrawAxis(mg_recbjet_lmis_pt, leg_recbjet_lmis_pt, ptMin, ptMax, 0.0, 0.005, "p_{T} [GeV]", "light - mistag rate (%)", true, false);
    26172592    leg_recbjet_lmis_pt->Draw();
    26182593    pave->Draw();
    … …  
    26252600
    26262601    mg_recbjet_lmis_eta->Draw("APE");
    2627     DrawAxis(mg_recbjet_lmis_eta, leg_recbjet_lmis_eta, etaMin, etaMax, 0.0, 0.5, " #eta ", "light - mistag rate (%)", false, false);
     2602    DrawAxis(mg_recbjet_lmis_eta, leg_recbjet_lmis_eta, etaMin, etaMax, 0.0, 0.005, " #eta ", "light - mistag rate (%)", false, false);
    26282603    leg_recbjet_lmis_eta->Draw();
    26292604    pave->Draw();
    … …  
    27432718
    27442719    mg_rectaujet_mis_pt->Draw("APE");
    2745     DrawAxis(mg_rectaujet_mis_pt, leg_rectaujet_mis_pt, ptMin, ptMax, 0.0, 5., "p_{T} [GeV]", "#tau - mistag(%)", true, false);
     2720    DrawAxis(mg_rectaujet_mis_pt, leg_rectaujet_mis_pt, ptMin, ptMax, 0.0, 10., "p_{T} [GeV]", "#tau - mistag(%)", true, false);
    27462721    leg_rectaujet_mis_pt->Draw();
    27472722    pave->Draw();
    … …  
    27752750       plots_trkpi_res_pt[k].at(bin).resolHist->Write();
    27762751       plots_trkele_res_pt[k].at(bin).resolHist->Write();
    2777        plots_trkmu_res_pt[k].at(bin).resolHist->Write();
     2752           plots_trkmu_res_pt[k].at(bin).resolHist->Write();
    27782753       plots_ecal_res_e[k].at(bin).resolHist->Write();
    27792754       plots_hcal_res_e[k].at(bin).resolHist->Write();
    … …  
    28052780
    28062781  fout->Write();
    2807 
    2808 
    2809 
    28102782
    28112783  cout << "** Exiting..." << endl;

  • external/fastjet/CMakeLists.txt

    ra7c9002 r837fa70  
    1 include_directories(BEFORE ${CMAKE_SOURCE_DIR}/external)
     1include_directories(${CMAKE_SOURCE_DIR}/external)
    22
    33file(GLOB sources *.cc plugins/*/*.cc contribs/*/*.cc tools/*.cc)

  • modules/CMakeLists.txt

    ra7c9002 r837fa70  
    1818
    1919add_library(modules OBJECT ${sources} FastJetDict.cxx ModulesDict.cxx)
    20 
    21 # install public headers
    22 install(FILES Delphes.h DESTINATION include/modules)

  • modules/StatusPidFilter.cc

    ra7c9002 r837fa70  
    104104    pass = kFALSE;
    105105
    106     // hard scattering particles (first condition for Py6, second for Py8)
     106    // status == 3
    107107    if(status == 3) pass = kTRUE;
    108     if(status > 20 && status < 30 ) pass = kTRUE;
    109108
    110109    // electrons, muons, taus and neutrinos
    … …  
    112111
    113112    // heavy quarks
    114     if(pdgCode == 4 ||pdgCode == 5 || pdgCode == 6) pass = kTRUE;
     113    if(pdgCode == 5 || pdgCode == 6) pass = kTRUE;
    115114
    116115    // Gauge bosons and other fundamental bosons
    117116    if(pdgCode > 22 && pdgCode < 43) pass = kTRUE;
    118117
    119     if(!pass || candidate->Momentum.Pt() < fPTMin) continue;
     118    if(!pass || candidate->Momentum.Pt() <= fPTMin) continue;
    120119
    121120    fOutputArray->Add(candidate);
Note: See TracChangeset for help on using the changeset viewer.