Fork me on GitHub

Changeset 3e5f05c in git


Ignore:
Timestamp:
02/26/20 10:57:43 (3 months 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.