Fork me on GitHub

Context Navigation

source: git/cards/CMS_PhaseII_0PU.tcl@ 181ff53

Visit:

ImprovedOutputFile Timing dual_readout llp

Last change on this file since 181ff53 was bd225dd, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago
udated track reso. for tilted
Property mode set to `100644`
File size: 102.4 KB

Rev	Line
[4f728cd]	1	#
	2	# Phase II - No Pile-Up
	3	#
	4	# Main authors: Michele Selvaggi (UCL)
	5	#
	6	# Released on:
	7	#
	8	# Version: v01
	9	#
	10	#
	11	#######################################
	12	# Order of execution of various modules
	13	#######################################
	14
	15	set ExecutionPath {
	16	ParticlePropagator
	17
	18	ChargedHadronTrackingEfficiency
	19	ElectronTrackingEfficiency
	20	MuonTrackingEfficiency
	21
	22	ChargedHadronMomentumSmearing
	23	ElectronMomentumSmearing
	24	MuonMomentumSmearing
	25
	26	TrackMerger
	27
	28	ECal
	29	HCal
	30
	31	TowerMerger
	32	EFlowMerger
	33
	34	PhotonIsolation
[576428a]	35	PhotonEfficiency
	36
	37	ElectronFilter
[4f728cd]	38	ElectronIsolation
[576428a]	39	ElectronEfficiency
	40
[4f728cd]	41	MuonIsolation
[2aaf039]	42	MuonLooseIdEfficiency
	43	MuonTightIdEfficiency
[576428a]	44
[4f728cd]	45	NeutrinoFilter
	46
	47	MissingET
	48	GenMissingET
	49
	50	GenJetFinder
	51	FastJetFinder
	52
	53	ScalarHT
	54
	55	JetEnergyScale
	56
	57	JetFlavorAssociation
	58
	59	BTaggingLoose
	60	BTaggingMedium
	61	BTaggingTight
	62
	63	TauTagging
	64
	65	TreeWriter
	66	}
	67
	68	#################################
	69	# Propagate particles in cylinder
	70	#################################
	71
	72	module ParticlePropagator ParticlePropagator {
	73	set InputArray Delphes/stableParticles
	74
	75	set OutputArray stableParticles
	76	set ChargedHadronOutputArray chargedHadrons
	77	set ElectronOutputArray electrons
	78	set MuonOutputArray muons
	79
	80	# radius of the magnetic field coverage, in m
	81	set Radius 1.29
	82	# half-length of the magnetic field coverage, in m
	83	set HalfLength 3.0
	84
	85	# magnetic field
	86	set Bz 3.8
	87	}
	88
	89
	90	####################################
	91	# Charged hadron tracking efficiency
	92	####################################
	93
	94	module Efficiency ChargedHadronTrackingEfficiency {
	95	## particles after propagation
	96	set InputArray ParticlePropagator/chargedHadrons
	97	set OutputArray chargedHadrons
	98	# tracking efficiency formula for charged hadrons
	99	set EfficiencyFormula {
	100	(pt <= 0.2) * (0.00) + \
	101	(abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
	102	(abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
	103	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
	104	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
	105	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
	106	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
	107	(abs(eta) > 4.0) * (0.00)
	108	}
	109	}
	110
	111
	112	#####################################
	113	# Electron tracking efficiency - ID
	114	####################################
	115
	116	module Efficiency ElectronTrackingEfficiency {
	117	set InputArray ParticlePropagator/electrons
	118	set OutputArray electrons
	119	# tracking efficiency formula for electrons
	120	set EfficiencyFormula {
	121	(pt <= 0.2) * (0.00) + \
	122	(abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
	123	(abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
	124	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
	125	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
	126	(abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
	127	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
	128	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
	129	(abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
	130	(abs(eta) > 4.0) * (0.00)
	131
	132	}
	133	}
	134
	135	##########################
	136	# Muon tracking efficiency
	137	##########################
	138
	139	module Efficiency MuonTrackingEfficiency {
	140	set InputArray ParticlePropagator/muons
	141	set OutputArray muons
	142	# tracking efficiency formula for muons
	143	set EfficiencyFormula {
	144	(pt <= 0.2) * (0.00) + \
[2aaf039]	145	(abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
	146	(abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
	147	(abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
	148	(abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
	149	(abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
	150	(abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
[4f728cd]	151	(abs(eta) > 4.0) * (0.00)
	152
	153	}
	154	}
	155
	156	########################################
	157	# Momentum resolution for charged tracks
	158	########################################
	159
	160	module MomentumSmearing ChargedHadronMomentumSmearing {
	161	## hadrons after having applied the tracking efficiency
	162	set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
	163	set OutputArray chargedHadrons
[bd225dd]	164	# resolution formula for charged hadrons ,
	165
	166	# from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
	167
	168	set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00467469) + \
	169	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004675 + (pt-1.000000)* 0.000056) + \
	170	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.005179 + (pt-10.000000)* 0.000064) + \
	171	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.010955*pt/100.000000) + \
	172	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00515885) + \
	173	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005159 + (pt-1.000000)* 0.000049) + \
	174	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005600 + (pt-10.000000)* 0.000062) + \
	175	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.011149*pt/100.000000) + \
	176	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00509775) + \
	177	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005098 + (pt-1.000000)* 0.000055) + \
	178	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005590 + (pt-10.000000)* 0.000063) + \
	179	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.011251*pt/100.000000) + \
	180	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00568785) + \
	181	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005688 + (pt-1.000000)* 0.000061) + \
	182	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.006237 + (pt-10.000000)* 0.000059) + \
	183	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.011563*pt/100.000000) + \
	184	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00668287) + \
	185	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.006683 + (pt-1.000000)* 0.000065) + \
	186	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007263 + (pt-10.000000)* 0.000067) + \
	187	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.013251*pt/100.000000) + \
	188	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01047734) + \
	189	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010477 + (pt-1.000000)* 0.000059) + \
	190	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.011005 + (pt-10.000000)* 0.000098) + \
	191	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.019785*pt/100.000000) + \
	192	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01430653) + \
	193	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014307 + (pt-1.000000)* 0.000038) + \
	194	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.014647 + (pt-10.000000)* 0.000087) + \
	195	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.022499*pt/100.000000) + \
	196	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01719240) + \
	197	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017192 + (pt-1.000000)* 0.000020) + \
	198	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.017372 + (pt-10.000000)* 0.000080) + \
	199	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.024604*pt/100.000000) + \
	200	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01783940) + \
	201	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.017839 + (pt-1.000000)* 0.000045) + \
	202	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.018244 + (pt-10.000000)* 0.000073) + \
	203	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.024839*pt/100.000000) + \
	204	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.01787758) + \
	205	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.017878 + (pt-1.000000)* 0.000154) + \
	206	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.019265 + (pt-10.000000)* 0.000141) + \
	207	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.031979*pt/100.000000) + \
	208	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01825549) + \
	209	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.018255 + (pt-1.000000)* 0.000270) + \
	210	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.020690 + (pt-10.000000)* 0.000212) + \
	211	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.039811*pt/100.000000) + \
	212	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01803308) + \
	213	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.018033 + (pt-1.000000)* 0.000220) + \
	214	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.020013 + (pt-10.000000)* 0.000247) + \
	215	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.042262*pt/100.000000) + \
	216	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.02156195) + \
	217	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.021562 + (pt-1.000000)* 0.000225) + \
	218	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.023591 + (pt-10.000000)* 0.000361) + \
	219	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.056075*pt/100.000000) + \
	220	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.02691276) + \
	221	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.026913 + (pt-1.000000)* 0.000357) + \
	222	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.030127 + (pt-10.000000)* 0.000483) + \
	223	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.073578*pt/100.000000) + \
	224	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
	225	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.032532 + (pt-1.000000)* 0.000382) + \
	226	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.035969 + (pt-10.000000)* 0.000673) + \
	227	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.096568*pt/100.000000) + \
	228	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
	229	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045517 + (pt-1.000000)* 0.002188) + \
	230	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.065207 + (pt-10.000000)* 0.003254) + \
	231	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.358089*pt/100.000000) + \
	232	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
	233	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.047248 + (pt-1.000000)* 0.003322) + \
	234	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.077142 + (pt-10.000000)* 0.004516) + \
	235	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.483608*pt/100.000000) + \
	236	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
	237	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.044930 + (pt-1.000000)* 0.002917) + \
	238	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.071183 + (pt-10.000000)* 0.004186) + \
	239	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.447915*pt/100.000000) + \
	240	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
	241	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.054642 + (pt-1.000000)* 0.003456) + \
	242	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.085748 + (pt-10.000000)* 0.004995) + \
	243	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.535262*pt/100.000000) + \
	244	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
	245	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.065512 + (pt-1.000000)* 0.006574) + \
	246	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.124674 + (pt-10.000000)* 0.009560) + \
	247	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.985040*pt/100.000000)
	248	}
	249
[4f728cd]	250	}
	251
	252	########################################
	253	# Momentum resolution for electrons
	254	########################################
	255
	256	module MomentumSmearing ElectronMomentumSmearing {
	257	## hadrons after having applied the tracking efficiency
	258	set InputArray ElectronTrackingEfficiency/electrons
	259	set OutputArray electrons
	260	# resolution formula for electrons
[bd225dd]	261	# from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
	262
	263	set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00467469) + \
	264	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004675 + (pt-1.000000)* 0.000056) + \
	265	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.005179 + (pt-10.000000)* 0.000064) + \
	266	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.010955*pt/100.000000) + \
	267	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00515885) + \
	268	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005159 + (pt-1.000000)* 0.000049) + \
	269	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005600 + (pt-10.000000)* 0.000062) + \
	270	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.011149*pt/100.000000) + \
	271	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00509775) + \
	272	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005098 + (pt-1.000000)* 0.000055) + \
	273	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005590 + (pt-10.000000)* 0.000063) + \
	274	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.011251*pt/100.000000) + \
	275	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00568785) + \
	276	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005688 + (pt-1.000000)* 0.000061) + \
	277	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.006237 + (pt-10.000000)* 0.000059) + \
	278	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.011563*pt/100.000000) + \
	279	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00668287) + \
	280	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.006683 + (pt-1.000000)* 0.000065) + \
	281	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007263 + (pt-10.000000)* 0.000067) + \
	282	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.013251*pt/100.000000) + \
	283	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01047734) + \
	284	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010477 + (pt-1.000000)* 0.000059) + \
	285	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.011005 + (pt-10.000000)* 0.000098) + \
	286	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.019785*pt/100.000000) + \
	287	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01430653) + \
	288	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014307 + (pt-1.000000)* 0.000038) + \
	289	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.014647 + (pt-10.000000)* 0.000087) + \
	290	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.022499*pt/100.000000) + \
	291	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01719240) + \
	292	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017192 + (pt-1.000000)* 0.000020) + \
	293	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.017372 + (pt-10.000000)* 0.000080) + \
	294	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.024604*pt/100.000000) + \
	295	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01783940) + \
	296	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.017839 + (pt-1.000000)* 0.000045) + \
	297	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.018244 + (pt-10.000000)* 0.000073) + \
	298	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.024839*pt/100.000000) + \
	299	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.01787758) + \
	300	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.017878 + (pt-1.000000)* 0.000154) + \
	301	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.019265 + (pt-10.000000)* 0.000141) + \
	302	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.031979*pt/100.000000) + \
	303	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01825549) + \
	304	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.018255 + (pt-1.000000)* 0.000270) + \
	305	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.020690 + (pt-10.000000)* 0.000212) + \
	306	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.039811*pt/100.000000) + \
	307	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01803308) + \
	308	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.018033 + (pt-1.000000)* 0.000220) + \
	309	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.020013 + (pt-10.000000)* 0.000247) + \
	310	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.042262*pt/100.000000) + \
	311	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.02156195) + \
	312	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.021562 + (pt-1.000000)* 0.000225) + \
	313	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.023591 + (pt-10.000000)* 0.000361) + \
	314	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.056075*pt/100.000000) + \
	315	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.02691276) + \
	316	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.026913 + (pt-1.000000)* 0.000357) + \
	317	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.030127 + (pt-10.000000)* 0.000483) + \
	318	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.073578*pt/100.000000) + \
	319	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
	320	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.032532 + (pt-1.000000)* 0.000382) + \
	321	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.035969 + (pt-10.000000)* 0.000673) + \
	322	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.096568*pt/100.000000) + \
	323	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
	324	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045517 + (pt-1.000000)* 0.002188) + \
	325	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.065207 + (pt-10.000000)* 0.003254) + \
	326	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.358089*pt/100.000000) + \
	327	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
	328	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.047248 + (pt-1.000000)* 0.003322) + \
	329	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.077142 + (pt-10.000000)* 0.004516) + \
	330	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.483608*pt/100.000000) + \
	331	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
	332	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.044930 + (pt-1.000000)* 0.002917) + \
	333	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.071183 + (pt-10.000000)* 0.004186) + \
	334	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.447915*pt/100.000000) + \
	335	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
	336	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.054642 + (pt-1.000000)* 0.003456) + \
	337	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.085748 + (pt-10.000000)* 0.004995) + \
	338	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.535262*pt/100.000000) + \
	339	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
	340	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.065512 + (pt-1.000000)* 0.006574) + \
	341	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.124674 + (pt-10.000000)* 0.009560) + \
	342	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.985040*pt/100.000000)
	343	}
	344
[4f728cd]	345	}
	346
	347
	348	###############################
	349	# Momentum resolution for muons
	350	###############################
	351
	352	module MomentumSmearing MuonMomentumSmearing {
	353	set InputArray MuonTrackingEfficiency/muons
	354	set OutputArray muons
	355	# resolution formula for muons
[a0236d5]	356
	357	# up to \|eta\| < 2.8 take measurement from tracking + muon chambers
[bd225dd]	358	# for \|eta\| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
	359	# http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
[a0236d5]	360
[4f728cd]	361	set ResolutionFormula {
[a0236d5]	362
	363	(pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00519 + \
	364	(pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00575 + \
	365	(pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00629 + \
	366	(pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00626 + \
	367	(pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00722 + \
	368	(pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.00968 + \
	369	(pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00983 + \
	370	(pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0102 + \
	371	(pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0126 + \
	372	(pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0178 + \
	373	(pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0248 + \
	374	(pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0347 + \
	375	(pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0359 + \
	376	(pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.0412 + \
	377	(pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0422 + \
	378	(pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0496 + \
	379
	380	(pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00533 + \
	381	(pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00592 + \
	382	(pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00602 + \
	383	(pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00614 + \
	384	(pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00736 + \
	385	(pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0104 + \
	386	(pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00968 + \
	387	(pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.00991 + \
	388	(pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0106 + \
	389	(pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0158 + \
	390	(pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0207 + \
	391	(pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0327 + \
	392	(pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0369 + \
	393	(pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.039 + \
	394	(pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0488 + \
	395	(pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0543 + \
	396
	397	(pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00559 + \
	398	(pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00586 + \
	399	(pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00614 + \
	400	(pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00654 + \
	401	(pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00759 + \
	402	(pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0103 + \
	403	(pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00963 + \
	404	(pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0102 + \
	405	(pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.011 + \
	406	(pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0142 + \
	407	(pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0197 + \
	408	(pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0315 + \
	409	(pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0409 + \
	410	(pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.0471 + \
	411	(pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0579 + \
	412	(pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0679 + \
	413
	414	(pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00656 + \
	415	(pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.007 + \
	416	(pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00733 + \
	417	(pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00752 + \
	418	(pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00851 + \
	419	(pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0105 + \
	420	(pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.0105 + \
	421	(pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0111 + \
	422	(pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0119 + \
	423	(pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0165 + \
	424	(pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0242 + \
	425	(pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0359 + \
	426	(pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0474 + \
	427	(pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.052 + \
	428	(pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0639 + \
	429	(pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0763 + \
	430
[bd225dd]	431	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
	432	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.032532 + (pt-1.000000)* 0.000382) + \
	433	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.035969 + (pt-10.000000)* 0.000673) + \
	434	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.096568*pt/100.000000) + \
	435	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
	436	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045517 + (pt-1.000000)* 0.002188) + \
	437	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.065207 + (pt-10.000000)* 0.003254) + \
	438	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.358089*pt/100.000000) + \
	439	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
	440	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.047248 + (pt-1.000000)* 0.003322) + \
	441	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.077142 + (pt-10.000000)* 0.004516) + \
	442	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.483608*pt/100.000000) + \
	443	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
	444	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.044930 + (pt-1.000000)* 0.002917) + \
	445	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.071183 + (pt-10.000000)* 0.004186) + \
	446	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.447915*pt/100.000000) + \
	447	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
	448	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.054642 + (pt-1.000000)* 0.003456) + \
	449	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.085748 + (pt-10.000000)* 0.004995) + \
	450	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.535262*pt/100.000000) + \
	451	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
	452	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.065512 + (pt-1.000000)* 0.006574) + \
	453	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.124674 + (pt-10.000000)* 0.009560) + \
	454	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.985040*pt/100.000000) + \
[a0236d5]	455
[bd225dd]	456	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00467469) + \
	457	(abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 5.000) * (0.004675 + (pt-1.000000)* 0.000056) + \
	458	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00515885) + \
	459	(abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 5.000) * (0.005159 + (pt-1.000000)* 0.000049) + \
	460	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00509775) + \
	461	(abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 5.000) * (0.005098 + (pt-1.000000)* 0.000055) + \
	462	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00568785) + \
	463	(abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 5.000) * (0.005688 + (pt-1.000000)* 0.000061) + \
	464	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00668287) + \
	465	(abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 5.000) * (0.006683 + (pt-1.000000)* 0.000065) + \
	466	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01047734) + \
	467	(abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 5.000) * (0.010477 + (pt-1.000000)* 0.000059) + \
	468	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01430653) + \
	469	(abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 5.000) * (0.014307 + (pt-1.000000)* 0.000038) + \
	470	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01719240) + \
	471	(abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 5.000) * (0.017192 + (pt-1.000000)* 0.000020) + \
	472	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01783940) + \
	473	(abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 5.000) * (0.017839 + (pt-1.000000)* 0.000045) + \
	474	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.01787758) + \
	475	(abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 5.000) * (0.017878 + (pt-1.000000)* 0.000154) + \
	476	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01825549) + \
	477	(abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 5.000) * (0.018255 + (pt-1.000000)* 0.000270) + \
	478	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01803308) + \
	479	(abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 5.000) * (0.018033 + (pt-1.000000)* 0.000220) + \
	480	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.02156195) + \
	481	(abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 5.000) * (0.021562 + (pt-1.000000)* 0.000225) + \
	482	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.02691276) + \
	483	(abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 5.000) * (0.026913 + (pt-1.000000)* 0.000357) + \
	484	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
	485	(abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 5.000) * (0.032532 + (pt-1.000000)* 0.000382) + \
	486	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
	487	(abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 5.000) * (0.045517 + (pt-1.000000)* 0.002188) + \
	488	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
	489	(abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 5.000) * (0.047248 + (pt-1.000000)* 0.003322) + \
	490	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
	491	(abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 5.000) * (0.044930 + (pt-1.000000)* 0.002917) + \
	492	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
	493	(abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 5.000) * (0.054642 + (pt-1.000000)* 0.003456) + \
	494	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
	495	(abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 5.000) * (0.065512 + (pt-1.000000)* 0.006574)
	496
[4f728cd]	497	}
	498	}
	499
	500
	501	##############
	502	# Track merger
	503	##############
	504
	505	module Merger TrackMerger {
	506	# add InputArray InputArray
	507	add InputArray ChargedHadronMomentumSmearing/chargedHadrons
	508	add InputArray ElectronMomentumSmearing/electrons
	509	add InputArray MuonMomentumSmearing/muons
	510	set OutputArray tracks
	511	}
	512
	513	#############
	514	# ECAL
	515	#############
	516
	517	module SimpleCalorimeter ECal {
	518	set ParticleInputArray ParticlePropagator/stableParticles
	519	set TrackInputArray TrackMerger/tracks
	520
	521	set TowerOutputArray ecalTowers
	522	set EFlowTrackOutputArray eflowTracks
	523	set EFlowTowerOutputArray eflowPhotons
	524
	525	set IsEcal true
	526
	527	set EnergyMin 0.5
	528	set EnergySignificanceMin 1.0
	529
	530	set SmearTowerCenter true
	531
	532	set pi [expr {acos(-1)}]
	533
	534	# lists of the edges of each tower in eta and phi
	535	# each list starts with the lower edge of the first tower
	536	# the list ends with the higher edged of the last tower
	537
	538	# assume 0.02 x 0.02 resolution in eta,phi in the barrel \|eta\| < 1.5
	539
	540	set PhiBins {}
[5303656]	541	for {set i -180} {$i <= 180} {incr i} {
	542	add PhiBins [expr {$i * $pi/180.0}]
[4f728cd]	543	}
	544
	545	# 0.02 unit in eta up to eta = 1.5 (barrel)
[5303656]	546	for {set i -85} {$i <= 86} {incr i} {
	547	set eta [expr {$i * 0.0174}]
[4f728cd]	548	add EtaPhiBins $eta $PhiBins
	549	}
	550
[5303656]	551	# assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < \|eta\| < 3.0 (HGCAL- ECAL)
[4f728cd]	552
	553	set PhiBins {}
[5303656]	554	for {set i -180} {$i <= 180} {incr i} {
	555	add PhiBins [expr {$i * $pi/180.0}]
[4f728cd]	556	}
	557
[5303656]	558	# 0.02 unit in eta up to eta = 3
	559	for {set i 1} {$i <= 84} {incr i} {
	560	set eta [expr { -2.958 + $i * 0.0174}]
[4f728cd]	561	add EtaPhiBins $eta $PhiBins
	562	}
	563
[5303656]	564	for {set i 1} {$i <= 84} {incr i} {
	565	set eta [expr { 1.4964 + $i * 0.0174}]
[4f728cd]	566	add EtaPhiBins $eta $PhiBins
	567	}
	568
	569	# take present CMS granularity for HF
	570
	571	# 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < \|eta\| < 5.0
	572	set PhiBins {}
	573	for {set i -18} {$i <= 18} {incr i} {
	574	add PhiBins [expr {$i * $pi/18.0}]
	575	}
	576
[5303656]	577	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} {
[4f728cd]	578	add EtaPhiBins $eta $PhiBins
	579	}
	580
	581
	582	add EnergyFraction {0} {0.0}
	583	# energy fractions for e, gamma and pi0
	584	add EnergyFraction {11} {1.0}
	585	add EnergyFraction {22} {1.0}
	586	add EnergyFraction {111} {1.0}
	587	# energy fractions for muon, neutrinos and neutralinos
	588	add EnergyFraction {12} {0.0}
	589	add EnergyFraction {13} {0.0}
	590	add EnergyFraction {14} {0.0}
	591	add EnergyFraction {16} {0.0}
	592	add EnergyFraction {1000022} {0.0}
	593	add EnergyFraction {1000023} {0.0}
	594	add EnergyFraction {1000025} {0.0}
	595	add EnergyFraction {1000035} {0.0}
	596	add EnergyFraction {1000045} {0.0}
	597	# energy fractions for K0short and Lambda
	598	add EnergyFraction {310} {0.3}
	599	add EnergyFraction {3122} {0.3}
	600
	601	# set ResolutionFormula {resolution formula as a function of eta and energy}
	602
	603	# for the ECAL barrel (\|eta\| < 1.5), see hep-ex/1306.2016.
	604	# for the endcaps (1.5 < \|eta\| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
	605
[5303656]	606	set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^20.005^2 + energy0.027^2 + 0.15^2) + \
[4f728cd]	607	(abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^20.006^2 + energy0.20^2) + \
	608	(abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^20.007^2 + energy0.21^2) + \
	609	(abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^20.008^2 + energy0.24^2) + \
[5303656]	610	(abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^20.08^2 + energy1.98^2)}
[4f728cd]	611
	612	}
	613
	614	#############
	615	# HCAL
	616	#############
	617
	618	module SimpleCalorimeter HCal {
	619	set ParticleInputArray ParticlePropagator/stableParticles
	620	set TrackInputArray ECal/eflowTracks
	621
	622	set TowerOutputArray hcalTowers
	623	set EFlowTrackOutputArray eflowTracks
	624	set EFlowTowerOutputArray eflowNeutralHadrons
	625
	626	set IsEcal false
	627
	628	set EnergyMin 1.0
	629	set EnergySignificanceMin 1.0
	630
	631	set SmearTowerCenter true
	632
	633	set pi [expr {acos(-1)}]
	634
	635	# lists of the edges of each tower in eta and phi
	636	# each list starts with the lower edge of the first tower
	637	# the list ends with the higher edged of the last tower
	638
	639	# assume 0.087 x 0.087 resolution in eta,phi in the barrel \|eta\| < 1.5
	640
	641	set PhiBins {}
	642	for {set i -36} {$i <= 36} {incr i} {
	643	add PhiBins [expr {$i * $pi/36.0}]
	644	}
	645	foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.65} {
	646	add EtaPhiBins $eta $PhiBins
	647	}
	648
[61c0925]	649	# assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < \|eta\| < 3.0 (HGCAL- HCAL)
	650
[4f728cd]	651	set PhiBins {}
[61c0925]	652	for {set i -180} {$i <= 180} {incr i} {
	653	add PhiBins [expr {$i * $pi/180.0}]
[4f728cd]	654	}
	655
[61c0925]	656	# 0.02 unit in eta up to eta = 3
	657	for {set i 1} {$i <= 84} {incr i} {
	658	set eta [expr { -2.958 + $i * 0.0174}]
[4f728cd]	659	add EtaPhiBins $eta $PhiBins
	660	}
	661
[61c0925]	662	for {set i 1} {$i <= 84} {incr i} {
	663	set eta [expr { 1.4964 + $i * 0.0174}]
[4f728cd]	664	add EtaPhiBins $eta $PhiBins
[61c0925]	665	}
[4f728cd]	666
	667	# take present CMS granularity for HF
	668
	669	# 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < \|eta\| < 5.0
	670	set PhiBins {}
	671	for {set i -18} {$i <= 18} {incr i} {
	672	add PhiBins [expr {$i * $pi/18.0}]
	673	}
	674
[61c0925]	675	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} {
[4f728cd]	676	add EtaPhiBins $eta $PhiBins
	677	}
	678
	679	# default energy fractions {abs(PDG code)} {Fecal Fhcal}
	680	add EnergyFraction {0} {1.0}
	681	# energy fractions for e, gamma and pi0
	682	add EnergyFraction {11} {0.0}
	683	add EnergyFraction {22} {0.0}
	684	add EnergyFraction {111} {0.0}
	685	# energy fractions for muon, neutrinos and neutralinos
	686	add EnergyFraction {12} {0.0}
	687	add EnergyFraction {13} {0.0}
	688	add EnergyFraction {14} {0.0}
	689	add EnergyFraction {16} {0.0}
	690	add EnergyFraction {1000022} {0.0}
	691	add EnergyFraction {1000023} {0.0}
	692	add EnergyFraction {1000025} {0.0}
	693	add EnergyFraction {1000035} {0.0}
	694	add EnergyFraction {1000045} {0.0}
	695	# energy fractions for K0short and Lambda
	696	add EnergyFraction {310} {0.7}
	697	add EnergyFraction {3122} {0.7}
	698
	699	# set ResolutionFormula {resolution formula as a function of eta and energy}
[5303656]	700	set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^20.05^2 + energy1.00^2) + \
	701	(abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^20.05^2 + energy1.00^2) + \
	702	(abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^20.11^2 + energy2.80^2)}
[4f728cd]	703
	704	}
	705
	706	#################
	707	# Electron filter
	708	#################
	709
	710	module PdgCodeFilter ElectronFilter {
	711	set InputArray HCal/eflowTracks
	712	set OutputArray electrons
	713	set Invert true
	714	add PdgCode {11}
	715	add PdgCode {-11}
	716	}
	717
	718	###################################################
	719	# Tower Merger (in case not using e-flow algorithm)
	720	###################################################
	721
	722	module Merger TowerMerger {
	723	# add InputArray InputArray
	724	add InputArray ECal/ecalTowers
	725	add InputArray HCal/hcalTowers
	726	set OutputArray towers
	727	}
	728
	729	####################
	730	# Energy flow merger
	731	####################
	732
	733	module Merger EFlowMerger {
	734	# add InputArray InputArray
	735	add InputArray HCal/eflowTracks
	736	add InputArray ECal/eflowPhotons
	737	add InputArray HCal/eflowNeutralHadrons
	738	set OutputArray eflow
	739	}
	740
	741
	742	###################
	743	# Missing ET merger
	744	###################
	745
	746	module Merger MissingET {
	747	# add InputArray InputArray
	748	add InputArray EFlowMerger/eflow
	749	set MomentumOutputArray momentum
	750	}
	751
	752
	753
	754	##################
	755	# Scalar HT merger
	756	##################
	757
	758	module Merger ScalarHT {
	759	# add InputArray InputArray
	760	add InputArray EFlowMerger/eflow
	761	set EnergyOutputArray energy
	762	}
	763
	764	#################
	765	# Neutrino Filter
	766	#################
	767
	768	module PdgCodeFilter NeutrinoFilter {
	769
	770	set InputArray Delphes/stableParticles
	771	set OutputArray filteredParticles
	772
	773	set PTMin 0.0
	774
	775	add PdgCode {12}
	776	add PdgCode {14}
	777	add PdgCode {16}
	778	add PdgCode {-12}
	779	add PdgCode {-14}
	780	add PdgCode {-16}
	781
	782	}
	783
	784
	785	#####################
	786	# MC truth jet finder
	787	#####################
	788
	789	module FastJetFinder GenJetFinder {
	790	set InputArray NeutrinoFilter/filteredParticles
	791
	792	set OutputArray jets
	793
	794	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	795	set JetAlgorithm 6
	796	set ParameterR 0.4
	797
	798	set JetPTMin 15.0
	799	}
	800
	801	#########################
	802	# Gen Missing ET merger
	803	########################
	804
	805	module Merger GenMissingET {
	806
	807	# add InputArray InputArray
	808	add InputArray NeutrinoFilter/filteredParticles
	809	set MomentumOutputArray momentum
	810	}
	811
	812
	813
	814	############
	815	# Jet finder
	816	############
	817
	818	module FastJetFinder FastJetFinder {
	819	# set InputArray TowerMerger/towers
	820	set InputArray EFlowMerger/eflow
	821
	822	set OutputArray jets
	823
	824	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	825	set JetAlgorithm 6
	826	set ParameterR 0.4
	827
	828	set JetPTMin 15.0
	829	}
	830
	831	##################
	832	# Jet Energy Scale
	833	##################
	834
	835	module EnergyScale JetEnergyScale {
	836	set InputArray FastJetFinder/jets
	837	set OutputArray jets
	838
	839	# scale formula for jets
	840	set ScaleFormula {1.00}
	841	}
	842
[576428a]	843	####################
	844	# Photon isolation #
	845	####################
	846
	847	module Isolation PhotonIsolation {
	848
	849	# particle for which calculate the isolation
	850	set CandidateInputArray ECal/eflowPhotons
	851
	852	# isolation collection
	853	set IsolationInputArray EFlowMerger/eflow
	854
	855	# output array
	856	set OutputArray photons
	857
	858	# isolation cone
	859	set DeltaRMax 0.3
	860
	861	# minimum pT
	862	set PTMin 1.0
	863
	864	# iso ratio to cut
	865	set PTRatioMax 9999.
	866	}
	867
	868
[4f728cd]	869
	870	#####################
	871	# Photon efficiency #
	872	#####################
	873
	874	module Efficiency PhotonEfficiency {
	875
	876	## input particles
[576428a]	877	set InputArray PhotonIsolation/photons
[4f728cd]	878	## output particles
	879	set OutputArray photons
	880	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	881	# efficiency formula for photons
	882	set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
	883	(abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
	884	(abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
	885	(abs(eta) > 4.0) * (0.00)}
	886
	887	}
	888
	889
[576428a]	890	######################
	891	# Electron isolation #
	892	######################
	893
	894	module Isolation ElectronIsolation {
[4f728cd]	895
[576428a]	896	set CandidateInputArray ElectronFilter/electrons
[4f728cd]	897
	898	# isolation collection
	899	set IsolationInputArray EFlowMerger/eflow
	900
[576428a]	901	set OutputArray electrons
[4f728cd]	902
[576428a]	903	set DeltaRMax 0.3
	904	set PTMin 1.0
[4f728cd]	905	set PTRatioMax 9999.
[576428a]	906
[4f728cd]	907	}
	908
	909
[576428a]	910
[4f728cd]	911	#######################
	912	# Electron efficiency #
	913	#######################
	914
	915	module Efficiency ElectronEfficiency {
	916
[576428a]	917	set InputArray ElectronIsolation/electrons
[4f728cd]	918	set OutputArray electrons
	919
	920	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	921	# efficiency formula for electrons
	922	set EfficiencyFormula {
	923	(pt <= 4.0) * (0.00) + \
	924	(abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
	925	(abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
	926	(abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
	927	(abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
	928	(abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
	929	(abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
	930	(abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
	931	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
	932	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
	933	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
	934	(abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
	935	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
	936	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
	937	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
	938	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
	939	(abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
	940	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	941	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	942	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	943	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	944	(abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
	945	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
	946	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
	947	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
	948	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
	949	(abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
	950	(abs(eta) > 4.0) * (0.00)
	951
	952	}
	953	}
	954
[576428a]	955	##################
	956	# Muon isolation #
	957	##################
[4f728cd]	958
[576428a]	959	module Isolation MuonIsolation {
	960	set CandidateInputArray MuonMomentumSmearing/muons
	961
[4f728cd]	962	# isolation collection
	963	set IsolationInputArray EFlowMerger/eflow
[576428a]	964
	965	set OutputArray muons
[4f728cd]	966
	967	set DeltaRMax 0.3
	968	set PTMin 1.0
	969	set PTRatioMax 9999.
	970
	971	}
	972
[2aaf039]	973	##################
	974	# Muon Loose Id #
	975	##################
[4f728cd]	976
[2aaf039]	977	module Efficiency MuonLooseIdEfficiency {
	978	set InputArray MuonIsolation/muons
	979	set OutputArray muons
	980	# tracking + LooseID efficiency formula for muons
	981	set EfficiencyFormula {
	982
	983	(pt <= 5.0) * (0.00) + \
	984	(pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.967 + \
	985	(pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.968 + \
	986	(pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
	987	(pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.982 + \
	988	(pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.972 + \
	989	(pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.974 + \
	990	(pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.969 + \
	991	(pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.988 + \
	992	(pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.989 + \
	993	(pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.992 + \
	994	(pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.979 + \
	995	(pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.976 + \
	996	(pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.882 + \
	997	(pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.883 + \
	998	(pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.851 + \
	999	(pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.828 + \
	1000
	1001	(pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.989 + \
	1002	(pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
	1003	(pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
	1004	(pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.995 + \
	1005	(pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.987 + \
	1006	(pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
	1007	(pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
	1008	(pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.99 + \
	1009	(pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.997 + \
	1010	(pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
	1011	(pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.989 + \
	1012	(pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.995 + \
	1013	(pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.888 + \
	1014	(pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.894 + \
	1015	(pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.886 + \
	1016	(pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.883 + \
	1017
	1018	(pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.981 + \
	1019	(pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.994 + \
	1020	(pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
	1021	(pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.997 + \
	1022	(pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.986 + \
	1023	(pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.991 + \
	1024	(pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.987 + \
	1025	(pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.992 + \
	1026	(pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.993 + \
	1027	(pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.999 + \
	1028	(pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.994 + \
	1029	(pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.997 + \
	1030	(pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.894 + \
	1031	(pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.879 + \
	1032	(pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.871 + \
	1033	(pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.864 + \
	1034
	1035	(pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.986 + \
	1036	(pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
	1037	(pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.996 + \
	1038	(pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.994 + \
	1039	(pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.991 + \
	1040	(pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
	1041	(pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.985 + \
	1042	(pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.993 + \
	1043	(pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.996 + \
	1044	(pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
	1045	(pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.99 + \
	1046	(pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.993 + \
	1047	(pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.905 + \
	1048	(pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.86 + \
	1049	(pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.845 + \
	1050	(pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.855 + \
	1051
	1052	(abs(eta) > 2.8) * (0.00)
	1053
	1054	}
[4f728cd]	1055	}
	1056
	1057
[2aaf039]	1058	##################
	1059	# Muon Tight Id #
	1060	##################
[4f728cd]	1061
[2aaf039]	1062	module Efficiency MuonTightIdEfficiency {
	1063	set InputArray MuonIsolation/muons
	1064	set OutputArray muons
	1065	# tracking + TightID efficiency formula for muons
	1066	set EfficiencyFormula {
	1067
	1068	(pt <= 5.0) * (0.00) + \
	1069	(pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.949 + \
	1070	(pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.894 + \
	1071	(pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.966 + \
	1072	(pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.964 + \
	1073	(pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.923 + \
	1074	(pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.94 + \
	1075	(pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.959 + \
	1076	(pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.972 + \
	1077	(pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.962 + \
	1078	(pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.972 + \
	1079	(pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.947 + \
	1080	(pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.911 + \
	1081	(pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.825 + \
	1082	(pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.821 + \
	1083	(pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.783 + \
	1084	(pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.746 + \
	1085
	1086	(pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
	1087	(pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
	1088	(pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
	1089	(pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.976 + \
	1090	(pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.954 + \
	1091	(pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.978 + \
	1092	(pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.975 + \
	1093	(pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.98 + \
	1094	(pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.982 + \
	1095	(pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
	1096	(pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.962 + \
	1097	(pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.946 + \
	1098	(pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.839 + \
	1099	(pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.844 + \
	1100	(pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.803 + \
	1101	(pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.791 + \
	1102
	1103	(pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.965 + \
	1104	(pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.938 + \
	1105	(pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.99 + \
	1106	(pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.977 + \
	1107	(pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.959 + \
	1108	(pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.976 + \
	1109	(pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.984 + \
	1110	(pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.981 + \
	1111	(pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.978 + \
	1112	(pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
	1113	(pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.969 + \
	1114	(pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.964 + \
	1115	(pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.863 + \
	1116	(pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.828 + \
	1117	(pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.807 + \
	1118	(pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.764 + \
	1119
	1120	(pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
	1121	(pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
	1122	(pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.991 + \
	1123	(pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.983 + \
	1124	(pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.957 + \
	1125	(pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.979 + \
	1126	(pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
	1127	(pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.983 + \
	1128	(pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.981 + \
	1129	(pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.984 + \
	1130	(pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.964 + \
	1131	(pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.956 + \
	1132	(pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.855 + \
	1133	(pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.817 + \
	1134	(pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.764 + \
	1135	(pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.752 + \
	1136
	1137	(abs(eta) > 2.8) * (0.00)
	1138
	1139	}
	1140	}
[4f728cd]	1141
	1142	########################
	1143	# Jet Flavor Association
	1144	########################
	1145
	1146	module JetFlavorAssociation JetFlavorAssociation {
	1147
	1148	set PartonInputArray Delphes/partons
	1149	set ParticleInputArray Delphes/allParticles
	1150	set ParticleLHEFInputArray Delphes/allParticlesLHEF
	1151	set JetInputArray JetEnergyScale/jets
	1152
	1153	set DeltaR 0.5
	1154	set PartonPTMin 10.0
	1155	set PartonEtaMax 4.0
	1156
	1157	}
	1158
	1159
	1160	#############
	1161	# b-tagging #
	1162	#############
	1163	module BTagging BTaggingLoose {
	1164
	1165	set JetInputArray JetEnergyScale/jets
	1166
	1167	set BitNumber 0
	1168
	1169	add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
	1170	(pt > 20.0) * (abs(eta) <= 3.4) * (0.1) + \
	1171	(abs(eta) > 3.4) * (0.000)
	1172	}
	1173
	1174
	1175	add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
	1176	(abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.341) + \
	1177	(abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.389) + \
	1178	(abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.416) + \
	1179	(abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.432) + \
	1180	(abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.447) + \
	1181	(abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.451) + \
	1182	(abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.451) + \
	1183	(abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.447) + \
	1184	(abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.447) + \
	1185	(abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.444) + \
	1186	(abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.430) + \
	1187	(abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.424) + \
	1188	(abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.419) + \
	1189	(abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.407) + \
	1190	(abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.384) + \
	1191	(abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.364) + \
	1192	(abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.353) + \
	1193	(abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.303) + \
	1194	(abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.273) + \
	1195	(abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.249) + \
	1196	(abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.233) + \
	1197	(abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.209) + \
	1198	(abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.194) + \
	1199	(abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.227) + \
	1200	(abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
	1201
	1202	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
	1203	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.236) + \
	1204	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.276) + \
	1205	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.302) + \
	1206	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.330) + \
	1207	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.338) + \
	1208	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.355) + \
	1209	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.353) + \
	1210	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.358) + \
	1211	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.350) + \
	1212	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.327) + \
	1213	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.313) + \
	1214	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.307) + \
	1215	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.322) + \
	1216	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.299) + \
	1217	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.274) + \
	1218	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.240) + \
	1219	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.203) + \
	1220	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.218) + \
	1221	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.203) + \
	1222	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.178) + \
	1223	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
	1224
	1225	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
	1226	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.179) + \
	1227	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.198) + \
	1228	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.219) + \
	1229	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.225) + \
	1230	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.232) + \
	1231	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.251) + \
	1232	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.267) + \
	1233	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.252) + \
	1234	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.243) + \
	1235	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.244) + \
	1236	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.242) + \
	1237	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.246) + \
	1238	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.217) + \
	1239	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.223)+ \
	1240	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.194) + \
	1241	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.187) + \
	1242	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
	1243
	1244	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
	1245	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.09) + \
	1246	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.10) + \
	1247	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.11) + \
	1248	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.11) + \
	1249	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.12) + \
	1250	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.13) + \
	1251	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.13) + \
	1252	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.13) + \
	1253	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.12) + \
	1254	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.12) + \
	1255	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.12) + \
	1256	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.12) + \
	1257	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.11) + \
	1258	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.11) + \
	1259	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.10) + \
	1260	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.09) + \
	1261	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
	1262	(abs(eta) > 3.4) * (0.000)
	1263	}
	1264
	1265	add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
	1266	(abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.712) + \
	1267	(abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.788) + \
	1268	(abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
	1269	(abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.856) + \
	1270	(abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.873) + \
	1271	(abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.886) + \
	1272	(abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.883) + \
	1273	(abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.882) + \
	1274	(abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.883) + \
	1275	(abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.883) + \
	1276	(abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.876) + \
	1277	(abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.872) + \
	1278	(abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.867) + \
	1279	(abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.853) + \
	1280	(abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.831) + \
	1281	(abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.819) + \
	1282	(abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.791) + \
	1283	(abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.740) + \
	1284	(abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.699) + \
	1285	(abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.631) + \
	1286	(abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.600) + \
	1287	(abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.567) + \
	1288	(abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.478) + \
	1289	(abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.470) + \
	1290	(abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
	1291
	1292	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
	1293	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.543) + \
	1294	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.627) + \
	1295	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.701) + \
	1296	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.741) + \
	1297	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.754) + \
	1298	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.796) + \
	1299	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.778) + \
	1300	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.798) + \
	1301	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.788) + \
	1302	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.775) + \
	1303	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.765) + \
	1304	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.752) + \
	1305	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.740) + \
	1306	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.716) + \
	1307	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.671) + \
	1308	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.637) + \
	1309	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.577) + \
	1310	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.572) + \
	1311	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.514) + \
	1312	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.516) + \
	1313	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
	1314
	1315	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
	1316	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.407) + \
	1317	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.479) + \
	1318	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.529) + \
	1319	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.578) + \
	1320	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.612) + \
	1321	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.673) + \
	1322	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.656) + \
	1323	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.658) + \
	1324	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.647) + \
	1325	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.629) + \
	1326	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.630) + \
	1327	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.615) + \
	1328	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.604) + \
	1329	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.532)+ \
	1330	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.532) + \
	1331	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.503) + \
	1332	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
	1333
	1334	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
	1335	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.20) + \
	1336	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.235) + \
	1337	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.26) + \
	1338	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.28) + \
	1339	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.31) + \
	1340	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.33) + \
	1341	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.32) + \
	1342	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.32) + \
	1343	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.33) + \
	1344	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.31) + \
	1345	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.32) + \
	1346	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.31) + \
	1347	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.30) + \
	1348	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.26) + \
	1349	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.27) + \
	1350	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.25) + \
	1351	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
	1352	(abs(eta) > 3.4) * (0.000)
	1353	}
	1354	}
	1355
	1356
	1357	#############
	1358	# b-tagging #
	1359	#############
	1360	module BTagging BTaggingMedium {
	1361
	1362	set JetInputArray JetEnergyScale/jets
	1363
	1364	set BitNumber 1
	1365
	1366	add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
	1367	(pt > 20.0) * (abs(eta) <= 3.4) * (0.01) + \
	1368	(abs(eta) > 3.4) * (0.000)
	1369	}
	1370
	1371	add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
	1372	(abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.144) + \
	1373	(abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.173) + \
	1374	(abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.188) + \
	1375	(abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.200) + \
	1376	(abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.209) + \
	1377	(abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.210) + \
	1378	(abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.206) + \
	1379	(abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.205) + \
	1380	(abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.198) + \
	1381	(abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.187) + \
	1382	(abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.175) + \
	1383	(abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.167) + \
	1384	(abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.157) + \
	1385	(abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.138) + \
	1386	(abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.118) + \
	1387	(abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.094) + \
	1388	(abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.083) + \
	1389	(abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.062) + \
	1390	(abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.049) + \
	1391	(abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.044) + \
	1392	(abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.036) + \
	1393	(abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.035) + \
	1394	(abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.026) + \
	1395	(abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.041) + \
	1396	(abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
	1397
	1398	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
	1399	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.070) + \
	1400	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.089) + \
	1401	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.104) + \
	1402	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.121) + \
	1403	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.121) + \
	1404	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.127) + \
	1405	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.127) + \
	1406	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.117) + \
	1407	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.117) + \
	1408	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.096) + \
	1409	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.099) + \
	1410	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.086) + \
	1411	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.075) + \
	1412	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.066) + \
	1413	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.056) + \
	1414	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.048) + \
	1415	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.031) + \
	1416	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.042) + \
	1417	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.030) + \
	1418	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.030) + \
	1419	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
	1420
	1421	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
	1422	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.041) + \
	1423	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.047) + \
	1424	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.059) + \
	1425	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.056) + \
	1426	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.068) + \
	1427	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.073) + \
	1428	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.080) + \
	1429	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.063) + \
	1430	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.065) + \
	1431	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.057) + \
	1432	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.049) + \
	1433	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.067) + \
	1434	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.065) + \
	1435	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.044)+ \
	1436	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.031) + \
	1437	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.026) + \
	1438	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
	1439
	1440	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
	1441	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.02) + \
	1442	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.02) + \
	1443	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.03) + \
	1444	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.03) + \
	1445	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.03) + \
	1446	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.035) + \
	1447	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.04) + \
	1448	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.03) + \
	1449	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.03) + \
	1450	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.025) + \
	1451	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.025) + \
	1452	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.035) + \
	1453	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.032) + \
	1454	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.022) + \
	1455	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.015) + \
	1456	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.013) + \
	1457	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
	1458	(abs(eta) > 3.4) * (0.000)
	1459	}
	1460
	1461
	1462	add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
	1463	(abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.528) + \
	1464	(abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.611) + \
	1465	(abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.666) + \
	1466	(abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.701) + \
	1467	(abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.734) + \
	1468	(abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.741) + \
	1469	(abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.746) + \
	1470	(abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.750) + \
	1471	(abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.742) + \
	1472	(abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.729) + \
	1473	(abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.711) + \
	1474	(abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.694) + \
	1475	(abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.684) + \
	1476	(abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.642) + \
	1477	(abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.602) + \
	1478	(abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.554) + \
	1479	(abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.488) + \
	1480	(abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.417) + \
	1481	(abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.343) + \
	1482	(abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.284) + \
	1483	(abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.250) + \
	1484	(abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.221) + \
	1485	(abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.164) + \
	1486	(abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.173) + \
	1487	(abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
	1488
	1489	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
	1490	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.322) + \
	1491	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.414) + \
	1492	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.488) + \
	1493	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.521) + \
	1494	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.554) + \
	1495	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.590) + \
	1496	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.571) + \
	1497	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.583) + \
	1498	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.575) + \
	1499	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.557) + \
	1500	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.536) + \
	1501	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.516) + \
	1502	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.483) + \
	1503	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.432) + \
	1504	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.366) + \
	1505	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.323) + \
	1506	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.241) + \
	1507	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.252) + \
	1508	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.237) + \
	1509	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.243) + \
	1510	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
	1511
	1512	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
	1513	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.200) + \
	1514	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.252) + \
	1515	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.310) + \
	1516	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.353) + \
	1517	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.377) + \
	1518	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.425) + \
	1519	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.406) + \
	1520	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.414) + \
	1521	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.412) + \
	1522	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.366) + \
	1523	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.360) + \
	1524	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.347) + \
	1525	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.310) + \
	1526	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.307)+ \
	1527	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.224) + \
	1528	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.194) + \
	1529	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
	1530
	1531	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
	1532	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.100) + \
	1533	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.12) + \
	1534	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.15) + \
	1535	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.17) + \
	1536	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.19) + \
	1537	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.21) + \
	1538	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.20) + \
	1539	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.21) + \
	1540	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.20) + \
	1541	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.18) + \
	1542	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.18) + \
	1543	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.17) + \
	1544	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.15) + \
	1545	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.15) + \
	1546	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.11) + \
	1547	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.10) + \
	1548	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
	1549	(abs(eta) > 3.4) * (0.000)
	1550	}
	1551	}
	1552
	1553
	1554
	1555	#############
	1556	# b-tagging #
	1557	#############
	1558	module BTagging BTaggingTight {
	1559
	1560	set JetInputArray JetEnergyScale/jets
	1561
	1562	set BitNumber 2
	1563
	1564	add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
	1565	(pt > 20.0) * (abs(eta) <= 3.4) * (0.001) + \
	1566	(abs(eta) > 3.4) * (0.000)
	1567	}
	1568
	1569
	1570	add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
	1571	(abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.038) + \
	1572	(abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.047) + \
	1573	(abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.046) + \
	1574	(abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.044) + \
	1575	(abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.041) + \
	1576	(abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.044) + \
	1577	(abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.041) + \
	1578	(abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.038) + \
	1579	(abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.034) + \
	1580	(abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.027) + \
	1581	(abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.022) + \
	1582	(abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.021) + \
	1583	(abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.016) + \
	1584	(abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.012) + \
	1585	(abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.009) + \
	1586	(abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.005) + \
	1587	(abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.004) + \
	1588	(abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.003) + \
	1589	(abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002) + \
	1590	(abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002) + \
	1591	(abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.001) + \
	1592	(abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002) + \
	1593	(abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.001) + \
	1594	(abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.006) + \
	1595	(abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
	1596
	1597	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
	1598	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.020) + \
	1599	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.0242) + \
	1600	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.0272) + \
	1601	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.0349) + \
	1602	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.0314) + \
	1603	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.0220) + \
	1604	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.028 ) + \
	1605	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0263) + \
	1606	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0245) + \
	1607	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.0150) + \
	1608	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0088) + \
	1609	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.0115) + \
	1610	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.0063) + \
	1611	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.0083) + \
	1612	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.0027) + \
	1613	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0052) + \
	1614	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0007) + \
	1615	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.0048) + \
	1616	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0013) + \
	1617	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.0013) + \
	1618	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
	1619
	1620	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
	1621	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.0092) + \
	1622	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.0094) + \
	1623	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.0141) + \
	1624	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.0108) + \
	1625	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.0161) + \
	1626	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.0137) + \
	1627	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.0135) + \
	1628	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.0176) + \
	1629	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.0149) + \
	1630	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.0137) + \
	1631	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.0114) + \
	1632	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.0079) + \
	1633	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.0112) + \
	1634	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.0081)+ \
	1635	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.0033) + \
	1636	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.0016) + \
	1637	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
	1638
[44771df]	1639	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
[4f728cd]	1640	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.005) + \
	1641	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.005) + \
	1642	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.005) + \
	1643	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.007) + \
	1644	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.005) + \
	1645	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.008) + \
	1646	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.006) + \
	1647	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.007) + \
	1648	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.008) + \
	1649	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.007) + \
	1650	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.006) + \
	1651	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.004) + \
	1652	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.005) + \
	1653	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.004) + \
	1654	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.001) + \
	1655	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.0005) + \
	1656	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
	1657	(abs(eta) > 3.4) * (0.000)
	1658	}
	1659
	1660
	1661
	1662	add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
	1663	(abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.349) + \
	1664	(abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.435) + \
	1665	(abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.485) + \
	1666	(abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.504) + \
	1667	(abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.519) + \
	1668	(abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.537) + \
	1669	(abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.541) + \
	1670	(abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.528) + \
	1671	(abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.514) + \
	1672	(abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.473) + \
	1673	(abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.447) + \
	1674	(abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.408) + \
	1675	(abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.373) + \
	1676	(abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.315) + \
	1677	(abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.267) + \
	1678	(abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.181) + \
	1679	(abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.138) + \
	1680	(abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.101) + \
	1681	(abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.067) + \
	1682	(abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.048) + \
	1683	(abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.032) + \
	1684	(abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.029) + \
	1685	(abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.022) + \
	1686	(abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.038) + \
	1687	(abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
	1688
	1689	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
	1690	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.186) + \
	1691	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.258) + \
	1692	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.306) + \
	1693	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.348) + \
	1694	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.359) + \
	1695	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.358) + \
	1696	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.367) + \
	1697	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.380) + \
	1698	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.359) + \
	1699	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.308) + \
	1700	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.238) + \
	1701	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.252) + \
	1702	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.229) + \
	1703	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.174) + \
	1704	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.125) + \
	1705	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.096) + \
	1706	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.066) + \
	1707	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.046) + \
	1708	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.076) + \
	1709	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.072) + \
	1710	(abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
	1711
	1712	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
	1713	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.096) + \
	1714	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.119) + \
	1715	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.156) + \
	1716	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.174) + \
	1717	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.192) + \
	1718	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.228) + \
	1719	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.218) + \
	1720	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.231) + \
	1721	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.224) + \
	1722	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.184) + \
	1723	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.164) + \
	1724	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.115) + \
	1725	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.123) + \
	1726	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.131)+ \
	1727	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.063) + \
	1728	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.039) + \
	1729	(abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
	1730
	1731	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
	1732	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.05) + \
	1733	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.06) + \
	1734	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.08) + \
	1735	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.09) + \
	1736	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.10) + \
	1737	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.11) + \
	1738	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.11) + \
	1739	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.11) + \
	1740	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.11) + \
	1741	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.09) + \
	1742	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.08) + \
	1743	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.06) + \
	1744	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.06) + \
	1745	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.06) + \
	1746	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.03) + \
	1747	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.02) + \
	1748	(abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
	1749	(abs(eta) > 3.4) * (0.000)
	1750	}
	1751	}
	1752
	1753
	1754
	1755	#############
	1756	# tau-tagging
	1757	#############
	1758
	1759
	1760	module TauTagging TauTagging {
	1761	set ParticleInputArray Delphes/allParticles
	1762	set PartonInputArray Delphes/partons
	1763	set JetInputArray JetEnergyScale/jets
	1764
	1765	set DeltaR 0.5
	1766
[5303656]	1767	set TauPTMin 20.0
[4f728cd]	1768
[5303656]	1769	set TauEtaMax 2.3
[4f728cd]	1770
	1771	# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
	1772
[5303656]	1773	add EfficiencyFormula {0} { (abs(eta) < 2.3) * ((( -0.00621816+0.00130097pt-2.19642e-5pt^2+1.49393e-7pt^3-4.58972e-10pt^4+5.27983e-13pt^5 )) (pt<250) + 0.0032*(pt>250)) + \
	1774	(abs(eta) > 2.3) * (0.000)
	1775	}
	1776	add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.970.77( (0.32 + 0.01pt - 0.000054ptpt )(pt<100)+0.78*(pt>100) ) + \
	1777	(abs(eta) > 2.3) * (0.000)
	1778	}
[4f728cd]	1779	}
	1780
	1781
	1782	##################
	1783	# ROOT tree writer
	1784	##################
	1785
	1786	module TreeWriter TreeWriter {
	1787	# add Branch InputArray BranchName BranchClass
[44771df]	1788	add Branch Delphes/allParticles Particle GenParticle
[4f728cd]	1789
	1790	add Branch GenJetFinder/jets GenJet Jet
	1791	add Branch GenMissingET/momentum GenMissingET MissingET
	1792
[44771df]	1793	add Branch HCal/eflowTracks EFlowTrack Track
	1794	add Branch ECal/eflowPhotons EFlowPhoton Tower
	1795	add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
[4f728cd]	1796
[2aaf039]	1797	add Branch PhotonEfficiency/photons Photon Photon
	1798	add Branch ElectronEfficiency/electrons Electron Electron
	1799	add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
	1800	add Branch MuonTightIdEfficiency/muons MuonTight Muon
[4f728cd]	1801
	1802	add Branch JetEnergyScale/jets Jet Jet
	1803
	1804	add Branch MissingET/momentum MissingET MissingET
	1805	add Branch ScalarHT/energy ScalarHT ScalarHT
	1806	}
	1807

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

Download in other formats: