Fork me on GitHub

Changeset 4f52352 in git


Ignore:
Timestamp:
Aug 31, 2020, 8:31:15 AM (4 years ago)
Author:
GitHub <noreply@…>
Branches:
master
Children:
fa7729a
Parents:
a643a7f (diff), f0f1443 (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.
git-author:
Michele Selvaggi <michele.selvaggi@…> (08/31/20 08:31:15)
git-committer:
GitHub <noreply@…> (08/31/20 08:31:15)
Message:

Merge pull request #77 from zarnecki/master

Code cleaning

Location:
cards/ILCgen
Files:
19 edited

Legend:

Unmodified
Added
Removed

  • cards/ILCgen/ILCgen_BTagging_N.tcl

    ra643a7f r4f52352  
    104104    source ILCgen/ILCgen_BTagging_50.tcl
    105105}
    106 

  • cards/ILCgen/ILCgen_BeamCalEfficiency.tcl

    ra643a7f r4f52352  
    22# Based on IDR Figure 8.8a and Moritz Hebermehl PhD Thesis Figure 4.10
    33# A.F.Zarnecki,  June 2020
    4 
    5 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    64
    75set EfficiencyFormula {

  • cards/ILCgen/ILCgen_BeamCalF_Binning.tcl

    ra643a7f r4f52352  
    22# corrected by A.F.Zarnecki
    33
    4 set pi [expr {acos(-1)}]
    5    
     4set pi [expr {acos(-1)} ]
     5
    66# BeamCal eta range 4.0 - 5.8 -> 2.099-0.347 =1.752deg
    77#

  • cards/ILCgen/ILCgen_BeamCalR_Binning.tcl

    ra643a7f r4f52352  
    22# corrected by A.F.Zarnecki
    33
    4 set pi [expr {acos(-1)}]
    5    
     4set pi [expr {acos(-1)} ]
     5
    66# BeamCal eta range 4.0 - 5.8 -> 2.099-0.347 =1.752deg
    77#

  • cards/ILCgen/ILCgen_BeamCal_EnergyFractions.tcl

    ra643a7f r4f52352  
    11#P.Sopicki
    2   # default energy fractions {abs(PDG code)} {fraction of energy deposited in BeamCal}
    32
    43  add EnergyFraction {0} {0.0}
    … …  
    2019  add EnergyFraction {310} {0.3}
    2120  add EnergyFraction {3122} {0.3}
    22 

  • cards/ILCgen/ILCgen_BeamCal_Resolution.tcl

    ra643a7f r4f52352  
    11# P.Sopicki, corrected by A.F.Zarnecki
    22# BeamCal resolution
    3 
    4 # set ResolutionFormula {resolution formula as a function of eta and energy}
    53
    64  set ResolutionFormula {
    … …  
    86  (abs(eta) > 4.8 && abs(eta) <= 5.8) * sqrt(energy^2*0.03^2 + energy*0.45^2)
    97  }
    10 

  • cards/ILCgen/ILCgen_CTagging_N.tcl

    ra643a7f r4f52352  
    104104    source ILCgen/ILCgen_CTagging_20.tcl
    105105}
    106 

  • cards/ILCgen/ILCgen_ChrgHadMomentumSmearing.tcl

    ra643a7f r4f52352  
    66
    77set ResolutionFormula {
    8 
    9 (abs(eta) <= 2.0) *
    10   exp( -6.82012  - 0.284133*abs(eta) + 0.64585*eta*eta
     8  (abs(eta) <= 2.0) *
     9    exp( -6.82012  - 0.284133*abs(eta) + 0.64585*eta*eta
    1110     + (0.829675 -0.188197*abs(eta) +  0.106226 * eta*eta  - 1.0) * log(pt)
    1211     + (0.0766417-0.00435667*abs(eta) + 0.00174804*eta*eta) * log(pt)*log(pt))
    13 
    14 +
    15 
    16 (abs(eta) > 2.0 && abs(eta) <= 3.0) *
    17   exp(  -4.54609
     12  +
     13  (abs(eta) > 2.0 && abs(eta) <= 3.0) *
     14    exp(  -4.54609
    1815     + ( 0.87526 - 1.0) * log(pt)
    1916     + (-0.0800672 + 0.0891669 * abs(eta)) * log(pt)*log(pt))
    20 
    21 }
    22 
     17  }

  • cards/ILCgen/ILCgen_ChrgHadTrackingEff.tcl

    ra643a7f r4f52352  
    44#  A.F.Zarnecki, June 12, 2020
    55
    6   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00)  +
    7                                    (abs(eta) <= 1.83)     *      (pt > 0.1 && pt <= 0.2)   * (0.70)  +
    8                                    (abs(eta) <= 1.83)     *      (pt > 0.2 && pt <= 0.3)   * (0.93)  +
    9                                    (abs(eta) <= 1.83)     *      (pt > 0.3 && pt <= 0.7)   * (0.995) +
    10                                    (abs(eta) <= 1.83)     *      (pt > 0.7)                * (1.00)  +
    11                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.1 && pt <= 0.2)   * (0.697) +
    12                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.2 && pt <= 0.3)   * (0.925) +
    13                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.3 && pt <= 0.7)   * (0.99)  +
    14                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.7)                * (0.995) +
    15                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.1 && pt <= 0.2)   * (0.665) +
    16                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.2 && pt <= 0.3)   * (0.884) +
    17                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.3 && pt <= 0.7)   * (0.945) +
    18                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.7)                * (0.95)  +
    19                 (abs(eta) > 3.00)                         *      (pt > 0.1)                * (0.00)  }
     6set EfficiencyFormula {
     7                                                                 (pt <= 0.1)   * (0.00)  +
     8                       (abs(eta) <= 1.83)     *      (pt > 0.1 && pt <= 0.2)   * (0.70)  +
     9                       (abs(eta) <= 1.83)     *      (pt > 0.2 && pt <= 0.3)   * (0.93)  +
     10                       (abs(eta) <= 1.83)     *      (pt > 0.3 && pt <= 0.7)   * (0.995) +
     11                       (abs(eta) <= 1.83)     *      (pt > 0.7)                * (1.00)  +
     12    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.1 && pt <= 0.2)   * (0.697) +
     13    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.2 && pt <= 0.3)   * (0.925) +
     14    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.3 && pt <= 0.7)   * (0.99)  +
     15    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.7)                * (0.995) +
     16    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.1 && pt <= 0.2)   * (0.665) +
     17    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.2 && pt <= 0.3)   * (0.884) +
     18    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.3 && pt <= 0.7)   * (0.945) +
     19    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.7)                * (0.95)  +
     20    (abs(eta) > 3.00)                         *      (pt > 0.1)                * (0.00)
     21  }

  • cards/ILCgen/ILCgen_ECAL_EnergyFractions.tcl

    ra643a7f r4f52352  
    11#P.Sopicki
    2 # default energy fractions {abs(PDG code)} {fraction in ECAL or LumiCal}
     2
    33  add EnergyFraction {0} {0.0}
    44  # energy fractions for e, gamma and pi0

  • cards/ILCgen/ILCgen_ECAL_Resolution.tcl

    ra643a7f r4f52352  
    22# ECAL and LumiCal resolution: same formula used at the moment
    33#
    4 # set ResolutionFormula {resolution formula as a function of eta and energy}
     4
    55set ResolutionFormula {
    66    (abs(eta) <= 3 )                 * sqrt(energy^2*0.01^2 + energy*0.17^2)+
    77    (abs(eta) > 3 && abs(eta) <= 4 ) * sqrt(energy^2*0.01^2 + energy*0.17^2)
    8  }
     8  }

  • cards/ILCgen/ILCgen_ElectronEfficiency.tcl

    ra643a7f r4f52352  
    1 
    2   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    31
    42# efficiency formula for electrons (ECAL and LumiCal)

  • cards/ILCgen/ILCgen_ElectronMomentumSmearing.tcl

    ra643a7f r4f52352  
    66
    77set ResolutionFormula {
    8 
    9 (abs(eta) <= 2.0) *
    10   exp( -6.82012  - 0.284133*abs(eta) + 0.64585*eta*eta
     8  (abs(eta) <= 2.0) *
     9    exp( -6.82012  - 0.284133*abs(eta) + 0.64585*eta*eta
    1110     + (0.829675 -0.188197*abs(eta) +  0.106226 * eta*eta  - 1.0) * log(pt)
    1211     + (0.0766417-0.00435667*abs(eta) + 0.00174804*eta*eta) * log(pt)*log(pt))
    13 
    14 +
    15 
    16 (abs(eta) > 2.0 && abs(eta) <= 3.0) *
    17   exp(  -4.54609
     12  +
     13  (abs(eta) > 2.0 && abs(eta) <= 3.0) *
     14    exp(  -4.54609
    1815     + ( 0.87526 - 1.0) * log(pt)
    1916     + (-0.0800672 + 0.0891669 * abs(eta)) * log(pt)*log(pt))
    20 
    21 }
    22 
     17  }

  • cards/ILCgen/ILCgen_ElectronTrackingEff.tcl

    ra643a7f r4f52352  
    44#  A.F.Zarnecki, June 12, 2020
    55
    6   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00)  +
    7                                    (abs(eta) <= 1.83)     *      (pt > 0.1 && pt <= 0.2)   * (0.70)  +
    8                                    (abs(eta) <= 1.83)     *      (pt > 0.2 && pt <= 0.3)   * (0.93)  +
    9                                    (abs(eta) <= 1.83)     *      (pt > 0.3 && pt <= 0.7)   * (0.995) +
    10                                    (abs(eta) <= 1.83)     *      (pt > 0.7)                * (1.00)  +
    11                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.1 && pt <= 0.2)   * (0.697) +
    12                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.2 && pt <= 0.3)   * (0.925) +
    13                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.3 && pt <= 0.7)   * (0.99)  +
    14                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.7)                * (0.995) +
    15                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.1 && pt <= 0.2)   * (0.665) +
    16                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.2 && pt <= 0.3)   * (0.884) +
    17                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.3 && pt <= 0.7)   * (0.945) +
    18                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.7)                * (0.95)  +
    19                 (abs(eta) > 3.00)                         *      (pt > 0.1)                * (0.00)  }
     6set EfficiencyFormula {
     7                                                                 (pt <= 0.1)   * (0.00)  +
     8                       (abs(eta) <= 1.83)     *      (pt > 0.1 && pt <= 0.2)   * (0.70)  +
     9                       (abs(eta) <= 1.83)     *      (pt > 0.2 && pt <= 0.3)   * (0.93)  +
     10                       (abs(eta) <= 1.83)     *      (pt > 0.3 && pt <= 0.7)   * (0.995) +
     11                       (abs(eta) <= 1.83)     *      (pt > 0.7)                * (1.00)  +
     12    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.1 && pt <= 0.2)   * (0.697) +
     13    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.2 && pt <= 0.3)   * (0.925) +
     14    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.3 && pt <= 0.7)   * (0.99)  +
     15    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.7)                * (0.995) +
     16    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.1 && pt <= 0.2)   * (0.665) +
     17    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.2 && pt <= 0.3)   * (0.884) +
     18    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.3 && pt <= 0.7)   * (0.945) +
     19    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.7)                * (0.95)  +
     20    (abs(eta) > 3.00)                         *      (pt > 0.1)                * (0.00)
     21  }

  • cards/ILCgen/ILCgen_HCAL_EnergyFractions.tcl

    ra643a7f r4f52352  
    11#P.Sopicki
    2 # default energy fractions {abs(PDG code)} {fraction in HCAL or LHCAL}
     2
    33  add EnergyFraction {0} {1.0}
    44  # energy fractions for e, gamma and pi0

  • cards/ILCgen/ILCgen_HCAL_Resolution.tcl

    ra643a7f r4f52352  
    22# HCAL and LHCAL resolution: same formula used at the moment
    33#
    4 # set ResolutionFormula {resolution formula as a function of eta and energy}
     4
    55set ResolutionFormula {
    66   (abs(eta) <= 2.8 )                 * sqrt(energy^2*0.017^2 + energy*0.45^2)+
    77   (abs(eta) > 2.8 && abs(eta)<=3.8 ) * sqrt(energy^2*0.017^2 + energy*0.45^2)
    8 }
    9 
     8  }

  • cards/ILCgen/ILCgen_MuonMomentumSmearing.tcl

    ra643a7f r4f52352  
    66
    77set ResolutionFormula {
    8 
    9 (abs(eta) <= 2.0) *
    10   exp( -6.82012  - 0.284133*abs(eta) + 0.64585*eta*eta
     8  (abs(eta) <= 2.0) *
     9    exp( -6.82012  - 0.284133*abs(eta) + 0.64585*eta*eta
    1110     + (0.829675 -0.188197*abs(eta) +  0.106226 * eta*eta  - 1.0) * log(pt)
    1211     + (0.0766417-0.00435667*abs(eta) + 0.00174804*eta*eta) * log(pt)*log(pt))
    13 
    14 +
    15 
    16 (abs(eta) > 2.0 && abs(eta) <= 3.0) *
    17   exp(  -4.54609
     12  +
     13  (abs(eta) > 2.0 && abs(eta) <= 3.0) *
     14    exp(  -4.54609
    1815     + ( 0.87526 - 1.0) * log(pt)
    1916     + (-0.0800672 + 0.0891669 * abs(eta)) * log(pt)*log(pt))
    20 
    21 }
    22 
     17  }

  • cards/ILCgen/ILCgen_MuonTrackingEff.tcl

    ra643a7f r4f52352  
    44#  A.F.Zarnecki, June 12, 2020
    55
    6   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00)  +
    7                                    (abs(eta) <= 1.83)     *      (pt > 0.1 && pt <= 0.2)   * (0.70)  +
    8                                    (abs(eta) <= 1.83)     *      (pt > 0.2 && pt <= 0.3)   * (0.93)  +
    9                                    (abs(eta) <= 1.83)     *      (pt > 0.3 && pt <= 0.7)   * (0.995) +
    10                                    (abs(eta) <= 1.83)     *      (pt > 0.7)                * (1.00)  +
    11                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.1 && pt <= 0.2)   * (0.697) +
    12                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.2 && pt <= 0.3)   * (0.925) +
    13                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.3 && pt <= 0.7)   * (0.99)  +
    14                 (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.7)                * (0.995) +
    15                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.1 && pt <= 0.2)   * (0.665) +
    16                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.2 && pt <= 0.3)   * (0.884) +
    17                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.3 && pt <= 0.7)   * (0.945) +
    18                 (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.7)                * (0.95)  +
    19                 (abs(eta) > 3.00)                         *      (pt > 0.1)                * (0.00)  }
     6set EfficiencyFormula {
     7                                                                 (pt <= 0.1)   * (0.00)  +
     8                       (abs(eta) <= 1.83)     *      (pt > 0.1 && pt <= 0.2)   * (0.70)  +
     9                       (abs(eta) <= 1.83)     *      (pt > 0.2 && pt <= 0.3)   * (0.93)  +
     10                       (abs(eta) <= 1.83)     *      (pt > 0.3 && pt <= 0.7)   * (0.995) +
     11                       (abs(eta) <= 1.83)     *      (pt > 0.7)                * (1.00)  +
     12    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.1 && pt <= 0.2)   * (0.697) +
     13    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.2 && pt <= 0.3)   * (0.925) +
     14    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.3 && pt <= 0.7)   * (0.99)  +
     15    (abs(eta) > 1.83 && abs(eta) <= 2.65)     *      (pt > 0.7)                * (0.995) +
     16    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.1 && pt <= 0.2)   * (0.665) +
     17    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.2 && pt <= 0.3)   * (0.884) +
     18    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.3 && pt <= 0.7)   * (0.945) +
     19    (abs(eta) > 2.65 && abs(eta) <= 3.00)     *      (pt > 0.7)                * (0.95)  +
     20    (abs(eta) > 3.00)                         *      (pt > 0.1)                * (0.00)
     21  }

  • cards/ILCgen/ILCgen_PhotonEfficiency.tcl

    ra643a7f r4f52352  
    1 
    2   # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    3 
    4   # efficiency formula for photons (ECAL + LumiCal)
     1# set EfficiencyFormula (efficiency formula as a function of eta and pt)
     2# efficiency formula for photons (ECAL + LumiCal)
    53
    64  set EfficiencyFormula {
Note: See TracChangeset for help on using the changeset viewer.