Fork me on GitHub

Changeset 05a7228 in git for cards


Ignore:
Timestamp:
Apr 23, 2019, 6:29:32 PM (5 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, master
Children:
0453894
Parents:
190cfa0
Message:

added tau tagging parametrisation based on DNN

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/CMS_PhaseII/CMS_PhaseII_200PU_v04.tcl

    r190cfa0 r05a7228  
    123123  BTaggingPUPPIAK8
    124124
    125   TauTagging
    126   TauTaggingAK8
    127   TauTaggingPUPPI
    128   TauTaggingPUPPIAK8
     125  TauTaggingCutBased
     126  TauTaggingDNNMedium
     127  TauTaggingDNNTight
     128
     129  TauTaggingAK8CutBased
     130  TauTaggingAK8DNNMedium
     131  TauTaggingAK8DNNTight
     132
     133  TauTaggingPUPPICutBased
     134  TauTaggingPUPPIDNNMedium
     135  TauTaggingPUPPIDNNTight
     136
     137  TauTaggingPUPPIAK8CutBased
     138  TauTaggingPUPPIAK8DNNMedium
     139  TauTaggingPUPPIAK8DNNTight
    129140
    130141  TreeWriter
    … …  
    626637  }
    627638
    628   # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
     639  # assume 0.07 x 0.07 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
    629640
    630641  set PhiBins {}
    … …  
    633644  }
    634645
    635   # 0.02 unit in eta up to eta = 3
     646  # 0.07 unit in eta up to eta = 3
    636647  for {set i 1} {$i <= 21} {incr i} {
    637648    set eta [expr { -2.958 + $i * 0.0696}]
    … …  
    39914002
    39924003
    3993 module TauTagging TauTagging {
     4004module TauTagging TauTaggingCutBased {
    39944005  set ParticleInputArray Delphes/allParticles
    39954006  set PartonInputArray Delphes/partons
    … …  
    40014012
    40024013  set TauEtaMax 2.3
     4014
     4015  set BitNumber 0
    40034016
    40044017  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    … …  
    40124025}
    40134026
    4014 module TauTagging TauTaggingAK8 {
     4027
     4028module TauTagging TauTaggingDNNMedium {
     4029  set ParticleInputArray Delphes/allParticles
     4030  set PartonInputArray Delphes/partons
     4031  set JetInputArray JetEnergyScale/jets
     4032
     4033  set DeltaR 0.5
     4034
     4035  set TauPTMin 20.0
     4036
     4037  set TauEtaMax 3.0
     4038
     4039  set BitNumber 1
     4040
     4041  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4042
     4043  add EfficiencyFormula {0} {
     4044                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4045                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4046                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4047                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4048                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.011) +
     4049                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.007) +
     4050                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4051                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.004) +
     4052                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4053                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.021) +
     4054                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4055                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4056                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.010) +
     4057                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.006) +
     4058                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4059                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.004) +
     4060                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4061                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4062                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.026) +
     4063                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.018) +
     4064                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.012) +
     4065                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.008) +
     4066                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4067                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.004) +
     4068                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4069                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4070                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.028) +
     4071                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.020) +
     4072                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.014) +
     4073                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.009) +
     4074                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.005) +
     4075                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.005) +
     4076                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4077                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4078                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.031) +
     4079                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.025) +
     4080                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.018) +
     4081                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.013) +
     4082                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.009) +
     4083                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.009)
     4084                                 
     4085                             }
     4086
     4087  add EfficiencyFormula {15} {
     4088                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.643) +
     4089                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.800) +
     4090                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.846) +
     4091                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.877) +
     4092                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.928) +
     4093                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.940) +
     4094                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.953) +
     4095                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.953) +
     4096                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.634) +
     4097                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.795) +
     4098                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.837) +
     4099                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.876) +
     4100                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.910) +
     4101                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.929) +
     4102                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.973) +
     4103                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.973) +
     4104                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.625) +
     4105                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.777) +
     4106                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.821) +
     4107                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.857) +
     4108                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.869) +
     4109                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.864) +
     4110                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.838) +
     4111                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.838) +
     4112                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.653) +
     4113                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.793) +
     4114                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.829) +
     4115                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.871) +
     4116                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.876) +
     4117                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.869) +
     4118                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.856) +
     4119                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.856) +
     4120                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.638) +
     4121                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.754) +
     4122                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.791) +
     4123                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.846) +
     4124                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.860) +
     4125                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.850) +
     4126                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.778) +
     4127                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.778)
     4128                           }
     4129}
     4130
     4131module TauTagging TauTaggingDNNTight {
     4132  set ParticleInputArray Delphes/allParticles
     4133  set PartonInputArray Delphes/partons
     4134  set JetInputArray JetEnergyScale/jets
     4135
     4136  set DeltaR 0.5
     4137
     4138  set TauPTMin 20.0
     4139
     4140  set TauEtaMax 3.0
     4141
     4142  set BitNumber 2
     4143
     4144  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4145  add EfficiencyFormula {0} {
     4146                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4147                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.003) +
     4148                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4149                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4150                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4151                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.001) +
     4152                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4153                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.001) +
     4154                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4155                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4156                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4157                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4158                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4159                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4160                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4161                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.001) +
     4162                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4163                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4164                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4165                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4166                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4167                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4168                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4169                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.001) +
     4170                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4171                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4172                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4173                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4174                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4175                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4176                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.002) +
     4177                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.001) +
     4178                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4179                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.001) +
     4180                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.002) +
     4181                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4182                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4183                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4184                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.003) +
     4185                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.002)
     4186                             }
     4187
     4188  add EfficiencyFormula {15} {
     4189                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.402) +
     4190                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.560) +
     4191                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.646) +
     4192                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.711) +
     4193                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.761) +
     4194                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.775) +
     4195                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.860) +
     4196                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.860) +
     4197                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.379) +
     4198                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.518) +
     4199                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.606) +
     4200                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.693) +
     4201                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.728) +
     4202                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.770) +
     4203                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.818) +
     4204                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.818) +
     4205                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.339) +
     4206                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.409) +
     4207                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.493) +
     4208                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.610) +
     4209                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.659) +
     4210                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.678) +
     4211                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.631) +
     4212                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.631) +
     4213                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.396) +
     4214                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.450) +
     4215                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.510) +
     4216                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.627) +
     4217                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.681) +
     4218                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.705) +
     4219                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.644) +
     4220                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.644) +
     4221                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.364) +
     4222                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.392) +
     4223                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.445) +
     4224                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.579) +
     4225                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.621) +
     4226                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.627) +
     4227                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.622) +
     4228                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.622)
     4229                           }
     4230}
     4231
     4232
     4233module TauTagging TauTaggingAK8CutBased {
    40154234  set ParticleInputArray Delphes/allParticles
    40164235  set PartonInputArray Delphes/partons
    … …  
    40224241
    40234242  set TauEtaMax 2.3
     4243
     4244  set BitNumber 0
    40244245
    40254246  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    … …  
    40334254}
    40344255
    4035 module TauTagging TauTaggingPUPPI {
     4256
     4257module TauTagging TauTaggingAK8DNNMedium {
     4258  set ParticleInputArray Delphes/allParticles
     4259  set PartonInputArray Delphes/partons
     4260  set JetInputArray JetEnergyScaleAK8/jets
     4261
     4262  set DeltaR 0.5
     4263
     4264  set TauPTMin 20.0
     4265
     4266  set TauEtaMax 3.0
     4267
     4268  set BitNumber 1
     4269
     4270  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4271  add EfficiencyFormula {0} {
     4272                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4273                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4274                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4275                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4276                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.011) +
     4277                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.007) +
     4278                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4279                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.004) +
     4280                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4281                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.021) +
     4282                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4283                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4284                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.010) +
     4285                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.006) +
     4286                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4287                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.004) +
     4288                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4289                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4290                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.026) +
     4291                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.018) +
     4292                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.012) +
     4293                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.008) +
     4294                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4295                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.004) +
     4296                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4297                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4298                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.028) +
     4299                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.020) +
     4300                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.014) +
     4301                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.009) +
     4302                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.005) +
     4303                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.005) +
     4304                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4305                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4306                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.031) +
     4307                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.025) +
     4308                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.018) +
     4309                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.013) +
     4310                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.009) +
     4311                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.009)
     4312                                 
     4313                             }
     4314
     4315  add EfficiencyFormula {15} {
     4316                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.643) +
     4317                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.800) +
     4318                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.846) +
     4319                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.877) +
     4320                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.928) +
     4321                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.940) +
     4322                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.953) +
     4323                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.953) +
     4324                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.634) +
     4325                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.795) +
     4326                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.837) +
     4327                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.876) +
     4328                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.910) +
     4329                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.929) +
     4330                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.973) +
     4331                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.973) +
     4332                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.625) +
     4333                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.777) +
     4334                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.821) +
     4335                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.857) +
     4336                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.869) +
     4337                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.864) +
     4338                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.838) +
     4339                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.838) +
     4340                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.653) +
     4341                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.793) +
     4342                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.829) +
     4343                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.871) +
     4344                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.876) +
     4345                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.869) +
     4346                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.856) +
     4347                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.856) +
     4348                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.638) +
     4349                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.754) +
     4350                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.791) +
     4351                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.846) +
     4352                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.860) +
     4353                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.850) +
     4354                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.778) +
     4355                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.778)
     4356                           }
     4357
     4358
     4359
     4360}
     4361
     4362module TauTagging TauTaggingAK8DNNTight {
     4363  set ParticleInputArray Delphes/allParticles
     4364  set PartonInputArray Delphes/partons
     4365  set JetInputArray JetEnergyScaleAK8/jets
     4366
     4367  set DeltaR 0.5
     4368
     4369  set TauPTMin 20.0
     4370
     4371  set TauEtaMax 3.0
     4372
     4373  set BitNumber 2
     4374
     4375  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4376  add EfficiencyFormula {0} {
     4377                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4378                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.003) +
     4379                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4380                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4381                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4382                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.001) +
     4383                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4384                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.001) +
     4385                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4386                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4387                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4388                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4389                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4390                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4391                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4392                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.001) +
     4393                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4394                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4395                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4396                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4397                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4398                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4399                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4400                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.001) +
     4401                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4402                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4403                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4404                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4405                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4406                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4407                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.002) +
     4408                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.001) +
     4409                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4410                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.001) +
     4411                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.002) +
     4412                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4413                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4414                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4415                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.003) +
     4416                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.002)
     4417                             }
     4418
     4419  add EfficiencyFormula {15} {
     4420                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.402) +
     4421                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.560) +
     4422                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.646) +
     4423                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.711) +
     4424                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.761) +
     4425                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.775) +
     4426                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.860) +
     4427                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.860) +
     4428                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.379) +
     4429                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.518) +
     4430                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.606) +
     4431                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.693) +
     4432                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.728) +
     4433                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.770) +
     4434                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.818) +
     4435                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.818) +
     4436                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.339) +
     4437                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.409) +
     4438                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.493) +
     4439                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.610) +
     4440                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.659) +
     4441                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.678) +
     4442                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.631) +
     4443                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.631) +
     4444                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.396) +
     4445                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.450) +
     4446                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.510) +
     4447                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.627) +
     4448                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.681) +
     4449                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.705) +
     4450                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.644) +
     4451                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.644) +
     4452                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.364) +
     4453                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.392) +
     4454                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.445) +
     4455                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.579) +
     4456                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.621) +
     4457                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.627) +
     4458                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.622) +
     4459                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.622)
     4460                           }
     4461}
     4462
     4463
     4464
     4465
     4466
     4467module TauTagging TauTaggingPUPPICutBased {
    40364468  set ParticleInputArray Delphes/allParticles
    40374469  set PartonInputArray Delphes/partons
    … …  
    40434475
    40444476  set TauEtaMax 2.3
     4477
     4478  set BitNumber 0
    40454479
    40464480  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    … …  
    40544488}
    40554489
    4056 module TauTagging TauTaggingPUPPIAK8 {
     4490
     4491module TauTagging TauTaggingPUPPIDNNMedium {
     4492  set ParticleInputArray Delphes/allParticles
     4493  set PartonInputArray Delphes/partons
     4494  set JetInputArray JetEnergyScalePUPPI/jets
     4495
     4496  set DeltaR 0.5
     4497
     4498  set TauPTMin 20.0
     4499
     4500  set TauEtaMax 3.0
     4501
     4502  set BitNumber 1
     4503
     4504  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4505  add EfficiencyFormula {0} {
     4506                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4507                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4508                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4509                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4510                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.011) +
     4511                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.007) +
     4512                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4513                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.004) +
     4514                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4515                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.021) +
     4516                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4517                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4518                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.010) +
     4519                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.006) +
     4520                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4521                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.004) +
     4522                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4523                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4524                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.026) +
     4525                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.018) +
     4526                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.012) +
     4527                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.008) +
     4528                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4529                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.004) +
     4530                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4531                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4532                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.028) +
     4533                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.020) +
     4534                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.014) +
     4535                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.009) +
     4536                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.005) +
     4537                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.005) +
     4538                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4539                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4540                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.031) +
     4541                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.025) +
     4542                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.018) +
     4543                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.013) +
     4544                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.009) +
     4545                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.009)
     4546                                 
     4547                             }
     4548
     4549  add EfficiencyFormula {15} {
     4550                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.643) +
     4551                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.800) +
     4552                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.846) +
     4553                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.877) +
     4554                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.928) +
     4555                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.940) +
     4556                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.953) +
     4557                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.953) +
     4558                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.634) +
     4559                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.795) +
     4560                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.837) +
     4561                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.876) +
     4562                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.910) +
     4563                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.929) +
     4564                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.973) +
     4565                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.973) +
     4566                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.625) +
     4567                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.777) +
     4568                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.821) +
     4569                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.857) +
     4570                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.869) +
     4571                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.864) +
     4572                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.838) +
     4573                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.838) +
     4574                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.653) +
     4575                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.793) +
     4576                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.829) +
     4577                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.871) +
     4578                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.876) +
     4579                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.869) +
     4580                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.856) +
     4581                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.856) +
     4582                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.638) +
     4583                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.754) +
     4584                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.791) +
     4585                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.846) +
     4586                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.860) +
     4587                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.850) +
     4588                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.778) +
     4589                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.778)
     4590                           }
     4591
     4592}
     4593
     4594
     4595
     4596module TauTagging TauTaggingPUPPIDNNTight {
     4597  set ParticleInputArray Delphes/allParticles
     4598  set PartonInputArray Delphes/partons
     4599  set JetInputArray JetEnergyScalePUPPI/jets
     4600
     4601  set DeltaR 0.5
     4602
     4603  set TauPTMin 20.0
     4604
     4605  set TauEtaMax 3.0
     4606
     4607  set BitNumber 2
     4608
     4609  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4610  add EfficiencyFormula {0} {
     4611                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4612                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.003) +
     4613                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4614                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4615                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4616                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.001) +
     4617                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4618                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.001) +
     4619                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4620                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4621                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4622                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4623                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4624                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4625                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4626                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.001) +
     4627                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4628                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4629                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4630                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4631                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4632                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4633                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4634                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.001) +
     4635                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4636                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4637                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4638                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4639                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4640                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4641                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.002) +
     4642                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.001) +
     4643                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4644                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.001) +
     4645                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.002) +
     4646                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4647                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4648                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4649                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.003) +
     4650                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.002)
     4651                             }
     4652
     4653  add EfficiencyFormula {15} {
     4654                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.402) +
     4655                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.560) +
     4656                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.646) +
     4657                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.711) +
     4658                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.761) +
     4659                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.775) +
     4660                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.860) +
     4661                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.860) +
     4662                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.379) +
     4663                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.518) +
     4664                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.606) +
     4665                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.693) +
     4666                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.728) +
     4667                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.770) +
     4668                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.818) +
     4669                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.818) +
     4670                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.339) +
     4671                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.409) +
     4672                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.493) +
     4673                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.610) +
     4674                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.659) +
     4675                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.678) +
     4676                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.631) +
     4677                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.631) +
     4678                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.396) +
     4679                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.450) +
     4680                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.510) +
     4681                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.627) +
     4682                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.681) +
     4683                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.705) +
     4684                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.644) +
     4685                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.644) +
     4686                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.364) +
     4687                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.392) +
     4688                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.445) +
     4689                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.579) +
     4690                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.621) +
     4691                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.627) +
     4692                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.622) +
     4693                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.622)
     4694                           }
     4695
     4696}
     4697
     4698
     4699
     4700module TauTagging TauTaggingPUPPIAK8CutBased {
    40574701  set ParticleInputArray Delphes/allParticles
    40584702  set PartonInputArray Delphes/partons
    … …  
    40654709  set TauEtaMax 2.3
    40664710
     4711  set BitNumber 0
     4712
    40674713  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    40684714
    4069   add EfficiencyFormula {0}  { (abs(eta) < 2.3) * ((( -0.00621816+0.00130097*pt-2.19642e-5*pt^2+1.49393e-7*pt^3-4.58972e-10*pt^4+5.27983e-13*pt^5 )) * (pt<250) + 0.0032*(pt>250)) + \
    4070                                (abs(eta) > 2.3) * (0.000)
     4715
     4716}
     4717
     4718module TauTagging TauTaggingPUPPIAK8DNNMedium {
     4719  set ParticleInputArray Delphes/allParticles
     4720  set PartonInputArray Delphes/partons
     4721  set JetInputArray JetEnergyScalePUPPIAK8/jets
     4722
     4723  set DeltaR 0.5
     4724
     4725  set TauPTMin 20.0
     4726
     4727  set TauEtaMax 3.0
     4728
     4729  set BitNumber 1
     4730
     4731  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4732  add EfficiencyFormula {0} {
     4733                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4734                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4735                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4736                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4737                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.011) +
     4738                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.007) +
     4739                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4740                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.004) +
     4741                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.005) +
     4742                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.021) +
     4743                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.023) +
     4744                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.015) +
     4745                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.010) +
     4746                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.006) +
     4747                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4748                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.004) +
     4749                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4750                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4751                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.026) +
     4752                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.018) +
     4753                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.012) +
     4754                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.008) +
     4755                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.004) +
     4756                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.004) +
     4757                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4758                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4759                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.028) +
     4760                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.020) +
     4761                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.014) +
     4762                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.009) +
     4763                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.005) +
     4764                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.005) +
     4765                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.006) +
     4766                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.022) +
     4767                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.031) +
     4768                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.025) +
     4769                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.018) +
     4770                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.013) +
     4771                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.009) +
     4772                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.009)
     4773                                 
    40714774                             }
    4072   add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.97*0.77*( (0.32 + 0.01*pt - 0.000054*pt*pt )*(pt<100)+0.78*(pt>100) ) + \
    4073                                (abs(eta) > 2.3) * (0.000)
     4775
     4776  add EfficiencyFormula {15} {
     4777                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.643) +
     4778                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.800) +
     4779                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.846) +
     4780                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.877) +
     4781                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.928) +
     4782                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.940) +
     4783                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.953) +
     4784                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.953) +
     4785                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.634) +
     4786                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.795) +
     4787                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.837) +
     4788                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.876) +
     4789                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.910) +
     4790                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.929) +
     4791                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.973) +
     4792                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.973) +
     4793                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.625) +
     4794                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.777) +
     4795                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.821) +
     4796                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.857) +
     4797                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.869) +
     4798                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.864) +
     4799                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.838) +
     4800                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.838) +
     4801                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.653) +
     4802                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.793) +
     4803                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.829) +
     4804                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.871) +
     4805                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.876) +
     4806                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.869) +
     4807                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.856) +
     4808                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.856) +
     4809                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.638) +
     4810                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.754) +
     4811                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.791) +
     4812                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.846) +
     4813                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.860) +
     4814                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.850) +
     4815                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.778) +
     4816                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.778)
     4817                           }
     4818
     4819
     4820}
     4821
     4822module TauTagging TauTaggingPUPPIAK8DNNTight {
     4823  set ParticleInputArray Delphes/allParticles
     4824  set PartonInputArray Delphes/partons
     4825  set JetInputArray JetEnergyScalePUPPIAK8/jets
     4826
     4827  set DeltaR 0.5
     4828
     4829  set TauPTMin 20.0
     4830
     4831  set TauEtaMax 3.0
     4832
     4833  set BitNumber 2
     4834
     4835  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
     4836  add EfficiencyFormula {0} {
     4837                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4838                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.003) +
     4839                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4840                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4841                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4842                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.001) +
     4843                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4844                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.001) +
     4845                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4846                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4847                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4848                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4849                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4850                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4851                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4852                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.001) +
     4853                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4854                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4855                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4856                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.002) +
     4857                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.002) +
     4858                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.002) +
     4859                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.001) +
     4860                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.001) +
     4861                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4862                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.002) +
     4863                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.003) +
     4864                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4865                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4866                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4867                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.002) +
     4868                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.001) +
     4869                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.0005) +
     4870                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.001) +
     4871                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.002) +
     4872                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.003) +
     4873                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.003) +
     4874                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.003) +
     4875                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.003) +
     4876                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.002)
    40744877                             }
     4878
     4879  add EfficiencyFormula {15} {
     4880                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  20.00 && pt <=  30.00) * (0.402) +
     4881                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  30.00 && pt <=  40.00) * (0.560) +
     4882                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  40.00 && pt <=  60.00) * (0.646) +
     4883                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  60.00 && pt <=  80.00) * (0.711) +
     4884                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt >  80.00 && pt <= 100.00) * (0.761) +
     4885                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 100.00 && pt <= 150.00) * (0.775) +
     4886                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 150.00 && pt <= 200.00) * (0.860) +
     4887                                  (abs(eta) > 0.00 && abs(eta) <= 0.50) * (pt > 200.00)                 * (0.860) +
     4888                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  20.00 && pt <=  30.00) * (0.379) +
     4889                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  30.00 && pt <=  40.00) * (0.518) +
     4890                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  40.00 && pt <=  60.00) * (0.606) +
     4891                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  60.00 && pt <=  80.00) * (0.693) +
     4892                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt >  80.00 && pt <= 100.00) * (0.728) +
     4893                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 100.00 && pt <= 150.00) * (0.770) +
     4894                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 150.00 && pt <= 200.00) * (0.818) +
     4895                                  (abs(eta) > 0.50 && abs(eta) <= 1.00) * (pt > 200.00)                 * (0.818) +
     4896                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  20.00 && pt <=  30.00) * (0.339) +
     4897                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  30.00 && pt <=  40.00) * (0.409) +
     4898                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  40.00 && pt <=  60.00) * (0.493) +
     4899                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  60.00 && pt <=  80.00) * (0.610) +
     4900                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt >  80.00 && pt <= 100.00) * (0.659) +
     4901                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 100.00 && pt <= 150.00) * (0.678) +
     4902                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 150.00 && pt <= 200.00) * (0.631) +
     4903                                  (abs(eta) > 1.00 && abs(eta) <= 1.60) * (pt > 200.00)                 * (0.631) +
     4904                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  20.00 && pt <=  30.00) * (0.396) +
     4905                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  30.00 && pt <=  40.00) * (0.450) +
     4906                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  40.00 && pt <=  60.00) * (0.510) +
     4907                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  60.00 && pt <=  80.00) * (0.627) +
     4908                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt >  80.00 && pt <= 100.00) * (0.681) +
     4909                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 100.00 && pt <= 150.00) * (0.705) +
     4910                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 150.00 && pt <= 200.00) * (0.644) +
     4911                                  (abs(eta) > 1.60 && abs(eta) <= 2.10) * (pt > 200.00)                 * (0.644) +
     4912                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  20.00 && pt <=  30.00) * (0.364) +
     4913                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  30.00 && pt <=  40.00) * (0.392) +
     4914                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  40.00 && pt <=  60.00) * (0.445) +
     4915                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  60.00 && pt <=  80.00) * (0.579) +
     4916                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt >  80.00 && pt <= 100.00) * (0.621) +
     4917                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 100.00 && pt <= 150.00) * (0.627) +
     4918                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 150.00 && pt <= 200.00) * (0.622) +
     4919                                  (abs(eta) > 2.10 && abs(eta) <= 3.00) * (pt > 200.00)                 * (0.622)
     4920                           }
     4921
    40754922}
    40764923
Note: See TracChangeset for help on using the changeset viewer.