Fork me on GitHub

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

cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_200PU.tcl

    r308ab1a r8f4a953  
     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
    … …  
    114132
    115133  # pre-generated minbias input file
    116   set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
    117   #set PileUpFile MinBias.pileup
     134  #set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
     135  set PileUpFile MinBias.pileup
    118136
    119137  # average expected pile up
    … …  
    126144  set TVertexSpread 800E-12
    127145
    128   # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
     146  # 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))}
    129147  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    130148
    … …  
    138156module ParticlePropagator ParticlePropagator {
    139157  set InputArray PileUpMerger/stableParticles
     158  #set InputArray Delphes/stableParticles
    140159
    141160  set OutputArray stableParticles
    … …  
    272291  source muonMomentumResolution.tcl
    273292}
    274 
    275 
    276 
    277293
    278294##############
    … …  
    526542}
    527543
     544########################
     545# Reco PU filter
     546########################
     547
     548module RecoPuFilter RecoPuFilter {
     549  set InputArray HCal/eflowTracks
     550  set OutputArray eflowTracks
     551}
    528552
    529553###################################################
    … …  
    538562}
    539563
    540 
    541564####################
    542565# Neutral eflow erger
    … …  
    550573}
    551574
    552 
    553 ####################
     575#####################
    554576# Energy flow merger
    555 ####################
     577#####################
    556578
    557579module Merger EFlowMerger {
    … …  
    563585}
    564586
    565 ##################################
    566 # Energy flow merger (all tracks)
    567 ##################################
    568 
    569 module Merger EFlowMergerAllTracks {
     587############################
     588# Energy flow merger no PU
     589############################
     590
     591module Merger EFlowMergerCHS {
    570592# add InputArray InputArray
    571   add InputArray TrackMerger/tracks
     593  add InputArray RecoPuFilter/eflowTracks
    572594  add InputArray PhotonEnergySmearing/eflowPhotons
    573595  add InputArray HCal/eflowNeutralHadrons
    … …  
    694716}
    695717
    696 
    697718#####################
    698719# MC truth jet finder
    … …  
    735756
    736757
    737 
    738 ############
     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##############
    739788# Jet finder
    740 ############
     789##############
    741790
    742791module FastJetFinder FastJetFinder {
    743792#  set InputArray TowerMerger/towers
    744   set InputArray EFlowMerger/eflow
     793  set InputArray EFlowMergerCHS/eflow
    745794
    746795  set OutputArray jets
     796
     797  set AreaAlgorithm 5
    747798
    748799  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    756807module FastJetFinder FastJetFinderAK8 {
    757808#  set InputArray TowerMerger/towers
    758   set InputArray EFlowMerger/eflow
     809  set InputArray EFlowMergerCHS/eflow
    759810
    760811  set OutputArray jets
     812
     813  set AreaAlgorithm 5
    761814
    762815  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    785838}
    786839
     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
    787866module FastJetFinder FastJetFinderPUPPI {
    788867#  set InputArray TowerMerger/towers
    … …  
    834913
    835914module EnergyScale JetEnergyScale {
    836   set InputArray FastJetFinder/jets
     915  set InputArray JetPileUpSubtractor/jets
    837916  set OutputArray jets
    838917
    … …  
    842921
    843922module EnergyScale JetEnergyScaleAK8 {
    844   set InputArray FastJetFinderAK8/jets
     923  set InputArray JetPileUpSubtractorAK8/jets
    845924  set OutputArray jets
    846925
    … …  
    908987}
    909988
     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}
    9101015
    9111016
    … …  
    9301035
    9311036
     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
    9321056######################
    9331057# Electron isolation #
    … …  
    9501074}
    9511075
     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}
    9521095
    9531096
    … …  
    9961139}
    9971140
     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
    9981186##################
    9991187# Muon isolation #
    … …  
    10141202}
    10151203
    1016 
    1017 ##################
    1018 # Muon Loose Id  #
    1019 ##################
     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#####################
    10201226
    10211227module Efficiency MuonLooseIdEfficiency {
    … …  
    10391245}
    10401246
     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}
    10411271
    10421272
    … …  
    12511481module StatusPidFilter GenParticleFilter {
    12521482
    1253     set InputArray  Delphes/allParticles
     1483    set InputArray Delphes/allParticles
    12541484    set OutputArray filteredParticles
    12551485    set PTMin 5.0
    … …  
    12791509  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    12801510  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
    12811516
    12821517  add Branch JetEnergyScale/jets Jet Jet
    … …  
    12901525  add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
    12911526  add Branch ScalarHT/energy ScalarHT ScalarHT
    1292 }
     1527
     1528}
Note: See TracChangeset for help on using the changeset viewer.