Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_IDEA.tcl@ 058667b

Visit:

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

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

Download in other formats: