Fork me on GitHub

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


Ignore:
Files:
2 added
4 edited

Legend:

Unmodified
Added
Removed

  • cards/FCC/FCChh.tcl

    r8dd735c r3cd4c6e  
    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

    r8dd735c r3cd4c6e  
    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/delphes_card_CircularEE.tcl

    r8dd735c r3cd4c6e  
    5252  BTagging
    5353  CTagging
    54   TauTagging
     54  TauTaggingTight
     55  TauTaggingLoose
    5556
    5657  BTaggingExclusive_N2
    … …  
    6263  CTaggingExclusive_N6
    6364
    64   TauTaggingExclusive_N2
    65   TauTaggingExclusive_N4
    66   TauTaggingExclusive_N6
     65  TauTaggingTightExclusive_N2
     66  TauTaggingTightExclusive_N4
     67  TauTaggingTightExclusive_N6
     68
     69  TauTaggingLooseExclusive_N2
     70  TauTaggingLooseExclusive_N4
     71  TauTaggingLooseExclusive_N6
    6772
    6873  ScalarHT
    … …  
    719724
    720725
    721 module TauTagging TauTagging {
     726module TauTagging TauTaggingTight {
    722727  set ParticleInputArray Delphes/allParticles
    723728  set PartonInputArray Delphes/partons
    724729  set JetInputArray JetEnergyScale/jets
    725730
     731  set BitNumber 0
     732
    726733  set DeltaR 0.5
    727734
    … …  
    740747
    741748
     749#############
     750# tau-tagging
     751#############
     752
     753
     754module 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
    742777###########
    743778# b-tagging
    … …  
    784819
    785820
    786 module TauTagging TauTaggingExclusive_N2 {
     821module TauTagging TauTaggingTightExclusive_N2 {
    787822  set ParticleInputArray Delphes/allParticles
    788823  set PartonInputArray Delphes/partons
    789824  set JetInputArray ExclusiveFastJetFinder_N2/jets
    790825
     826  set BitNumber 0
     827
    791828  set DeltaR 0.5
    792829
    … …  
    804841
    805842
     843#############
     844# tau-tagging
     845#############
     846
     847
     848module 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
    806873###########
    807874# b-tagging
    … …  
    848915
    849916
    850 module TauTagging TauTaggingExclusive_N4 {
     917module TauTagging TauTaggingTightExclusive_N4 {
    851918  set ParticleInputArray Delphes/allParticles
    852919  set PartonInputArray Delphes/partons
    853920  set JetInputArray ExclusiveFastJetFinder_N4/jets
     921 
     922  set BitNumber 0
    854923
    855924  set DeltaR 0.5
    … …  
    866935  add EfficiencyFormula {15} {0.4}
    867936}
     937
     938
     939
     940#############
     941# tau-tagging
     942#############
     943
     944
     945module 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
    868968###########
    869969# b-tagging
    … …  
    9101010
    9111011
    912 module TauTagging TauTaggingExclusive_N6 {
     1012module TauTagging TauTaggingTightExclusive_N6 {
    9131013  set ParticleInputArray Delphes/allParticles
    9141014  set PartonInputArray Delphes/partons
    9151015  set JetInputArray ExclusiveFastJetFinder_N6/jets
    9161016
     1017  set BitNumber 0
     1018
    9171019  set DeltaR 0.5
    9181020
    … …  
    9271029  # efficiency formula for tau-jets
    9281030  add EfficiencyFormula {15} {0.4}
     1031}
     1032
     1033
     1034#############
     1035# tau-tagging
     1036#############
     1037
     1038
     1039module 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}
    9291058}
    9301059

  • modules/TrackCovariance.cc

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