Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_IDEA.tcl@ 0b8551f

Visit:

Last change on this file since 0b8551f was ebf40fd, checked in by Franco BEDESCHI <bed@…>, 3 years ago
Major update to handle highly displaced tracks
Property mode set to `100644`
File size: 29.2 KB

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

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

Download in other formats: