Fork me on GitHub

Changeset 667a02a in git for cards/delphes_card_HLLHC.tcl


Ignore:
Timestamp:
Jun 8, 2018, 3:23:13 PM (6 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
e57c062
Parents:
001ee95 (diff), 17d0ab8 (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:

Merge branch 'master' of github.com:delphes/delphes

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_HLLHC.tcl

    r001ee95 r667a02a  
    1414
    1515  ParticlePropagator
     16  TrackMergerProp
     17
     18  DenseProp
     19  DenseMergeTracks
     20  DenseTrackFilter
    1621
    1722  ChargedHadronTrackingEfficiency
    … …  
    2328  MuonMomentumSmearing
    2429
    25   DenseProp
    26   DenseMergeTracks
    27   TrackMergerProp 
    2830  TrackMerger
    2931
    … …  
    5355  GenMissingET
    5456
    55   FastJetFinder
    56   FatJetFinder
     57  GenJetFinder02
     58  GenJetFinder04
     59  GenJetFinder08
     60
     61  FastJetFinder02
     62  FastJetFinder04
     63  FastJetFinder08
     64
     65  CaloJetFinder02
     66  CaloJetFinder04
     67  CaloJetFinder08
     68
     69  TrackJetFinder02
     70  TrackJetFinder04
     71  TrackJetFinder08
    5772
    5873  JetEnergyScale
    … …  
    7186
    7287
    73 #################################
    74 # Propagate particles in cylinder
    75 #################################
     88#####################################
     89# Track propagation to calorimeters
     90#####################################
    7691
    7792module ParticlePropagator ParticlePropagator {
    … …  
    92107}
    93108
     109##############
     110# Track merger
     111##############
     112
     113module Merger TrackMergerProp {
     114# add InputArray InputArray
     115  add InputArray ParticlePropagator/chargedHadrons
     116  add InputArray ParticlePropagator/electrons
     117  add InputArray ParticlePropagator/muons
     118  set OutputArray tracks
     119}
     120
     121####################################
     122# Track propagation to pseudo-pixel
     123####################################
     124
     125module ParticlePropagator DenseProp {
     126
     127  set InputArray TrackMergerProp/tracks
     128
     129  # radius of the magnetic field coverage, in m
     130  set Radius 0.3
     131  set RadiusMax 1.2
     132  # half-length of the magnetic field coverage, in m
     133  set HalfLength 0.7
     134  set HalfLengthMax 3.25
     135
     136  # magnetic field
     137  set Bz 3.0
     138}
     139
     140#####################
     141# Dense Track merger
     142#####################
     143
     144module Merger DenseMergeTracks {
     145# add InputArray InputArray
     146  add InputArray DenseProp/chargedHadrons
     147  add InputArray DenseProp/electrons
     148  add InputArray DenseProp/muons
     149  set OutputArray tracks
     150}
     151
     152
     153######################
     154#   Dense Track Filter
     155######################
     156
     157module DenseTrackFilter DenseTrackFilter {
     158
     159  set TrackInputArray DenseMergeTracks/tracks
     160
     161  set TrackOutputArray tracks
     162  set ChargedHadronOutputArray chargedHadrons
     163  set ElectronOutputArray electrons
     164  set MuonOutputArray muons
     165
     166  set EtaPhiRes 0.003
     167  set EtaMax 4.0
     168
     169  set pi [expr {acos(-1)}]
     170
     171  set nbins_phi [expr {$pi/$EtaPhiRes} ]
     172  set nbins_phi [expr {int($nbins_phi)} ]
     173
     174  set PhiBins {}
     175  for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
     176    add PhiBins [expr {$i * $pi/$nbins_phi}]
     177  }
     178
     179  set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
     180  set nbins_eta [expr {int($nbins_eta)} ]
     181
     182  for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
     183    set eta [expr {$i * $EtaPhiRes}]
     184    add EtaPhiBins $eta $PhiBins
     185  }
     186}
     187
     188
     189
    94190####################################
    95191# Charged hadron tracking efficiency
    … …  
    97193
    98194module Efficiency ChargedHadronTrackingEfficiency {
    99   set InputArray ParticlePropagator/chargedHadrons
     195  set InputArray DenseTrackFilter/chargedHadrons
    100196  set OutputArray chargedHadrons
    101197
    … …  
    122218
    123219module Efficiency ElectronTrackingEfficiency {
    124   set InputArray ParticlePropagator/electrons
     220  set InputArray DenseTrackFilter/electrons
    125221  set OutputArray electrons
    126222
    … …  
    145241
    146242module Efficiency MuonTrackingEfficiency {
    147   set InputArray ParticlePropagator/muons
     243  set InputArray DenseTrackFilter/muons
    148244  set OutputArray muons
    149245
    … …  
    443539}
    444540
     541
     542
    445543##############
    446544# Track merger
    447545##############
    448546
    449 module Merger TrackMergerProp {
     547module Merger TrackMerger {
    450548# add InputArray InputArray
    451549  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
    … …  
    453551  add InputArray MuonMomentumSmearing/muons
    454552  set OutputArray tracks
    455 }
    456 
    457 
    458 #################################
    459 # Dense Track propagation
    460 #################################
    461 
    462 module ParticlePropagator DenseProp {
    463   set InputArray Delphes/stableParticles
    464 
    465   set OutputArray stableParticles
    466   set NeutralOutputArray neutralParticles
    467   set ChargedHadronOutputArray chargedHadrons
    468   set ElectronOutputArray electrons
    469   set MuonOutputArray muons
    470 
    471   # radius of the magnetic field coverage, in m
    472   set Radius 0.3
    473   set RadiusMax 1.2
    474   # half-length of the magnetic field coverage, in m
    475   set HalfLength 0.7
    476   set HalfLengthMax 3.25
    477 
    478   # magnetic field
    479   set Bz 3.0
    480 }
    481 
    482 
    483 ####################
    484 # Dense Track merger
    485 ###################
    486 
    487 module Merger DenseMergeTracks {
    488 # add InputArray InputArray
    489   add InputArray DenseProp/chargedHadrons
    490   add InputArray DenseProp/electrons
    491   add InputArray DenseProp/muons
    492   set OutputArray tracks
    493 }
    494 
    495 
    496 ######################
    497 #   Dense Track Filter
    498 ######################
    499 
    500 module DenseTrackFilter TrackMerger {
    501   set TrackInputArray TrackMergerProp/tracks
    502   set DenseChargedInputArray DenseMergeTracks/tracks
    503 
    504   set TrackOutputArray tracks
    505 
    506   set EtaPhiRes 0.003
    507   set EtaMax 4.0
    508 
    509   set pi [expr {acos(-1)}]
    510 
    511   set nbins_phi [expr {$pi/$EtaPhiRes} ]
    512   set nbins_phi [expr {int($nbins_phi)} ]
    513 
    514   set PhiBins {}
    515   for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
    516     add PhiBins [expr {$i * $pi/$nbins_phi}]
    517   }
    518 
    519   set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
    520   set nbins_eta [expr {int($nbins_eta)} ]
    521 
    522   for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
    523     set eta [expr {$i * $EtaPhiRes}]
    524     add EtaPhiBins $eta $PhiBins
    525   }
    526553}
    527554
    … …  
    841868
    842869
    843 ############
    844 # Jet finder
    845 ############
    846 
    847 module FastJetFinder FastJetFinder {
    848   set InputArray Calorimeter/towers
    849 #  set InputArray EFlowMerger/eflow
    850 
    851   set OutputArray jets
    852 
     870
     871#####################
     872# MC truth jet finder
     873#####################
     874
     875# TBC: is jet radius fine?
     876
     877module FastJetFinder GenJetFinder02 {
     878  set InputArray NeutrinoFilter/filteredParticles
     879
     880  set OutputArray jets
     881
     882  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     883  set JetAlgorithm 6
     884  set ParameterR 0.2
     885
     886  set ComputeNsubjettiness 1
     887  set Beta 1.0
     888  set AxisMode 4
     889
     890  set ComputeSoftDrop 1
     891  set BetaSoftDrop 0.0
     892  set SymmetryCutSoftDrop 0.1
     893  set R0SoftDrop 0.2
     894
     895  set JetPTMin 25.0
     896}
     897
     898
     899#####################
     900# MC truth jet finder
     901#####################
     902
     903# TBC: is jet radius fine?
     904
     905module FastJetFinder GenJetFinder04 {
     906  set InputArray NeutrinoFilter/filteredParticles
     907
     908  set OutputArray jets
     909
     910  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    853911  set JetAlgorithm 6
    854912  set ParameterR 0.4
    … …  
    863921  set R0SoftDrop 0.4
    864922
    865   set JetPTMin 30.0
    866 }
    867 
    868 ##################
    869 # Fat Jet finder
    870 ##################
    871 
    872 module FastJetFinder FatJetFinder {
    873   set InputArray EFlowMerger/eflow
     923  set JetPTMin 25.0
     924}
     925#####################
     926# MC truth jet finder
     927#####################
     928
     929# TBC: is jet radius fine?
     930
     931module FastJetFinder GenJetFinder08 {
     932  set InputArray NeutrinoFilter/filteredParticles
    874933
    875934  set OutputArray jets
    … …  
    888947  set R0SoftDrop 0.8
    889948
    890   set JetPTMin 200.0
    891 }
    892 
     949  set JetPTMin 25.0
     950}
     951
     952#####################
     953# MC truth jet finder
     954#####################
     955
     956# TBC: is jet radius fine?
     957
     958module FastJetFinder GenJetFinder15 {
     959  set InputArray NeutrinoFilter/filteredParticles
     960
     961  set OutputArray jets
     962
     963  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     964  set JetAlgorithm 6
     965  set ParameterR 1.5
     966
     967  set ComputeNsubjettiness 1
     968  set Beta 1.0
     969  set AxisMode 4
     970
     971  set ComputeSoftDrop 1
     972  set BetaSoftDrop 0.0
     973  set SymmetryCutSoftDrop 0.1
     974  set R0SoftDrop 1.5
     975
     976  set JetPTMin 25.0
     977}
     978
     979
     980##################
     981# Fast Jet finder
     982##################
     983
     984module FastJetFinder FastJetFinder02 {
     985  set InputArray EFlowMerger/eflow
     986
     987  set OutputArray jets
     988
     989  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     990  set JetAlgorithm 6
     991  set ParameterR 0.2
     992
     993  set ComputeNsubjettiness 1
     994  set Beta 1.0
     995  set AxisMode 4
     996
     997  set ComputeSoftDrop 1
     998  set BetaSoftDrop 0.0
     999  set SymmetryCutSoftDrop 0.1
     1000  set R0SoftDrop 0.2
     1001
     1002  set JetPTMin 25.0
     1003}
     1004
     1005##################
     1006# Fast Jet finder
     1007##################
     1008
     1009module FastJetFinder FastJetFinder04 {
     1010  set InputArray EFlowMerger/eflow
     1011
     1012  set OutputArray jets
     1013
     1014  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1015  set JetAlgorithm 6
     1016  set ParameterR 0.4
     1017
     1018  set ComputeNsubjettiness 1
     1019  set Beta 1.0
     1020  set AxisMode 4
     1021
     1022  set ComputeSoftDrop 1
     1023  set BetaSoftDrop 0.0
     1024  set SymmetryCutSoftDrop 0.1
     1025  set R0SoftDrop 0.4
     1026
     1027  set JetPTMin 25.0
     1028}
     1029
     1030
     1031##################
     1032# Fast Jet finder
     1033##################
     1034
     1035module FastJetFinder FastJetFinder08 {
     1036  set InputArray EFlowMerger/eflow
     1037
     1038  set OutputArray jets
     1039
     1040  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1041  set JetAlgorithm 6
     1042  set ParameterR 0.8
     1043
     1044  set ComputeNsubjettiness 1
     1045  set Beta 1.0
     1046  set AxisMode 4
     1047
     1048  set ComputeSoftDrop 1
     1049  set BetaSoftDrop 0.0
     1050  set SymmetryCutSoftDrop 0.1
     1051  set R0SoftDrop 0.8
     1052
     1053  set JetPTMin 25.0
     1054}
     1055
     1056##################
     1057# Fast Jet finder
     1058##################
     1059
     1060module FastJetFinder FastJetFinder15 {
     1061  set InputArray EFlowMerger/eflow
     1062
     1063  set OutputArray jets
     1064
     1065  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1066  set JetAlgorithm 6
     1067  set ParameterR 1.5
     1068
     1069  set ComputeNsubjettiness 1
     1070  set Beta 1.0
     1071  set AxisMode 4
     1072
     1073  set ComputeSoftDrop 1
     1074  set BetaSoftDrop 0.0
     1075  set SymmetryCutSoftDrop 0.1
     1076  set R0SoftDrop 1.5
     1077
     1078  set JetPTMin 25.0
     1079}
     1080
     1081
     1082##################
     1083# Fast Jet finder
     1084##################
     1085
     1086module FastJetFinder CaloJetFinder02 {
     1087  set InputArray Calorimeter/towers
     1088
     1089  set OutputArray jets
     1090
     1091  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1092  set JetAlgorithm 6
     1093  set ParameterR 0.2
     1094
     1095  set ComputeNsubjettiness 1
     1096  set Beta 1.0
     1097  set AxisMode 4
     1098
     1099  set ComputeSoftDrop 1
     1100  set BetaSoftDrop 0.0
     1101  set SymmetryCutSoftDrop 0.1
     1102  set R0SoftDrop 0.2
     1103
     1104  set JetPTMin 25.0
     1105}
     1106
     1107##################
     1108# Fast Jet finder
     1109##################
     1110
     1111module FastJetFinder CaloJetFinder04 {
     1112  set InputArray Calorimeter/towers
     1113
     1114  set OutputArray jets
     1115
     1116  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1117  set JetAlgorithm 6
     1118  set ParameterR 0.4
     1119
     1120  set ComputeNsubjettiness 1
     1121  set Beta 1.0
     1122  set AxisMode 4
     1123
     1124  set ComputeSoftDrop 1
     1125  set BetaSoftDrop 0.0
     1126  set SymmetryCutSoftDrop 0.1
     1127  set R0SoftDrop 0.4
     1128
     1129  set JetPTMin 25.0
     1130}
     1131
     1132
     1133##################
     1134# Fast Jet finder
     1135##################
     1136
     1137module FastJetFinder CaloJetFinder08 {
     1138  set InputArray Calorimeter/towers
     1139
     1140  set OutputArray jets
     1141
     1142  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1143  set JetAlgorithm 6
     1144  set ParameterR 0.8
     1145
     1146  set ComputeNsubjettiness 1
     1147  set Beta 1.0
     1148  set AxisMode 4
     1149
     1150  set ComputeSoftDrop 1
     1151  set BetaSoftDrop 0.0
     1152  set SymmetryCutSoftDrop 0.1
     1153  set R0SoftDrop 0.8
     1154
     1155  set JetPTMin 25.0
     1156}
     1157
     1158##################
     1159# Fast Jet finder
     1160##################
     1161
     1162module FastJetFinder CaloJetFinder15 {
     1163  set InputArray Calorimeter/towers
     1164
     1165  set OutputArray jets
     1166
     1167  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1168  set JetAlgorithm 6
     1169  set ParameterR 1.5
     1170
     1171  set ComputeNsubjettiness 1
     1172  set Beta 1.0
     1173  set AxisMode 4
     1174
     1175  set ComputeSoftDrop 1
     1176  set BetaSoftDrop 0.0
     1177  set SymmetryCutSoftDrop 0.1
     1178  set R0SoftDrop 1.5
     1179
     1180  set JetPTMin 25.0
     1181}
     1182
     1183
     1184##################
     1185# Fast Jet finder
     1186##################
     1187
     1188module FastJetFinder TrackJetFinder02 {
     1189  set InputArray TrackMerger/tracks
     1190
     1191  set OutputArray jets
     1192
     1193  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1194  set JetAlgorithm 6
     1195  set ParameterR 0.2
     1196
     1197  set ComputeNsubjettiness 1
     1198  set Beta 1.0
     1199  set AxisMode 4
     1200
     1201  set ComputeSoftDrop 1
     1202  set BetaSoftDrop 0.0
     1203  set SymmetryCutSoftDrop 0.1
     1204  set R0SoftDrop 0.2
     1205
     1206  set JetPTMin 25.0
     1207}
     1208
     1209##################
     1210# Fast Jet finder
     1211##################
     1212
     1213module FastJetFinder TrackJetFinder04 {
     1214  set InputArray TrackMerger/tracks
     1215
     1216  set OutputArray jets
     1217
     1218  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1219  set JetAlgorithm 6
     1220  set ParameterR 0.4
     1221
     1222  set ComputeNsubjettiness 1
     1223  set Beta 1.0
     1224  set AxisMode 4
     1225
     1226  set ComputeSoftDrop 1
     1227  set BetaSoftDrop 0.0
     1228  set SymmetryCutSoftDrop 0.1
     1229  set R0SoftDrop 0.4
     1230
     1231  set JetPTMin 25.0
     1232}
     1233
     1234
     1235##################
     1236# Fast Jet finder
     1237##################
     1238
     1239module FastJetFinder TrackJetFinder08 {
     1240  set InputArray TrackMerger/tracks
     1241
     1242  set OutputArray jets
     1243
     1244  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1245  set JetAlgorithm 6
     1246  set ParameterR 0.8
     1247
     1248  set ComputeNsubjettiness 1
     1249  set Beta 1.0
     1250  set AxisMode 4
     1251
     1252  set ComputeSoftDrop 1
     1253  set BetaSoftDrop 0.0
     1254  set SymmetryCutSoftDrop 0.1
     1255  set R0SoftDrop 0.8
     1256
     1257  set JetPTMin 25.0
     1258}
     1259
     1260##################
     1261# Fast Jet finder
     1262##################
     1263
     1264module FastJetFinder TrackJetFinder15 {
     1265  set InputArray TrackMerger/tracks
     1266
     1267  set OutputArray jets
     1268
     1269  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     1270  set JetAlgorithm 6
     1271  set ParameterR 1.5
     1272
     1273  set ComputeNsubjettiness 1
     1274  set Beta 1.0
     1275  set AxisMode 4
     1276
     1277  set ComputeSoftDrop 1
     1278  set BetaSoftDrop 0.0
     1279  set SymmetryCutSoftDrop 0.1
     1280  set R0SoftDrop 1.5
     1281
     1282  set JetPTMin 25.0
     1283}
    8931284
    8941285
    … …  
    8981289
    8991290module EnergyScale JetEnergyScale {
    900   set InputArray FastJetFinder/jets
     1291  set InputArray CaloJetFinder04/jets
    9011292  set OutputArray jets
    9021293
    … …  
    12581649  add Branch UniqueObjectFinder/jets Jet Jet
    12591650
    1260   add Branch FatJetFinder/jets FatJet Jet
     1651  add Branch GenJetFinder02/jets GenJet02 Jet
     1652  add Branch GenJetFinder04/jets GenJet04 Jet
     1653  add Branch GenJetFinder08/jets GenJet08 Jet
     1654  add Branch GenJetFinder15/jets GenJet15 Jet
     1655
     1656  add Branch FastJetFinder02/jets ParticleFlowJet02 Jet
     1657  add Branch FastJetFinder04/jets ParticleFlowJet04 Jet
     1658  add Branch FastJetFinder08/jets ParticleFlowJet08 Jet
     1659  add Branch FastJetFinder15/jets ParticleFlowJet15 Jet
     1660
     1661  add Branch CaloJetFinder02/jets CaloJet02 Jet
     1662  add Branch CaloJetFinder04/jets CaloJet04 Jet
     1663  add Branch CaloJetFinder08/jets CaloJet08 Jet
     1664  add Branch CaloJetFinder15/jets CaloJet15 Jet
     1665
     1666  add Branch TrackJetFinder02/jets TrackJet02 Jet
     1667  add Branch TrackJetFinder04/jets TrackJet04 Jet
     1668  add Branch TrackJetFinder08/jets TrackJet08 Jet
     1669  add Branch TrackJetFinder15/jets TrackJet15 Jet
    12611670
    12621671  add Branch MissingET/momentum MissingET MissingET
Note: See TracChangeset for help on using the changeset viewer.