Fork me on GitHub

Ignore:
Timestamp:
Sep 29, 2015, 2:08:10 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
a98c7ef
Parents:
d870fc5 (diff), 06ec139 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.
Message:

Merge remote-tracking branch 'upstream/master'

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_FCC_basic.tcl

    rd870fc5 rd77b51d  
    1717  AngularSmearing
    1818  ImpactParameterSmearing
    19  
     19
    2020  ECal
    2121  HCal
    22  
     22
    2323  TowerMerger
    2424  EFlowMerger
    2525
    2626  MissingET
    27  
     27
    2828  NeutrinoFilter
    2929  GenJetFinder
    … …  
    3131
    3232  JetEnergyScale
    33  
     33
    3434  TrackCountingBTagging
    3535  TauTagging
    … …  
    7272
    7373  # tracking efficiency formula for charged hadrons
    74   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    75                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    76                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    77                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    78                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.90) + \
     74  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     75                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     76                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     77                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     78                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.90) +
    7979                         (abs(eta) > 4.0)                                                  * (0.00)}
    80                          
     80
    8181}
    8282
    … …  
    9292
    9393  # tracking efficiency formula for electrons
    94   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    95                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    96                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    97                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    98                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) + \
     94  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     95                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     96                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     97                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     98                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) +
    9999                         (abs(eta) > 4.0)                                                  * (0.00)}
    100100}
    … …  
    111111
    112112  # tracking efficiency formula for muons
    113   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    114                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    115                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    116                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    117                          (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) + \
     113  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     114                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     115                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     116                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     117                         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0)                  * (0.98) +
    118118                         (abs(eta) > 4.0)                                                  * (0.00)}
    119119}
    … …  
    130130
    131131  # resolution formula for charged hadrons
    132   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) + \
     132  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*2.e-5) +
    133133                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*3.e-5)}
    134134
    … …  
    147147
    148148  # resolution formula for electrons
    149   set ResolutionFormula {                  (abs(eta) <= 4.0) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.007) + \
    150                                            (abs(eta) <= 4.0) * (energy > 2.0e1)                    * sqrt(energy^2*0.005^2 + energy*0.02^2) + \
     149  set ResolutionFormula {                  (abs(eta) <= 4.0) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.007) +
     150                                           (abs(eta) <= 4.0) * (energy > 2.0e1)                    * sqrt(energy^2*0.005^2 + energy*0.02^2) +
    151151                                           (abs(eta) > 4.0 && abs(eta) <= 6.0)                     * sqrt(energy^2*0.05^2 + energy*1.00^2)}
    152152
    … …  
    164164
    165165  # resolution formula for muons
    166   set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*5.e-6) + \
     166  set ResolutionFormula {    (abs(eta) <= 1.5)                   * (pt > 0.1) * (0.01 + pt*5.e-6) +
    167167                             (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1) * (0.02 + pt*1.e-5)}
    168168
    … …  
    210210
    211211  # absolute impact parameter smearing formula (in mm) as a function of pt and eta
    212   set ResolutionFormula {(pt > 0.1  && pt <= 5.0)   * (0.010) + \
     212  set ResolutionFormula {(pt > 0.1  && pt <= 5.0)   * (0.010) +
    213213                         (pt > 5.0)                 * (0.005)}
    214214
    … …  
    225225  set TowerOutputArray ecalTowers
    226226  set EFlowTowerOutputArray eflowPhotons
    227  
     227
     228  set IsEcal true
     229 
    228230  set EnergyMin 0.5
    229231  set EnergySignificanceMin 1.0
    230  
     232
     233  set SmearTowerCenter true
     234
    231235  set pi [expr {acos(-1)}]
    232236
    … …  
    234238  # each list starts with the lower edge of the first tower
    235239  # the list ends with the higher edged of the last tower
    236    
     240
    237241  # 0.5 degree towers
    238242  set PhiBins {}
    … …  
    246250    add EtaPhiBins $eta $PhiBins
    247251  }
    248  
     252
    249253  # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
    250  
     254
    251255  add EnergyFraction {0} {0.0}
    252256  # energy fractions for e, gamma and pi0
    … …  
    269273
    270274  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    271  
     275
    272276  # This is the CMS ECAL resolution, extended up eta = 6.0
    273   set ResolutionFormula { (abs(eta) <= 3.0)                   * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2)  + \
     277  set ResolutionFormula { (abs(eta) <= 3.0)                   * sqrt(energy^2*0.003^2 + energy*0.029^2 + 0.125^2)  +
    274278                          (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    275279
    276  
     280
    277281}
    278282
    … …  
    287291  set TowerOutputArray hcalTowers
    288292  set EFlowTowerOutputArray eflowNeutralHadrons
    289  
     293
     294  set IsEcal false
     295 
    290296  set EnergyMin 1.0
    291297  set EnergySignificanceMin 1.0
    292  
     298
     299  set SmearTowerCenter true
     300
    293301  set pi [expr {acos(-1)}]
    294302
    … …  
    296304  # each list starts with the lower edge of the first tower
    297305  # the list ends with the higher edged of the last tower
    298  
    299  
     306
     307
    300308  # 5 degree towers
    301309  set PhiBins {}
    … …  
    309317    add EtaPhiBins $eta $PhiBins
    310318  }
    311  
    312    
     319
     320
    313321  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    314322  add EnergyFraction {0} {1.0}
    … …  
    332340
    333341  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    334  
     342
    335343  # This is the ATLAS HCAL resolution, extended up eta = 6.0
    336   set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
    337                              (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
     344  set ResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
     345                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
    338346                             (abs(eta) > 3.2 && abs(eta) <= 6.0) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
    339347
    … …  
    373381  set MomentumOutputArray momentum
    374382}
    375                    
     383
    376384
    377385##################
    … …  
    390398
    391399module PdgCodeFilter NeutrinoFilter {
    392  
     400
    393401  set InputArray Delphes/stableParticles
    394402  set OutputArray filteredParticles
    395403
    396404  set PTMin 0.0
    397  
     405
    398406  add PdgCode {12}
    399407  add PdgCode {14}
    … …  
    427435
    428436module FastJetFinder FastJetFinder {
    429 #  set InputArray Calorimeter/towers
     437#  set InputArray TowerMerger/towers
    430438  set InputArray EFlowMerger/eflow
    431439
    … …  
    463471  # maximum distance between jet and track
    464472  set DeltaR 0.3
    465  
    466   # minimum pt of tracks 
     473
     474  # minimum pt of tracks
    467475  set TrackPTMin 1.0
    468  
     476
    469477  # minimum transverse impact parameter (in mm)
    470478  set TrackIPMax 2.0
    471  
     479
    472480  # minimum ip significance for the track to be counted
    473481  set SigMin 6.5
    474  
     482
    475483  # minimum number of tracks (high efficiency n=2, high purity n=3)
    476484  set Ntracks 3
    … …  
    510518  add Branch Delphes/allParticles Particle GenParticle
    511519  add Branch GenJetFinder/jets GenJet Jet
    512  
     520
    513521  add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
    514522  add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
    … …  
    521529  add Branch ScalarHT/energy ScalarHT ScalarHT
    522530}
    523  
Note: See TracChangeset for help on using the changeset viewer.