Fork me on GitHub

Changes in / [3cd4c6e:8dd735c] in git


Ignore:
Files:
2 deleted
4 edited

Legend:

Unmodified
Added
Removed
  • cards/FCC/FCChh.tcl

    r3cd4c6e r8dd735c  
    44#  Main authors:  Michele Selvaggi (CERN)
    55#
    6 #  Released on: October 14th, 2020
    7 #
    8 #  - fix muon resolution at high pT
    9 #  - updated btagging, tau tagging and photon ID
    10 #
     6#  Released on: Dec. 1st, 2017
    117#
    128#  Configuration: FCC-hh baseline detector
     
    10531049  (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0)  * (0.70) +
    10541050  (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
    1055   (abs(eta) <= 2.5) * (pt > 10.0)             * (0.90) +
     1051  (abs(eta) <= 2.5) * (pt > 10.0)             * (0.95) +
    10561052   
    10571053  (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0)  * (0.60) +
     
    11421138  (abs(eta) < 2.5)                   * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \
    11431139  (abs(eta) < 2.5)                   * (pt > 15000.0)               * (0.00) + \
    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.) + \
     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.) + \
    11461142  (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0)               * (0.000) + \
    11471143  (abs(eta) > 4.0) * (0.00)}
     
    11501146
    11511147  (pt <= 10.0)                       * (0.00) +
    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.) + \
     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.) + \
    11541150  (abs(eta) < 2.5)                   * (pt > 15000.0)               * (0.000) + \
    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.) + \
     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.) + \
    11571153  (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0)               * (0.000) + \
    11581154  (abs(eta) > 4.0) * (0.00)}
     
    11611157
    11621158  (pt <= 10.0)                                                       * (0.00) +
    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.) +
     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.) +
    11651161  (abs(eta) < 2.5)                    * (pt > 15000.0)               * (0.000) +
    11661162  (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500)      * (0.64) +
     
    12211217
    12221218  (pt <= 10.0)                                                       * (0.00) +
    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.) + \
     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.) + \
    12251221  (abs(eta) < 2.5)                   * (pt > 34000.0)                * (0.000) + \
    12261222  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0)    * (0.0075) + \
     
    12321228
    12331229  (pt <= 10.0)                                                       * (0.00) +
    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.) + \
     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.) + \
    12361232  (abs(eta) < 2.5)                   * (pt > 34000.0)                * (0.000) + \
    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.) + \
     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.) + \
    12391235  (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0)                * (0.00) + \
    12401236  (abs(eta) > 4.0)                   * (0.00)}
  • cards/FCC/muonMomentumResolutionVsP.tcl

    r3cd4c6e r8dd735c  
    1212set ResolutionFormula {
    1313 
    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     ) +
     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) +
    2928 
    3029( abs(eta) > 1.35 && abs(eta) < 2.00) *
  • cards/delphes_card_CircularEE.tcl

    r3cd4c6e r8dd735c  
    5252  BTagging
    5353  CTagging
    54   TauTaggingTight
    55   TauTaggingLoose
     54  TauTagging
    5655
    5756  BTaggingExclusive_N2
     
    6362  CTaggingExclusive_N6
    6463
    65   TauTaggingTightExclusive_N2
    66   TauTaggingTightExclusive_N4
    67   TauTaggingTightExclusive_N6
    68 
    69   TauTaggingLooseExclusive_N2
    70   TauTaggingLooseExclusive_N4
    71   TauTaggingLooseExclusive_N6
     64  TauTaggingExclusive_N2
     65  TauTaggingExclusive_N4
     66  TauTaggingExclusive_N6
    7267
    7368  ScalarHT
     
    724719
    725720
    726 module TauTagging TauTaggingTight {
     721module TauTagging TauTagging {
    727722  set ParticleInputArray Delphes/allParticles
    728723  set PartonInputArray Delphes/partons
    729724  set JetInputArray JetEnergyScale/jets
    730725
    731   set BitNumber 0
    732 
    733726  set DeltaR 0.5
    734727
     
    747740
    748741
    749 #############
    750 # tau-tagging
    751 #############
    752 
    753 
    754 module TauTagging TauTaggingLoose {
    755   set ParticleInputArray Delphes/allParticles
    756   set PartonInputArray Delphes/partons
    757   set JetInputArray JetEnergyScale/jets
    758 
    759   set BitNumber 1
    760 
    761   set DeltaR 0.5
    762 
    763   set TauPTMin 1.0
    764 
    765   set TauEtaMax 3.0
    766 
    767   # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    768 
    769   # default efficiency formula (misidentification rate)
    770   add EfficiencyFormula {0} {0.1}
    771   # efficiency formula for tau-jets
    772   add EfficiencyFormula {15} {0.85}
    773 }
    774 
    775 
    776 
    777742###########
    778743# b-tagging
     
    819784
    820785
    821 module TauTagging TauTaggingTightExclusive_N2 {
     786module TauTagging TauTaggingExclusive_N2 {
    822787  set ParticleInputArray Delphes/allParticles
    823788  set PartonInputArray Delphes/partons
    824789  set JetInputArray ExclusiveFastJetFinder_N2/jets
    825790
    826   set BitNumber 0
    827 
    828791  set DeltaR 0.5
    829792
     
    841804
    842805
    843 #############
    844 # tau-tagging
    845 #############
    846 
    847 
    848 module TauTagging TauTaggingLooseExclusive_N2 {
    849   set ParticleInputArray Delphes/allParticles
    850   set PartonInputArray Delphes/partons
    851   set JetInputArray ExclusiveFastJetFinder_N2/jets
    852 
    853   set BitNumber 1
    854 
    855   set DeltaR 0.5
    856 
    857   set TauPTMin 1.0
    858 
    859   set TauEtaMax 3.0
    860 
    861   # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    862 
    863   # default efficiency formula (misidentification rate)
    864   add EfficiencyFormula {0} {0.1}
    865   # efficiency formula for tau-jets
    866   add EfficiencyFormula {15} {0.85}
    867 }
    868 
    869 
    870 
    871 
    872 
    873806###########
    874807# b-tagging
     
    915848
    916849
    917 module TauTagging TauTaggingTightExclusive_N4 {
     850module TauTagging TauTaggingExclusive_N4 {
    918851  set ParticleInputArray Delphes/allParticles
    919852  set PartonInputArray Delphes/partons
    920853  set JetInputArray ExclusiveFastJetFinder_N4/jets
    921  
    922   set BitNumber 0
    923854
    924855  set DeltaR 0.5
     
    935866  add EfficiencyFormula {15} {0.4}
    936867}
    937 
    938 
    939 
    940 #############
    941 # tau-tagging
    942 #############
    943 
    944 
    945 module TauTagging TauTaggingLooseExclusive_N4 {
    946   set ParticleInputArray Delphes/allParticles
    947   set PartonInputArray Delphes/partons
    948   set JetInputArray ExclusiveFastJetFinder_N4/jets
    949  
    950   set BitNumber 1
    951 
    952   set DeltaR 0.5
    953 
    954   set TauPTMin 1.0
    955 
    956   set TauEtaMax 3.0
    957 
    958   # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    959 
    960   # default efficiency formula (misidentification rate)
    961   add EfficiencyFormula {0} {0.1}
    962   # efficiency formula for tau-jets
    963   add EfficiencyFormula {15} {0.85}
    964 }
    965 
    966 
    967 
    968868###########
    969869# b-tagging
     
    1010910
    1011911
    1012 module TauTagging TauTaggingTightExclusive_N6 {
     912module TauTagging TauTaggingExclusive_N6 {
    1013913  set ParticleInputArray Delphes/allParticles
    1014914  set PartonInputArray Delphes/partons
    1015915  set JetInputArray ExclusiveFastJetFinder_N6/jets
    1016916
    1017   set BitNumber 0
    1018 
    1019917  set DeltaR 0.5
    1020918
     
    1029927  # efficiency formula for tau-jets
    1030928  add EfficiencyFormula {15} {0.4}
    1031 }
    1032 
    1033 
    1034 #############
    1035 # tau-tagging
    1036 #############
    1037 
    1038 
    1039 module TauTagging TauTaggingLooseExclusive_N6 {
    1040   set ParticleInputArray Delphes/allParticles
    1041   set PartonInputArray Delphes/partons
    1042   set JetInputArray ExclusiveFastJetFinder_N6/jets
    1043 
    1044   set BitNumber 1
    1045 
    1046   set DeltaR 0.5
    1047 
    1048   set TauPTMin 1.0
    1049 
    1050   set TauEtaMax 3.0
    1051 
    1052   # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    1053 
    1054   # default efficiency formula (misidentification rate)
    1055   add EfficiencyFormula {0} {0.1}
    1056   # efficiency formula for tau-jets
    1057   add EfficiencyFormula {15} {0.85}
    1058929}
    1059930
  • modules/TrackCovariance.cc

    r3cd4c6e r8dd735c  
    133133    candidate->D0       = track.GetObsPar()[0]*1e03;
    134134    candidate->Phi      = track.GetObsPar()[1];
    135 
    136     // inverse of curvature
    137     candidate->C        = track.GetObsPar()[2]*1e-03;
     135    candidate->C        = track.GetObsPar()[2]*1e03;
    138136    candidate->DZ       = track.GetObsPar()[3]*1e03;
    139137    candidate->CtgTheta = track.GetObsPar()[4];
     
    148146    dpt       = 2 * TMath::Sqrt( track.GetCov()(2, 2))*pt*pt / (0.2998*fBz);
    149147    dp        = TMath::Sqrt((1.+ct*ct)*dpt*dpt + 4*pt*pt*ct*ct*dct*dct/(1.+ct*ct)/(1.+ct*ct));
    150     dC        = TMath::Sqrt(track.GetCov()(2, 2))*1e-03;
     148    dC        = TMath::Sqrt(track.GetCov()(2, 2))*1e03;
    151149
    152150    candidate->ErrorD0 = dd0;
Note: See TracChangeset for help on using the changeset viewer.