Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_HLLHC.tcl@ b4786d3

Visit:

Last change on this file since b4786d3 was e57c062, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago
fixed large R jet collections
Property mode set to `100644`
File size: 66.6 KB

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

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

Download in other formats: