Fork me on GitHub

Changes in / [837fa70:a7c9002] in git


Ignore:
Files:
1 added
13 edited

Legend:

Unmodified
Added
Removed

  • cards/FCC/FCChh.tcl

    r837fa70 ra7c9002  
    1515
    1616set ExecutionPath {
     17
    1718  ParticlePropagator
    1819
    … …  
    3233  Calorimeter
    3334  EFlowMerger
    34 
     35  EFlowFilter
     36 
    3537  PhotonEfficiency
    3638  PhotonIsolation
    … …  
    6264
    6365  UniqueObjectFinder
     66 
     67  GenParticleFilter
    6468
    6569  TreeWriter
    … …  
    410414}
    411415
     416######################
     417# EFlowFilter
     418######################
     419
     420module 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
    412430
    413431###################
    … …  
    573591module Isolation PhotonIsolation {
    574592  set CandidateInputArray PhotonEfficiency/photons
    575   set IsolationInputArray EFlowMerger/eflow
     593  set IsolationInputArray EFlowFilter/eflow
    576594
    577595  set OutputArray photons
    … …  
    593611module Isolation ElectronIsolation {
    594612  set CandidateInputArray ElectronFilter/electrons
    595   set IsolationInputArray EFlowMerger/eflow
     613  set IsolationInputArray EFlowFilter/eflow
    596614
    597615  set OutputArray electrons
    … …  
    613631module Isolation MuonIsolation {
    614632  set CandidateInputArray MuonMomentumSmearing/muons
    615   set IsolationInputArray EFlowMerger/eflow
     633  set IsolationInputArray EFlowFilter/eflow
    616634
    617635  set OutputArray muons
    … …  
    637655
    638656  (pt <= 10.0)                       * (0.00) +
    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.) + \
     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) + \
    643663  (abs(eta) > 4.0) * (0.00)}
    644664
    … …  
    646666
    647667  (pt <= 10.0)                       * (0.00) +
    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.) + \
     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) + \
    652674  (abs(eta) > 4.0) * (0.00)}
    653675
    654676  add EfficiencyFormula {5} {
    655677
    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)}
     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)}
    662686
    663687}
    … …  
    712736  add EfficiencyFormula {0} {
    713737
    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.) + \
     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) + \
    719745  (abs(eta) > 4.0)                   * (0.00)}
    720746
    721747  add EfficiencyFormula {11} {
    722748
    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.) + \
     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) + \
    728756  (abs(eta) > 4.0)                   * (0.00)}
    729757
    730758  add EfficiencyFormula {15} {
    731759
    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)}
     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)}
    738768
    739769}
    … …  
    753783}
    754784
     785######################
     786# GenParticleFilter
     787######################
     788
     789# store only interesting particles
     790module StatusPidFilter GenParticleFilter {
     791  set InputArray Delphes/allParticles
     792  set OutputArray filteredParticles
     793  set PTMin 0.0
     794
     795}
    755796
    756797##################
    … …  
    761802# add Branch InputArray BranchName BranchClass
    762803  add Branch Delphes/allParticles Particle GenParticle
     804  add Branch GenParticleFilter/filteredParticles FilteredParticle GenParticle
    763805
    764806  add Branch GenJetFinder/jets GenJet Jet

  • cards/delphes_card_ATLAS.tcl

    r837fa70 ra7c9002  
    2121  Calorimeter
    2222  EFlowMerger
    23 
     23  EFlowFilter
     24 
    2425  PhotonEfficiency
    2526  PhotonIsolation
    … …  
    405406}
    406407
     408######################
     409# EFlowFilter
     410######################
     411
     412module 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
    407422###################
    408423# Photon efficiency
    … …  
    428443module Isolation PhotonIsolation {
    429444  set CandidateInputArray PhotonEfficiency/photons
    430   set IsolationInputArray EFlowMerger/eflow
     445  set IsolationInputArray EFlowFilter/eflow
    431446
    432447  set OutputArray photons
    … …  
    439454}
    440455
    441 #################
    442 # Electron filter
    443 #################
    444 
    445 module PdgCodeFilter ElectronFilter {
    446   set InputArray HCal/eflowTracks
    447   set OutputArray electrons
    448   set Invert true
    449   add PdgCode {11}
    450   add PdgCode {-11}
    451 }
    452456
    453457#####################
    … …  
    474478module Isolation ElectronIsolation {
    475479  set CandidateInputArray ElectronEfficiency/electrons
    476   set IsolationInputArray EFlowMerger/eflow
     480  set IsolationInputArray EFlowFilter/eflow
    477481
    478482  set OutputArray electrons
    … …  
    508512module Isolation MuonIsolation {
    509513  set CandidateInputArray MuonEfficiency/muons
    510   set IsolationInputArray EFlowMerger/eflow
     514  set IsolationInputArray EFlowFilter/eflow
    511515
    512516  set OutputArray muons

  • cards/delphes_card_ATLAS_PileUp.tcl

    r837fa70 ra7c9002  
    2323  EFlowMergerAllTracks
    2424  EFlowMerger
    25 
     25  EFlowFilter
     26 
    2627  NeutrinoFilter
    2728  GenJetFinder
    … …  
    372373}
    373374
     375######################
     376# EFlowFilter
     377######################
     378
     379module 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
    374390#############
    375391# Rho pile-up
    … …  
    508524module Isolation PhotonIsolation {
    509525  set CandidateInputArray PhotonEfficiency/photons
    510   set IsolationInputArray EFlowMergerAllTracks/eflow
     526  set IsolationInputArray EFlowFilter/eflow
    511527  set RhoInputArray Rho/rho
    512528
    … …  
    543559module Isolation ElectronIsolation {
    544560  set CandidateInputArray ElectronEfficiency/electrons
    545   set IsolationInputArray EFlowMergerAllTracks/eflow
     561  set IsolationInputArray EFlowFilter/eflow
    546562  set RhoInputArray Rho/rho
    547563
    … …  
    578594module Isolation MuonIsolation {
    579595  set CandidateInputArray MuonEfficiency/muons
    580   set IsolationInputArray EFlowMergerAllTracks/eflow
     596  set IsolationInputArray EFlowFilter/eflow
    581597  set RhoInputArray Rho/rho
    582598

  • cards/delphes_card_CMS.tcl

    r837fa70 ra7c9002  
    2121  Calorimeter
    2222  EFlowMerger
    23 
     23  EFlowFilter
     24 
    2425  PhotonEfficiency
    2526  PhotonIsolation
    … …  
    317318
    318319  set EnergyMin 1.0
    319   set EnergySignificanceMin 2.0
     320  set EnergySignificanceMin 1.0
    320321
    321322  set SmearTowerCenter true
    … …  
    418419}
    419420
     421######################
     422# EFlowFilter
     423######################
     424
     425module 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
    420436###################
    421437# Photon efficiency
    … …  
    441457module Isolation PhotonIsolation {
    442458  set CandidateInputArray PhotonEfficiency/photons
    443   set IsolationInputArray EFlowMerger/eflow
     459  set IsolationInputArray EFlowFilter/eflow
    444460
    445461  set OutputArray photons
    … …  
    476492module Isolation ElectronIsolation {
    477493  set CandidateInputArray ElectronEfficiency/electrons
    478   set IsolationInputArray EFlowMerger/eflow
     494  set IsolationInputArray EFlowFilter/eflow
    479495
    480496  set OutputArray electrons
    … …  
    510526module Isolation MuonIsolation {
    511527  set CandidateInputArray MuonEfficiency/muons
    512   set IsolationInputArray EFlowMerger/eflow
     528  set IsolationInputArray EFlowFilter/eflow
    513529
    514530  set OutputArray muons

  • cards/delphes_card_CMS_PileUp.tcl

    r837fa70 ra7c9002  
    2626  EFlowMergerAllTracks
    2727  EFlowMerger
    28 
     28  EFlowFilter
     29 
    2930  NeutrinoFilter
    3031  GenJetFinder
    … …  
    436437}
    437438
     439######################
     440# EFlowFilter
     441######################
     442
     443module 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}
    438452
    439453##########################
    … …  
    653667module Isolation PhotonIsolation {
    654668  set CandidateInputArray PhotonEfficiency/photons
    655   set IsolationInputArray EFlowMergerAllTracks/eflow
     669  set IsolationInputArray EFlowFilter/eflow
    656670  set RhoInputArray Rho/rho
    657671
    … …  
    688702module Isolation ElectronIsolation {
    689703  set CandidateInputArray ElectronEfficiency/electrons
    690   set IsolationInputArray EFlowMergerAllTracks/eflow
     704  set IsolationInputArray EFlowFilter/eflow
    691705  set RhoInputArray Rho/rho
    692706
    … …  
    724738module Isolation MuonIsolation {
    725739  set CandidateInputArray MuonEfficiency/muons
    726   set IsolationInputArray EFlowMergerAllTracks/eflow
     740  set IsolationInputArray EFlowFilter/eflow
    727741  set RhoInputArray Rho/rho
    728742

  • cards/delphes_card_ILD.tcl

    r837fa70 ra7c9002  
    2323  Calorimeter
    2424  EFlowMerger
    25 
     25  EFlowFilter
     26 
    2627  PhotonEfficiency
    2728  PhotonIsolation
    … …  
    350351}
    351352
     353######################
     354# EFlowFilter
     355######################
     356
     357module 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
    352367
    353368###################
    … …  
    491506module Isolation PhotonIsolation {
    492507  set CandidateInputArray PhotonEfficiency/photons
    493   set IsolationInputArray EFlowMerger/eflow
     508  set IsolationInputArray EFlowFilter/eflow
    494509
    495510  set OutputArray photons
    … …  
    525540module Isolation ElectronIsolation {
    526541  set CandidateInputArray ElectronEfficiency/electrons
    527   set IsolationInputArray EFlowMerger/eflow
     542  set IsolationInputArray EFlowFilter/eflow
    528543
    529544  set OutputArray electrons
    … …  
    561576module Isolation MuonIsolation {
    562577  set CandidateInputArray MuonEfficiency/muons
    563   set IsolationInputArray EFlowMerger/eflow
     578  set IsolationInputArray EFlowFilter/eflow
    564579
    565580  set OutputArray muons

  • classes/DelphesHepMCReader.cc

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

  • classes/DelphesLHEFReader.cc

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

  • examples/Example1.py

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

  • examples/Validation.cpp

    r837fa70 ra7c9002  
    250250      eta = TMath::Abs(genMomentum.Eta());
    251251
     252
    252253      if(eta > etamax || eta < etamin ) continue;
    253254
     255      //cout<<"b parton: "<<pt<<endl;
    254256      if (particle->PID == pdgID && genMomentum.Pt() > ptmin && genMomentum.Pt() < ptmax )
    255257      {
    … …  
    279281            if( !(jet->BTag & (1 << 0)) ) continue;
    280282          }
    281        
    282           if(TMath::Abs(pdgID) == 1)
     283
     284            if(TMath::Abs(pdgID) == 1)
    283285          {
    284286            Jet *jet = (Jet *)recoObj;
    … …  
    299301          }
    300302        }
    301 
    302303        histGenPt->Fill(pt);
    303         if(deltaR < 0.3) { histRecoPt->Fill(pt); }
     304        if(deltaR < 0.3 && bestRecoMomentum.Pt()> 0.20*pt ) { histRecoPt->Fill(pt);}
    304305
    305306      }
    … …  
    13521353  TClonesArray *branchCaloJet = treeReaderJet->UseBranch("CaloJet");
    13531354  TClonesArray *branchJet = treeReaderJet->UseBranch("Jet");
    1354  
     1355
    13551356  TClonesArray *branchParticleBJet = treeReaderBJet->UseBranch("Particle");
    13561357  TClonesArray *branchPFBJet = treeReaderBJet->UseBranch("Jet");
    … …  
    14041405
    14051406  TPaveText *pave = new TPaveText(0.0, 0.89, 0.94, 0.94,"NDC");
    1406   pave->SetTextAlign(kHAlignRight);
     1407  pave->SetTextAlign(30);
    14071408  pave->SetTextFont(132);
    14081409  pave->SetBorderSize(0);
    … …  
    14391440  const int n_ptbins = ptVals.size();
    14401441
    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
    19081911  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];
    19091913
    19101914  TCanvas *c_pfele_res_e[n_etabins];
    … …  
    19221926     gr_pfele_res_e[k] = EresGraph(&plots_pfele_res_e[k]);
    19231927
     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
    19241932     s_etaMin = Form("%.1f",etaVals.at(k));
    19251933     s_etaMax = Form("%.1f",etaVals.at(k+1));
    … …  
    19301938
    19311939     addResoGraph(mg_pfele_res_e[k], &gr_ecal_res_e[k], leg_pfele_res_e[k], markerStyles.at(0), colors.at(0), "ECAL");
    1932      addResoGraph(mg_pfele_res_e[k], &gr_trkele_res_pt[k], leg_pfele_res_e[k], markerStyles.at(1), colors.at(1), "Track");
     1940     addResoGraph(mg_pfele_res_e[k], &gr_trkele_res_e[k], leg_pfele_res_e[k], markerStyles.at(1), colors.at(1), "Track");
    19331941     addResoGraph(mg_pfele_res_e[k], &gr_pfele_res_e[k], leg_pfele_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    19341942
    … …  
    19611969     gr_pfele_res_eta[k] = EresGraphVsEta(&plots_pfele_res_eta[k]);
    19621970
     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
    19631975     s_e = Form("e^{ #pm}, E = %.0f GeV",ptVals.at(k));
    19641976     if(ptVals.at(k) >= 1000.) s_e = Form("e^{ #pm}, E = %.0f TeV",ptVals.at(k)/1000.);
    … …  
    19691981
    19701982     addResoGraph(mg_pfele_res_eta[k], &gr_ecal_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(0), colors.at(0), "ECAL");
    1971      addResoGraph(mg_pfele_res_eta[k], &gr_trkele_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
     1983     addResoGraph(mg_pfele_res_eta[k], &gr_trkele_res_eeta[k], leg_pfele_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
    19721984     addResoGraph(mg_pfele_res_eta[k], &gr_pfele_res_eta[k], leg_pfele_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    19731985
    … …  
    20012013  TGraphErrors gr_pfpi_res_eta[n_ptbins];
    20022014
     2015  TGraphErrors gr_trkpi_res_e[n_etabins];
     2016  TGraphErrors gr_trkpi_res_eeta[n_ptbins];
     2017
    20032018  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];
    20042020
    20052021  TCanvas *c_pfpi_res_e[n_etabins];
    … …  
    20172033     gr_pfpi_res_e[k] = EresGraph(&plots_pfpi_res_e[k]);
    20182034
     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
    20192040     s_etaMin = Form("%.1f",etaVals.at(k));
    20202041     s_etaMax = Form("%.1f",etaVals.at(k+1));
    … …  
    20252046
    20262047     addResoGraph(mg_pfpi_res_e[k], &gr_hcal_res_e[k], leg_pfpi_res_e[k], markerStyles.at(0), colors.at(0), "HCAL");
    2027      addResoGraph(mg_pfpi_res_e[k], &gr_trkpi_res_pt[k], leg_pfpi_res_e[k], markerStyles.at(1), colors.at(1), "Track");
     2048     addResoGraph(mg_pfpi_res_e[k], &gr_trkpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(1), colors.at(1), "Track");
    20282049     addResoGraph(mg_pfpi_res_e[k], &gr_pfpi_res_e[k], leg_pfpi_res_e[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    20292050
    … …  
    20562077     gr_pfpi_res_eta[k] = EresGraphVsEta(&plots_pfpi_res_eta[k]);
    20572078
     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
    20582084     s_e = Form("#pi^{ #pm}, E = %.0f GeV",ptVals.at(k));
    20592085     if(ptVals.at(k) >= 1000.) s_e = Form("#pi^{ #pm}, E = %.0f TeV",ptVals.at(k)/1000.);
    … …  
    20632089
    20642090     addResoGraph(mg_pfpi_res_eta[k], &gr_hcal_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(0), colors.at(0), "HCAL");
    2065      addResoGraph(mg_pfpi_res_eta[k], &gr_trkpi_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
     2091     addResoGraph(mg_pfpi_res_eta[k], &gr_trkpi_res_eeta[k], leg_pfpi_res_eta[k], markerStyles.at(1), colors.at(1), "Track");
    20662092     addResoGraph(mg_pfpi_res_eta[k], &gr_pfpi_res_eta[k], leg_pfpi_res_eta[k], markerStyles.at(2), colors.at(2), "Particle-flow");
    20672093
    … …  
    24182444    c_recpho_eff_eta->Print(figPath+"img_recpho_eff_eta.png","png");
    24192445
     2446
    24202447    /////////////////////////////////////////
    24212448    // B-jets  Efficiency/ mistag rates   ///
    … …  
    24822509    c_recbjet_eff_eta->Print(figPath+"img_recbjet_eff_eta.png","png");
    24832510
    2484 
    2485 
    24862511    // ------ c - mistag  ------
    24872512
    … …  
    25602585    for (k = 0; k < etaVals.size()-1; k++)
    25612586    {
    2562        
     2587
    25632588       h_recbjet_lmis_pt = GetEffPt<Jet>(branchJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
    25642589       gr_recbjet_lmis_pt[k] = TGraphErrors(h_recbjet_lmis_pt);
    … …  
    25892614
    25902615    mg_recbjet_lmis_pt->Draw("APE");
    2591     DrawAxis(mg_recbjet_lmis_pt, leg_recbjet_lmis_pt, ptMin, ptMax, 0.0, 0.005, "p_{T} [GeV]", "light - mistag rate (%)", true, false);
     2616    DrawAxis(mg_recbjet_lmis_pt, leg_recbjet_lmis_pt, ptMin, ptMax, 0.0, 0.5, "p_{T} [GeV]", "light - mistag rate (%)", true, false);
    25922617    leg_recbjet_lmis_pt->Draw();
    25932618    pave->Draw();
    … …  
    26002625
    26012626    mg_recbjet_lmis_eta->Draw("APE");
    2602     DrawAxis(mg_recbjet_lmis_eta, leg_recbjet_lmis_eta, etaMin, etaMax, 0.0, 0.005, " #eta ", "light - mistag rate (%)", false, false);
     2627    DrawAxis(mg_recbjet_lmis_eta, leg_recbjet_lmis_eta, etaMin, etaMax, 0.0, 0.5, " #eta ", "light - mistag rate (%)", false, false);
    26032628    leg_recbjet_lmis_eta->Draw();
    26042629    pave->Draw();
    … …  
    27182743
    27192744    mg_rectaujet_mis_pt->Draw("APE");
    2720     DrawAxis(mg_rectaujet_mis_pt, leg_rectaujet_mis_pt, ptMin, ptMax, 0.0, 10., "p_{T} [GeV]", "#tau - mistag(%)", true, false);
     2745    DrawAxis(mg_rectaujet_mis_pt, leg_rectaujet_mis_pt, ptMin, ptMax, 0.0, 5., "p_{T} [GeV]", "#tau - mistag(%)", true, false);
    27212746    leg_rectaujet_mis_pt->Draw();
    27222747    pave->Draw();
    … …  
    27502775       plots_trkpi_res_pt[k].at(bin).resolHist->Write();
    27512776       plots_trkele_res_pt[k].at(bin).resolHist->Write();
    2752            plots_trkmu_res_pt[k].at(bin).resolHist->Write();
     2777       plots_trkmu_res_pt[k].at(bin).resolHist->Write();
    27532778       plots_ecal_res_e[k].at(bin).resolHist->Write();
    27542779       plots_hcal_res_e[k].at(bin).resolHist->Write();
    … …  
    27802805
    27812806  fout->Write();
     2807
     2808
     2809
    27822810
    27832811  cout << "** Exiting..." << endl;

  • external/fastjet/CMakeLists.txt

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

  • modules/CMakeLists.txt

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

  • modules/StatusPidFilter.cc

    r837fa70 ra7c9002  
    104104    pass = kFALSE;
    105105
    106     // status == 3
     106    // hard scattering particles (first condition for Py6, second for Py8)
    107107    if(status == 3) pass = kTRUE;
     108    if(status > 20 && status < 30 ) pass = kTRUE;
    108109
    109110    // electrons, muons, taus and neutrinos
    … …  
    111112
    112113    // heavy quarks
    113     if(pdgCode == 5 || pdgCode == 6) pass = kTRUE;
     114    if(pdgCode == 4 ||pdgCode == 5 || pdgCode == 6) pass = kTRUE;
    114115
    115116    // Gauge bosons and other fundamental bosons
    116117    if(pdgCode > 22 && pdgCode < 43) pass = kTRUE;
    117118
    118     if(!pass || candidate->Momentum.Pt() <= fPTMin) continue;
     119    if(!pass || candidate->Momentum.Pt() < fPTMin) continue;
    119120
    120121    fOutputArray->Add(candidate);
Note: See TracChangeset for help on using the changeset viewer.