Fork me on GitHub

Changeset 3e5f05c in git


Ignore:
Timestamp:
Feb 26, 2020, 10:57:43 AM (5 years ago)
Author:
michele <michele.selvaggi@…>
Branches:
ImprovedOutputFile, master
Children:
95e4e44
Parents:
31f3bca
Message:

updated IDA card with geometry in delphes card

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_IDEAtrkCov.tcl

    r31f3bca r3e5f05c  
    11####################################################################                                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
    55# 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#
    1011#######################################
    1112# Order of execution of various modules
    … …  
    2223  TrackSmearing
    2324
    24   TrackMerger 
     25  TrackMerger
    2526  Calorimeter
    2627  EFlowMerger
    … …  
    2829  PhotonEfficiency
    2930  PhotonIsolation
     31
     32  MuonFilter
    3033
    3134  ElectronFilter
    … …  
    4144  GenJetFinder
    4245  GenMissingET
    43  
     46
    4447  FastJetFinder
    4548
    … …  
    97100}
    98101
    99 #       (pt <= 0.1)                                     * (0.00) + 
     102#       (pt <= 0.1)                                     * (0.00) +
    100103#       (abs(eta) <= 3.0)               * (pt > 0.1)    * (1.00) +
    101104#       (abs(eta) > 3)                                  * (0.00)
    … …  
    118121        (energy < 0.5 && energy >= 0.3) * (abs(eta) <= 3.0)    * (0.65) +
    119122        (energy < 0.3) * (abs(eta) <= 3.0)                     * (0.06)
    120     } 
     123    }
    121124}
    122125
    … …  
    160163    set OutputArray tracks
    161164
    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
    163173    set Bz 2.0
    164 
    165174}
    166175
    … …  
    172181module Merger TrackMerger {
    173182# add InputArray InputArray
    174   add InputArray TrackCovariance/tracks
     183  add InputArray TrackSmearing/tracks
    175184  set OutputArray tracks
    176185}
    177186
    178187
    179 #############                                                                                                                         
    180 # Calorimeter                                                                                                                                           
    181 #############                                                                                                                                           
     188#############
     189# Calorimeter
     190#############
    182191module DualReadoutCalorimeter Calorimeter {
    183192  set ParticleInputArray ParticlePropagator/stableParticles
    … …  
    201210    set pi [expr {acos(-1)}]
    202211
    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:
    211220    set PhiBins {}
    212221    for {set i -120} {$i <= 120} {incr i} {
    213222        add PhiBins [expr {$i * $pi/120}]
    214223    }
    215     #deta=0.02 units for |eta| <= 0.88                                                             
     224    #deta=0.02 units for |eta| <= 0.88
    216225    for {set i -44} {$i < 45} {incr i} {
    217226        set eta [expr {$i * 0.02}]
    … …  
    219228    }
    220229
    221     #endcaps:                                                                                       
     230    #endcaps:
    222231    set PhiBins {}
    223232    for {set i -120} {$i <= 120} {incr i} {
    224233        add PhiBins [expr {$i* $pi/120}]
    225234    }
    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
    228237    for {set i 1} {$i <=106} {incr i} {
    229238        set eta [expr {-3.00 + $i * 0.02}]
    230239        add EtaPhiBins $eta $PhiBins
    231240    }
    232     #same for 0.88 to 3.0                                                                           
     241    #same for 0.88 to 3.0
    233242    for  {set i 1} {$i <=106} {incr i} {
    234243        set eta [expr {0.88 + $i * 0.02}]
    … …  
    236245    }
    237246
    238     # default energy fractions {abs(PDG code)} {Fecal Fhcal}                                                                                 
     247    # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    239248    add EnergyFraction {0} {0.0 1.0}
    240     # energy fractions for e, gamma and pi0                                                                                                           
     249    # energy fractions for e, gamma and pi0
    241250    add EnergyFraction {11} {1.0 0.0}
    242251    add EnergyFraction {22} {1.0 0.0}
    243252    add EnergyFraction {111} {1.0 0.0}
    244     # energy fractions for muon, neutrinos and neutralinos                                                                             
     253    # energy fractions for muon, neutrinos and neutralinos
    245254    add EnergyFraction {12} {0.0 0.0}
    246255    add EnergyFraction {13} {0.0 0.0}
    … …  
    252261    add EnergyFraction {1000035} {0.0 0.0}
    253262    add EnergyFraction {1000045} {0.0 0.0}
    254     # energy fractions for K0short and Lambda                                                                                                           
     263    # energy fractions for K0short and Lambda
    255264    add EnergyFraction {310} {0.3 0.7}
    256265    add EnergyFraction {3122} {0.3 0.7}
    257266
    258267
    259     # set ECalResolutionFormula {resolution formula as a function of eta and energy}                               
     268    # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    260269    set ECalResolutionFormula {
    261270    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.11^2)+
    … …  
    263272    }
    264273
    265     # set HCalResolutionFormula {resolution formula as a function of eta and energy}                                               
     274    # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    266275    set HCalResolutionFormula {
    267276    (abs(eta) <= 0.88 )                     * sqrt(energy^2*0.01^2 + energy*0.30^2)+
    … …  
    329338}
    330339
     340#################
     341# Muon filter
     342#################
     343
     344module 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
    331353#####################
    332354# Electron efficiency
    … …  
    340362
    341363  # efficiency formula for electrons
    342   set EfficiencyFormula {         
     364  set EfficiencyFormula {
    343365        (energy < 2.0)                                         * (0.000)+
    344366        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
    … …  
    370392
    371393module Efficiency MuonEfficiency {
    372   set InputArray MuonMomentumSmearing/muons
     394  set InputArray MuonFilter/muons
    373395  set OutputArray muons
    374396
    … …  
    376398
    377399  # efficiency formula for muons
    378   set EfficiencyFormula {                                   
     400  set EfficiencyFormula {
    379401        (energy < 2.0)                                         * (0.000)+
    380402        (energy >= 2.0) * (abs(eta) <= 0.88)                   * (0.99) +
    … …  
    524546
    525547  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    526  
     548
    527549  # default efficiency formula (misidentification rate)
    528550  add EfficiencyFormula {0} {0.01}
    … …  
    580602module TreeWriter TreeWriter {
    581603    # add Branch InputArray BranchName BranchClass
    582    
     604
    583605    add Branch Delphes/allParticles Particle GenParticle
    584606
    585607    add Branch TrackMerger/tracks Track Track
    586608    add Branch Calorimeter/towers Tower Tower
    587    
     609
    588610    add Branch Calorimeter/eflowTracks EFlowTrack Track
    589611    add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
    … …  
    593615    add Branch PhotonEfficiency/photons PhotonEff Photon
    594616    add Branch PhotonIsolation/photons PhotonIso Photon
    595    
     617
    596618    add Branch GenJetFinder/jets GenJet Jet
    597619    add Branch GenMissingET/momentum GenMissingET MissingET
    598    
     620
    599621    add Branch UniqueObjectFinder/jets Jet Jet
    600622    add Branch UniqueObjectFinder/electrons Electron Electron
    601623    add Branch UniqueObjectFinder/photons Photon Photon
    602624    add Branch UniqueObjectFinder/muons Muon Muon
    603    
    604     add Branch JetEnergyScale/jets AntiKtJet Jet 
    605    
     625
     626    add Branch JetEnergyScale/jets AntiKtJet Jet
     627
    606628    add Branch MissingET/momentum MissingET MissingET
    607629    add Branch ScalarHT/energy ScalarHT ScalarHT
    608630}
    609 
Note: See TracChangeset for help on using the changeset viewer.