Changeset 05a7228 in git
- Timestamp:
- Apr 23, 2019, 6:29:32 PM (6 years ago)
- Branches:
- ImprovedOutputFile, Timing, master
- Children:
- 0453894
- Parents:
- 190cfa0
- File:
-
- 1 edited
-
cards/CMS_PhaseII/CMS_PhaseII_200PU_v04.tcl (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
cards/CMS_PhaseII/CMS_PhaseII_200PU_v04.tcl
r190cfa0 r05a7228 123 123 BTaggingPUPPIAK8 124 124 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 129 140 130 141 TreeWriter … … 626 637 } 627 638 628 # assume 0.0 2 x 0.02resolution 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) 629 640 630 641 set PhiBins {} … … 633 644 } 634 645 635 # 0.0 2unit in eta up to eta = 3646 # 0.07 unit in eta up to eta = 3 636 647 for {set i 1} {$i <= 21} {incr i} { 637 648 set eta [expr { -2.958 + $i * 0.0696}] … … 3991 4002 3992 4003 3993 module TauTagging TauTagging {4004 module TauTagging TauTaggingCutBased { 3994 4005 set ParticleInputArray Delphes/allParticles 3995 4006 set PartonInputArray Delphes/partons … … 4001 4012 4002 4013 set TauEtaMax 2.3 4014 4015 set BitNumber 0 4003 4016 4004 4017 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt} … … 4012 4025 } 4013 4026 4014 module TauTagging TauTaggingAK8 { 4027 4028 module 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 4131 module 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 4233 module TauTagging TauTaggingAK8CutBased { 4015 4234 set ParticleInputArray Delphes/allParticles 4016 4235 set PartonInputArray Delphes/partons … … 4022 4241 4023 4242 set TauEtaMax 2.3 4243 4244 set BitNumber 0 4024 4245 4025 4246 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt} … … 4033 4254 } 4034 4255 4035 module TauTagging TauTaggingPUPPI { 4256 4257 module 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 4362 module 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 4467 module TauTagging TauTaggingPUPPICutBased { 4036 4468 set ParticleInputArray Delphes/allParticles 4037 4469 set PartonInputArray Delphes/partons … … 4043 4475 4044 4476 set TauEtaMax 2.3 4477 4478 set BitNumber 0 4045 4479 4046 4480 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt} … … 4054 4488 } 4055 4489 4056 module TauTagging TauTaggingPUPPIAK8 { 4490 4491 module 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 4596 module 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 4700 module TauTagging TauTaggingPUPPIAK8CutBased { 4057 4701 set ParticleInputArray Delphes/allParticles 4058 4702 set PartonInputArray Delphes/partons … … 4065 4709 set TauEtaMax 2.3 4066 4710 4711 set BitNumber 0 4712 4067 4713 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt} 4068 4714 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 4718 module 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 4071 4774 } 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 4822 module 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) 4074 4877 } 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 4075 4922 } 4076 4923
Note:
See TracChangeset
for help on using the changeset viewer.
