Changeset 3e5f05c in git for cards/delphes_card_IDEAtrkCov.tcl
- Timestamp:
- Feb 26, 2020, 10:57:43 AM (5 years ago)
- Branches:
- ImprovedOutputFile, master
- Children:
- 95e4e44
- Parents:
- 31f3bca
- File:
-
- 1 edited
-
cards/delphes_card_IDEAtrkCov.tcl (modified) (20 diffs)
Legend:
- Unmodified
- Added
- Removed
-
cards/delphes_card_IDEAtrkCov.tcl
r31f3bca r3e5f05c 1 1 #################################################################### l 2 # FCC-ee IDEA detector model 3 # 4 # Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello 2 # FCC-ee IDEA detector model 3 # 4 # Authors: Elisa Fontanesi, Lorenzo Pezzotti, Massimiliano Antonello, Michele Selvaggi 5 5 # email: efontane@bo.infn.it, 6 # lorenzo.pezzotti01@universitadipavia.it, 7 # m.antonello@uninsubria.it, 8 ##################################################################### 9 # 6 # lorenzo.pezzotti01@universitadipavia.it, 7 # m.antonello@uninsubria.it, 8 # michele.selvaggi@cern.ch 9 ##################################################################### 10 # 10 11 ####################################### 11 12 # Order of execution of various modules … … 22 23 TrackSmearing 23 24 24 TrackMerger 25 TrackMerger 25 26 Calorimeter 26 27 EFlowMerger … … 28 29 PhotonEfficiency 29 30 PhotonIsolation 31 32 MuonFilter 30 33 31 34 ElectronFilter … … 41 44 GenJetFinder 42 45 GenMissingET 43 46 44 47 FastJetFinder 45 48 … … 97 100 } 98 101 99 # (pt <= 0.1) * (0.00) + 102 # (pt <= 0.1) * (0.00) + 100 103 # (abs(eta) <= 3.0) * (pt > 0.1) * (1.00) + 101 104 # (abs(eta) > 3) * (0.00) … … 118 121 (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0) * (0.65) + 119 122 (energy < 0.3) * (abs(eta) <= 3.0) * (0.06) 120 } 123 } 121 124 } 122 125 … … 160 163 set OutputArray tracks 161 164 162 set DetectorGeometry GeoIDEA_BASE.txt 165 ## uses https://raw.githubusercontent.com/selvaggi/FastTrackCovariance/master/GeoIDEA_BASE.txt 166 set DetectorGeometry { 167 168 1 PIPE -100 100 0.015 0.0012 0.35276 0 0 0 0 0 0 169 1 VTXLOW -0.12 0.12 0.017 0.00028 0.0937 2 0 1.5708 3e-006 3e-006 1 170 171 } 172 163 173 set Bz 2.0 164 165 174 } 166 175 … … 172 181 module Merger TrackMerger { 173 182 # add InputArray InputArray 174 add InputArray Track Covariance/tracks183 add InputArray TrackSmearing/tracks 175 184 set OutputArray tracks 176 185 } 177 186 178 187 179 ############# 180 # Calorimeter 181 ############# 188 ############# 189 # Calorimeter 190 ############# 182 191 module DualReadoutCalorimeter Calorimeter { 183 192 set ParticleInputArray ParticlePropagator/stableParticles … … 201 210 set pi [expr {acos(-1)}] 202 211 203 # Lists of the edges of each tower in eta and phi; 204 # each list starts with the lower edge of the first tower; 205 # the list ends with the higher edged of the last tower. 206 # Barrel: deta=0.02 towers up to |eta| <= 0.88 ( up to 45°) 207 # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad) 208 # Cell size: about 6 cm x 6 cm 209 210 #barrel: 212 # Lists of the edges of each tower in eta and phi; 213 # each list starts with the lower edge of the first tower; 214 # the list ends with the higher edged of the last tower. 215 # Barrel: deta=0.02 towers up to |eta| <= 0.88 ( up to 45°) 216 # Endcaps: deta=0.02 towers up to |eta| <= 3.0 (8.6° = 100 mrad) 217 # Cell size: about 6 cm x 6 cm 218 219 #barrel: 211 220 set PhiBins {} 212 221 for {set i -120} {$i <= 120} {incr i} { 213 222 add PhiBins [expr {$i * $pi/120}] 214 223 } 215 #deta=0.02 units for |eta| <= 0.88 224 #deta=0.02 units for |eta| <= 0.88 216 225 for {set i -44} {$i < 45} {incr i} { 217 226 set eta [expr {$i * 0.02}] … … 219 228 } 220 229 221 #endcaps: 230 #endcaps: 222 231 set PhiBins {} 223 232 for {set i -120} {$i <= 120} {incr i} { 224 233 add PhiBins [expr {$i* $pi/120}] 225 234 } 226 #deta=0.02 units for 0.88 < |eta| <= 3.0 227 #first, from -3.0 to -0.88 235 #deta=0.02 units for 0.88 < |eta| <= 3.0 236 #first, from -3.0 to -0.88 228 237 for {set i 1} {$i <=106} {incr i} { 229 238 set eta [expr {-3.00 + $i * 0.02}] 230 239 add EtaPhiBins $eta $PhiBins 231 240 } 232 #same for 0.88 to 3.0 241 #same for 0.88 to 3.0 233 242 for {set i 1} {$i <=106} {incr i} { 234 243 set eta [expr {0.88 + $i * 0.02}] … … 236 245 } 237 246 238 # default energy fractions {abs(PDG code)} {Fecal Fhcal} 247 # default energy fractions {abs(PDG code)} {Fecal Fhcal} 239 248 add EnergyFraction {0} {0.0 1.0} 240 # energy fractions for e, gamma and pi0 249 # energy fractions for e, gamma and pi0 241 250 add EnergyFraction {11} {1.0 0.0} 242 251 add EnergyFraction {22} {1.0 0.0} 243 252 add EnergyFraction {111} {1.0 0.0} 244 # energy fractions for muon, neutrinos and neutralinos 253 # energy fractions for muon, neutrinos and neutralinos 245 254 add EnergyFraction {12} {0.0 0.0} 246 255 add EnergyFraction {13} {0.0 0.0} … … 252 261 add EnergyFraction {1000035} {0.0 0.0} 253 262 add EnergyFraction {1000045} {0.0 0.0} 254 # energy fractions for K0short and Lambda 263 # energy fractions for K0short and Lambda 255 264 add EnergyFraction {310} {0.3 0.7} 256 265 add EnergyFraction {3122} {0.3 0.7} 257 266 258 267 259 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 268 # set ECalResolutionFormula {resolution formula as a function of eta and energy} 260 269 set ECalResolutionFormula { 261 270 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.11^2)+ … … 263 272 } 264 273 265 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 274 # set HCalResolutionFormula {resolution formula as a function of eta and energy} 266 275 set HCalResolutionFormula { 267 276 (abs(eta) <= 0.88 ) * sqrt(energy^2*0.01^2 + energy*0.30^2)+ … … 329 338 } 330 339 340 ################# 341 # Muon filter 342 ################# 343 344 module PdgCodeFilter MuonFilter { 345 set InputArray Calorimeter/eflowTracks 346 set OutputArray muons 347 set Invert true 348 add PdgCode {13} 349 add PdgCode {-13} 350 } 351 352 331 353 ##################### 332 354 # Electron efficiency … … 340 362 341 363 # efficiency formula for electrons 342 set EfficiencyFormula { 364 set EfficiencyFormula { 343 365 (energy < 2.0) * (0.000)+ 344 366 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) + … … 370 392 371 393 module Efficiency MuonEfficiency { 372 set InputArray Muon MomentumSmearing/muons394 set InputArray MuonFilter/muons 373 395 set OutputArray muons 374 396 … … 376 398 377 399 # efficiency formula for muons 378 set EfficiencyFormula { 400 set EfficiencyFormula { 379 401 (energy < 2.0) * (0.000)+ 380 402 (energy >= 2.0) * (abs(eta) <= 0.88) * (0.99) + … … 524 546 525 547 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt} 526 548 527 549 # default efficiency formula (misidentification rate) 528 550 add EfficiencyFormula {0} {0.01} … … 580 602 module TreeWriter TreeWriter { 581 603 # add Branch InputArray BranchName BranchClass 582 604 583 605 add Branch Delphes/allParticles Particle GenParticle 584 606 585 607 add Branch TrackMerger/tracks Track Track 586 608 add Branch Calorimeter/towers Tower Tower 587 609 588 610 add Branch Calorimeter/eflowTracks EFlowTrack Track 589 611 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower … … 593 615 add Branch PhotonEfficiency/photons PhotonEff Photon 594 616 add Branch PhotonIsolation/photons PhotonIso Photon 595 617 596 618 add Branch GenJetFinder/jets GenJet Jet 597 619 add Branch GenMissingET/momentum GenMissingET MissingET 598 620 599 621 add Branch UniqueObjectFinder/jets Jet Jet 600 622 add Branch UniqueObjectFinder/electrons Electron Electron 601 623 add Branch UniqueObjectFinder/photons Photon Photon 602 624 add Branch UniqueObjectFinder/muons Muon Muon 603 604 add Branch JetEnergyScale/jets AntiKtJet Jet 605 625 626 add Branch JetEnergyScale/jets AntiKtJet Jet 627 606 628 add Branch MissingET/momentum MissingET MissingET 607 629 add Branch ScalarHT/energy ScalarHT ScalarHT 608 630 } 609
Note:
See TracChangeset
for help on using the changeset viewer.
