Fork me on GitHub

Ignore:
Timestamp:
Aug 25, 2016, 2:04:58 PM (8 years ago)
Author:
Alexandre Mertens <alexandre.mertens@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
629e819
Parents:
7bb13cd (diff), 1408174 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.
Message:

dev_01

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl

    r7bb13cd r70b9632  
    3333  PhotonEnergySmearing
    3434  ElectronFilter
     35
    3536  TrackPileUpSubtractor
     37  RecoPuFilter
    3638
    3739  TowerMerger
    3840  NeutralEFlowMerger
    39   EFlowMergerAllTracks
     41
    4042  EFlowMerger
     43  EFlowMergerCHS
     44  Rho
    4145
    4246  LeptonFilterNoLep
    … …  
    4852
    4953  PhotonIsolation
     54  PhotonIsolationCHS
    5055  PhotonEfficiency
     56  PhotonEfficiencyCHS
    5157
    5258  ElectronIsolation
     59  ElectronIsolationCHS
     60
    5361  ElectronEfficiency
     62  ElectronEfficiencyCHS
    5463
    5564  MuonIsolation
     65  MuonIsolationCHS
     66
    5667  MuonLooseIdEfficiency
    5768  MuonTightIdEfficiency
     69
     70  MuonLooseIdEfficiencyCHS
     71  MuonTightIdEfficiencyCHS
    5872
    5973  NeutrinoFilter
    … …  
    6882  FastJetFinder
    6983  FastJetFinderAK8
     84  JetPileUpSubtractor
     85  JetPileUpSubtractorAK8
    7086  FastJetFinderPUPPI
    7187  FastJetFinderPUPPIAK8
    … …  
    121137
    122138  # maximum spread in the beam direction in m
    123   set ZVertexSpread 0
     139  set ZVertexSpread 0.25
    124140
    125141  # maximum spread in time in s
    126   set TVertexSpread 0
     142  set TVertexSpread 800E-12
    127143
    128144  # 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
     145  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    130146
    131147}
    … …  
    273289  source muonMomentumResolution.tcl
    274290}
    275 
    276 
    277 
    278291
    279292##############
    … …  
    527540}
    528541
     542########################
     543# Reco PU filter
     544########################
     545
     546module RecoPuFilter RecoPuFilter {
     547  set InputArray HCal/eflowTracks
     548  set OutputArray eflowTracks
     549}
    529550
    530551###################################################
    … …  
    539560}
    540561
    541 
    542562####################
    543563# Neutral eflow erger
    … …  
    551571}
    552572
    553 
    554 ####################
     573#####################
    555574# Energy flow merger
    556 ####################
     575#####################
    557576
    558577module Merger EFlowMerger {
    … …  
    564583}
    565584
    566 ##################################
    567 # Energy flow merger (all tracks)
    568 ##################################
    569 
    570 module Merger EFlowMergerAllTracks {
     585############################
     586# Energy flow merger no PU
     587############################
     588
     589module Merger EFlowMergerCHS {
    571590# add InputArray InputArray
    572   add InputArray TrackMerger/tracks
     591  add InputArray RecoPuFilter/eflowTracks
    573592  add InputArray PhotonEnergySmearing/eflowPhotons
    574593  add InputArray HCal/eflowNeutralHadrons
    … …  
    695714}
    696715
    697 
    698716#####################
    699717# MC truth jet finder
    … …  
    735753}
    736754
    737 ############
     755
     756#############
     757# Rho pile-up
     758#############
     759
     760module FastJetFinder Rho {
     761#  set InputArray Calorimeter/towers
     762  set InputArray EFlowMergerCHS/eflow
     763
     764  set ComputeRho true
     765  set RhoOutputArray rho
     766
     767  # 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
     768  set AreaAlgorithm 5
     769
     770  # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     771  set JetAlgorithm 4
     772  set ParameterR 0.4
     773  set GhostEtaMax 5.0
     774
     775  add RhoEtaRange -5.0 -4.0
     776  add RhoEtaRange -4.0 -1.5
     777  add RhoEtaRange -1.5 1.5
     778  add RhoEtaRange 1.5 4.0
     779  add RhoEtaRange 4.0 5.0
     780
     781  set JetPTMin 0.0
     782}
     783
     784
     785##############
    738786# Jet finder
    739 ############
     787##############
    740788
    741789module FastJetFinder FastJetFinder {
    742790#  set InputArray TowerMerger/towers
    743   set InputArray EFlowMerger/eflow
     791  set InputArray EFlowMergerCHS/eflow
    744792
    745793  set OutputArray jets
     794
     795  set AreaAlgorithm 5
    746796
    747797  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    755805module FastJetFinder FastJetFinderAK8 {
    756806#  set InputArray TowerMerger/towers
    757   set InputArray EFlowMerger/eflow
     807  set InputArray EFlowMergerCHS/eflow
    758808
    759809  set OutputArray jets
     810
     811  set AreaAlgorithm 5
    760812
    761813  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    … …  
    784836}
    785837
     838###########################
     839# Jet Pile-Up Subtraction
     840###########################
     841
     842module JetPileUpSubtractor JetPileUpSubtractor {
     843  set JetInputArray FastJetFinder/jets
     844  set RhoInputArray Rho/rho
     845
     846  set OutputArray jets
     847
     848  set JetPTMin 15.0
     849}
     850
     851##############################
     852# Jet Pile-Up Subtraction AK8
     853##############################
     854
     855module JetPileUpSubtractor JetPileUpSubtractorAK8 {
     856  set JetInputArray FastJetFinderAK8/jets
     857  set RhoInputArray Rho/rho
     858
     859  set OutputArray jets
     860
     861  set JetPTMin 15.0
     862}
     863
    786864module FastJetFinder FastJetFinderPUPPI {
    787865#  set InputArray TowerMerger/towers
    … …  
    833911
    834912module EnergyScale JetEnergyScale {
    835   set InputArray FastJetFinder/jets
     913  set InputArray JetPileUpSubtractor/jets
    836914  set OutputArray jets
    837915
    … …  
    841919
    842920module EnergyScale JetEnergyScaleAK8 {
    843   set InputArray FastJetFinderAK8/jets
     921  set InputArray JetPileUpSubtractorAK8/jets
    844922  set OutputArray jets
    845923
    … …  
    907985}
    908986
     987
     988########################
     989# Photon isolation CHS #
     990########################
     991
     992module Isolation PhotonIsolationCHS {
     993
     994  # particle for which calculate the isolation
     995  set CandidateInputArray PhotonFilter/photons
     996
     997  # isolation collection
     998  set IsolationInputArray EFlowMerger/eflow
     999
     1000  # output array
     1001  set OutputArray photons
     1002
     1003  # isolation cone
     1004  set DeltaRMax 0.3
     1005
     1006  # minimum pT
     1007  set PTMin     1.0
     1008
     1009  # iso ratio to cut
     1010  set PTRatioMax 9999.
     1011
     1012}
    9091013
    9101014
    … …  
    9291033
    9301034
     1035#####################
     1036# Photon efficiency #
     1037#####################
     1038
     1039module Efficiency PhotonEfficiencyCHS {
     1040
     1041  ## input particles
     1042  set InputArray PhotonIsolationCHS/photons
     1043  ## output particles
     1044  set OutputArray photons
     1045  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     1046  # efficiency formula for photons
     1047  set EfficiencyFormula {                      (pt <= 10.0) * (0.00) + \
     1048                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.9635) + \
     1049         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0)  * (0.9624) + \
     1050         (abs(eta) > 4.0)                                   * (0.00)}
     1051
     1052}
     1053
    9311054######################
    9321055# Electron isolation #
    … …  
    9491072}
    9501073
     1074
     1075##########################
     1076# Electron isolation CHS #
     1077##########################
     1078
     1079module Isolation ElectronIsolationCHS {
     1080
     1081  set CandidateInputArray ElectronFilter/electrons
     1082
     1083  # isolation collection
     1084  set IsolationInputArray EFlowMerger/eflow
     1085
     1086  set OutputArray electrons
     1087
     1088  set DeltaRMax 0.3
     1089  set PTMin 1.0
     1090  set PTRatioMax 9999.
     1091
     1092}
    9511093
    9521094
    … …  
    9951137}
    9961138
     1139###########################
     1140# Electron efficiency CHS #
     1141###########################
     1142
     1143module Efficiency ElectronEfficiencyCHS {
     1144
     1145  set InputArray ElectronIsolationCHS/electrons
     1146  set OutputArray electrons
     1147
     1148  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     1149  # efficiency formula for electrons
     1150  set EfficiencyFormula {
     1151                                      (pt <= 4.0)  * (0.00) + \
     1152                         (abs(eta) <= 1.45 ) * (pt >  4.0 && pt <= 6.0)   * (0.50) + \
     1153                         (abs(eta) <= 1.45 ) * (pt >  6.0 && pt <= 8.0)   * (0.70) + \
     1154                         (abs(eta) <= 1.45 ) * (pt >  8.0 && pt <= 10.0)  * (0.85) + \
     1155                         (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0)  * (0.94) + \
     1156                         (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0)  * (0.97) + \
     1157                         (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0)  * (0.98) + \
     1158                         (abs(eta) <= 1.45 ) * (pt > 70.0 )  * (1.0) + \
     1159                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt >  4.0 && pt <= 10.0)   * (0.35) + \
     1160                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0)   * (0.40) + \
     1161                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0)   * (0.45) + \
     1162                         (abs(eta) > 1.45  && abs(eta) <= 1.55) * (pt > 70.0 )  * (0.55) + \
     1163                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt >  4.0 && pt <= 10.0)  * (0.75) + \
     1164                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0)  * (0.85) + \
     1165                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0)  * (0.95) + \
     1166                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0)  * (0.95) + \
     1167                         (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 )  * (1.0) + \
     1168                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt >  4.0 && pt <= 10.0)  * (0.65) + \
     1169                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0)  * (0.75) + \
     1170                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0)  * (0.90) + \
     1171                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0)  * (0.90) + \
     1172                         (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 )  * (0.90) + \
     1173                         (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
     1174                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
     1175                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
     1176                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
     1177                                          (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
     1178                                          (abs(eta) > 4.0) * (0.00)
     1179
     1180  }
     1181}
     1182
     1183
    9971184##################
    9981185# Muon isolation #
    … …  
    10131200}
    10141201
    1015 
    1016 ##################
    1017 # Muon Loose Id  #
    1018 ##################
     1202######################
     1203# Muon isolation CHS #
     1204######################
     1205
     1206module Isolation MuonIsolationCHS {
     1207  set CandidateInputArray MuonMomentumSmearing/muons
     1208
     1209  # isolation collection
     1210  set IsolationInputArray EFlowMerger/eflow
     1211
     1212  set OutputArray muons
     1213
     1214  set DeltaRMax 0.3
     1215  set PTMin 1.0
     1216  set PTRatioMax 9999.
     1217
     1218}
     1219
     1220
     1221#####################
     1222# Muon Loose Id     #
     1223#####################
    10191224
    10201225module Efficiency MuonLooseIdEfficiency {
    … …  
    10381243}
    10391244
     1245
     1246#####################
     1247# Muon Loose Id CHS #
     1248#####################
     1249
     1250module Efficiency MuonLooseIdEfficiencyCHS {
     1251    set InputArray MuonIsolationCHS/muons
     1252    set OutputArray muons
     1253    # tracking + LooseID efficiency formula for muons
     1254    source muonLooseId.tcl
     1255
     1256}
     1257
     1258
     1259######################
     1260# Muon Tight Id  CHS #
     1261######################
     1262
     1263module Efficiency MuonTightIdEfficiencyCHS {
     1264    set InputArray MuonIsolationCHS/muons
     1265    set OutputArray muons
     1266    # tracking + TightID efficiency formula for muons
     1267    source muonTightId.tcl
     1268}
    10401269
    10411270
    … …  
    12501479module StatusPidFilter GenParticleFilter {
    12511480
    1252     set InputArray  Delphes/allParticles
     1481    set InputArray Delphes/allParticles
    12531482    set OutputArray filteredParticles
    12541483    set PTMin 5.0
    … …  
    12781507  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    12791508  add Branch MuonTightIdEfficiency/muons MuonTight Muon
     1509
     1510  add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
     1511  add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
     1512  add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
     1513  add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon
    12801514
    12811515  add Branch JetEnergyScale/jets Jet Jet
    … …  
    12891523  add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
    12901524  add Branch ScalarHT/energy ScalarHT ScalarHT
    1291 }
     1525
     1526}
Note: See TracChangeset for help on using the changeset viewer.