Fork me on GitHub

Changes in / [fa7729a:1a8ef08] in git


Ignore:
Files:
20 edited

Legend:

Unmodified
Added
Removed

  • cards/ILCgen/ILCgen_BTagging_N.tcl

    rfa7729a r1a8ef08  
    104104    source ILCgen/ILCgen_BTagging_50.tcl
    105105}
     106

  • cards/ILCgen/ILCgen_BeamCalEfficiency.tcl

    rfa7729a r1a8ef08  
    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}
    46
    57set EfficiencyFormula {

  • cards/ILCgen/ILCgen_BeamCalF_Binning.tcl

    rfa7729a r1a8ef08  
    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

    rfa7729a r1a8ef08  
    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

    rfa7729a r1a8ef08  
    11#P.Sopicki
     2  # default energy fractions {abs(PDG code)} {fraction of energy deposited in BeamCal}
    23
    34  add EnergyFraction {0} {0.0}
    … …  
    1920  add EnergyFraction {310} {0.3}
    2021  add EnergyFraction {3122} {0.3}
     22

  • cards/ILCgen/ILCgen_BeamCal_Resolution.tcl

    rfa7729a r1a8ef08  
    11# P.Sopicki, corrected by A.F.Zarnecki
    22# BeamCal resolution
     3
     4# set ResolutionFormula {resolution formula as a function of eta and energy}
    35
    46  set ResolutionFormula {
    … …  
    68  (abs(eta) > 4.8 && abs(eta) <= 5.8) * sqrt(energy^2*0.03^2 + energy*0.45^2)
    79  }
     10

  • cards/ILCgen/ILCgen_CTagging_N.tcl

    rfa7729a r1a8ef08  
    104104    source ILCgen/ILCgen_CTagging_20.tcl
    105105}
     106

  • cards/ILCgen/ILCgen_ChrgHadMomentumSmearing.tcl

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

  • cards/ILCgen/ILCgen_ChrgHadTrackingEff.tcl

    rfa7729a r1a8ef08  
    44#  A.F.Zarnecki, June 12, 2020
    55
    6 set 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   }
     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)  }

  • cards/ILCgen/ILCgen_ECAL_EnergyFractions.tcl

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

  • cards/ILCgen/ILCgen_ECAL_Resolution.tcl

    rfa7729a r1a8ef08  
    22# ECAL and LumiCal resolution: same formula used at the moment
    33#
    4 
     4# set ResolutionFormula {resolution formula as a function of eta and energy}
    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

    rfa7729a r1a8ef08  
     1
     2  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
    13
    24# efficiency formula for electrons (ECAL and LumiCal)

  • cards/ILCgen/ILCgen_ElectronMomentumSmearing.tcl

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

  • cards/ILCgen/ILCgen_ElectronTrackingEff.tcl

    rfa7729a r1a8ef08  
    44#  A.F.Zarnecki, June 12, 2020
    55
    6 set 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   }
     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)  }

  • cards/ILCgen/ILCgen_HCAL_EnergyFractions.tcl

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

  • cards/ILCgen/ILCgen_HCAL_Resolution.tcl

    rfa7729a r1a8ef08  
    22# HCAL and LHCAL resolution: same formula used at the moment
    33#
    4 
     4# set ResolutionFormula {resolution formula as a function of eta and energy}
    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   }
     8}
     9

  • cards/ILCgen/ILCgen_MuonMomentumSmearing.tcl

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

  • cards/ILCgen/ILCgen_MuonTrackingEff.tcl

    rfa7729a r1a8ef08  
    44#  A.F.Zarnecki, June 12, 2020
    55
    6 set 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   }
     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)  }

  • cards/ILCgen/ILCgen_PhotonEfficiency.tcl

    rfa7729a r1a8ef08  
    1 # set EfficiencyFormula (efficiency formula as a function of eta and pt)
    2 # efficiency formula for photons (ECAL + LumiCal)
     1
     2  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     3
     4  # efficiency formula for photons (ECAL + LumiCal)
    35
    46  set EfficiencyFormula {

  • readers/DelphesPythia8.cpp

    rfa7729a r1a8ef08  
    296296
    297297    // jet matching
    298 #if PYTHIA_VERSION_INTEGER < 8300
    299298    matching = combined->getHook(*pythia);
    300299    if(!matching)
    … …  
    303302    }
    304303    pythia->setUserHooksPtr(matching);
    305 #endif
    306304
    307305    if(pythia == NULL)
Note: See TracChangeset for help on using the changeset viewer.