Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_HLLHC.tcl@ 001ee95

Visit:

ImprovedOutputFile Timing dual_readout llp

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

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

Download in other formats: