← Back to branch summary

~maddevelopers/mg5amcnlo/laser_beam_option

« back to all changes in this revision

Viewing changes to Template/Cards/delphes_card_ATLAS.dat

  • Committer: olivier Mattelaer
  • Date: 2013-03-05 18:52:27 UTC
  • mfrom: (239.1.31 mg5_1_5_8)
  • Revision ID: olivier.mattelaer@uclouvain.be-20130305185227-h1358ux804s2wc26
Pass to 1.5.8

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/Cards/delphes_card_ATLAS.dat
1
 
DETECTOR CARD      # DO NOT REMOVE THIS IS A TAG!
2
 
 
3
 
# Detector extension, in pseudorapidity units
4
 
CEN_max_tracker    2.5                // Maximum tracker coverage          
5
 
CEN_max_calo_cen   3.2                // central calorimeter coverage
6
 
CEN_max_calo_fwd   4.9                // forward calorimeter pseudorapidity coverage
7
 
CEN_max_mu         2.7                // muon chambers pseudorapidity coverage
8
 
  
9
 
# Energy resolution for electron/photon
10
 
# \sigma/E = C + N/E + S/\sqrt{E}, E in GeV
11
 
ELG_Scen          0.10                // S term for central ECAL  #http://arxiv.org/pdf/physics/0608012v1
12
 
ELG_Ncen          0.0                 // N term for central ECAL  #http://arxiv.org/pdf/physics/0608012v1
13
 
ELG_Ccen          0.0017              // C term for central ECAL  #http://arxiv.org/pdf/physics/0608012v1
14
 
ELG_Sfwd          0.30407             // S term for forward FCAL  #http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
15
 
ELG_Nfwd          1.1533              // N term for central FCAL  #http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
16
 
ELG_Cfwd          0.0313              // C term for forward FCAL  #http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
17
 
ELG_Szdc          0.70              // S term for ZDC
18
 
ELG_Nzdc          0.0               // N term for ZDC
19
 
ELG_Czdc          0.08              // C term for ZDC
20
 
 
21
 
 
22
 
# Energy resolution for hadrons in ecal/hcal/hf
23
 
# \sigma/E = C + N/E + S/\sqrt{E}, E in GeV
24
 
HAD_Shcal         0.58               // S term for central HCAL    #http://arxiv.org/pdf/hep-ex/0004009v1
25
 
HAD_Nhcal         1.7                // N term for central HCAL    #http://arxiv.org/pdf/hep-ex/0004009v1
26
 
HAD_Chcal         0.025              // C term for central HCAL    #http://arxiv.org/pdf/hep-ex/0004009v1
27
 
HAD_Shf           0.9681             // S term for FCAL            #http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
28
 
HAD_Nhf           0.                 // N term for FCAL            #http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
29
 
HAD_Chf           0.0524             // C term for FCAL            #http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
30
 
HAD_Szdc          1.38              // S term for ZDC
31
 
HAD_Nzdc          0.                // N term for ZDC
32
 
HAD_Czdc          0.13              // C term for ZDC
33
 
 
34
 
# Time resolution for ZDC/RP220/RP420
35
 
ZDC_T_resolution   0                // in s
36
 
RP220_T_resolution 0                // in s
37
 
RP420_T_resolution 0                // in s
38
 
 
39
 
# Muon smearing
40
 
MU_SmearPt        0.01              // transverse momentum Pt in GeV
41
 
 
42
 
# Tracking efficiencies
43
 
TRACK_ptmin       0.9               // minimal pt needed to reach the calorimeter in GeV
44
 
TRACK_eff          90               // efficiency associated to the tracking (%)
45
 
 
46
 
# Calorimetric towers
47
 
TOWER_number         38
48
 
TOWER_eta_edges 0.0 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 1.10 1.20 1.30 1.40 1.50 1.60 1.70 1.80 1.90 2.0 2.10 2.20 2.30 2.40 2.50 2.60 2.80 3.00 3.20 3.30 3.50 3.70 3.90 4.10 4.30 4.50 4.70 4.90 
49
 
## list of the edges of each tower in eta for eta>0 assuming a symmetric detector in eta<0
50
 
### the list starts with the lower edge of the most central tower
51
 
### the list ends with the higher edged of the most forward tower
52
 
### there should be NTOWER+1 values
53
 
TOWER_dphi  10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 20 20 20 10 20 20 20 20 20 20 20 20
54
 
### list of the tower size in phi (in degrees), assuming that all towers are similar in phi for a given eta value
55
 
### the list starts with the phi-size of the most central tower (eta=0)
56
 
### the list ends with the phi-size of the most forward tower
57
 
### there should be NTOWER values
58
 
##
59
 
#
60
 
 
61
 
# Thresholds for reconstructed objetcs, Pt in GeV
62
 
PTCUT_elec       10.0
63
 
PTCUT_muon       10.0
64
 
PTCUT_jet        20.0
65
 
PTCUT_gamma      10.0
66
 
PTCUT_taujet     10.0
67
 
 
68
 
# Thresholds for reconstructed objects in ZDC, E in GeV
69
 
ZDC_gamma_E      20
70
 
ZDC_n_E          50
71
 
 
72
 
# Charged lepton isolation. Pt and Et in GeV
73
 
ISOL_PT          2.0      //minimal pt of tracks for isolation criteria
74
 
ISOL_Cone        0.5      //Cone  for isolation criteria
75
 
ISOL_Calo_ET     1E99     //minimal tower transverse energy for isolation criteria. 1E99 means "off"
76
 
ISOL_Calo_Cone   0.4      //Cone for calorimetric isolation
77
 
ISOL_Calo_Grid   3        //Grid size (N x N) for calorimetric isolation
78
 
 
79
 
# General jet variable
80
 
JET_coneradius   0.7            // generic jet radius ; not for tau's !!!
81
 
JET_jetalgo      1              // 1 for Cone algorithm, 2 for MidPoint algorithm, 3 for SIScone algorithm, 4 for kt algorithm
82
 
JET_seed         1.0            // minimum seed to start jet reconstruction, in GeV
83
 
JET_Eflow        0              // Perfect energy assumed in the tracker coverage
84
 
 
85
 
# Tagging definition
86
 
BTAG_b           50             // b-tag efficiency (%)
87
 
BTAG_mistag_c    10             // mistagging (%)
88
 
BTAG_mistag_l    1              // mistagging (%)
89
 
 
90
 
# FLAGS
91
 
FLAG_bfield      1                       // 1 to run the bfield propagation else 0
92
 
FLAG_vfd         1                       // 1 to run the very forward detectors else 0
93
 
FLAG_RP          1                       // 1 to run the very forward detectors else 0
94
 
FLAG_trigger     1                       // 1 to run the trigger selection else 0
95
 
FLAG_frog        0                       // 1 to run the FROG event display
96
 
FLAG_lhco        0                       //1 to run the LHCO
97
 
 
98
 
# In case BField propagation allowed
99
 
TRACK_radius      115                   // radius of the BField coverage, in cm #chapter 1 you gave me
100
 
TRACK_length      351                   // length of the BField coverage, in cm #also
101
 
TRACK_bfield_x    0                     // X component of the BField, in T
102
 
TRACK_bfield_y    0                     // Y component of the BField, in T
103
 
TRACK_bfield_z    2.0                   // Z component of the BField, in T #also
104
 
 
105
 
# Very forward detector extension, in pseudorapidity
106
 
# if allowed (WARNING: not expected for ATLAS!)
107
 
VFD_min_calo_vfd  5.2                   // very forward calorimeter (if any) like CASTOR
108
 
VFD_max_calo_vfd  6.6
109
 
VFD_min_zdc       8.3                   // zero-degree neutral calorimeter
110
 
VFD_s_zdc         140                   // distance of the Zero Degree Calorimeter, from the IP, in [m]
111
 
 
112
 
RP_220_s          220                   // distance of the RP to the IP, in meters
113
 
RP_220_x          0.002                 // distance of the RP to the beam, in meters
114
 
RP_420_s          420                   // distance of the RP to the IP, in meters
115
 
RP_420_x          0.004                 // distance of the RP to the beam, in meters
116
 
RP_beam1Card      DELPHESDIR/data/LHCB1IR1.tfs
117
 
RP_beam2Card      DELPHESDIR/data/LHCB2IR1.tfs
118
 
RP_IP_name        IP1
119
 
RP_offsetEl_x     0.097                 // horizontal separation between both beam, in meters
120
 
RP_offsetEl_y     0                     // vertical separation between both beam, in meters
121
 
RP_offsetEl_s     120                   // distance of beam separation point, from IP
122
 
RP_cross_x        0                     // IP offset in horizontal plane, in micrometers
123
 
RP_cross_y        +500                  // IP offset in vertical plane, in micrometers
124
 
RP_cross_ang_x    142.5                 // half-crossing angle in horizontal plane, in microrad
125
 
RP_cross_ang_y    0                     // half-crossing angle in vertical plane, in microrad
126
 
 
127
 
 
128
 
# In case FROG event display allowed
129
 
NEvents_Frog      100
130
 
# Number of events to process
131
 
NEvents           -1                    // -1 means 'all'
132
 
 
133
 
# input PDG tables
134
 
PdgTableFilename  DELPHESDIR/data/particle.tbl     // table with particle pid,mass,charge,...
 
1
#######################################
 
2
# Order of execution of various modules
 
3
#######################################
 
4
 
 
5
set ExecutionPath {
 
6
  ParticlePropagator
 
7
 
 
8
  ChargedHadronTrackingEfficiency
 
9
  ElectronTrackingEfficiency
 
10
  MuonTrackingEfficiency
 
11
 
 
12
  ChargedHadronMomentumSmearing
 
13
  ElectronEnergySmearing
 
14
  MuonMomentumSmearing
 
15
 
 
16
  TrackMerger
 
17
  Calorimeter
 
18
  EFlowMerger
 
19
 
 
20
  PhotonEfficiency
 
21
  PhotonIsolation
 
22
  
 
23
  ElectronEfficiency
 
24
  ElectronIsolation
 
25
  
 
26
  MuonEfficiency
 
27
  MuonIsolation
 
28
  
 
29
  MissingET
 
30
  
 
31
  FastJetFinder
 
32
  BTagging
 
33
  TauTagging
 
34
  
 
35
  UniqueObjectFinder
 
36
  
 
37
  ScalarHT
 
38
  
 
39
  TreeWriter
 
40
}
 
41
 
 
42
#################################
 
43
# Propagate particles in cylinder
 
44
#################################
 
45
 
 
46
module ParticlePropagator ParticlePropagator {
 
47
  set InputArray Delphes/stableParticles
 
48
 
 
49
  set OutputArray stableParticles
 
50
  set ChargedHadronOutputArray chargedHadrons
 
51
  set ElectronOutputArray electrons
 
52
  set MuonOutputArray muons
 
53
 
 
54
  # radius of the magnetic field coverage, in m
 
55
  set Radius 1.15
 
56
  # half-length of the magnetic field coverage, in m
 
57
  set HalfLength 3.51
 
58
 
 
59
  # magnetic field
 
60
  set Bz 2.0
 
61
}
 
62
 
 
63
####################################
 
64
# Charged hadron tracking efficiency
 
65
####################################
 
66
 
 
67
module Efficiency ChargedHadronTrackingEfficiency {
 
68
  set InputArray ParticlePropagator/chargedHadrons
 
69
  set OutputArray chargedHadrons
 
70
 
 
71
  # add EfficiencyFormula {efficiency formula as a function of eta and pt}
 
72
 
 
73
  # tracking efficiency formula for charged hadrons
 
74
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
 
75
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
 
76
                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
 
77
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
 
78
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
 
79
                         (abs(eta) > 2.5)                                                  * (0.00)}
 
80
}
 
81
 
 
82
##############################
 
83
# Electron tracking efficiency
 
84
##############################
 
85
 
 
86
module Efficiency ElectronTrackingEfficiency {
 
87
  set InputArray ParticlePropagator/electrons
 
88
  set OutputArray electrons
 
89
 
 
90
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
 
91
              
 
92
  # tracking efficiency formula for electrons
 
93
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
 
94
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
 
95
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
 
96
                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
 
97
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
 
98
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
 
99
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
 
100
                         (abs(eta) > 2.5)                                                  * (0.00)}
 
101
}
 
102
 
 
103
##############################
 
104
# Muon tracking efficiency
 
105
##############################
 
106
 
 
107
module Efficiency MuonTrackingEfficiency {
 
108
  set InputArray ParticlePropagator/muons
 
109
  set OutputArray muons
 
110
 
 
111
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
 
112
              
 
113
  # tracking efficiency formula for muons
 
114
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
 
115
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
 
116
                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
 
117
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
 
118
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
 
119
                         (abs(eta) > 2.5)                                                  * (0.00)}
 
120
}
 
121
 
 
122
########################################
 
123
# Momentum resolution for charged tracks
 
124
########################################
 
125
 
 
126
module MomentumSmearing ChargedHadronMomentumSmearing {
 
127
  set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
 
128
  set OutputArray chargedHadrons
 
129
 
 
130
  # set ResolutionFormula {resolution formula as a function of eta and pt}
 
131
 
 
132
  # resolution formula for charged hadrons
 
133
  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
 
134
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
 
135
                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
 
136
                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
 
137
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
 
138
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
 
139
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
 
140
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
 
141
}
 
142
 
 
143
#################################
 
144
# Energy resolution for electrons
 
145
#################################
 
146
 
 
147
module EnergySmearing ElectronEnergySmearing {
 
148
  set InputArray ElectronTrackingEfficiency/electrons
 
149
  set OutputArray electrons
 
150
 
 
151
  # set ResolutionFormula {resolution formula as a function of eta and energy}
 
152
 
 
153
  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
 
154
                                           (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
 
155
                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
 
156
                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
 
157
 
 
158
}
 
159
 
 
160
###############################
 
161
# Momentum resolution for muons
 
162
###############################
 
163
 
 
164
module MomentumSmearing MuonMomentumSmearing {
 
165
  set InputArray MuonTrackingEfficiency/muons
 
166
  set OutputArray muons
 
167
 
 
168
  # set ResolutionFormula {resolution formula as a function of eta and pt}
 
169
 
 
170
  # resolution formula for muons
 
171
  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
 
172
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
 
173
                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
 
174
                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
 
175
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
 
176
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.03) + \
 
177
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
 
178
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
 
179
}
 
180
 
 
181
##############
 
182
# Track merger
 
183
##############
 
184
 
 
185
module Merger TrackMerger {
 
186
# add InputArray InputArray
 
187
  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
 
188
  add InputArray ElectronEnergySmearing/electrons
 
189
  set OutputArray tracks
 
190
}
 
191
 
 
192
#############
 
193
# Calorimeter
 
194
#############
 
195
 
 
196
module Calorimeter Calorimeter {
 
197
  set ParticleInputArray ParticlePropagator/stableParticles
 
198
  set TrackInputArray TrackMerger/tracks
 
199
 
 
200
  set TowerOutputArray towers
 
201
  set PhotonOutputArray photons
 
202
 
 
203
  set EFlowTrackOutputArray eflowTracks
 
204
  set EFlowTowerOutputArray eflowTowers
 
205
 
 
206
  set pi [expr {acos(-1)}]
 
207
  
 
208
  # lists of the edges of each tower in eta and phi
 
209
  # each list starts with the lower edge of the first tower
 
210
  # the list ends with the higher edged of the last tower
 
211
 
 
212
  # 10 degrees towers
 
213
  set PhiBins {}
 
214
  for {set i -18} {$i <= 18} {incr i} { 
 
215
    add PhiBins [expr {$i * $pi/18.0}]
 
216
  }
 
217
  foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} { 
 
218
    add EtaPhiBins $eta $PhiBins
 
219
  }
 
220
 
 
221
  # 20 degrees towers
 
222
  set PhiBins {}
 
223
  for {set i -9} {$i <= 9} {incr i} { 
 
224
    add PhiBins [expr {$i * $pi/9.0}]
 
225
  }  
 
226
  foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.3 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} { 
 
227
    add EtaPhiBins $eta $PhiBins
 
228
  }
 
229
 
 
230
  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
 
231
  add EnergyFraction {0} {0.0 1.0}
 
232
  # energy fractions for e, gamma and pi0
 
233
  add EnergyFraction {11} {1.0 0.0}
 
234
  add EnergyFraction {22} {1.0 0.0}
 
235
  add EnergyFraction {111} {1.0 0.0}
 
236
  # energy fractions for muon and neutrinos
 
237
  add EnergyFraction {12} {0.0 0.0}
 
238
  add EnergyFraction {13} {0.0 0.0}
 
239
  add EnergyFraction {14} {0.0 0.0}
 
240
  add EnergyFraction {16} {0.0 0.0}
 
241
  # energy fractions for K0short and Lambda
 
242
  add EnergyFraction {310} {0.3 0.7}
 
243
  add EnergyFraction {3122} {0.3 0.7}
 
244
  
 
245
  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
 
246
  # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
 
247
  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
 
248
  set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
 
249
                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*2.085^2)}
 
250
 
 
251
  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
 
252
  # http://arxiv.org/pdf/hep-ex/0004009v1
 
253
  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
 
254
  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
 
255
                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
 
256
                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
 
257
}
 
258
 
 
259
####################
 
260
# Energy flow merger
 
261
####################
 
262
 
 
263
module Merger EFlowMerger {
 
264
# add InputArray InputArray
 
265
  add InputArray Calorimeter/eflowTracks
 
266
  add InputArray Calorimeter/eflowTowers
 
267
  add InputArray MuonMomentumSmearing/muons
 
268
  set OutputArray eflow
 
269
}
 
270
 
 
271
###################
 
272
# Photon efficiency
 
273
###################
 
274
 
 
275
module Efficiency PhotonEfficiency {
 
276
  set InputArray Calorimeter/photons
 
277
  set OutputArray photons
 
278
 
 
279
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
 
280
 
 
281
  # efficiency formula for photons
 
282
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
 
283
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
 
284
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
 
285
                         (abs(eta) > 2.5)                                   * (0.00)}
 
286
}
 
287
 
 
288
##################
 
289
# Photon isolation
 
290
##################
 
291
 
 
292
module Isolation PhotonIsolation {
 
293
  set CandidateInputArray PhotonEfficiency/photons
 
294
  set IsolationInputArray EFlowMerger/eflow
 
295
  
 
296
  set OutputArray photons
 
297
 
 
298
  set DeltaRMax 0.5
 
299
 
 
300
  set PTMin 0.5
 
301
 
 
302
  set PTRatioMax 0.1
 
303
}
 
304
 
 
305
#####################
 
306
# Electron efficiency
 
307
#####################
 
308
 
 
309
module Efficiency ElectronEfficiency {
 
310
  set InputArray ElectronEnergySmearing/electrons
 
311
  set OutputArray electrons
 
312
 
 
313
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
 
314
 
 
315
  # efficiency formula for electrons
 
316
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
 
317
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
 
318
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
 
319
                         (abs(eta) > 2.5)                                   * (0.00)}
 
320
}
 
321
 
 
322
####################
 
323
# Electron isolation
 
324
####################
 
325
 
 
326
module Isolation ElectronIsolation {
 
327
  set CandidateInputArray ElectronEfficiency/electrons
 
328
  set IsolationInputArray EFlowMerger/eflow
 
329
  
 
330
  set OutputArray electrons
 
331
 
 
332
  set DeltaRMax 0.5
 
333
 
 
334
  set PTMin 0.5
 
335
 
 
336
  set PTRatioMax 0.1
 
337
}
 
338
 
 
339
#################
 
340
# Muon efficiency
 
341
#################
 
342
 
 
343
module Efficiency MuonEfficiency {
 
344
  set InputArray MuonMomentumSmearing/muons
 
345
  set OutputArray muons
 
346
 
 
347
  # set EfficiencyFormula {efficiency as a function of eta and pt}
 
348
 
 
349
  # efficiency formula for muons
 
350
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
 
351
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
 
352
                         (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) + \
 
353
                         (abs(eta) > 2.7)                                   * (0.00)}
 
354
}
 
355
 
 
356
################
 
357
# Muon isolation
 
358
################
 
359
 
 
360
module Isolation MuonIsolation {
 
361
  set CandidateInputArray MuonEfficiency/muons
 
362
  set IsolationInputArray EFlowMerger/eflow
 
363
  
 
364
  set OutputArray muons
 
365
 
 
366
  set DeltaRMax 0.5
 
367
 
 
368
  set PTMin 0.5
 
369
 
 
370
  set PTRatioMax 0.1
 
371
}
 
372
 
 
373
###################
 
374
# Missing ET merger
 
375
###################
 
376
 
 
377
module Merger MissingET {
 
378
# add InputArray InputArray
 
379
  add InputArray EFlowMerger/eflow
 
380
  set MomentumOutputArray momentum
 
381
}
 
382
 
 
383
##################
 
384
# Scalar HT merger
 
385
##################
 
386
 
 
387
module Merger ScalarHT {
 
388
# add InputArray InputArray
 
389
  add InputArray UniqueObjectFinder/jets
 
390
  add InputArray UniqueObjectFinder/electrons
 
391
  add InputArray UniqueObjectFinder/photons
 
392
  add InputArray MuonIsolation/muons
 
393
  set EnergyOutputArray energy
 
394
}
 
395
 
 
396
############
 
397
# Jet finder
 
398
############
 
399
 
 
400
module FastJetFinder FastJetFinder {
 
401
#  set InputArray Calorimeter/towers
 
402
  set InputArray EFlowMerger/eflow
 
403
  
 
404
  set OutputArray jets
 
405
 
 
406
  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
 
407
  set JetAlgorithm 6
 
408
  set ParameterR 0.7
 
409
 
 
410
  set ConeRadius 0.5
 
411
  set SeedThreshold 1.0
 
412
  set ConeAreaFraction 1.0
 
413
  set AdjacencyCut 2.0
 
414
  set OverlapThreshold 0.75
 
415
 
 
416
  set MaxIterations 100
 
417
  set MaxPairSize 2
 
418
  set Iratch 1
 
419
  
 
420
  set JetPTMin 20.0
 
421
}
 
422
 
 
423
###########
 
424
# b-tagging
 
425
###########
 
426
 
 
427
module BTagging BTagging {
 
428
  set PartonInputArray Delphes/partons
 
429
  set JetInputArray FastJetFinder/jets
 
430
 
 
431
  set DeltaR 0.5
 
432
 
 
433
  set PartonPTMin 1.0
 
434
 
 
435
  set PartonEtaMax 2.5
 
436
 
 
437
  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
 
438
  # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
 
439
  # gluon's PDG code has the lowest priority
 
440
 
 
441
  # default efficiency formula (misidentification rate)
 
442
  add EfficiencyFormula {0} {0.001}
 
443
  # efficiency formula for c-jets (misidentification rate)
 
444
  add EfficiencyFormula {4} {0.1}
 
445
  # efficiency formula for b-jets
 
446
  add EfficiencyFormula {5} {0.4}
 
447
}
 
448
 
 
449
module TauTagging TauTagging {
 
450
  set ParticleInputArray Delphes/allParticles
 
451
  set PartonInputArray Delphes/partons
 
452
  set JetInputArray FastJetFinder/jets
 
453
 
 
454
  set DeltaR 0.5
 
455
 
 
456
  set TauPTMin 1.0
 
457
 
 
458
  set TauEtaMax 2.5
 
459
 
 
460
  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
 
461
 
 
462
  # default efficiency formula (misidentification rate)
 
463
  add EfficiencyFormula {0} {0.001}
 
464
  # efficiency formula for tau-jets
 
465
  add EfficiencyFormula {15} {0.4}
 
466
}
 
467
 
 
468
#####################################################
 
469
# Find uniquely identified photons/electrons/tau/jets
 
470
#####################################################
 
471
 
 
472
module UniqueObjectFinder UniqueObjectFinder {
 
473
# earlier arrays take precedence over later ones
 
474
# add InputArray InputArray OutputArray
 
475
  add InputArray PhotonIsolation/photons photons
 
476
  add InputArray ElectronIsolation/electrons electrons
 
477
  add InputArray FastJetFinder/jets jets
 
478
}
 
479
 
 
480
##################
 
481
# ROOT tree writer
 
482
##################
 
483
 
 
484
module TreeWriter TreeWriter {
 
485
# add Branch InputArray BranchName BranchClass
 
486
  add Branch Delphes/allParticles Particle GenParticle
 
487
  add Branch TrackMerger/tracks Track Track
 
488
  add Branch Calorimeter/towers Tower Tower
 
489
  add Branch Calorimeter/eflowTracks EFlowTrack Track
 
490
  add Branch Calorimeter/eflowTowers EFlowTower Tower
 
491
  add Branch MuonMomentumSmearing/muons EFlowMuon Muon
 
492
  add Branch UniqueObjectFinder/jets Jet Jet
 
493
  add Branch UniqueObjectFinder/electrons Electron Electron
 
494
  add Branch UniqueObjectFinder/photons Photon Photon
 
495
  add Branch MuonIsolation/muons Muon Muon
 
496
  add Branch MissingET/momentum MissingET MissingET
 
497
  add Branch ScalarHT/energy ScalarHT ScalarHT
 
498
}
135
499