Fork me on GitHub

Context Navigation

source: git/cards/delphes_card_IDEA.tcl@ 83e77ee

Visit:

Last change on this file since 83e77ee was 83e77ee, checked in by michele <michele.selvaggi@…>, 3 years ago
added neutral time of flight
Property mode set to `100644`
File size: 29.1 KB

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

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

Download in other formats: