Fork me on GitHub

Changeset 740d430 in git


Ignore:
Timestamp:
Jun 6, 2016, 4:43:16 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
479f808, 4a0d9d5
Parents:
e2339af
Message:

fixed photon isolation with puppi

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl

    re2339af r740d430  
     1set MaxEvents 10
     2
    13#
    24#  Phase II - Pile-Up
    35#
    46#  Main authors: Michele Selvaggi (UCL)
    5 #                               
    6 #  Released on: 
     7#
     8#  Released on:
    79#
    810#  Version: v01
    … …  
    1416
    1517set ExecutionPath {
    16  
    17   PileUpMerger
    18  
     18
     19  PileUpMerger
    1920  ParticlePropagator
    2021
    … …  
    2829
    2930  TrackMerger
    30  
     31
    3132  ECal
    3233  HCal
    … …  
    3435  ElectronFilter
    3536  TrackPileUpSubtractor
    36  
     37
    3738  TowerMerger
    3839  NeutralEFlowMerger
    … …  
    4546  RunPUPPI
    4647
     48  PhotonFilter
    4749
    4850  PhotonIsolation
    4951  PhotonEfficiency
    50  
     52
    5153  ElectronIsolation
    5254  ElectronEfficiency
    53  
     55
    5456  MuonIsolation
    5557  MuonLooseIdEfficiency
    5658  MuonTightIdEfficiency
    57  
     59
    5860  NeutrinoFilter
    5961
    … …  
    6769
    6870  ScalarHT
    69  
     71
    7072  JetEnergyScale
    7173
    … …  
    7779
    7880  TauTagging
    79  
     81
    8082  GenParticleFilter
    81  
     83
    8284  TreeWriter
    8385}
    … …  
    9698
    9799  # pre-generated minbias input file
    98   set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup 
     100  set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
    99101
    100102  # average expected pile up
    101103  set MeanPileUp 200
    102  
     104
    103105  # maximum spread in the beam direction in m
    104106  set ZVertexSpread 0.25
    … …  
    109111  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
    110112  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    111  
     113
    112114}
    113115
    … …  
    165167  set OutputArray electrons
    166168  # tracking efficiency formula for electrons
    167   set EfficiencyFormula { 
     169  set EfficiencyFormula {
    168170      (pt <= 0.2) * (0.00) + \
    169171          (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
    … …  
    197199          (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
    198200          (abs(eta) > 4.0) * (0.00)
    199          
     201
    200202  }
    201203}
    … …  
    210212  set InputArray  ChargedHadronTrackingEfficiency/chargedHadrons
    211213  set OutputArray chargedHadrons
    212   # resolution formula for charged hadrons , 
    213  
     214  # resolution formula for charged hadrons ,
     215
    214216  # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
    215217  source trackMomentumResolution.tcl
    … …  
    227229
    228230  # resolution formula for electrons
    229  
     231
    230232  # taking something flat in energy for now, ECAL will take over at high energy anyway.
    231233  # inferred from hep-ex/1306.2016 and 1502.02701
    232   set ResolutionFormula { 
    233  
     234  set ResolutionFormula {
     235
    234236                        (abs(eta) <= 1.5)  * (1+0.64*abs(eta)^2)*(energy*0.028) +
    235237    (abs(eta) > 1.5  && abs(eta) <= 1.75)  * (energy*0.037) +
    236238    (abs(eta) > 1.75  && abs(eta) <= 2.15) * (energy*0.038) +
    237239    (abs(eta) > 2.15  && abs(eta) <= 3.00) * (energy*0.044) +
    238     (abs(eta) > 3.00  && abs(eta) <= 4.00) * (energy*0.10)}   
     240    (abs(eta) > 3.00  && abs(eta) <= 4.00) * (energy*0.10)}
    239241
    240242}
    … …  
    248250  set OutputArray muons
    249251  # resolution formula for muons
    250  
    251   # up to |eta| < 2.8 take measurement from tracking + muon chambers 
    252   # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from 
     252
     253  # up to |eta| < 2.8 take measurement from tracking + muon chambers
     254  # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
    253255  # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
    254   source muonMomentumResolution.tcl 
     256  source muonMomentumResolution.tcl
    255257}
    256258
    … …  
    282284  set EFlowTowerOutputArray eflowPhotons
    283285
    284   set IsEcal true 
    285  
     286  set IsEcal true
     287
    286288  set EnergyMin 0.5
    287289  set EnergySignificanceMin 1.0
    … …  
    296298
    297299  # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
    298  
     300
    299301  set PhiBins {}
    300302  for {set i -180} {$i <= 180} {incr i} {
    … …  
    309311
    310312  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
    311  
     313
    312314  set PhiBins {}
    313315  for {set i -180} {$i <= 180} {incr i} {
    … …  
    326328  }
    327329
    328   # take present CMS granularity for HF 
    329  
     330  # take present CMS granularity for HF
     331
    330332  # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
    331333  set PhiBins {}
    … …  
    333335    add PhiBins [expr {$i * $pi/18.0}]
    334336  }
    335  
     337
    336338  foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
    337339    add EtaPhiBins $eta $PhiBins
    … …  
    359361
    360362  # set ResolutionFormula {resolution formula as a function of eta and energy}
    361  
     363
    362364  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
    363365  # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL  see LHCC-P-008, Fig. 3.39, p.117
    364366
    365   set ResolutionFormula {  (abs(eta) <= 1.50)                    * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) + 
     367  set ResolutionFormula {  (abs(eta) <= 1.50)                    * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
    366368                           (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
    367369                           (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
    … …  
    383385  set EFlowTowerOutputArray eflowNeutralHadrons
    384386
    385   set IsEcal false 
    386  
     387  set IsEcal false
     388
    387389  set EnergyMin 1.0
    388390  set EnergySignificanceMin 1.0
    … …  
    395397  # each list starts with the lower edge of the first tower
    396398  # the list ends with the higher edged of the last tower
    397  
     399
    398400  # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
    399401
    … …  
    407409
    408410  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
    409  
     411
    410412  set PhiBins {}
    411413  for {set i -180} {$i <= 180} {incr i} {
    … …  
    424426  }
    425427
    426   # take present CMS granularity for HF 
    427  
     428  # take present CMS granularity for HF
     429
    428430  # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
    429431  set PhiBins {}
    … …  
    431433    add PhiBins [expr {$i * $pi/18.0}]
    432434  }
    433  
     435
    434436  foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
    435437    add EtaPhiBins $eta $PhiBins
    … …  
    476478}
    477479
     480
    478481##########################
    479482# Track pile-up subtractor
    … …  
    573576  set UseExp            false
    574577  set UseNoLep          false
    575  
    576   ## define puppi algorithm parameters (more than one for the same eta region is possible)                                                                                     
     578
     579  ## define puppi algorithm parameters (more than one for the same eta region is possible)
    577580  add EtaMinBin           0.0   1.5   4.0
    578581  add EtaMaxBin           1.5   4.0   10.0
    … …  
    585588  add ApplyCHS            true  true  true
    586589  add UseCharged          true  true  false
    587   add ApplyLowPUCorr      true  true  true 
     590  add ApplyLowPUCorr      true  true  true
    588591  add MetricId            5     5     5
    589592  add CombId              0     0     0
    … …  
    593596  set OutputArrayTracks   puppiTracks
    594597  set OutputArrayNeutrals puppiNeutrals
    595 } 
     598}
    596599
    597600module Merger RunPUPPI {
    … …  
    620623
    621624###################
    622 # Ger PileUp Missing ET 
     625# Ger PileUp Missing ET
    623626###################
    624627
    … …  
    667670module FastJetFinder GenJetFinder {
    668671  set InputArray NeutrinoFilter/filteredParticles
    669  
     672
    670673  set OutputArray jets
    671674
    … …  
    720723
    721724
     725#################
     726# Photon filter
     727#################
     728
     729module PdgCodeFilter PhotonFilter {
     730  set InputArray ECal/eflowPhotons
     731  set OutputArray photons
     732  set Invert true
     733  set PTMin 5.0
     734  add PdgCode {22}
     735}
     736
     737
    722738####################
    723739# Photon isolation #
    … …  
    725741
    726742module Isolation PhotonIsolation {
    727  
     743
    728744  # particle for which calculate the isolation
    729   set CandidateInputArray ECal/eflowPhotons
    730  
     745  set CandidateInputArray PhotonFilter/photons
     746
    731747  # isolation collection
    732   set IsolationInputArray EFlowMerger/eflow
    733  
     748  set IsolationInputArray RunPUPPI/PuppiParticles
     749
    734750  # output array
    735751  set OutputArray photons
    736  
     752
     753  # veto isolation cand. based on proximity to input cand.
     754  set DeltaRMin 0.01
     755  set UseMiniCone true
     756
    737757  # isolation cone
    738758  set DeltaRMax 0.3
    739  
    740   # minimum pT 
     759
     760  # minimum pT
    741761  set PTMin     1.0
    742  
     762
    743763  # iso ratio to cut
    744764  set PTRatioMax 9999.
     765
    745766}
    746767
    … …  
    752773
    753774module Efficiency PhotonEfficiency {
    754  
     775
    755776  ## input particles
    756   set InputArray PhotonIsolation/photons 
     777  set InputArray PhotonIsolation/photons
    757778  ## output particles
    758779  set OutputArray photons
    … …  
    772793
    773794module Isolation ElectronIsolation {
    774  
     795
    775796  set CandidateInputArray ElectronFilter/electrons
    776  
     797
    777798  # isolation collection
    778   set IsolationInputArray EFlowMerger/eflow
    779  
     799  set IsolationInputArray RunPUPPI/PuppiParticles
     800  #set IsolationInputArray EFlowMerger/eflow
     801
    780802  set OutputArray electrons
    781  
     803
    782804  set DeltaRMax 0.3
    783805  set PTMin 1.0
    … …  
    793815
    794816module Efficiency ElectronEfficiency {
    795  
     817
    796818  set InputArray ElectronIsolation/electrons
    797819  set OutputArray electrons
    798  
     820
    799821  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    800822  # efficiency formula for electrons
    … …  
    804826                         (abs(eta) <= 1.45 ) * (pt >  6.0 && pt <= 8.0)   * (0.70) + \
    805827                         (abs(eta) <= 1.45 ) * (pt >  8.0 && pt <= 10.0)  * (0.85) + \
    806                          (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0)  * (0.94) + \                                                     
    807                          (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0)  * (0.97) + \                         
    808                          (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0)  * (0.98) + \         
    809                          (abs(eta) <= 1.45 ) * (pt > 70.0 )  * (1.0) + \                                                                                                 
     828                         (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0)  * (0.94) + \
     829                         (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0)  * (0.97) + \
     830                         (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0)  * (0.98) + \
     831                         (abs(eta) <= 1.45 ) * (pt > 70.0 )  * (1.0) + \
    810832                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt >  4.0 && pt <= 10.0)   * (0.35) + \
    811                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0)   * (0.40) + \   
    812                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0)   * (0.45) + \                                 
    813                          (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 70.0 )  * (0.55) + \   
     833                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0)   * (0.40) + \
     834                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0)   * (0.45) + \
     835                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 70.0 )  * (0.55) + \
    814836                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt >  4.0 && pt <= 10.0)  * (0.75) + \
    815                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0)  * (0.85) + \                                                     
    816                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0)  * (0.95) + \                         
    817                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0)  * (0.95) + \         
    818                          (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 )  * (1.0) + \   
     837                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0)  * (0.85) + \
     838                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0)  * (0.95) + \
     839                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0)  * (0.95) + \
     840                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 )  * (1.0) + \
    819841                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt >  4.0 && pt <= 10.0)  * (0.65) + \
    820                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0)  * (0.75) + \                                                     
    821                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0)  * (0.90) + \                         
    822                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0)  * (0.90) + \         
    823                          (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 )  * (0.90) + \                                                                           
     842                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0)  * (0.75) + \
     843                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0)  * (0.90) + \
     844                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0)  * (0.90) + \
     845                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 )  * (0.90) + \
    824846                         (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
    825847                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
    … …  
    828850                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
    829851                                          (abs(eta) > 4.0) * (0.00)
    830  
     852
    831853  }
    832854}
    … …  
    838860module Isolation MuonIsolation {
    839861  set CandidateInputArray MuonMomentumSmearing/muons
    840  
     862
    841863  # isolation collection
    842   set IsolationInputArray EFlowMerger/eflow
    843  
     864  set IsolationInputArray RunPUPPI/PuppiParticles
     865
    844866  set OutputArray muons
    845  
     867
    846868  set DeltaRMax 0.3
    847869  set PTMin 1.0
    … …  
    916938
    917939  set BitNumber 1
    918  
     940
    919941  source btagMedium.tcl
    920942}
    921943
    922  
    923944
    924945#############
    … …  
    931952  set BitNumber 2
    932953
    933   source btagTight.tcl
    934 }
    935 
    936 
    937 
     954  source btagTight.tcl
     955}
    938956
    939957
    … …  
    961979  add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.97*0.77*( (0.32 + 0.01*pt - 0.000054*pt*pt )*(pt<100)+0.78*(pt>100) ) + \
    962980                               (abs(eta) > 2.3) * (0.000)
    963                              } 
     981                             }
    964982}
    965983
    … …  
    970988
    971989module StatusPidFilter GenParticleFilter {
    972  
     990
    973991    set InputArray  Delphes/allParticles
    974992    set OutputArray filteredParticles
    975993    set PTMin 5.0
    976    
     994
    977995}
    978996
    … …  
    9861004  add Branch GenParticleFilter/filteredParticles Particle GenParticle
    9871005  add Branch PileUpMerger/vertices Vertex Vertex
    988  
     1006
    9891007  add Branch GenJetFinder/jets GenJet Jet
    9901008  add Branch GenMissingET/momentum GenMissingET MissingET
    … …  
    9931011#  add Branch ECal/eflowPhotons EFlowPhoton Tower
    9941012#  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
    995  
     1013
    9961014  add Branch PhotonEfficiency/photons Photon Photon
    9971015  add Branch ElectronEfficiency/electrons Electron Electron
    9981016  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    9991017  add Branch MuonTightIdEfficiency/muons MuonTight Muon
    1000    
     1018
    10011019  add Branch JetEnergyScale/jets Jet Jet
    1002  
     1020
    10031021  add Branch MissingET/momentum MissingET MissingET
    10041022  add Branch PuppiMissingET/momentum PuppiMissingET MissingET
Note: See TracChangeset for help on using the changeset viewer.