Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_HLLHC.tcl@ 0c5bf12

Visit:

ImprovedOutputFile Timing dual_readout llp

Last change on this file since 0c5bf12 was fce5dd5, checked in by Michele Selvaggi (michele.selvaggi@…>, 7 years ago
changed header
Property mode set to `100644`
File size: 58.3 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) + \
	135	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90) + \
	136	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
	137	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85) + \
	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
	174	set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
	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
	267	set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
	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) + \
	364	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244/5.*pt/100.000000) + \
	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) + \
	368	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172/5.*pt/100.000000) + \
	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) + \
	372	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077/5.*pt/100.000000) + \
	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) + \
	376	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177/5.*pt/100.000000) + \
	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) + \
	380	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429/5.*pt/100.000000) + \
	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) + \
	384	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867/5.*pt/100.000000) + \
	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) + \
	388	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800/5.*pt/100.000000) + \
	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) + \
	392	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161/5.*pt/100.000000) + \
	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) + \
	396	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039/5.*pt/100.000000) + \
	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) + \
	400	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098/5.*pt/100.000000) + \
	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) + \
	404	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289/5.*pt/100.000000) + \
	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) + \
	408	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524/5.*pt/100.000000) + \
	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) + \
	412	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477/5.*pt/100.000000) + \
	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) + \
	416	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392/5.*pt/100.000000) + \
	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) + \
	420	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445/5.*pt/100.000000) + \
	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) + \
	424	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053/5.*pt/100.000000) + \
	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) + \
	428	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710/5.*pt/100.000000) + \
	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) + \
	432	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629/5.*pt/100.000000) + \
	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) + \
	436	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345/5.*pt/100.000000) + \
	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) + \
	440	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617/5.*pt/100.000000))
	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
	622	set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^20.009^2 + energy0.12^2 + 0.45^2) +
	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) + \
	625	(abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^20.008^2 + energy0.24^2) + \
	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}
	697	set ResolutionFormula { (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) +
	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
	737	add InputArray ECal/ecalTowers
	738	add InputArray HCal/hcalTowers
	739	set OutputArray towers
	740	}
	741
	742	####################
	743	# Energy flow merger
	744	####################
	745
	746	module Merger EFlowMerger {
	747	# add InputArray InputArray
	748	add InputArray HCal/eflowTracks
	749	add InputArray ECal/eflowPhotons
	750	add InputArray HCal/eflowNeutralHadrons
	751	set OutputArray eflow
	752	}
	753
	754	######################
	755	# EFlowFilter
	756	######################
	757
	758	module PdgCodeFilter EFlowFilter {
	759	set InputArray EFlowMerger/eflow
	760	set OutputArray eflow
	761
	762	add PdgCode {11}
	763	add PdgCode {-11}
	764	add PdgCode {13}
	765	add PdgCode {-13}
	766	}
	767
	768
	769	###################
	770	# Missing ET merger
	771	###################
	772
	773	module Merger MissingET {
	774	# add InputArray InputArray
	775	add InputArray EFlowMerger/eflow
	776	set MomentumOutputArray momentum
	777	}
	778
	779
	780
	781	##################
	782	# Scalar HT merger
	783	##################
	784
	785	module Merger ScalarHT {
	786	# add InputArray InputArray
	787	add InputArray EFlowMerger/eflow
	788	set EnergyOutputArray energy
	789	}
	790
	791	#################
	792	# Neutrino Filter
	793	#################
	794
	795	module PdgCodeFilter NeutrinoFilter {
	796
	797	set InputArray Delphes/stableParticles
	798	set OutputArray filteredParticles
	799
	800	set PTMin 0.0
	801
	802	add PdgCode {12}
	803	add PdgCode {14}
	804	add PdgCode {16}
	805	add PdgCode {-12}
	806	add PdgCode {-14}
	807	add PdgCode {-16}
	808
	809	}
	810
	811
	812	#####################
	813	# MC truth jet finder
	814	#####################
	815
	816	module FastJetFinder GenJetFinder {
	817	set InputArray NeutrinoFilter/filteredParticles
	818
	819	set OutputArray jets
	820
	821	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	822	set JetAlgorithm 6
	823	set ParameterR 0.4
	824
	825	set JetPTMin 15.0
	826	}
	827
	828	#########################
	829	# Gen Missing ET merger
	830	########################
	831
	832	module Merger GenMissingET {
	833
	834	# add InputArray InputArray
	835	add InputArray NeutrinoFilter/filteredParticles
	836	set MomentumOutputArray momentum
	837	}
	838
	839
	840	############
	841	# Jet finder
	842	############
	843
	844	module FastJetFinder FastJetFinder {
	845	# set InputArray Calorimeter/towers
	846	set InputArray EFlowMerger/eflow
	847
	848	set OutputArray jets
	849
	850	set JetAlgorithm 6
	851	set ParameterR 0.4
	852
	853	set ComputeNsubjettiness 1
	854	set Beta 1.0
	855	set AxisMode 4
	856
	857	set ComputeSoftDrop 1
	858	set BetaSoftDrop 0.0
	859	set SymmetryCutSoftDrop 0.1
	860	set R0SoftDrop 0.4
	861
	862	set JetPTMin 30.0
	863	}
	864
	865	##################
	866	# Fat Jet finder
	867	##################
	868
	869	module FastJetFinder FatJetFinder {
	870	set InputArray EFlowMerger/eflow
	871
	872	set OutputArray jets
	873
	874	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	875	set JetAlgorithm 6
	876	set ParameterR 0.8
	877
	878	set ComputeNsubjettiness 1
	879	set Beta 1.0
	880	set AxisMode 4
	881
	882	set ComputeSoftDrop 1
	883	set BetaSoftDrop 0.0
	884	set SymmetryCutSoftDrop 0.1
	885	set R0SoftDrop 0.8
	886
	887	set JetPTMin 30.0
	888	}
	889
	890
	891
	892	##################
	893	# Jet Energy Scale
	894	##################
	895
	896	module EnergyScale JetEnergyScale {
	897	set InputArray FastJetFinder/jets
	898	set OutputArray jets
	899
	900	# scale formula for jets
	901	set ScaleFormula {1.00}
	902	}
	903
	904
	905	########################
	906	# Jet Flavor Association
	907	########################
	908
	909	module JetFlavorAssociation JetFlavorAssociation {
	910
	911	set PartonInputArray Delphes/partons
	912	set ParticleInputArray Delphes/allParticles
	913	set ParticleLHEFInputArray Delphes/allParticlesLHEF
	914	set JetInputArray JetEnergyScale/jets
	915
	916	set DeltaR 0.5
	917	set PartonPTMin 5.0
	918	set PartonEtaMax 4.0
	919
	920	}
	921
	922	##################
	923	# Photon isolation
	924	##################
	925
	926	# TBC: check values for iso cuts
	927
	928	module Isolation PhotonIsolation {
	929	set CandidateInputArray ECal/eflowPhotons
	930	set IsolationInputArray EFlowFilter/eflow
	931
	932	set OutputArray photons
	933
	934	set DeltaRMax 0.3
	935	set PTMin 0.0
	936	set PTRatioMax 0.1
	937	}
	938
	939
	940
	941	###################
	942	# Photon efficiency
	943	###################
	944
	945	module Efficiency PhotonEfficiency {
	946	set InputArray PhotonIsolation/photons
	947	set OutputArray photons
	948
	949	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	950
	951	set EfficiencyFormula {
	952	(pt <= 1.0) * (0.00) + \
	953	(abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
	954	(abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
	955	(abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
	956
	957	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
	958	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
	959	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
	960
	961	(abs(eta) > 4.0) * (0.00)}
	962
	963	}
	964
	965
	966
	967	####################
	968	# Electron isolation
	969	####################
	970
	971	# TBC: check values for iso cuts
	972
	973	module Isolation ElectronIsolation {
	974	set CandidateInputArray ElectronFilter/electrons
	975	set IsolationInputArray EFlowFilter/eflow
	976
	977	set OutputArray electrons
	978
	979	set DeltaRMax 0.3
	980
	981	set PTMin 0.0
	982
	983	set PTRatioMax 0.1
	984	}
	985
	986
	987	#######################
	988	# Electron efficiency #
	989	#######################
	990
	991	module Efficiency ElectronEfficiency {
	992
	993	set InputArray ElectronIsolation/electrons
	994	set OutputArray electrons
	995
	996	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	997	# efficiency formula for electrons (taken from CMS PhaseII card)
	998	# efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
	999	set EfficiencyFormula {
	1000	(pt <= 2.0) * (0.00) +
	1001	(abs(eta) <= 1.45 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1002	(abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
	1003	(abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
	1004	(abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
	1005	(abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
	1006	(abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
	1007	(abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
	1008	(abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
	1009	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1010	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
	1011	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
	1012	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
	1013	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
	1014	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1015	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
	1016	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
	1017	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
	1018	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
	1019	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
	1020	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1021	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	1022	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	1023	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	1024	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	1025	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
	1026	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1027	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1028	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	1029	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	1030	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	1031	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	1032	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
	1033	(abs(eta) > 4.0) * (0.00)
	1034
	1035	}
	1036
	1037	}
	1038
	1039
	1040
	1041	################
	1042	# Muon isolation
	1043	################
	1044
	1045	# TBC: check values for iso cuts
	1046
	1047	module Isolation MuonIsolation {
	1048	set CandidateInputArray MuonMomentumSmearing/muons
	1049	set IsolationInputArray EFlowFilter/eflow
	1050
	1051	set OutputArray muons
	1052
	1053	set DeltaRMax 0.3
	1054
	1055	set PTMin 0.0
	1056
	1057	set PTRatioMax 0.2
	1058	}
	1059
	1060	####################
	1061	# Muon Efficiency #
	1062	####################
	1063
	1064	module Efficiency MuonEfficiency {
	1065	set InputArray MuonIsolation/muons
	1066	set OutputArray muons
	1067	# tracking + TightID efficiency formula for muons
	1068	set EfficiencyFormula {
	1069	(pt <= 2.0) * (0.00) +
	1070	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
	1071	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
	1072	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
	1073	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.68) +
	1074	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
	1075	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.91) +
	1076	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
	1077	(abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1078	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
	1079	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.47) +
	1080	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
	1081	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.69) +
	1082	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
	1083	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.93) +
	1084	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.94) +
	1085	(abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1086	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.16) +
	1087	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.48) +
	1088	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
	1089	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.66) +
	1090	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.79) +
	1091	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
	1092	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.95) +
	1093	(abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1094	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.24) +
	1095	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
	1096	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.51) +
	1097	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.71) +
	1098	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.77) +
	1099	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
	1100	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.92) +
	1101	(abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
	1102	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.23) +
	1103	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.35) +
	1104	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.43) +
	1105	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
	1106	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
	1107	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
	1108	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
	1109	(abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.81) +
	1110	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
	1111	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	1112	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	1113	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	1114	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	1115	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
	1116	(abs(eta) > 4.0) * (0.00)
	1117	}
	1118	}
	1119
	1120
	1121
	1122	###########
	1123	# b-tagging
	1124	###########
	1125
	1126	module BTagging BTagging {
	1127	set JetInputArray JetEnergyScale/jets
	1128
	1129	set BitNumber 0
	1130
	1131	add EfficiencyFormula {0} {
	1132
	1133	(pt <= 10.0) * (0.00) +
	1134	(abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
	1135	(abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.01)*(1.0 - pt/5000.) + \
	1136	(abs(eta) < 2.5) * (pt > 5000.0) * (0.00) + \
	1137	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
	1138	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.0075)*(1.0 - pt/5000.) + \
	1139	(abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
	1140	(abs(eta) > 4.0) * (0.00)}
	1141
	1142	add EfficiencyFormula {4} {
	1143
	1144	(pt <= 10.0) * (0.00) +
	1145	(abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.10) + \
	1146	(abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.10)*(1.0 - pt/5000.) + \
	1147	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1148	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.06) + \
	1149	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.06)*(1.0 - pt/5000.) + \
	1150	(abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
	1151	(abs(eta) > 4.0) * (0.00)}
	1152
	1153	add EfficiencyFormula {5} {
	1154
	1155	(pt <= 10.0) * (0.00) +
	1156	(abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.75) +
	1157	(abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.75)*(1.0 - pt/5000.) +
	1158	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) +
	1159	(abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.60) +
	1160	(abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.60)*(1.0 - pt/5000.) +
	1161	(abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) +
	1162	(abs(eta) >= 4.0) * (0.00)}
	1163
	1164	}
	1165
	1166
	1167	#############
	1168	# tau-tagging
	1169	#############
	1170
	1171
	1172	module TauTagging TauTagging {
	1173	set ParticleInputArray Delphes/allParticles
	1174	set PartonInputArray Delphes/partons
	1175	set JetInputArray JetEnergyScale/jets
	1176
	1177	set DeltaR 0.5
	1178
	1179	set TauPTMin 1.0
	1180
	1181	set TauEtaMax 4.0
	1182
	1183	# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
	1184	add EfficiencyFormula {0} {
	1185
	1186	(pt <= 10.0) * (0.00) +
	1187	(abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.01) + \
	1188	(abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.01) *(1. - pt/5000.) + \
	1189	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1190	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.0075) + \
	1191	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.0075)*(1. - pt/5000.) + \
	1192	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
	1193	(abs(eta) > 4.0) * (0.00)}
	1194
	1195	add EfficiencyFormula {11} {
	1196
	1197	(pt <= 10.0) * (0.00) +
	1198	(abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.005) + \
	1199	(abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.005) *(1. - pt/5000.) + \
	1200	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1201	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.00375) + \
	1202	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.00375)*(1. - pt/5000.) + \
	1203	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
	1204	(abs(eta) > 4.0) * (0.00)}
	1205
	1206	add EfficiencyFormula {15} {
	1207
	1208	(pt <= 10.0) * (0.00) +
	1209	(abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.6) + \
	1210	(abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.6) *(1. - pt/5000.) + \
	1211	(abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
	1212	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.45) + \
	1213	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.45)*(1. - pt/5000.) + \
	1214	(abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
	1215	(abs(eta) > 4.0) * (0.00)}
	1216
	1217	}
	1218
	1219	#####################################################
	1220	# Find uniquely identified photons/electrons/tau/jets
	1221	#####################################################
	1222
	1223	module UniqueObjectFinder UniqueObjectFinder {
	1224	# earlier arrays take precedence over later ones
	1225	# add InputArray InputArray OutputArray
	1226	add InputArray PhotonEfficiency/photons photons
	1227	add InputArray ElectronEfficiency/electrons electrons
	1228	add InputArray MuonEfficiency/muons muons
	1229	add InputArray JetEnergyScale/jets jets
	1230
	1231	}
	1232
	1233
	1234	##################
	1235	# ROOT tree writer
	1236	##################
	1237
	1238	module TreeWriter TreeWriter {
	1239	# add Branch InputArray BranchName BranchClass
	1240	add Branch Delphes/allParticles Particle GenParticle
	1241
	1242	add Branch GenJetFinder/jets GenJet Jet
	1243	add Branch GenMissingET/momentum GenMissingET MissingET
	1244
	1245	add Branch TrackMerger/tracks Track Track
	1246	add Branch Calorimeter/towers Tower Tower
	1247
	1248	add Branch HCal/eflowTracks EFlowTrack Track
	1249	add Branch ECal/eflowPhotons EFlowPhoton Tower
	1250	add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
	1251
	1252	add Branch UniqueObjectFinder/photons Photon Photon
	1253	add Branch UniqueObjectFinder/electrons Electron Electron
	1254	add Branch UniqueObjectFinder/muons Muon Muon
	1255	add Branch UniqueObjectFinder/jets Jet Jet
	1256
	1257	add Branch FatJetFinder/jets FatJet Jet
	1258
	1259	add Branch MissingET/momentum MissingET MissingET
	1260	add Branch ScalarHT/energy ScalarHT ScalarHT
	1261	}

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

Download in other formats: