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_ATLAS_PileUp.tcl

    rd870fc5 rd77b51d  
    2222  EFlowMergerAllTracks
    2323  EFlowMerger
    24  
     24
    2525  NeutrinoFilter
    2626  GenJetFinder
    … …  
    4242
    4343  MissingET
     44 
     45  JetFlavorAssociation
    4446
    4547  BTagging
    … …  
    7981  set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
    8082
    81   #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
    82   #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
    83   #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) + \
     83  #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
     84  #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) +
     85  #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) +
    8486  #                               (abs(t) >  1.0e-09) * (abs(z) > 0.15)  * (0.00)}
    8587
    … …  
    119121
    120122  # tracking efficiency formula for charged hadrons
    121   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    122                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    123                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    124                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    125                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
     123  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     124                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     125                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     126                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     127                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    126128                         (abs(eta) > 2.5)                                                  * (0.00)}
    127129}
    … …  
    138140
    139141  # tracking efficiency formula for electrons
    140   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    141                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
    142                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
    143                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
    144                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
    145                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
    146                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
     142  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     143                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     144                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     145                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     146                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     147                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     148                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    147149                         (abs(eta) > 2.5)                                                  * (0.00)}
    148150}
    … …  
    159161
    160162  # tracking efficiency formula for muons
    161   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    162                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    163                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    164                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    165                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
     163  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     164                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     165                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     166                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     167                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    166168                         (abs(eta) > 2.5)                                                  * (0.00)}
    167169}
    … …  
    178180
    179181  # resolution formula for charged hadrons
    180   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    181                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    182                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    183                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    184                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    185                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    186                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
     182  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     183                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     184                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     185                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     186                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     187                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     188                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    187189                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    188190}
    … …  
    198200  # set ResolutionFormula {resolution formula as a function of eta and energy}
    199201
    200   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
    201                                            (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
    202                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
     202  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) +
     203                                           (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
     204                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) +
    203205                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    204206
    … …  
    216218
    217219  # resolution formula for muons
    218   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    219                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) + \
    220                                            (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
    221                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) + \
    222                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
    223                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) + \
    224                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
     220  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     221                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) +
     222                                           (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) +
     223                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) +
     224                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) +
     225                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) +
     226                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) +
    225227                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.10)}
    226228}
    … …  
    249251  set PhotonOutputArray photons
    250252
    251   set ECalEnergyMin 0.50
    252   set HCalEnergyMin 1.00
    253 
    254   set ECalEnergySignificanceMin 1.0
    255   set HCalEnergySignificanceMin 1.0
    256 
    257253  set EFlowTrackOutputArray eflowTracks
    258254  set EFlowPhotonOutputArray eflowPhotons
    259255  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     256
     257  set ECalEnergyMin 0.5
     258  set HCalEnergyMin 1.0
     259
     260  set ECalEnergySignificanceMin 1.0
     261  set HCalEnergySignificanceMin 1.0
     262
     263  set SmearTowerCenter true
    260264
    261265  set pi [expr {acos(-1)}]
    … …  
    307311  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    308312  # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
    309   set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
     313  set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
    310314                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
    311315
    … …  
    313317  # http://arxiv.org/pdf/hep-ex/0004009v1
    314318  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    315   set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
    316                              (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
     319  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
     320                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
    317321                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
    318322}
    … …  
    374378
    375379module FastJetGridMedianEstimator Rho {
    376  
     380
    377381  set InputArray Calorimeter/towers
    378382  set RhoOutputArray rho
    379383
    380   # etamin etamax gridsize_eta gridsize_phi
    381  
    382   add GridRange 0.0 2.5 0.5 0.5
    383   add GridRange 2.5 5.0 0.5 0.5
     384  # add GridRange rapmin rapmax drap dphi
     385  # rapmin - the minimum rapidity extent of the grid
     386  # rapmax - the maximum rapidity extent of the grid
     387  # drap - the grid spacing in rapidity
     388  # dphi - the grid spacing in azimuth
     389
     390  add GridRange -5.0 -2.5 1.0 1.0
     391  add GridRange -2.5 2.5 0.5 0.5
     392  add GridRange 2.5 5.0 1.0 1.0
    384393
    385394}
    … …  
    391400
    392401module PdgCodeFilter NeutrinoFilter {
    393  
     402
    394403  set InputArray Delphes/stableParticles
    395404  set OutputArray filteredParticles
    396405
    397406  set PTMin 0.0
    398  
     407
    399408  add PdgCode {12}
    400409  add PdgCode {14}
    … …  
    471480
    472481module Efficiency PhotonEfficiency {
    473   set InputArray Calorimeter/photons
     482  set InputArray Calorimeter/eflowPhotons
    474483  set OutputArray photons
    475484
    … …  
    477486
    478487  # efficiency formula for photons
    479   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    480                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    481                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     488  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     489                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     490                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    482491                         (abs(eta) > 2.5)                                   * (0.00)}
    483492}
    … …  
    512521
    513522  # efficiency formula for electrons
    514   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    515                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    516                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     523  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     524                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     525                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    517526                         (abs(eta) > 2.5)                                   * (0.00)}
    518527}
    … …  
    547556
    548557  # efficiency formula for muons
    549   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    550                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    551                          (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) + \
     558  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     559                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     560                         (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) +
    552561                         (abs(eta) > 2.7)                                   * (0.00)}
    553562}
    … …  
    595604}
    596605
     606########################
     607# Jet Flavor Association
     608########################
     609
     610module JetFlavorAssociation JetFlavorAssociation {
     611 
     612  set PartonInputArray Delphes/partons
     613  set ParticleInputArray Delphes/allParticles
     614  set ParticleLHEFInputArray Delphes/allParticlesLHEF
     615  set JetInputArray JetEnergyScale/jets
     616 
     617  set DeltaR 0.5
     618  set PartonPTMin 1.0
     619  set PartonEtaMax 2.5
     620
     621}
     622
    597623###########
    598624# b-tagging
    … …  
    600626
    601627module BTagging BTagging {
    602   set PartonInputArray Delphes/partons
    603628  set JetInputArray JetEnergyScale/jets
    604629
    605630  set BitNumber 0
    606 
    607   set DeltaR 0.5
    608 
    609   set PartonPTMin 1.0
    610 
    611   set PartonEtaMax 2.5
    612631
    613632  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    … …  
    619638
    620639  # efficiency formula for c-jets (misidentification rate)
    621   add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
    622                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
    623                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
     640  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) +
     641                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
     642                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
    624643                              (abs(eta) > 2.5)                                  * (0.000)}
    625644
    626645  # efficiency formula for b-jets
    627   add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
    628                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
    629                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
     646  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) +
     647                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
     648                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
    630649                              (abs(eta) > 2.5)                                  * (0.000)}
    631650}
     651
     652#############
     653# tau-tagging
     654#############
    632655
    633656module TauTagging TauTagging {
    … …  
    693716
    694717}
    695 
    696 
    697 
Note: See TracChangeset for help on using the changeset viewer.