Fork me on GitHub

Changeset 25291dd in git


Ignore:
Timestamp:
Oct 17, 2020, 4:05:03 AM (4 years ago)
Author:
Stephen Sekula <sekula@…>
Children:
44f4a4e
Parents:
4e5ed02 (diff), be1a3be (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'

Files:
23 added
8 deleted
28 edited

Legend:

Unmodified
Added
Removed
  • Makefile

    r4e5ed02 r25291dd  
    8989DISPLAYLIB = libDelphesDisplay.lib
    9090
    91 VERSION = $(shell cat VERSION)
     91VERSION = x.y.z
    9292DISTDIR = Delphes-$(VERSION)
    9393DISTTAR = $(DISTDIR).tar.gz
     
    174174tmp/examples/Example1.$(ObjSuf): \
    175175        examples/Example1.cpp \
    176         classes/DelphesClasses.h \
    177         external/ExRootAnalysis/ExRootResult.h \
    178         external/ExRootAnalysis/ExRootTreeBranch.h \
    179         external/ExRootAnalysis/ExRootTreeReader.h \
    180         external/ExRootAnalysis/ExRootTreeWriter.h \
    181         external/ExRootAnalysis/ExRootUtilities.h
    182 DelphesValidation$(ExeSuf): \
    183         tmp/validation/DelphesValidation.$(ObjSuf)
    184 
    185 tmp/validation/DelphesValidation.$(ObjSuf): \
    186         validation/DelphesValidation.cpp \
    187176        classes/DelphesClasses.h \
    188177        external/ExRootAnalysis/ExRootResult.h \
     
    199188        stdhep2pileup$(ExeSuf) \
    200189        CaloGrid$(ExeSuf) \
    201         Example1$(ExeSuf) \
    202         DelphesValidation$(ExeSuf)
     190        Example1$(ExeSuf)
    203191
    204192EXECUTABLE_OBJ +=  \
     
    210198        tmp/converters/stdhep2pileup.$(ObjSuf) \
    211199        tmp/examples/CaloGrid.$(ObjSuf) \
    212         tmp/examples/Example1.$(ObjSuf) \
    213         tmp/validation/DelphesValidation.$(ObjSuf)
     200        tmp/examples/Example1.$(ObjSuf)
    214201
    215202DelphesHepMC$(ExeSuf): \
     
    24342421        @echo ">> Building $(DISTTAR)"
    24352422        @mkdir -p $(DISTDIR)
    2436         @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb VERSION Makefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR)
     2423        @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb Makefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR)
    24372424        @find $(DISTDIR) -depth -name .\* -exec rm -rf {} \;
    24382425        @tar -czf $(DISTTAR) $(DISTDIR)
  • cards/FCC/FCChh.tcl

    r4e5ed02 r25291dd  
    44#  Main authors:  Michele Selvaggi (CERN)
    55#
    6 #  Released on: Dec. 1st, 2017
     6#  Released on: October 14th, 2020
     7#
     8#  - fix muon resolution at high pT
     9#  - updated btagging, tau tagging and photon ID
     10#
    711#
    812#  Configuration: FCC-hh baseline detector
     
    10491053  (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0)  * (0.70) +
    10501054  (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
    1051   (abs(eta) <= 2.5) * (pt > 10.0)             * (0.95) +
     1055  (abs(eta) <= 2.5) * (pt > 10.0)             * (0.90) +
    10521056   
    10531057  (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0)  * (0.60) +
     
    11381142  (abs(eta) < 2.5)                   * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \
    11391143  (abs(eta) < 2.5)                   * (pt > 15000.0)               * (0.00) + \
    1140   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.0075) + \
    1141   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.0075)*(1.0 - pt/15000.) + \
     1144  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.01) + \
     1145  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \
    11421146  (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0)               * (0.000) + \
    11431147  (abs(eta) > 4.0) * (0.00)}
     
    11461150
    11471151  (pt <= 10.0)                       * (0.00) +
    1148   (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 500)      * (0.05) + \
    1149   (abs(eta) < 2.5)                   * (pt > 500.0 && pt < 15000.0) * (0.05)*(1.0 - pt/15000.) + \
     1152  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 500)      * (0.15) + \
     1153  (abs(eta) < 2.5)                   * (pt > 500.0 && pt < 15000.0) * (0.15)*(1.0 - pt/15000.) + \
    11501154  (abs(eta) < 2.5)                   * (pt > 15000.0)               * (0.000) + \
    1151   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.03) + \
    1152   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.03)*(1.0 - pt/15000.) + \
     1155  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.10) + \
     1156  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.10)*(1.0 - pt/15000.) + \
    11531157  (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0)               * (0.000) + \
    11541158  (abs(eta) > 4.0) * (0.00)}
     
    11571161
    11581162  (pt <= 10.0)                                                       * (0.00) +
    1159   (abs(eta) < 2.5)                    * (pt > 10.0 && pt < 500)      * (0.85) +
    1160   (abs(eta) < 2.5)                    * (pt > 500.0 && pt < 15000.0) * (0.85)*(1.0 - pt/15000.) +
     1163  (abs(eta) < 2.5)                    * (pt > 10.0 && pt < 500)      * (0.82) +
     1164  (abs(eta) < 2.5)                    * (pt > 500.0 && pt < 15000.0) * (0.82)*(1.0 - pt/15000.) +
    11611165  (abs(eta) < 2.5)                    * (pt > 15000.0)               * (0.000) +
    11621166  (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.64) +
     
    12171221
    12181222  (pt <= 10.0)                                                       * (0.00) +
    1219   (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0)    * (0.01) + \
    1220   (abs(eta) < 2.5)                   * (pt > 5000.0 && pt < 34000.0) * (0.01)  *(8./9. - pt/30000.) + \
     1223  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0)    * (0.02) + \
     1224  (abs(eta) < 2.5)                   * (pt > 5000.0 && pt < 34000.0) * (0.02)  *(8./9. - pt/30000.) + \
    12211225  (abs(eta) < 2.5)                   * (pt > 34000.0)                * (0.000) + \
    12221226  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0)    * (0.0075) + \
     
    12281232
    12291233  (pt <= 10.0)                                                       * (0.00) +
    1230   (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0)    * (0.005) + \
    1231   (abs(eta) < 2.5)                   * (pt > 5000.0 && pt < 34000.0) * (0.005)  *(8./9. - pt/30000.) + \
     1234  (abs(eta) < 2.5)                   * (pt > 10.0 && pt < 5000.0)    * (0.001) + \
     1235  (abs(eta) < 2.5)                   * (pt > 5000.0 && pt < 34000.0) * (0.001)  *(8./9. - pt/30000.) + \
    12321236  (abs(eta) < 2.5)                   * (pt > 34000.0)                * (0.000) + \
    1233   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0)    * (0.00375) + \
    1234   (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
     1237  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0)    * (0.001) + \
     1238  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.001)*(8./9. - pt/30000.) + \
    12351239  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0)                * (0.00) + \
    12361240  (abs(eta) > 4.0)                   * (0.00)}
  • cards/FCC/muonMomentumResolutionVsP.tcl

    r4e5ed02 r25291dd  
    1212set ResolutionFormula {
    1313 
    14 ( abs(eta) < 1.35 ) *
    15 
    16 (sqrt(0.0000364164 + (
    17    9.06262e-8 *pt^2* cosh(
    18      eta)^2 *(2.82074e-7/sin(2*atan(exp(-abs(eta))))^2 + (
    19       504.525 *(1/400000000 + (0.117945* 1/cosh(eta)^2)/(
    20          pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
    21       sin(2*atan(exp(-abs(eta))))^2) *sin(2*atan(exp(-abs(eta))))^2)/(
    22    0.00516429/sin(2*atan(exp(-abs(eta))))^2 + (
    23     96868.8 *(1/400000000 + (0.117945 * 1/cosh(eta)^2)/(
    24        pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
    25     sin(2*atan(exp(-abs(eta))))^2))
    26 
    27 ) +
     14
     15( abs(eta) < 1.35 ) *
     16
     17    (sqrt(0.0000364164 + (
     18                          3*9.06262e-8 *pt^2* cosh(
     19         eta)^2 *(2.82074e-7/sin(2*atan(exp(-abs(eta))))^2 + (
     20          504.525 *(1/400000000 + (0.117945* 1/cosh(eta)^2)/(
     21             pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
     22          sin(2*atan(exp(-abs(eta))))^2) *sin(2*atan(exp(-abs(eta))))^2)/(
     23       0.00516429/sin(2*atan(exp(-abs(eta))))^2 + (
     24        96868.8 *(1/400000000 + 5*(0.117945 * 1/cosh(eta)^2)/(
     25           pt^2 *sin(2*atan(exp(-abs(eta))))^2)))/
     26        sin(2*atan(exp(-abs(eta))))^2))
     27
     28    ) +
    2829 
    2930( abs(eta) > 1.35 && abs(eta) < 2.00) *
  • cards/ILCgen/ILCgen_BTagging_N.tcl

    r4e5ed02 r25291dd  
    104104    source ILCgen/ILCgen_BTagging_50.tcl
    105105}
    106 
  • cards/ILCgen/ILCgen_BeamCalEfficiency.tcl

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    104104    source ILCgen/ILCgen_CTagging_20.tcl
    105105}
    106 
  • cards/ILCgen/ILCgen_ChrgHadMomentumSmearing.tcl

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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

    r4e5ed02 r25291dd  
    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 {
  • doc/genMakefile.tcl

    r4e5ed02 r25291dd  
    282282DISPLAYLIB = libDelphesDisplay.lib
    283283
    284 VERSION = $(shell cat VERSION)
     284VERSION = x.y.z
    285285DISTDIR = Delphes-$(VERSION)
    286286DISTTAR = $(DISTDIR).tar.gz
     
    290290}
    291291
    292 executableDeps {converters/*.cpp} {examples/*.cpp} {validation/*.cpp}
     292executableDeps {converters/*.cpp} {examples/*.cpp}
    293293
    294294executableDeps {readers/DelphesHepMC.cpp} {readers/DelphesLHEF.cpp} {readers/DelphesSTDHEP.cpp} {readers/DelphesROOT.cpp}
     
    434434        @echo ">> Building $(DISTTAR)"
    435435        @mkdir -p $(DISTDIR)
    436         @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb VERSION Makefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR)
     436        @cp -a AUTHORS CHANGELOG CMakeLists.txt COPYING DelphesEnv.sh LICENSE NOTICE README README_4LHCb Makefile MinBias.pileup configure cards classes converters display doc examples external modules python readers validation $(DISTDIR)
    437437        @find $(DISTDIR) -depth -name .\* -exec rm -rf {} \;
    438438        @tar -czf $(DISTTAR) $(DISTDIR)
  • modules/DualReadoutCalorimeter.cc

    r4e5ed02 r25291dd  
    645645    //cout<<"Creating tower with Pt, Eta, Phi, Energy: "<<pt<<","<<eta<<","<<phi<<","<<neutralEnergy<<endl;
    646646    tower->Momentum.SetPtEtaPhiE(pt, eta, phi, neutralEnergy);
    647     tower->Eem = neutralEnergy;
    648     tower->Ehad = 0.0;
    649     tower->PID = 22;
    650 
     647
     648    // if no hadronic energy, use ECAL resolution
     649    if (fHCalTowerEnergy <= fHCalEnergyMin)
     650    {
     651      tower->Eem = neutralEnergy;
     652      tower->Ehad = 0.0;
     653      tower->PID = 22;
     654    }
     655
     656    // if hadronic fraction > 0, use HCAL resolution
     657    else
     658    {
     659      tower->Eem = 0;
     660      tower->Ehad = neutralEnergy;
     661      tower->PID = 130;
     662    }
     663
     664    fEFlowPhotonOutputArray->Add(tower);
    651665   
    652 
    653     fEFlowPhotonOutputArray->Add(tower);
    654 
    655666
    656667    //clone tracks
  • modules/TrackCovariance.cc

    r4e5ed02 r25291dd  
    102102  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
    103103  {
    104     const TLorentzVector &candidatePosition = candidate->InitialPosition;
     104
     105    // converting to meters
     106    const TLorentzVector &candidatePosition = candidate->InitialPosition*1e-03;
    105107    const TLorentzVector &candidateMomentum = candidate->Momentum;
    106108
     
    113115
    114116    candidate->Momentum.SetVectM(track.GetObsP(), mass);
    115     candidate->InitialPosition.SetXYZT(track.GetObsX().X(),track.GetObsX().Y(),track.GetObsX().Z(),candidatePosition.T());
    116117   
     118    // converting back to mm
     119    candidate->InitialPosition.SetXYZT(track.GetObsX().X()*1e03,track.GetObsX().Y()*1e03,track.GetObsX().Z()*1e03,candidatePosition.T()*1e03);
     120
    117121    // save full covariance 5x5 matrix internally (D0, phi, Curvature, dz, ctg(theta))
    118122    candidate->TrackCovariance = track.GetCov();
     
    123127    ct = track.GetObsPar()[4];
    124128
    125     candidate->Xd = track.GetObsX().X();
    126     candidate->Yd = track.GetObsX().Y();
    127     candidate->Zd = track.GetObsX().Z();
     129    candidate->Xd = track.GetObsX().X()*1e03;
     130    candidate->Yd = track.GetObsX().Y()*1e03;
     131    candidate->Zd = track.GetObsX().Z()*1e03;
    128132   
    129     candidate->D0 = track.GetObsPar()[0];
    130     candidate->Phi = track.GetObsPar()[1];
    131     candidate->C = track.GetObsPar()[2];
    132     candidate->DZ = track.GetObsPar()[3];
     133    candidate->D0       = track.GetObsPar()[0]*1e03;
     134    candidate->Phi      = track.GetObsPar()[1];
     135    candidate->C        = track.GetObsPar()[2]*1e03;
     136    candidate->DZ       = track.GetObsPar()[3]*1e03;
    133137    candidate->CtgTheta = track.GetObsPar()[4];
    134     candidate->P  = track.GetObsP().Mag();
    135     candidate->PT = pt;
    136     candidate->Charge = q;
     138    candidate->P        = track.GetObsP().Mag();
     139    candidate->PT       = pt;
     140    candidate->Charge   = q;
    137141
    138     dd0       = TMath::Sqrt(track.GetCov()(0, 0));
    139     ddz       = TMath::Sqrt(track.GetCov()(3, 3));
     142    dd0       = TMath::Sqrt(track.GetCov()(0, 0))*1e03;
     143    ddz       = TMath::Sqrt(track.GetCov()(3, 3))*1e03;
    140144    dphi      = TMath::Sqrt(track.GetCov()(1, 1));
    141145    dct       = TMath::Sqrt(track.GetCov()(4, 4));
    142146    dpt       = 2 * TMath::Sqrt( track.GetCov()(2, 2))*pt*pt / (0.2998*fBz);
    143147    dp        = TMath::Sqrt((1.+ct*ct)*dpt*dpt + 4*pt*pt*ct*ct*dct*dct/(1.+ct*ct)/(1.+ct*ct));
    144     dC        = TMath::Sqrt(track.GetCov()(2, 2));
     148    dC        = TMath::Sqrt(track.GetCov()(2, 2))*1e03;
    145149
    146150    candidate->ErrorD0 = dd0;
  • readers/DelphesPythia8.cpp

    r4e5ed02 r25291dd  
    296296
    297297    // jet matching
     298#if PYTHIA_VERSION_INTEGER < 8300
    298299    matching = combined->getHook(*pythia);
    299300    if(!matching)
     
    302303    }
    303304    pythia->setUserHooksPtr(matching);
     305#endif
    304306
    305307    if(pythia == NULL)
  • validation/DelphesValidation.cpp

    r4e5ed02 r25291dd  
    3232#include "TGraph.h"
    3333#include "TGraphErrors.h"
     34#include "TF1.h"
    3435#include "TH1.h"
    3536#include "TH2.h"
  • validation/validation.sh

    r4e5ed02 r25291dd  
    3838outputrootdir=report_${cardbase%.*}/root
    3939cardlabel=${cardbase%.*}
    40 version=$(cat VERSION)
     40version=x.y.z
    4141outpdf=$mainoutputdir/${output%.*}.pdf
    4242cardsdir=validation/cards
Note: See TracChangeset for help on using the changeset viewer.