Fork me on GitHub

Changeset 8cbe720 in git for cards/CMS_PhaseII


Ignore:
Timestamp:
Aug 22, 2016, 7:27:26 PM (8 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
1270747
Parents:
8f4a953
Message:

added CHS+rho corrections

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl

    r8f4a953 r8cbe720  
     1set MaxEvents 100
     2
    13#
    24#  Phase II - No Pile-Up
    … …  
    3335  PhotonEnergySmearing
    3436  ElectronFilter
     37
    3538  TrackPileUpSubtractor
     39  RecoPuFilter
    3640
    3741  TowerMerger
    3842  NeutralEFlowMerger
    39   EFlowMergerAllTracks
     43
    4044  EFlowMerger
     45  EFlowMergerCHS
     46  Rho
    4147
    4248  LeptonFilterNoLep
    … …  
    4854
    4955  PhotonIsolation
     56  PhotonIsolationCHS
    5057  PhotonEfficiency
     58  PhotonEfficiencyCHS
    5159
    5260  ElectronIsolation
     61  ElectronIsolationCHS
     62
    5363  ElectronEfficiency
     64  ElectronEfficiencyCHS
    5465
    5566  MuonIsolation
     67  MuonIsolationCHS
     68
    5669  MuonLooseIdEfficiency
    5770  MuonTightIdEfficiency
     71
     72  MuonLooseIdEfficiencyCHS
     73  MuonTightIdEfficiencyCHS
    5874
    5975  NeutrinoFilter
    … …  
    6884  FastJetFinder
    6985  FastJetFinderAK8
     86  JetPileUpSubtractor
     87  JetPileUpSubtractorAK8
    7088  FastJetFinderPUPPI
    7189  FastJetFinderPUPPIAK8
    … …  
    121139
    122140  # maximum spread in the beam direction in m
    123   set ZVertexSpread 0
     141  set ZVertexSpread 0.25
    124142
    125143  # maximum spread in time in s
    126   set TVertexSpread 0
     144  set TVertexSpread 800E-12
    127145
    128146  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) - {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    129   set VertexDistributionFormula 1
     147  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    130148
    131149}
    … …  
    273291  source muonMomentumResolution.tcl
    274292}
    275 
    276 
    277 
    278293
    279294##############
    … …  
    527542}
    528543
     544########################
     545# Reco PU filter
     546########################
     547
     548module RecoPuFilter RecoPuFilter {
     549  set InputArray HCal/eflowTracks
     550  set OutputArray eflowTracks
     551}
    529552
    530553###################################################
    … …  
    539562}
    540563
    541 
    542564####################
    543565# Neutral eflow erger
    … …  
    551573}
    552574
    553 
    554 ####################
     575#####################
    555576# Energy flow merger
    556 ####################
     577#####################
    557578
    558579module Merger EFlowMerger {
    … …  
    564585}
    565586
    566 ##################################
    567 # Energy flow merger (all tracks)
    568 ##################################
    569 
    570 module Merger EFlowMergerAllTracks {
     587############################
     588# Energy flow merger no PU
     589############################
     590
     591module Merger EFlowMergerCHS {
    571592# add InputArray InputArray
    572   add InputArray TrackMerger/tracks
     593  add InputArray RecoPuFilter/eflowTracks
    573594  add InputArray PhotonEnergySmearing/eflowPhotons
    574595  add InputArray HCal/eflowNeutralHadrons
    … …  
    695716}
    696717
    697 
    698718#####################
    699719# MC truth jet finder
    … …  
    735755}
    736756
    737 ############
     757
     758#############
     759# Rho pile-up
     760#############
     761
     762module FastJetFinder Rho {
     763#  set InputArray Calorimeter/towers
     764  set InputArray EFlowMergerCHS/eflow
     765
     766  set ComputeRho true
     767  set RhoOutputArray rho
     768
     769  # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
     770  set AreaAlgorithm 5
     771
     772  # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     773  set JetAlgorithm 4
     774  set ParameterR 0.4
     775  set GhostEtaMax 5.0
     776
     777  add RhoEtaRange -5.0 -4.0
     778  add RhoEtaRange -4.0 -1.5
     779  add RhoEtaRange -1.5 1.5
     780  add RhoEtaRange 1.5 4.0
     781  add RhoEtaRange 4.0 5.0
     782
     783  set JetPTMin 0.0
     784}
     785
     786
     787##############
    738788# Jet finder
    739 ############
     789##############
    740790
    741791module FastJetFinder FastJetFinder {
    742792#  set InputArray TowerMerger/towers
    743   set InputArray EFlowMerger/eflow
     793  set InputArray EFlowMergerCHS/eflow
    744794
    745795  set OutputArray jets
     796
     797  set AreaAlgorithm 5
    746798
    747799  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    755807module FastJetFinder FastJetFinderAK8 {
    756808#  set InputArray TowerMerger/towers
    757   set InputArray EFlowMerger/eflow
     809  set InputArray EFlowMergerCHS/eflow
    758810
    759811  set OutputArray jets
     812
     813  set AreaAlgorithm 5
    760814
    761815  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    784838}
    785839
     840###########################
     841# Jet Pile-Up Subtraction
     842###########################
     843
     844module JetPileUpSubtractor JetPileUpSubtractor {
     845  set JetInputArray FastJetFinder/jets
     846  set RhoInputArray Rho/rho
     847
     848  set OutputArray jets
     849
     850  set JetPTMin 15.0
     851}
     852
     853##############################
     854# Jet Pile-Up Subtraction AK8
     855##############################
     856
     857module JetPileUpSubtractor JetPileUpSubtractorAK8 {
     858  set JetInputArray FastJetFinderAK8/jets
     859  set RhoInputArray Rho/rho
     860
     861  set OutputArray jets
     862
     863  set JetPTMin 15.0
     864}
     865
    786866module FastJetFinder FastJetFinderPUPPI {
    787867#  set InputArray TowerMerger/towers
    … …  
    833913
    834914module EnergyScale JetEnergyScale {
    835   set InputArray FastJetFinder/jets
     915  set InputArray JetPileUpSubtractor/jets
    836916  set OutputArray jets
    837917
    … …  
    841921
    842922module EnergyScale JetEnergyScaleAK8 {
    843   set InputArray FastJetFinderAK8/jets
     923  set InputArray JetPileUpSubtractorAK8/jets
    844924  set OutputArray jets
    845925
    … …  
    907987}
    908988
     989
     990########################
     991# Photon isolation CHS #
     992########################
     993
     994module Isolation PhotonIsolationCHS {
     995
     996  # particle for which calculate the isolation
     997  set CandidateInputArray PhotonFilter/photons
     998
     999  # isolation collection
     1000  set IsolationInputArray EFlowMerger/eflow
     1001
     1002  # output array
     1003  set OutputArray photons
     1004
     1005  # isolation cone
     1006  set DeltaRMax 0.3
     1007
     1008  # minimum pT
     1009  set PTMin     1.0
     1010
     1011  # iso ratio to cut
     1012  set PTRatioMax 9999.
     1013
     1014}
    9091015
    9101016
    … …  
    9291035
    9301036
     1037#####################
     1038# Photon efficiency #
     1039#####################
     1040
     1041module Efficiency PhotonEfficiencyCHS {
     1042
     1043  ## input particles
     1044  set InputArray PhotonIsolationCHS/photons
     1045  ## output particles
     1046  set OutputArray photons
     1047  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     1048  # efficiency formula for photons
     1049  set EfficiencyFormula {                      (pt <= 10.0) * (0.00) + \
     1050                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.9635) + \
     1051         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0)  * (0.9624) + \
     1052         (abs(eta) > 4.0)                                   * (0.00)}
     1053
     1054}
     1055
    9311056######################
    9321057# Electron isolation #
    … …  
    9491074}
    9501075
     1076
     1077##########################
     1078# Electron isolation CHS #
     1079##########################
     1080
     1081module Isolation ElectronIsolationCHS {
     1082
     1083  set CandidateInputArray ElectronFilter/electrons
     1084
     1085  # isolation collection
     1086  set IsolationInputArray EFlowMerger/eflow
     1087
     1088  set OutputArray electrons
     1089
     1090  set DeltaRMax 0.3
     1091  set PTMin 1.0
     1092  set PTRatioMax 9999.
     1093
     1094}
    9511095
    9521096
    … …  
    9951139}
    9961140
     1141###########################
     1142# Electron efficiency CHS #
     1143###########################
     1144
     1145module Efficiency ElectronEfficiencyCHS {
     1146
     1147  set InputArray ElectronIsolationCHS/electrons
     1148  set OutputArray electrons
     1149
     1150  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     1151  # efficiency formula for electrons
     1152  set EfficiencyFormula {
     1153                                      (pt <= 4.0)  * (0.00) + \
     1154                         (abs(eta) <= 1.45 ) * (pt >  4.0 && pt <= 6.0)   * (0.50) + \
     1155                         (abs(eta) <= 1.45 ) * (pt >  6.0 && pt <= 8.0)   * (0.70) + \
     1156                         (abs(eta) <= 1.45 ) * (pt >  8.0 && pt <= 10.0)  * (0.85) + \
     1157                         (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0)  * (0.94) + \
     1158                         (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0)  * (0.97) + \
     1159                         (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0)  * (0.98) + \
     1160                         (abs(eta) <= 1.45 ) * (pt > 70.0 )  * (1.0) + \
     1161                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt >  4.0 && pt <= 10.0)   * (0.35) + \
     1162                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0)   * (0.40) + \
     1163                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0)   * (0.45) + \
     1164                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 70.0 )  * (0.55) + \
     1165                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt >  4.0 && pt <= 10.0)  * (0.75) + \
     1166                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0)  * (0.85) + \
     1167                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0)  * (0.95) + \
     1168                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0)  * (0.95) + \
     1169                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 )  * (1.0) + \
     1170                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt >  4.0 && pt <= 10.0)  * (0.65) + \
     1171                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0)  * (0.75) + \
     1172                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0)  * (0.90) + \
     1173                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0)  * (0.90) + \
     1174                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 )  * (0.90) + \
     1175                         (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
     1176                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
     1177                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
     1178                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
     1179                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
     1180                                          (abs(eta) > 4.0) * (0.00)
     1181
     1182  }
     1183}
     1184
     1185
    9971186##################
    9981187# Muon isolation #
    … …  
    10131202}
    10141203
    1015 
    1016 ##################
    1017 # Muon Loose Id  #
    1018 ##################
     1204######################
     1205# Muon isolation CHS #
     1206######################
     1207
     1208module Isolation MuonIsolationCHS {
     1209  set CandidateInputArray MuonMomentumSmearing/muons
     1210
     1211  # isolation collection
     1212  set IsolationInputArray EFlowMerger/eflow
     1213
     1214  set OutputArray muons
     1215
     1216  set DeltaRMax 0.3
     1217  set PTMin 1.0
     1218  set PTRatioMax 9999.
     1219
     1220}
     1221
     1222
     1223#####################
     1224# Muon Loose Id     #
     1225#####################
    10191226
    10201227module Efficiency MuonLooseIdEfficiency {
    … …  
    10381245}
    10391246
     1247
     1248#####################
     1249# Muon Loose Id CHS #
     1250#####################
     1251
     1252module Efficiency MuonLooseIdEfficiencyCHS {
     1253    set InputArray MuonIsolationCHS/muons
     1254    set OutputArray muons
     1255    # tracking + LooseID efficiency formula for muons
     1256    source muonLooseId.tcl
     1257
     1258}
     1259
     1260
     1261######################
     1262# Muon Tight Id  CHS #
     1263######################
     1264
     1265module Efficiency MuonTightIdEfficiencyCHS {
     1266    set InputArray MuonIsolationCHS/muons
     1267    set OutputArray muons
     1268    # tracking + TightID efficiency formula for muons
     1269    source muonTightId.tcl
     1270}
    10401271
    10411272
    … …  
    12501481module StatusPidFilter GenParticleFilter {
    12511482
    1252     set InputArray  Delphes/allParticles
     1483    set InputArray Delphes/allParticles
    12531484    set OutputArray filteredParticles
    12541485    set PTMin 5.0
    … …  
    12781509  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    12791510  add Branch MuonTightIdEfficiency/muons MuonTight Muon
     1511
     1512  add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
     1513  add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
     1514  add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
     1515  add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon
    12801516
    12811517  add Branch JetEnergyScale/jets Jet Jet
    … …  
    12891525  add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
    12901526  add Branch ScalarHT/energy ScalarHT ScalarHT
    1291 }
     1527
     1528}
Note: See TracChangeset for help on using the changeset viewer.