Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_HLLHC.tcl@ 8a58fff

Visit:

ImprovedOutputFile Timing dual_readout llp

Last change on this file since 8a58fff was 8a58fff, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago
updated HLLHC with new DTF
Property mode set to `100644`
File size: 58.9 KB

Rev	Line
[8562e4e]	1	#
[fce5dd5]	2	# Beta card for HL-LHC and HE-LHC studies
[8562e4e]	3	#
	4	# Main authors: Michele Selvaggi (CERN)
	5	#
[fce5dd5]	6	# Released on: Dec. 6th, 2017
[8562e4e]	7	#
	8	#
	9	#######################################
	10	# Order of execution of various modules
	11	#######################################
	12
	13	set ExecutionPath {
	14
	15	ParticlePropagator
[8a58fff]	16	TrackMergerProp
	17
	18	DenseProp
	19	DenseMergeTracks
	20	DenseTrackFilter
[8562e4e]	21
	22	ChargedHadronTrackingEfficiency
	23	ElectronTrackingEfficiency
	24	MuonTrackingEfficiency
	25
	26	ChargedHadronMomentumSmearing
	27	ElectronMomentumSmearing
	28	MuonMomentumSmearing
	29
	30	TrackMerger
	31
	32	ECal
	33	HCal
	34
	35	Calorimeter
	36	EFlowMerger
	37	EFlowFilter
	38
	39	PhotonIsolation
	40	PhotonEfficiency
	41
	42	ElectronFilter
	43	ElectronIsolation
	44	ElectronEfficiency
	45
	46	ChargedHadronFilter
	47
	48	MuonIsolation
	49	MuonEfficiency
	50
	51	MissingET
	52
	53	NeutrinoFilter
	54	GenJetFinder
	55	GenMissingET
	56
	57	FastJetFinder
	58	FatJetFinder
	59
	60	JetEnergyScale
	61
	62	JetFlavorAssociation
	63
	64	BTagging
	65	TauTagging
	66
	67	UniqueObjectFinder
	68
	69	ScalarHT
	70
	71	TreeWriter
	72	}
	73
	74
[8a58fff]	75	#####################################
	76	# Track propagation to calorimeters
	77	#####################################
[8562e4e]	78
	79	module ParticlePropagator ParticlePropagator {
	80	set InputArray Delphes/stableParticles
	81
	82	set OutputArray stableParticles
	83	set ChargedHadronOutputArray chargedHadrons
	84	set ElectronOutputArray electrons
	85	set MuonOutputArray muons
	86
	87	# radius of the magnetic field coverage, in m
	88	set Radius 1.2
	89	# half-length of the magnetic field coverage, in m
	90	set HalfLength 3.25
	91
	92	# magnetic field
	93	set Bz 3.0
	94	}
	95
[8a58fff]	96	##############
	97	# Track merger
	98	##############
	99
	100	module Merger TrackMergerProp {
	101	# add InputArray InputArray
	102	add InputArray ParticlePropagator/chargedHadrons
	103	add InputArray ParticlePropagator/electrons
	104	add InputArray ParticlePropagator/muons
	105	set OutputArray tracks
	106	}
	107
	108	####################################
	109	# Track propagation to pseudo-pixel
	110	####################################
	111
	112	module ParticlePropagator DenseProp {
	113
	114	set InputArray TrackMergerProp/tracks
	115
	116	# radius of the magnetic field coverage, in m
	117	set Radius 0.3
	118	set RadiusMax 1.2
	119	# half-length of the magnetic field coverage, in m
	120	set HalfLength 0.7
	121	set HalfLengthMax 3.25
	122
	123	# magnetic field
	124	set Bz 3.0
	125	}
	126
	127	#####################
	128	# Dense Track merger
	129	#####################
	130
	131	module Merger DenseMergeTracks {
	132	# add InputArray InputArray
	133	add InputArray DenseProp/chargedHadrons
	134	add InputArray DenseProp/electrons
	135	add InputArray DenseProp/muons
	136	set OutputArray tracks
	137	}
	138
	139
	140	######################
	141	# Dense Track Filter
	142	######################
	143
	144	module DenseTrackFilter DenseTrackFilter {
	145
	146	set TrackInputArray DenseMergeTracks/tracks
	147
	148	set TrackOutputArray tracks
	149	set ChargedHadronOutputArray chargedHadrons
	150	set ElectronOutputArray electrons
	151	set MuonOutputArray muons
	152
	153	set EtaPhiRes 0.003
	154	set EtaMax 4.0
	155
	156	set pi [expr {acos(-1)}]
	157
	158	set nbins_phi [expr {$pi/$EtaPhiRes} ]
	159	set nbins_phi [expr {int($nbins_phi)} ]
	160
	161	set PhiBins {}
	162	for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
	163	add PhiBins [expr {$i * $pi/$nbins_phi}]
	164	}
	165
	166	set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
	167	set nbins_eta [expr {int($nbins_eta)} ]
	168
	169	for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
	170	set eta [expr {$i * $EtaPhiRes}]
	171	add EtaPhiBins $eta $PhiBins
	172	}
	173	}
	174
	175
	176
[8562e4e]	177	####################################
	178	# Charged hadron tracking efficiency
	179	####################################
	180
	181	module Efficiency ChargedHadronTrackingEfficiency {
[8a58fff]	182	set InputArray DenseTrackFilter/chargedHadrons
[8562e4e]	183	set OutputArray chargedHadrons
	184
	185	# TBC (which eta_max ? which pT min?)
	186
	187	# tracking efficiency formula for charged hadrons
	188
	189	set EfficiencyFormula {
	190	(pt <= 0.5) * (0.00) + \
	191	(abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
	192	(abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
	193	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
	194	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90) + \
	195	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
	196	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85) + \
	197	(abs(eta) > 4.0) * (0.00)
	198	}
	199
	200	}
	201
	202	##############################
	203	# Electron tracking efficiency
	204	##############################
	205
	206	module Efficiency ElectronTrackingEfficiency {
[8a58fff]	207	set InputArray DenseTrackFilter/electrons
[8562e4e]	208	set OutputArray electrons
	209
	210	# TBC (which eta_max ?)
	211	# putting same as charged hadrons for now...
	212
	213	set EfficiencyFormula {
	214	(pt <= 0.5) * (0.00) + \
	215	(abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
	216	(abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
	217	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
[b2ab97c]	218	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90)/(-0.033*log10(pt) + 1.033) + \
[8562e4e]	219	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
[b2ab97c]	220	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85)/(-0.033*log10(pt) + 1.033) + \
[8562e4e]	221	(abs(eta) > 4.0) * (0.00)
	222	}
	223
	224	}
	225	##########################
	226	# Muon tracking efficiency
	227	##########################
	228
	229	module Efficiency MuonTrackingEfficiency {
[8a58fff]	230	set InputArray DenseTrackFilter/muons
[8562e4e]	231	set OutputArray muons
	232
	233	set EfficiencyFormula {
	234	(pt <= 0.5) * (0.00) + \
	235	(abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 1.00) + \
	236	(abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
	237	(abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.5 && pt <= 1.0) * (pt*1.00) + \
	238	(abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
	239	(abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.95) + \
	240	(abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
	241	(abs(eta) > 4.0) * (0.00)
	242
	243	}
	244
	245	}
	246
	247	########################################
	248	# Momentum resolution for charged tracks
	249	########################################
	250
	251	module MomentumSmearing ChargedHadronMomentumSmearing {
	252	set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
	253	set OutputArray chargedHadrons
	254
	255	# --- CMS resolution Phase II ---
	256
[72c141f]	257	set ResolutionFormula { 2((abs(eta) >= 0.0000 && abs(eta) < 0.2000) (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
[8562e4e]	258	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
	259	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
	260	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
	261	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
	262	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
	263	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
	264	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
	265	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
	266	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
	267	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
	268	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
	269	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
	270	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
	271	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
	272	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
	273	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
	274	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
	275	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
	276	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
	277	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
	278	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
	279	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
	280	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
	281	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
	282	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
	283	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
	284	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
	285	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
	286	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
	287	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
	288	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
	289	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
	290	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
	291	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
	292	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
	293	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
	294	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
	295	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
	296	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
	297	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
	298	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
	299	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
	300	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
	301	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
	302	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
	303	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
	304	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
	305	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
	306	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
	307	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
	308	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
	309	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
	310	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
	311	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
	312	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
	313	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
	314	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
	315	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
	316	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
	317	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
	318	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
	319	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
	320	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
	321	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
	322	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
	323	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
	324	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
	325	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
	326	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
	327	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
	328	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
	329	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
	330	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
	331	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
	332	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
	333	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
	334	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
	335	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
	336	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
	337	}
	338
	339	}
	340
	341
	342	###################################
	343	# Momentum resolution for electrons
	344	###################################
	345
	346	module MomentumSmearing ElectronMomentumSmearing {
	347	set InputArray ElectronTrackingEfficiency/electrons
	348	set OutputArray electrons
	349
[72c141f]	350	set ResolutionFormula { 2((abs(eta) >= 0.0000 && abs(eta) < 0.2000) (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
[8562e4e]	351	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
	352	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
	353	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
	354	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
	355	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
	356	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
	357	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
	358	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
	359	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
	360	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
	361	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
	362	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
	363	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
	364	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
	365	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
	366	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
	367	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
	368	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
	369	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
	370	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
	371	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
	372	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
	373	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
	374	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
	375	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
	376	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
	377	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
	378	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
	379	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
	380	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
	381	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
	382	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
	383	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
	384	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
	385	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
	386	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
	387	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
	388	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
	389	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
	390	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
	391	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
	392	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
	393	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
	394	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
	395	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
	396	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
	397	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
	398	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
	399	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
	400	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
	401	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
	402	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
	403	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
	404	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
	405	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
	406	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
	407	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
	408	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
	409	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
	410	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
	411	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
	412	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
	413	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
	414	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
	415	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
	416	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
	417	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
	418	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
	419	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
	420	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
	421	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
	422	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
	423	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
	424	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
	425	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
	426	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
	427	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
	428	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
	429	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
	430	}
	431	}
	432
	433
	434	###############################
	435	# Momentum resolution for muons
	436	###############################
	437
	438	module MomentumSmearing MuonMomentumSmearing {
	439	set InputArray MuonTrackingEfficiency/muons
	440	set OutputArray muons
	441
	442	# Equal to Tk resolution up to pT = 100 GeV, then tk slope/5.
	443
	444	set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
	445	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
	446	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
[97f1b2d]	447	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244/5.pt/100.000000 + 4./5.0.009244) + \
[8562e4e]	448	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
	449	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
	450	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
[97f1b2d]	451	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172/5.pt/100.000000 + 4./5.0.009172) + \
[8562e4e]	452	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
	453	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
	454	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
[97f1b2d]	455	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077/5.pt/100.000000 + 4./5.0.009077) + \
[8562e4e]	456	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
	457	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
	458	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
[97f1b2d]	459	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177/5.pt/100.000000 + 4./5.0.009177) + \
[8562e4e]	460	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
	461	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
	462	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
[97f1b2d]	463	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429/5.pt/100.000000 + 4./5.0.010429) + \
[8562e4e]	464	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
	465	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
	466	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
[97f1b2d]	467	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867/5.pt/100.000000 + 4./5.0.013867) + \
[8562e4e]	468	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
	469	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
	470	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
[97f1b2d]	471	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800/5.pt/100.000000 + 4./5.0.016800) + \
[8562e4e]	472	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
	473	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
	474	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
[97f1b2d]	475	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161/5.pt/100.000000 + 4./5.0.018161) + \
[8562e4e]	476	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
	477	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
	478	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
[97f1b2d]	479	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039/5.pt/100.000000 + 4./5.0.018039) + \
[8562e4e]	480	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
	481	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
	482	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
[97f1b2d]	483	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098/5.pt/100.000000 + 4./5.0.020098) + \
[8562e4e]	484	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
	485	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
	486	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
[97f1b2d]	487	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289/5.pt/100.000000 + 4./5.0.020289) + \
[8562e4e]	488	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
	489	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
	490	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
[97f1b2d]	491	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524/5.pt/100.000000 + 4./5.0.021524) + \
[8562e4e]	492	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
	493	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
	494	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
[97f1b2d]	495	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477/5.pt/100.000000 + 4./5.0.024477) + \
[8562e4e]	496	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
	497	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
	498	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
[97f1b2d]	499	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392/5.pt/100.000000 + 4./5.0.030392) + \
[8562e4e]	500	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
	501	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
	502	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
[97f1b2d]	503	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445/5.pt/100.000000 + 4./5.0.027445) + \
[8562e4e]	504	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
	505	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
	506	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
[97f1b2d]	507	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053/5.pt/100.000000 + 4./5.0.025053) + \
[8562e4e]	508	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
	509	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
	510	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
[97f1b2d]	511	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710/5.pt/100.000000 + 4./5.0.030710) + \
[8562e4e]	512	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
	513	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
	514	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
[97f1b2d]	515	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629/5.pt/100.000000 + 4./5.0.040629) + \
[8562e4e]	516	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
	517	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
	518	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
[97f1b2d]	519	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345/5.pt/100.000000 + 4./5.0.052345) + \
[8562e4e]	520	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
	521	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
	522	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
[97f1b2d]	523	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617/5.pt/100.000000 + 4./5.0.069617))
[8562e4e]	524	}
	525
	526	}
	527
[8a58fff]	528
	529
[8562e4e]	530	##############
	531	# Track merger
	532	##############
	533
[8a58fff]	534	module Merger TrackMerger {
[8562e4e]	535	# add InputArray InputArray
	536	add InputArray ChargedHadronMomentumSmearing/chargedHadrons
	537	add InputArray ElectronMomentumSmearing/electrons
	538	add InputArray MuonMomentumSmearing/muons
	539	set OutputArray tracks
	540	}
	541
	542
	543	#############
	544	# ECAL
	545	#############
	546
	547	module SimpleCalorimeter ECal {
	548	set ParticleInputArray ParticlePropagator/stableParticles
	549	set TrackInputArray TrackMerger/tracks
	550
	551	set TowerOutputArray ecalTowers
	552	set EFlowTrackOutputArray eflowTracks
	553	set EFlowTowerOutputArray eflowPhotons
	554
	555	set IsEcal true
	556
	557	set EnergyMin 0.5
	558	set EnergySignificanceMin 1.0
	559
	560	set SmearTowerCenter true
	561
	562	set pi [expr {acos(-1)}]
	563
	564	# lists of the edges of each tower in eta and phi
	565	# each list starts with the lower edge of the first tower
	566	# the list ends with the higher edged of the last tower
	567
	568	# assume 0.02 x 0.02 resolution in eta,phi in the barrel \|eta\| < 1.5
	569
	570	set PhiBins {}
	571	for {set i -180} {$i <= 180} {incr i} {
	572	add PhiBins [expr {$i * $pi/180.0}]
	573	}
	574
	575	# 0.02 unit in eta up to eta = 1.5 (barrel)
	576	for {set i -85} {$i <= 86} {incr i} {
	577	set eta [expr {$i * 0.0174}]
	578	add EtaPhiBins $eta $PhiBins
	579	}
	580
	581	# assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < \|eta\| < 3.0 (HGCAL- ECAL)
	582
	583	set PhiBins {}
	584	for {set i -180} {$i <= 180} {incr i} {
	585	add PhiBins [expr {$i * $pi/180.0}]
	586	}
	587
	588	# 0.02 unit in eta up to eta = 3
	589	for {set i 1} {$i <= 84} {incr i} {
	590	set eta [expr { -2.958 + $i * 0.0174}]
	591	add EtaPhiBins $eta $PhiBins
	592	}
	593
	594	for {set i 1} {$i <= 84} {incr i} {
	595	set eta [expr { 1.4964 + $i * 0.0174}]
	596	add EtaPhiBins $eta $PhiBins
	597	}
	598
	599	# take present CMS granularity for HF
	600
	601	# 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < \|eta\| < 5.0
	602	set PhiBins {}
	603	for {set i -18} {$i <= 18} {incr i} {
	604	add PhiBins [expr {$i * $pi/18.0}]
	605	}
	606
	607	foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
	608	add EtaPhiBins $eta $PhiBins
	609	}
	610
	611
	612	add EnergyFraction {0} {0.0}
	613	# energy fractions for e, gamma and pi0
	614	add EnergyFraction {11} {1.0}
	615	add EnergyFraction {22} {1.0}
	616	add EnergyFraction {111} {1.0}
	617	# energy fractions for muon, neutrinos and neutralinos
	618	add EnergyFraction {12} {0.0}
	619	add EnergyFraction {13} {0.0}
	620	add EnergyFraction {14} {0.0}
	621	add EnergyFraction {16} {0.0}
	622	add EnergyFraction {1000022} {0.0}
	623	add EnergyFraction {1000023} {0.0}
	624	add EnergyFraction {1000025} {0.0}
	625	add EnergyFraction {1000035} {0.0}
	626	add EnergyFraction {1000045} {0.0}
	627	# energy fractions for K0short and Lambda
	628	add EnergyFraction {310} {0.3}
	629	add EnergyFraction {3122} {0.3}
	630
	631
	632	# ECAL is CMS
	633	# for the ECAL barrel (\|eta\| < 1.5), see hep-ex/1306.2016 and 1502.02701
	634	# for the endcaps (1.5 < \|eta\| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
	635
[72c141f]	636	set ResolutionFormula { 2((abs(eta) <= 1.50) sqrt(energy^20.009^2 + energy0.12^2 + 0.45^2) +
[8562e4e]	637	(abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^20.006^2 + energy0.20^2) + \
	638	(abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^20.007^2 + energy0.21^2) + \
[72c141f]	639	(abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^20.008^2 + energy0.24^2)) + \
[8562e4e]	640	(abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^20.08^2 + energy1.98^2)}
	641
	642
	643	}
	644
	645	#############
	646	# HCAL
	647	#############
	648	module SimpleCalorimeter HCal {
	649	set ParticleInputArray ParticlePropagator/stableParticles
	650	set TrackInputArray ECal/eflowTracks
	651
	652	set TowerOutputArray hcalTowers
	653	set EFlowTrackOutputArray eflowTracks
	654	set EFlowTowerOutputArray eflowNeutralHadrons
	655
	656	set IsEcal false
	657
	658	set EnergyMin 1.0
	659	set EnergySignificanceMin 2.0
	660
	661	set SmearTowerCenter true
	662
	663	set pi [expr {acos(-1)}]
	664
	665	# lists of the edges of each tower in eta and phi
	666	# each list starts with the lower edge of the first tower
	667	# the list ends with the higher edged of the last tower
	668
	669	# 10 degrees towers
	670	set PhiBins {}
	671	for {set i -18} {$i <= 18} {incr i} {
	672	add PhiBins [expr {$i * $pi/18.0}]
	673	}
	674	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} {
	675	add EtaPhiBins $eta $PhiBins
	676	}
	677
	678	# 20 degrees towers
	679	set PhiBins {}
	680	for {set i -9} {$i <= 9} {incr i} {
	681	add PhiBins [expr {$i * $pi/9.0}]
	682	}
	683	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.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
	684	add EtaPhiBins $eta $PhiBins
	685	}
	686
	687	# default energy fractions {abs(PDG code)} {Fecal Fhcal}
	688	add EnergyFraction {0} {1.0}
	689	# energy fractions for e, gamma and pi0
	690	add EnergyFraction {11} {0.0}
	691	add EnergyFraction {22} {0.0}
	692	add EnergyFraction {111} {0.0}
	693	# energy fractions for muon, neutrinos and neutralinos
	694	add EnergyFraction {12} {0.0}
	695	add EnergyFraction {13} {0.0}
	696	add EnergyFraction {14} {0.0}
	697	add EnergyFraction {16} {0.0}
	698	add EnergyFraction {1000022} {0.0}
	699	add EnergyFraction {1000023} {0.0}
	700	add EnergyFraction {1000025} {0.0}
	701	add EnergyFraction {1000035} {0.0}
	702	add EnergyFraction {1000045} {0.0}
	703	# energy fractions for K0short and Lambda
	704	add EnergyFraction {310} {0.7}
	705	add EnergyFraction {3122} {0.7}
	706
	707	# HCAL is ALTAS
	708	# http://arxiv.org/pdf/hep-ex/0004009v1
	709	# http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
	710	# set HCalResolutionFormula {resolution formula as a function of eta and energy}
[72c141f]	711	set ResolutionFormula { 2((abs(eta) <= 1.7) sqrt(energy^20.0302^2 + energy0.5205^2 + 1.59^2) +
	712	(abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^20.0500^2 + energy0.706^2)) +
[8562e4e]	713	(abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^20.09420^2 + energy1.00^2)}
	714	}
	715
	716
	717	#################
	718	# Electron filter
	719	#################
	720
	721	module PdgCodeFilter ElectronFilter {
	722	set InputArray HCal/eflowTracks
	723	set OutputArray electrons
	724	set Invert true
	725	add PdgCode {11}
	726	add PdgCode {-11}
	727	}
	728
	729
	730	######################
	731	# ChargedHadronFilter
	732	######################
	733
	734	module PdgCodeFilter ChargedHadronFilter {
	735	set InputArray HCal/eflowTracks
	736	set OutputArray chargedHadrons
	737
	738	add PdgCode {11}
	739	add PdgCode {-11}
	740	add PdgCode {13}
	741	add PdgCode {-13}
	742	}
	743
	744
	745	###################################################
	746	# Tower Merger (in case not using e-flow algorithm)
	747	###################################################
	748
	749	module Merger Calorimeter {
	750	# add InputArray InputArray
	751	add InputArray ECal/ecalTowers
	752	add InputArray HCal/hcalTowers
	753	set OutputArray towers
	754	}
	755
	756	####################
	757	# Energy flow merger
	758	####################
	759
	760	module Merger EFlowMerger {
	761	# add InputArray InputArray
	762	add InputArray HCal/eflowTracks
	763	add InputArray ECal/eflowPhotons
	764	add InputArray HCal/eflowNeutralHadrons
	765	set OutputArray eflow
	766	}
	767
	768	######################
	769	# EFlowFilter
	770	######################
	771
	772	module PdgCodeFilter EFlowFilter {
	773	set InputArray EFlowMerger/eflow
	774	set OutputArray eflow
	775
	776	add PdgCode {11}
	777	add PdgCode {-11}
	778	add PdgCode {13}
	779	add PdgCode {-13}
	780	}
	781
	782
	783	###################
	784	# Missing ET merger
	785	###################
	786
	787	module Merger MissingET {
	788	# add InputArray InputArray
	789	add InputArray EFlowMerger/eflow
	790	set MomentumOutputArray momentum
	791	}
	792
	793
	794
	795	##################
	796	# Scalar HT merger
	797	##################
	798
	799	module Merger ScalarHT {
	800	# add InputArray InputArray
	801	add InputArray EFlowMerger/eflow
	802	set EnergyOutputArray energy
	803	}
	804
	805	#################
	806	# Neutrino Filter
	807	#################
	808
	809	module PdgCodeFilter NeutrinoFilter {
	810
	811	set InputArray Delphes/stableParticles
	812	set OutputArray filteredParticles
	813
	814	set PTMin 0.0
	815
	816	add PdgCode {12}
	817	add PdgCode {14}
	818	add PdgCode {16}
	819	add PdgCode {-12}
	820	add PdgCode {-14}
	821	add PdgCode {-16}
	822
	823	}
	824
	825
	826	#####################
	827	# MC truth jet finder
	828	#####################
	829
	830	module FastJetFinder GenJetFinder {
	831	set InputArray NeutrinoFilter/filteredParticles
	832
	833	set OutputArray jets
	834
	835	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	836	set JetAlgorithm 6
	837	set ParameterR 0.4
	838
	839	set JetPTMin 15.0
	840	}
	841
	842	#########################
	843	# Gen Missing ET merger
	844	########################
	845
	846	module Merger GenMissingET {
	847
	848	# add InputArray InputArray
	849	add InputArray NeutrinoFilter/filteredParticles
	850	set MomentumOutputArray momentum
	851	}
	852
	853
	854	############
	855	# Jet finder
	856	############
	857
	858	module FastJetFinder FastJetFinder {
[08fe3d2]	859	set InputArray Calorimeter/towers
	860	# set InputArray EFlowMerger/eflow
[8562e4e]	861
	862	set OutputArray jets
	863
	864	set JetAlgorithm 6
	865	set ParameterR 0.4
	866
	867	set ComputeNsubjettiness 1
	868	set Beta 1.0
	869	set AxisMode 4
	870
	871	set ComputeSoftDrop 1
	872	set BetaSoftDrop 0.0
	873	set SymmetryCutSoftDrop 0.1
	874	set R0SoftDrop 0.4
	875
	876	set JetPTMin 30.0
	877	}
	878
	879	##################
	880	# Fat Jet finder
	881	##################
	882
	883	module FastJetFinder FatJetFinder {
	884	set InputArray EFlowMerger/eflow
	885
	886	set OutputArray jets
	887
	888	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	889	set JetAlgorithm 6
	890	set ParameterR 0.8
	891
	892	set ComputeNsubjettiness 1
	893	set Beta 1.0
	894	set AxisMode 4
	895
	896	set ComputeSoftDrop 1
	897	set BetaSoftDrop 0.0
	898	set SymmetryCutSoftDrop 0.1
	899	set R0SoftDrop 0.8
	900
[08fe3d2]	901	set JetPTMin 200.0
[8562e4e]	902	}
	903
	904
	905
	906	##################
	907	# Jet Energy Scale
	908	##################
	909
	910	module EnergyScale JetEnergyScale {
	911	set InputArray FastJetFinder/jets
	912	set OutputArray jets
	913
	914	# scale formula for jets
	915	set ScaleFormula {1.00}
	916	}
	917
	918
	919	########################
	920	# Jet Flavor Association
	921	########################
	922
	923	module JetFlavorAssociation JetFlavorAssociation {
	924
	925	set PartonInputArray Delphes/partons
	926	set ParticleInputArray Delphes/allParticles
	927	set ParticleLHEFInputArray Delphes/allParticlesLHEF
	928	set JetInputArray JetEnergyScale/jets
	929
	930	set DeltaR 0.5
	931	set PartonPTMin 5.0
	932	set PartonEtaMax 4.0
	933
	934	}
	935
	936	##################
	937	# Photon isolation
	938	##################
	939
	940	# TBC: check values for iso cuts
	941
	942	module Isolation PhotonIsolation {
	943	set CandidateInputArray ECal/eflowPhotons
	944	set IsolationInputArray EFlowFilter/eflow
	945
	946	set OutputArray photons
	947
	948	set DeltaRMax 0.3
	949	set PTMin 0.0
	950	set PTRatioMax 0.1
	951	}
	952
	953
	954
	955	###################
	956	# Photon efficiency
	957	###################
	958
	959	module Efficiency PhotonEfficiency {
	960	set InputArray PhotonIsolation/photons
	961	set OutputArray photons
	962
	963	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	964
	965	set EfficiencyFormula {
	966	(pt <= 1.0) * (0.00) + \
	967	(abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
	968	(abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
	969	(abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
	970
	971	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
	972	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
	973	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
	974
	975	(abs(eta) > 4.0) * (0.00)}
	976
	977	}
	978
	979
	980
	981	####################
	982	# Electron isolation
	983	####################
	984
	985	# TBC: check values for iso cuts
	986
	987	module Isolation ElectronIsolation {
	988	set CandidateInputArray ElectronFilter/electrons
	989	set IsolationInputArray EFlowFilter/eflow
	990
	991	set OutputArray electrons
	992
	993	set DeltaRMax 0.3
	994
	995	set PTMin 0.0
	996
	997	set PTRatioMax 0.1
	998	}
	999
	1000
	1001	#######################
	1002	# Electron efficiency #
	1003	#######################
	1004
	1005	module Efficiency ElectronEfficiency {
	1006
	1007	set InputArray ElectronIsolation/electrons
	1008	set OutputArray electrons
	1009
	1010	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	1011	# efficiency formula for electrons (taken from CMS PhaseII card)
	1012	# efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
	1013	set EfficiencyFormula {
	1014	(pt <= 2.0) * (0.00) +
	1015	(abs(eta) <= 1.45 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1016	(abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
	1017	(abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
	1018	(abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
	1019	(abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
	1020	(abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
	1021	(abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
	1022	(abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
	1023	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1024	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
	1025	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
	1026	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
	1027	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
	1028	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1029	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
	1030	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
	1031	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
	1032	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
	1033	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
	1034	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1035	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	1036	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	1037	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	1038	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	1039	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
	1040	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1041	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1042	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	1043	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	1044	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	1045	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	1046	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
	1047	(abs(eta) > 4.0) * (0.00)
	1048
	1049	}
	1050
	1051	}
	1052
	1053
	1054
	1055	################
	1056	# Muon isolation
	1057	################
	1058
	1059	# TBC: check values for iso cuts
	1060
	1061	module Isolation MuonIsolation {
	1062	set CandidateInputArray MuonMomentumSmearing/muons
	1063	set IsolationInputArray EFlowFilter/eflow
	1064
	1065	set OutputArray muons
	1066
	1067	set DeltaRMax 0.3
	1068
	1069	set PTMin 0.0
	1070
	1071	set PTRatioMax 0.2
	1072	}
	1073
	1074	####################
	1075	# Muon Efficiency #
	1076	####################
	1077
	1078	module Efficiency MuonEfficiency {
	1079	set InputArray MuonIsolation/muons
	1080	set OutputArray muons
	1081	# tracking + TightID efficiency formula for muons
	1082	set EfficiencyFormula {
	1083	(pt <= 2.0) * (0.00) +
	1084	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
	1085	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
	1086	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
	1087	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.68) +
	1088	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
	1089	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.91) +
	1090	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
	1091	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1092	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
	1093	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.47) +
	1094	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
	1095	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.69) +
	1096	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
	1097	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.93) +
	1098	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.94) +
	1099	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1100	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.16) +
	1101	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.48) +
	1102	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
	1103	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.66) +
	1104	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.79) +
	1105	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
	1106	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.95) +
	1107	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1108	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.24) +
	1109	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
	1110	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.51) +
	1111	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.71) +
	1112	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.77) +
	1113	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
	1114	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.92) +
	1115	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1116	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.23) +
	1117	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.35) +
	1118	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.43) +
	1119	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
	1120	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
	1121	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
	1122	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
	1123	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.81) +
	1124	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1125	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	1126	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	1127	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	1128	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	1129	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
	1130	(abs(eta) > 4.0) * (0.00)
	1131	}
	1132	}
	1133
	1134
	1135
	1136	###########
	1137	# b-tagging
	1138	###########
	1139
	1140	module BTagging BTagging {
	1141	set JetInputArray JetEnergyScale/jets
	1142
	1143	set BitNumber 0
	1144
	1145	add EfficiencyFormula {0} {
	1146
	1147	(pt <= 10.0) * (0.00) +
	1148	(abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
	1149	(abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.01)*(1.0 - pt/5000.) + \
	1150	(abs(eta) < 2.5) * (pt > 5000.0) * (0.00) + \
	1151	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
	1152	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.0075)*(1.0 - pt/5000.) + \
	1153	(abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
	1154	(abs(eta) > 4.0) * (0.00)}
	1155
	1156	add EfficiencyFormula {4} {
	1157
	1158	(pt <= 10.0) * (0.00) +
	1159	(abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.10) + \
	1160	(abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.10)*(1.0 - pt/5000.) + \
	1161	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1162	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.06) + \
	1163	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.06)*(1.0 - pt/5000.) + \
	1164	(abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
	1165	(abs(eta) > 4.0) * (0.00)}
	1166
	1167	add EfficiencyFormula {5} {
	1168
	1169	(pt <= 10.0) * (0.00) +
	1170	(abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.75) +
	1171	(abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.75)*(1.0 - pt/5000.) +
	1172	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) +
	1173	(abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.60) +
	1174	(abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.60)*(1.0 - pt/5000.) +
	1175	(abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) +
	1176	(abs(eta) >= 4.0) * (0.00)}
	1177
	1178	}
	1179
	1180
	1181	#############
	1182	# tau-tagging
	1183	#############
	1184
	1185
	1186	module TauTagging TauTagging {
	1187	set ParticleInputArray Delphes/allParticles
	1188	set PartonInputArray Delphes/partons
	1189	set JetInputArray JetEnergyScale/jets
	1190
	1191	set DeltaR 0.5
	1192
	1193	set TauPTMin 1.0
	1194
	1195	set TauEtaMax 4.0
	1196
	1197	# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
	1198	add EfficiencyFormula {0} {
	1199
	1200	(pt <= 10.0) * (0.00) +
	1201	(abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.01) + \
	1202	(abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.01) *(1. - pt/5000.) + \
	1203	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1204	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.0075) + \
	1205	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.0075)*(1. - pt/5000.) + \
	1206	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
	1207	(abs(eta) > 4.0) * (0.00)}
	1208
	1209	add EfficiencyFormula {11} {
	1210
	1211	(pt <= 10.0) * (0.00) +
	1212	(abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.005) + \
	1213	(abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.005) *(1. - pt/5000.) + \
	1214	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1215	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.00375) + \
	1216	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.00375)*(1. - pt/5000.) + \
	1217	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
	1218	(abs(eta) > 4.0) * (0.00)}
	1219
	1220	add EfficiencyFormula {15} {
	1221
	1222	(pt <= 10.0) * (0.00) +
	1223	(abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.6) + \
	1224	(abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.6) *(1. - pt/5000.) + \
	1225	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1226	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.45) + \
	1227	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.45)*(1. - pt/5000.) + \
	1228	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
	1229	(abs(eta) > 4.0) * (0.00)}
	1230
	1231	}
	1232
	1233	#####################################################
	1234	# Find uniquely identified photons/electrons/tau/jets
	1235	#####################################################
	1236
	1237	module UniqueObjectFinder UniqueObjectFinder {
	1238	# earlier arrays take precedence over later ones
	1239	# add InputArray InputArray OutputArray
	1240	add InputArray PhotonEfficiency/photons photons
	1241	add InputArray ElectronEfficiency/electrons electrons
	1242	add InputArray MuonEfficiency/muons muons
	1243	add InputArray JetEnergyScale/jets jets
	1244
	1245	}
	1246
	1247
	1248	##################
	1249	# ROOT tree writer
	1250	##################
	1251
	1252	module TreeWriter TreeWriter {
	1253	# add Branch InputArray BranchName BranchClass
	1254	add Branch Delphes/allParticles Particle GenParticle
	1255
	1256	add Branch GenJetFinder/jets GenJet Jet
	1257	add Branch GenMissingET/momentum GenMissingET MissingET
	1258
	1259	add Branch TrackMerger/tracks Track Track
	1260	add Branch Calorimeter/towers Tower Tower
	1261
	1262	add Branch HCal/eflowTracks EFlowTrack Track
	1263	add Branch ECal/eflowPhotons EFlowPhoton Tower
	1264	add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
	1265
	1266	add Branch UniqueObjectFinder/photons Photon Photon
	1267	add Branch UniqueObjectFinder/electrons Electron Electron
	1268	add Branch UniqueObjectFinder/muons Muon Muon
	1269	add Branch UniqueObjectFinder/jets Jet Jet
	1270
	1271	add Branch FatJetFinder/jets FatJet Jet
	1272
	1273	add Branch MissingET/momentum MissingET MissingET
	1274	add Branch ScalarHT/energy ScalarHT ScalarHT
	1275	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats: