Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_IDEA.tcl@ 952bbbc

Visit:

Last change on this file since 952bbbc was a95da74, checked in by michele <michele.selvaggi@…>, 4 years ago
added v0 of Ionisation cluster counting
Property mode set to `100644`
File size: 27.3 KB

Rev	Line
[7c9f3da]	1	#################################################################### l
[a0db751]	2	# FCC-ee IDEA detector model
	3	#
[b8a6aa3]	4	# Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi
[7c9f3da]	5	# email: efontane@bo.infn.it,
[a0db751]	6	# lorenzo.pezzotti01@universitadipavia.it,
	7	# m.antonello@uninsubria.it,
[b8a6aa3]	8	# michele.selvaggi@cern.ch
[a0db751]	9	#####################################################################
[b8a6aa3]	10
	11	set B 2.0
	12
[ef8a06d]	13	#######################################
	14	# Order of execution of various modules
	15	#######################################
	16
	17	set ExecutionPath {
	18	ParticlePropagator
	19
	20	ChargedHadronTrackingEfficiency
	21	ElectronTrackingEfficiency
	22	MuonTrackingEfficiency
	23
[a0db751]	24	TrackMergerPre
	25	TrackSmearing
[a95da74]	26	ClusterCounting
[ef8a06d]	27
[a0db751]	28	TrackMerger
[ef8a06d]	29	Calorimeter
	30	EFlowMerger
	31
	32	PhotonEfficiency
	33	PhotonIsolation
	34
[b8a6aa3]	35	MuonFilter
	36
[ef8a06d]	37	ElectronFilter
	38	ElectronEfficiency
	39	ElectronIsolation
	40
	41	MuonEfficiency
	42	MuonIsolation
	43
	44	MissingET
	45
	46	NeutrinoFilter
	47	GenJetFinder
	48	GenMissingET
[a0db751]	49
[ef8a06d]	50	FastJetFinder
	51
	52	JetEnergyScale
	53
	54	JetFlavorAssociation
	55
	56	BTagging
	57	TauTagging
	58
	59	UniqueObjectFinder
	60
	61	ScalarHT
	62	TreeWriter
	63	}
	64
	65	#################################
	66	# Propagate particles in cylinder
	67	#################################
	68
	69	module ParticlePropagator ParticlePropagator {
	70	set InputArray Delphes/stableParticles
	71
	72	set OutputArray stableParticles
	73	set ChargedHadronOutputArray chargedHadrons
	74	set ElectronOutputArray electrons
	75	set MuonOutputArray muons
	76
[7c9f3da]	77	# inner radius of the solenoid, in m
	78	set Radius 2.25
	79
	80	# half-length: z of the solenoid, in m
	81	set HalfLength 2.5
[ef8a06d]	82
[7c9f3da]	83	# magnetic field, in T
[b8a6aa3]	84	set Bz $B
[ef8a06d]	85	}
	86
	87	####################################
	88	# Charged hadron tracking efficiency
	89	####################################
	90
	91	module Efficiency ChargedHadronTrackingEfficiency {
[7c9f3da]	92	set InputArray ParticlePropagator/chargedHadrons
	93	set OutputArray chargedHadrons
	94	# We use only one efficiency, we set only 0 effincency out of eta bounds:
	95
	96	set EfficiencyFormula {
	97	(abs(eta) > 3.0) * (0.000) +
	98	(energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
	99	(energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
	100	(energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
	101	}
[ef8a06d]	102	}
	103
[a0db751]	104	# (pt <= 0.1) * (0.00) +
[7c9f3da]	105	# (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) +
	106	# (abs(eta) > 3) * (0.00)
	107
	108
	109
[ef8a06d]	110	##############################
	111	# Electron tracking efficiency
	112	##############################
	113
	114	module Efficiency ElectronTrackingEfficiency {
[7c9f3da]	115	set InputArray ParticlePropagator/electrons
	116	set OutputArray electrons
[ef8a06d]	117
	118
[7c9f3da]	119	# Current full simulation with CLICdet provides for electrons:
	120	set EfficiencyFormula {
	121	(abs(eta) > 3.0) * (0.000) +
	122	(energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
	123	(energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
	124	(energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
[a0db751]	125	}
[ef8a06d]	126	}
	127
[7c9f3da]	128
[ef8a06d]	129	##########################
	130	# Muon tracking efficiency
	131	##########################
	132
	133	module Efficiency MuonTrackingEfficiency {
[7c9f3da]	134	set InputArray ParticlePropagator/muons
	135	set OutputArray muons
	136
	137	# Current full simulation with CLICdet provides for muons:
	138	set EfficiencyFormula {
	139	(abs(eta) > 3.0) * (0.000) +
	140	(energy >= 0.5) * (abs(eta) <= 3.0) * (0.997) +
	141	(energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) +
	142	(energy < 0.3) * (abs(eta) <= 3.0) * (0.06)
	143	}
[ef8a06d]	144	}
	145
[a0db751]	146	##############
	147	# Track merger
	148	##############
[7c9f3da]	149
[a0db751]	150	module Merger TrackMergerPre {
	151	# add InputArray InputArray
	152	add InputArray ChargedHadronTrackingEfficiency/chargedHadrons
	153	add InputArray ElectronTrackingEfficiency/electrons
	154	add InputArray MuonTrackingEfficiency/muons
	155	set OutputArray tracks
	156	}
[ef8a06d]	157
	158
[a0db751]	159	########################################
	160	# Smearing for charged tracks
	161	########################################
[ef8a06d]	162
[a0db751]	163	module TrackCovariance TrackSmearing {
[b8a6aa3]	164
	165	set InputArray TrackMergerPre/tracks
	166	set OutputArray tracks
	167
[a0db751]	168	## minimum number of hits to accept a track
	169	set NMinHits 6
	170
[a95da74]	171	## magnetic field
	172	set Bz $B
	173
[a0db751]	174	## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt
	175	set DetectorGeometry {
	176
	177
	178	# Layer type 1 = R (barrel) or 2 = z (forward/backward)
	179	# Layer label
	180	# Minimum dimension z for barrel or R for forward
	181	# Maximum dimension z for barrel or R for forward
	182	# R/z location of layer
	183	# Thickness (meters)
	184	# Radiation length (meters)
	185	# Number of measurements in layers (1D or 2D)
	186	# Stereo angle (rad) - 0(pi/2) = axial(z) layer - Upper side
	187	# Stereo angle (rad) - 0(pi/2) = axial(z) layer - Lower side
	188	# Resolution Upper side (meters) - 0 = no measurement
	189	# Resolution Lower side (meters) - 0 = no measurement
	190	# measurement flag = T, scattering only = F
	191
[b8a6aa3]	192	# barrel name zmin zmax r w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
[a0db751]	193
	194	# barrel name zmin zmax r w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
	195
[b8a6aa3]	196	1 PIPE -100 100 0.015 0.001655 0.2805 0 0 0 0 0 0
[a0db751]	197	1 VTXLOW -0.12 0.12 0.017 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
	198	1 VTXLOW -0.16 0.16 0.023 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
	199	1 VTXLOW -0.16 0.16 0.031 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1
	200	1 VTXHIGH -1 1 0.32 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
	201	1 VTXHIGH -1.05 1.05 0.34 0.00047 0.0937 2 0 1.5708 7e-006 7e-006 1
	202
	203	# endcap name rmin rmax z w (m) X0 n_meas th_up (rad) th_down (rad) reso_up (m) reso_down (m) flag
	204
	205	2 VTXDSK 0.141 0.3 -0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	206	2 VTXDSK 0.138 0.3 -0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	207	2 VTXDSK 0.065 0.3 -0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	208	2 VTXDSK 0.062 0.3 -0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	209	2 VTXDSK 0.062 0.3 0.4 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	210	2 VTXDSK 0.065 0.3 0.42 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	211	2 VTXDSK 0.138 0.3 0.9 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	212	2 VTXDSK 0.141 0.3 0.92 0.00028 0.0937 2 0 1.5708 7e-006 7e-006 1
	213
	214	1 DCHCANI -2.125 2.125 0.345 0.0002 0.237223 0 0 0 0 0 0
	215	1 DCH -2 2 0.36 0.0147748 1400 1 0.0203738 0 0.0001 0 1
	216	1 DCH -2 2 0.374775 0.0147748 1400 1 -0.0212097 0 0.0001 0 1
	217	1 DCH -2 2 0.38955 0.0147748 1400 1 0.0220456 0 0.0001 0 1
	218	1 DCH -2 2 0.404324 0.0147748 1400 1 -0.0228814 0 0.0001 0 1
	219	1 DCH -2 2 0.419099 0.0147748 1400 1 0.0237172 0 0.0001 0 1
	220	1 DCH -2 2 0.433874 0.0147748 1400 1 -0.024553 0 0.0001 0 1
	221	1 DCH -2 2 0.448649 0.0147748 1400 1 0.0253888 0 0.0001 0 1
	222	1 DCH -2 2 0.463423 0.0147748 1400 1 -0.0262245 0 0.0001 0 1
	223	1 DCH -2 2 0.478198 0.0147748 1400 1 0.0270602 0 0.0001 0 1
	224	1 DCH -2 2 0.492973 0.0147748 1400 1 -0.0278958 0 0.0001 0 1
	225	1 DCH -2 2 0.507748 0.0147748 1400 1 0.0287314 0 0.0001 0 1
	226	1 DCH -2 2 0.522523 0.0147748 1400 1 -0.029567 0 0.0001 0 1
	227	1 DCH -2 2 0.537297 0.0147748 1400 1 0.0304025 0 0.0001 0 1
	228	1 DCH -2 2 0.552072 0.0147748 1400 1 -0.031238 0 0.0001 0 1
	229	1 DCH -2 2 0.566847 0.0147748 1400 1 0.0320734 0 0.0001 0 1
	230	1 DCH -2 2 0.581622 0.0147748 1400 1 -0.0329088 0 0.0001 0 1
	231	1 DCH -2 2 0.596396 0.0147748 1400 1 0.0337442 0 0.0001 0 1
	232	1 DCH -2 2 0.611171 0.0147748 1400 1 -0.0345795 0 0.0001 0 1
	233	1 DCH -2 2 0.625946 0.0147748 1400 1 0.0354147 0 0.0001 0 1
	234	1 DCH -2 2 0.640721 0.0147748 1400 1 -0.0362499 0 0.0001 0 1
	235	1 DCH -2 2 0.655495 0.0147748 1400 1 0.0370851 0 0.0001 0 1
	236	1 DCH -2 2 0.67027 0.0147748 1400 1 -0.0379202 0 0.0001 0 1
	237	1 DCH -2 2 0.685045 0.0147748 1400 1 0.0387552 0 0.0001 0 1
	238	1 DCH -2 2 0.69982 0.0147748 1400 1 -0.0395902 0 0.0001 0 1
	239	1 DCH -2 2 0.714595 0.0147748 1400 1 0.0404252 0 0.0001 0 1
	240	1 DCH -2 2 0.729369 0.0147748 1400 1 -0.04126 0 0.0001 0 1
	241	1 DCH -2 2 0.744144 0.0147748 1400 1 0.0420949 0 0.0001 0 1
	242	1 DCH -2 2 0.758919 0.0147748 1400 1 -0.0429296 0 0.0001 0 1
	243	1 DCH -2 2 0.773694 0.0147748 1400 1 0.0437643 0 0.0001 0 1
	244	1 DCH -2 2 0.788468 0.0147748 1400 1 -0.044599 0 0.0001 0 1
	245	1 DCH -2 2 0.803243 0.0147748 1400 1 0.0454336 0 0.0001 0 1
	246	1 DCH -2 2 0.818018 0.0147748 1400 1 -0.0462681 0 0.0001 0 1
	247	1 DCH -2 2 0.832793 0.0147748 1400 1 0.0471025 0 0.0001 0 1
	248	1 DCH -2 2 0.847568 0.0147748 1400 1 -0.0479369 0 0.0001 0 1
	249	1 DCH -2 2 0.862342 0.0147748 1400 1 0.0487713 0 0.0001 0 1
	250	1 DCH -2 2 0.877117 0.0147748 1400 1 -0.0496055 0 0.0001 0 1
	251	1 DCH -2 2 0.891892 0.0147748 1400 1 0.0504397 0 0.0001 0 1
	252	1 DCH -2 2 0.906667 0.0147748 1400 1 -0.0512738 0 0.0001 0 1
	253	1 DCH -2 2 0.921441 0.0147748 1400 1 0.0521079 0 0.0001 0 1
	254	1 DCH -2 2 0.936216 0.0147748 1400 1 -0.0529418 0 0.0001 0 1
	255	1 DCH -2 2 0.950991 0.0147748 1400 1 0.0537757 0 0.0001 0 1
	256	1 DCH -2 2 0.965766 0.0147748 1400 1 -0.0546095 0 0.0001 0 1
	257	1 DCH -2 2 0.980541 0.0147748 1400 1 0.0554433 0 0.0001 0 1
	258	1 DCH -2 2 0.995315 0.0147748 1400 1 -0.056277 0 0.0001 0 1
	259	1 DCH -2 2 1.01009 0.0147748 1400 1 0.0571106 0 0.0001 0 1
	260	1 DCH -2 2 1.02486 0.0147748 1400 1 -0.0579441 0 0.0001 0 1
	261	1 DCH -2 2 1.03964 0.0147748 1400 1 0.0587775 0 0.0001 0 1
	262	1 DCH -2 2 1.05441 0.0147748 1400 1 -0.0596108 0 0.0001 0 1
	263	1 DCH -2 2 1.06919 0.0147748 1400 1 0.0604441 0 0.0001 0 1
	264	1 DCH -2 2 1.08396 0.0147748 1400 1 -0.0612773 0 0.0001 0 1
	265	1 DCH -2 2 1.09874 0.0147748 1400 1 0.0621104 0 0.0001 0 1
	266	1 DCH -2 2 1.11351 0.0147748 1400 1 -0.0629434 0 0.0001 0 1
	267	1 DCH -2 2 1.12829 0.0147748 1400 1 0.0637763 0 0.0001 0 1
	268	1 DCH -2 2 1.14306 0.0147748 1400 1 -0.0646092 0 0.0001 0 1
	269	1 DCH -2 2 1.15784 0.0147748 1400 1 0.0654419 0 0.0001 0 1
	270	1 DCH -2 2 1.17261 0.0147748 1400 1 -0.0662746 0 0.0001 0 1
	271	1 DCH -2 2 1.18739 0.0147748 1400 1 0.0671071 0 0.0001 0 1
	272	1 DCH -2 2 1.20216 0.0147748 1400 1 -0.0679396 0 0.0001 0 1
	273	1 DCH -2 2 1.21694 0.0147748 1400 1 0.068772 0 0.0001 0 1
	274	1 DCH -2 2 1.23171 0.0147748 1400 1 -0.0696042 0 0.0001 0 1
	275	1 DCH -2 2 1.24649 0.0147748 1400 1 0.0704364 0 0.0001 0 1
	276	1 DCH -2 2 1.26126 0.0147748 1400 1 -0.0712685 0 0.0001 0 1
	277	1 DCH -2 2 1.27604 0.0147748 1400 1 0.0721005 0 0.0001 0 1
	278	1 DCH -2 2 1.29081 0.0147748 1400 1 -0.0729324 0 0.0001 0 1
	279	1 DCH -2 2 1.30559 0.0147748 1400 1 0.0737642 0 0.0001 0 1
	280	1 DCH -2 2 1.32036 0.0147748 1400 1 -0.0745958 0 0.0001 0 1
	281	1 DCH -2 2 1.33514 0.0147748 1400 1 0.0754274 0 0.0001 0 1
	282	1 DCH -2 2 1.34991 0.0147748 1400 1 -0.0762589 0 0.0001 0 1
	283	1 DCH -2 2 1.36468 0.0147748 1400 1 0.0770903 0 0.0001 0 1
	284	1 DCH -2 2 1.37946 0.0147748 1400 1 -0.0779215 0 0.0001 0 1
	285	1 DCH -2 2 1.39423 0.0147748 1400 1 0.0787527 0 0.0001 0 1
	286	1 DCH -2 2 1.40901 0.0147748 1400 1 -0.0795837 0 0.0001 0 1
	287	1 DCH -2 2 1.42378 0.0147748 1400 1 0.0804147 0 0.0001 0 1
	288	1 DCH -2 2 1.43856 0.0147748 1400 1 -0.0812455 0 0.0001 0 1
	289	1 DCH -2 2 1.45333 0.0147748 1400 1 0.0820762 0 0.0001 0 1
	290	1 DCH -2 2 1.46811 0.0147748 1400 1 -0.0829068 0 0.0001 0 1
	291	1 DCH -2 2 1.48288 0.0147748 1400 1 0.0837373 0 0.0001 0 1
	292	1 DCH -2 2 1.49766 0.0147748 1400 1 -0.0845677 0 0.0001 0 1
	293	1 DCH -2 2 1.51243 0.0147748 1400 1 0.0853979 0 0.0001 0 1
	294	1 DCH -2 2 1.52721 0.0147748 1400 1 -0.086228 0 0.0001 0 1
	295	1 DCH -2 2 1.54198 0.0147748 1400 1 0.087058 0 0.0001 0 1
	296	1 DCH -2 2 1.55676 0.0147748 1400 1 -0.0878879 0 0.0001 0 1
	297	1 DCH -2 2 1.57153 0.0147748 1400 1 0.0887177 0 0.0001 0 1
	298	1 DCH -2 2 1.58631 0.0147748 1400 1 -0.0895474 0 0.0001 0 1
	299	1 DCH -2 2 1.60108 0.0147748 1400 1 0.0903769 0 0.0001 0 1
	300	1 DCH -2 2 1.61586 0.0147748 1400 1 -0.0912063 0 0.0001 0 1
	301	1 DCH -2 2 1.63063 0.0147748 1400 1 0.0920356 0 0.0001 0 1
	302	1 DCH -2 2 1.64541 0.0147748 1400 1 -0.0928647 0 0.0001 0 1
	303	1 DCH -2 2 1.66018 0.0147748 1400 1 0.0936937 0 0.0001 0 1
	304	1 DCH -2 2 1.67495 0.0147748 1400 1 -0.0945226 0 0.0001 0 1
	305	1 DCH -2 2 1.68973 0.0147748 1400 1 0.0953514 0 0.0001 0 1
	306	1 DCH -2 2 1.7045 0.0147748 1400 1 -0.09618 0 0.0001 0 1
	307	1 DCH -2 2 1.71928 0.0147748 1400 1 0.0970085 0 0.0001 0 1
	308	1 DCH -2 2 1.73405 0.0147748 1400 1 -0.0978369 0 0.0001 0 1
	309	1 DCH -2 2 1.74883 0.0147748 1400 1 0.0986651 0 0.0001 0 1
	310	1 DCH -2 2 1.7636 0.0147748 1400 1 -0.0994932 0 0.0001 0 1
	311	1 DCH -2 2 1.77838 0.0147748 1400 1 0.100321 0 0.0001 0 1
	312	1 DCH -2 2 1.79315 0.0147748 1400 1 -0.101149 0 0.0001 0 1
	313	1 DCH -2 2 1.80793 0.0147748 1400 1 0.101977 0 0.0001 0 1
	314	1 DCH -2 2 1.8227 0.0147748 1400 1 -0.102804 0 0.0001 0 1
	315	1 DCH -2 2 1.83748 0.0147748 1400 1 0.103632 0 0.0001 0 1
	316	1 DCH -2 2 1.85225 0.0147748 1400 1 -0.104459 0 0.0001 0 1
	317	1 DCH -2 2 1.86703 0.0147748 1400 1 0.105286 0 0.0001 0 1
	318	1 DCH -2 2 1.8818 0.0147748 1400 1 -0.106113 0 0.0001 0 1
	319	1 DCH -2 2 1.89658 0.0147748 1400 1 0.10694 0 0.0001 0 1
	320	1 DCH -2 2 1.91135 0.0147748 1400 1 -0.107766 0 0.0001 0 1
	321	1 DCH -2 2 1.92613 0.0147748 1400 1 0.108593 0 0.0001 0 1
	322	1 DCH -2 2 1.9409 0.0147748 1400 1 -0.109419 0 0.0001 0 1
	323	1 DCH -2 2 1.95568 0.0147748 1400 1 0.110246 0 0.0001 0 1
	324	1 DCH -2 2 1.97045 0.0147748 1400 1 -0.111072 0 0.0001 0 1
	325	1 DCH -2 2 1.98523 0.0147748 1400 1 0.111898 0 0.0001 0 1
	326	1 DCH -2 2 2 0.0147748 1400 1 -0.112723 0 0.0001 0 1
	327	1 DCHCANO -2.125 2.125 2.02 0.02 1.667 0 0 0 0 0 0
	328	1 BSILWRP -2.35 2.35 2.04 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
	329	1 BSILWRP -2.35 2.35 2.06 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
	330	1 MAG -2.5 2.5 2.25 0.05 0.0658 0 0 0 0 0 0
	331	1 BPRESH -2.55 2.55 2.45 0.02 1 2 0 1.5708 7e-005 0.01 1
	332	2 DCHWALL 0.345 2.02 2.125 0.25 5.55 0 0 0 0 0 0
	333	2 DCHWALL 0.345 2.02 -2.125 0.25 5.55 0 0 0 0 0 0
	334	2 FSILWRP 0.354 2.02 -2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
	335	2 FSILWRP 0.35 2.02 -2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
	336	2 FSILWRP 0.35 2.02 2.3 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
	337	2 FSILWRP 0.354 2.02 2.32 0.00047 0.0937 2 0 1.5708 7e-006 9e-005 1
	338	2 FRAD 0.38 2.09 2.49 0.0043 0.005612 0 0 0 0 0 0
	339	2 FRAD 0.38 2.09 -2.49 0.0043 0.005612 0 0 0 0 0 0
	340	2 FPRESH 0.39 2.43 -2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
	341	2 FPRESH 0.39 2.43 2.55 0.02 1 2 0 1.5708 7e-005 0.01 1
[7c9f3da]	342	}
[ef8a06d]	343
	344	}
	345
[a95da74]	346	###################
	347	# Cluster Counting
	348	###################
	349
	350	module ClusterCounting ClusterCounting {
	351
	352	add InputArray TrackSmearing/tracks
	353	set OutputArray tracks
	354
	355	set Bz $B
	356
	357	set Rmin 0.35
	358	set Rmax 2.0
	359	set Zmin -2.0
	360	set Zmax 2.0
	361
	362	# gas mix option: 0
	363	# 0: Helium 90% - Isobutane 10%
	364	# 1: Helium 100%
	365	# 2: Argon 50% - Ethane 50%
	366	# 3: Argon 100%
	367
	368	set GasOption 0
	369
	370	}
	371
	372
	373
[ef8a06d]	374	##############
	375	# Track merger
	376	##############
	377
	378	module Merger TrackMerger {
	379	# add InputArray InputArray
[a95da74]	380	add InputArray ClusterCounting/tracks
[ef8a06d]	381	set OutputArray tracks
	382	}
	383
	384
[a0db751]	385	#############
	386	# Calorimeter
	387	#############
[ef8a06d]	388	module DualReadoutCalorimeter Calorimeter {
	389	set ParticleInputArray ParticlePropagator/stableParticles
	390	set TrackInputArray TrackMerger/tracks
	391
	392	set TowerOutputArray towers
	393	set PhotonOutputArray photons
	394
	395	set EFlowTrackOutputArray eflowTracks
	396	set EFlowPhotonOutputArray eflowPhotons
	397	set EFlowNeutralHadronOutputArray eflowNeutralHadrons
	398
	399	set ECalEnergyMin 0.5
[7c9f3da]	400	set HCalEnergyMin 0.5
[a1b19ea]	401	set EnergyMin 0.5
[ef8a06d]	402	set ECalEnergySignificanceMin 1.0
	403	set HCalEnergySignificanceMin 1.0
	404	set EnergySignificanceMin 1.0
	405
	406	set SmearTowerCenter true
[7c9f3da]	407	set pi [expr {acos(-1)}]
	408
[a0db751]	409	# Lists of the edges of each tower in eta and phi;
	410	# each list starts with the lower edge of the first tower;
	411	# the list ends with the higher edged of the last tower.
	412	# Barrel: deta=0.02 towers up to \|eta\| <= 0.88 ( up to 45Â°)
	413	# Endcaps: deta=0.02 towers up to \|eta\| <= 3.0 (8.6Â° = 100 mrad)
	414	# Cell size: about 6 cm x 6 cm
[7c9f3da]	415
[a0db751]	416	#barrel:
[7c9f3da]	417	set PhiBins {}
	418	for {set i -120} {$i <= 120} {incr i} {
	419	add PhiBins [expr {$i * $pi/120}]
	420	}
[a0db751]	421	#deta=0.02 units for \|eta\| <= 0.88
[7c9f3da]	422	for {set i -44} {$i < 45} {incr i} {
	423	set eta [expr {$i * 0.02}]
	424	add EtaPhiBins $eta $PhiBins
	425	}
	426
[a0db751]	427	#endcaps:
[7c9f3da]	428	set PhiBins {}
	429	for {set i -120} {$i <= 120} {incr i} {
	430	add PhiBins [expr {$i* $pi/120}]
	431	}
[a0db751]	432	#deta=0.02 units for 0.88 < \|eta\| <= 3.0
	433	#first, from -3.0 to -0.88
[7c9f3da]	434	for {set i 1} {$i <=106} {incr i} {
	435	set eta [expr {-3.00 + $i * 0.02}]
	436	add EtaPhiBins $eta $PhiBins
	437	}
[a0db751]	438	#same for 0.88 to 3.0
[7c9f3da]	439	for {set i 1} {$i <=106} {incr i} {
	440	set eta [expr {0.88 + $i * 0.02}]
	441	add EtaPhiBins $eta $PhiBins
	442	}
	443
[a0db751]	444	# default energy fractions {abs(PDG code)} {Fecal Fhcal}
[7c9f3da]	445	add EnergyFraction {0} {0.0 1.0}
[a0db751]	446	# energy fractions for e, gamma and pi0
[7c9f3da]	447	add EnergyFraction {11} {1.0 0.0}
	448	add EnergyFraction {22} {1.0 0.0}
	449	add EnergyFraction {111} {1.0 0.0}
[a0db751]	450	# energy fractions for muon, neutrinos and neutralinos
[7c9f3da]	451	add EnergyFraction {12} {0.0 0.0}
	452	add EnergyFraction {13} {0.0 0.0}
	453	add EnergyFraction {14} {0.0 0.0}
	454	add EnergyFraction {16} {0.0 0.0}
	455	add EnergyFraction {1000022} {0.0 0.0}
	456	add EnergyFraction {1000023} {0.0 0.0}
	457	add EnergyFraction {1000025} {0.0 0.0}
	458	add EnergyFraction {1000035} {0.0 0.0}
	459	add EnergyFraction {1000045} {0.0 0.0}
[a0db751]	460	# energy fractions for K0short and Lambda
[7c9f3da]	461	add EnergyFraction {310} {0.3 0.7}
	462	add EnergyFraction {3122} {0.3 0.7}
	463
	464
[a0db751]	465	# set ECalResolutionFormula {resolution formula as a function of eta and energy}
[7c9f3da]	466	set ECalResolutionFormula {
	467	(abs(eta) <= 0.88 ) * sqrt(energy^20.01^2 + energy0.11^2)+
	468	(abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^20.01^2 + energy0.11^2)
	469	}
	470
[a0db751]	471	# set HCalResolutionFormula {resolution formula as a function of eta and energy}
[7c9f3da]	472	set HCalResolutionFormula {
	473	(abs(eta) <= 0.88 ) * sqrt(energy^20.01^2 + energy0.30^2)+
	474	(abs(eta) > 0.88 && abs(eta) <= 3.0) * sqrt(energy^20.01^2 + energy0.30^2)
	475	}
[ef8a06d]	476	}
	477
	478	####################
	479	# Energy flow merger
	480	####################
	481
	482	module Merger EFlowMerger {
	483	# add InputArray InputArray
	484	add InputArray Calorimeter/eflowTracks
	485	add InputArray Calorimeter/eflowPhotons
	486	add InputArray Calorimeter/eflowNeutralHadrons
	487	set OutputArray eflow
	488	}
	489
	490	###################
	491	# Photon efficiency
	492	###################
	493
	494	module Efficiency PhotonEfficiency {
	495	set InputArray Calorimeter/eflowPhotons
	496	set OutputArray photons
	497
	498	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	499	# efficiency formula for photons
[7c9f3da]	500	set EfficiencyFormula {
	501	(energy < 2.0) * (0.000)+
	502	(energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
	503	(energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
	504	(abs(eta) > 3.0) * (0.000)
	505	}
[ef8a06d]	506	}
	507
	508	##################
	509	# Photon isolation
	510	##################
	511
	512	module Isolation PhotonIsolation {
	513	set CandidateInputArray PhotonEfficiency/photons
	514	set IsolationInputArray EFlowMerger/eflow
	515
	516	set OutputArray photons
	517
	518	set DeltaRMax 0.5
	519
	520	set PTMin 0.5
	521
[7c9f3da]	522	set PTRatioMax 999.
[ef8a06d]	523	}
	524
	525	#################
	526	# Electron filter
	527	#################
	528
	529	module PdgCodeFilter ElectronFilter {
	530	set InputArray Calorimeter/eflowTracks
	531	set OutputArray electrons
	532	set Invert true
	533	add PdgCode {11}
	534	add PdgCode {-11}
	535	}
	536
[a0db751]	537	#################
	538	# Muon filter
	539	#################
	540
	541	module PdgCodeFilter MuonFilter {
	542	set InputArray Calorimeter/eflowTracks
	543	set OutputArray muons
	544	set Invert true
	545	add PdgCode {13}
	546	add PdgCode {-13}
	547	}
	548
	549
[ef8a06d]	550	#####################
	551	# Electron efficiency
	552	#####################
	553
	554	module Efficiency ElectronEfficiency {
	555	set InputArray ElectronFilter/electrons
	556	set OutputArray electrons
	557
	558	# set EfficiencyFormula {efficiency formula as a function of eta and pt}
	559
	560	# efficiency formula for electrons
[a0db751]	561	set EfficiencyFormula {
[7c9f3da]	562	(energy < 2.0) * (0.000)+
	563	(energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
	564	(energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
	565	(abs(eta) > 3.0) * (0.000)
	566	}
[ef8a06d]	567	}
	568
	569	####################
	570	# Electron isolation
	571	####################
	572
	573	module Isolation ElectronIsolation {
	574	set CandidateInputArray ElectronEfficiency/electrons
	575	set IsolationInputArray EFlowMerger/eflow
	576
	577	set OutputArray electrons
	578
	579	set DeltaRMax 0.5
	580
	581	set PTMin 0.5
	582
	583	set PTRatioMax 0.12
	584	}
	585
	586	#################
	587	# Muon efficiency
	588	#################
	589
	590	module Efficiency MuonEfficiency {
[a0db751]	591	set InputArray MuonFilter/muons
[ef8a06d]	592	set OutputArray muons
	593
	594	# set EfficiencyFormula {efficiency as a function of eta and pt}
	595
	596	# efficiency formula for muons
[a0db751]	597	set EfficiencyFormula {
[7c9f3da]	598	(energy < 2.0) * (0.000)+
	599	(energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) +
	600	(energy >= 2.0) * (abs(eta) >0.88 && abs(eta) <= 3.0) * (0.99) +
	601	(abs(eta) > 3.0) * (0.000)
	602	}
[ef8a06d]	603	}
	604
	605	################
	606	# Muon isolation
	607	################
	608
	609	module Isolation MuonIsolation {
	610	set CandidateInputArray MuonEfficiency/muons
	611	set IsolationInputArray EFlowMerger/eflow
	612
	613	set OutputArray muons
	614
	615	set DeltaRMax 0.5
	616
	617	set PTMin 0.5
	618
	619	set PTRatioMax 0.25
	620	}
	621
	622	###################
	623	# Missing ET merger
	624	###################
	625
	626	module Merger MissingET {
	627	# add InputArray InputArray
	628	add InputArray EFlowMerger/eflow
	629	set MomentumOutputArray momentum
	630	}
	631
	632	##################
	633	# Scalar HT merger
	634	##################
	635
	636	module Merger ScalarHT {
	637	# add InputArray InputArray
	638	add InputArray UniqueObjectFinder/jets
	639	add InputArray UniqueObjectFinder/electrons
	640	add InputArray UniqueObjectFinder/photons
	641	add InputArray UniqueObjectFinder/muons
	642	set EnergyOutputArray energy
	643	}
	644
	645	#####################
	646	# Neutrino Filter
	647	#####################
	648
	649	module PdgCodeFilter NeutrinoFilter {
	650
	651	set InputArray Delphes/stableParticles
	652	set OutputArray filteredParticles
	653
	654	set PTMin 0.0
	655
	656	add PdgCode {12}
	657	add PdgCode {14}
	658	add PdgCode {16}
	659	add PdgCode {-12}
	660	add PdgCode {-14}
	661	add PdgCode {-16}
	662	}
	663
	664
	665	#####################
	666	# MC truth jet finder
	667	#####################
	668
	669	module FastJetFinder GenJetFinder {
	670	set InputArray NeutrinoFilter/filteredParticles
	671
	672	set OutputArray jets
	673
	674	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	675	set JetAlgorithm 6
[7c9f3da]	676	set ParameterR 0.4
	677	set JetPTMin 1.0
[ef8a06d]	678	}
	679
[7c9f3da]	680
[ef8a06d]	681	#########################
	682	# Gen Missing ET merger
	683	########################
	684
	685	module Merger GenMissingET {
	686	# add InputArray InputArray
	687	add InputArray NeutrinoFilter/filteredParticles
	688	set MomentumOutputArray momentum
	689	}
	690
	691	############
	692	# Jet finder
	693	############
	694
	695	module FastJetFinder FastJetFinder {
	696	# set InputArray Calorimeter/towers
	697	set InputArray EFlowMerger/eflow
	698
	699	set OutputArray jets
	700
	701	# algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
	702	set JetAlgorithm 6
[7c9f3da]	703	set ParameterR 0.4
	704	set JetPTMin 1.0
[ef8a06d]	705	}
	706
	707	##################
	708	# Jet Energy Scale
	709	##################
	710
	711	module EnergyScale JetEnergyScale {
	712	set InputArray FastJetFinder/jets
	713	set OutputArray jets
	714
	715	# scale formula for jets
[7c9f3da]	716	set ScaleFormula {1.08}
[ef8a06d]	717	}
	718
	719	########################
	720	# Jet Flavor Association
	721	########################
	722
	723	module JetFlavorAssociation JetFlavorAssociation {
	724
	725	set PartonInputArray Delphes/partons
	726	set ParticleInputArray Delphes/allParticles
	727	set ParticleLHEFInputArray Delphes/allParticlesLHEF
	728	set JetInputArray JetEnergyScale/jets
	729
	730	set DeltaR 0.5
	731	set PartonPTMin 1.0
[7c9f3da]	732	set PartonEtaMax 3.0
[ef8a06d]	733	}
	734
	735	###########
	736	# b-tagging
	737	###########
	738
	739	module BTagging BTagging {
	740	set JetInputArray JetEnergyScale/jets
	741
	742	set BitNumber 0
	743
	744	# add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
[a0db751]	745
[ef8a06d]	746	# default efficiency formula (misidentification rate)
[7c9f3da]	747	add EfficiencyFormula {0} {0.01}
[ef8a06d]	748
	749	# efficiency formula for c-jets (misidentification rate)
[7c9f3da]	750	add EfficiencyFormula {4} {0.10}
[ef8a06d]	751
	752	# efficiency formula for b-jets
[7c9f3da]	753	add EfficiencyFormula {5} {0.80}
[ef8a06d]	754	}
	755
	756	#############
	757	# tau-tagging
	758	#############
	759
	760	module TauTagging TauTagging {
	761	set ParticleInputArray Delphes/allParticles
	762	set PartonInputArray Delphes/partons
	763	set JetInputArray JetEnergyScale/jets
	764
	765	set DeltaR 0.5
	766	set TauPTMin 1.0
[7c9f3da]	767	set TauEtaMax 3.0
[ef8a06d]	768
	769	# default efficiency formula (misidentification rate)
[7c9f3da]	770	add EfficiencyFormula {0} {0.001}
[ef8a06d]	771	# efficiency formula for tau-jets
	772	add EfficiencyFormula {15} {0.6}
	773	}
	774
[7c9f3da]	775
[ef8a06d]	776	#####################################################
	777	# Find uniquely identified photons/electrons/tau/jets
	778	#####################################################
	779
	780	module UniqueObjectFinder UniqueObjectFinder {
	781	# earlier arrays take precedence over later ones
	782	# add InputArray InputArray OutputArray
	783	add InputArray PhotonIsolation/photons photons
	784	add InputArray ElectronIsolation/electrons electrons
	785	add InputArray MuonIsolation/muons muons
	786	add InputArray JetEnergyScale/jets jets
	787	}
	788
[7c9f3da]	789
	790
[ef8a06d]	791	##################
	792	# ROOT tree writer
	793	##################
	794
[7c9f3da]	795	# Tracks, towers and eflow objects are not stored by default in the output.
	796	# If needed (for jet constituent or other studies), uncomment the relevant
[ef8a06d]	797	# "add Branch ..." lines.
	798
	799	module TreeWriter TreeWriter {
[7c9f3da]	800	# add Branch InputArray BranchName BranchClass
[a0db751]	801
[7c9f3da]	802	add Branch Delphes/allParticles Particle GenParticle
	803
	804	add Branch TrackMerger/tracks Track Track
	805	add Branch Calorimeter/towers Tower Tower
[a0db751]	806
[7c9f3da]	807	add Branch Calorimeter/eflowTracks EFlowTrack Track
	808	add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
	809	add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
	810
	811	add Branch Calorimeter/photons CaloPhoton Photon
	812	add Branch PhotonEfficiency/photons PhotonEff Photon
	813	add Branch PhotonIsolation/photons PhotonIso Photon
[a0db751]	814
[7c9f3da]	815	add Branch GenJetFinder/jets GenJet Jet
	816	add Branch GenMissingET/momentum GenMissingET MissingET
[a0db751]	817
[7c9f3da]	818	add Branch UniqueObjectFinder/jets Jet Jet
	819	add Branch UniqueObjectFinder/electrons Electron Electron
	820	add Branch UniqueObjectFinder/photons Photon Photon
	821	add Branch UniqueObjectFinder/muons Muon Muon
[a0db751]	822
	823	add Branch JetEnergyScale/jets AntiKtJet Jet
	824
[7c9f3da]	825	add Branch MissingET/momentum MissingET MissingET
	826	add Branch ScalarHT/energy ScalarHT ScalarHT
[b8a6aa3]	827
	828	# add Info InfoName InfoValue
	829	add Info Bz $B
[ef8a06d]	830	}

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

Download in other formats: