Fork me on GitHub

Changes in / [fbad4c7:792092a] in git


Ignore:
Files:
2 deleted
5 edited

Legend:

Unmodified
Added
Removed
  • Makefile

    rfbad4c7 r792092a  
    379379        modules/TaggingParticlesSkimmer.h \
    380380        modules/PileUpJetID.h \
    381         modules/PhotonID.h \
    382381        modules/ConstituentFilter.h \
    383382        modules/StatusPidFilter.h \
     
    820819        external/ExRootAnalysis/ExRootFilter.h \
    821820        external/ExRootAnalysis/ExRootClassifier.h
    822 tmp/modules/PhotonID.$(ObjSuf): \
    823         modules/PhotonID.$(SrcSuf) \
    824         modules/PhotonID.h \
    825         classes/DelphesClasses.h \
    826         classes/DelphesFactory.h \
    827         classes/DelphesFormula.h \
    828         external/ExRootAnalysis/ExRootResult.h \
    829         external/ExRootAnalysis/ExRootFilter.h \
    830         external/ExRootAnalysis/ExRootClassifier.h
    831821tmp/modules/PileUpJetID.$(ObjSuf): \
    832822        modules/PileUpJetID.$(SrcSuf) \
     
    10761066        tmp/modules/PdgCodeFilter.$(ObjSuf) \
    10771067        tmp/modules/PhotonConversions.$(ObjSuf) \
    1078         tmp/modules/PhotonID.$(ObjSuf) \
    10791068        tmp/modules/PileUpJetID.$(ObjSuf) \
    10801069        tmp/modules/PileUpMerger.$(ObjSuf) \
     
    19411930        @touch $@
    19421931
    1943 modules/PhotonID.h: \
    1944         classes/DelphesModule.h
    1945         @touch $@
    1946 
    19471932external/fastjet/tools/Pruner.hh: \
    19481933        external/fastjet/ClusterSequence.hh \
  • cards/CMS_PhaseII/CMS_PhaseII_200PU_v03.tcl

    rfbad4c7 r792092a  
    22#  Main authors: Michele Selvaggi (CERN)
    33#
    4 #  Released on: May 2018
     4#  Released on: Jun 26 - 2017
    55#
    6 #  Version: v04
     6#  Version: v02
    77#
    88#
     
    5454  EFlowFilterCHS
    5555
    56   GenParticleFilter
    5756  PhotonFilter
     57
    5858  PhotonCloner
    59 
    6059  PhotonIsolation
    6160  PhotonIsolationCHS
     61  PhotonEfficiency
    6262  PhotonEfficiencyCHS
    6363
    64   PhotonLooseID
    65   PhotonTightID
    66  
    6764  ElectronCloner
    6865  ElectronIsolation
     
    124121  TauTaggingPUPPI
    125122  TauTaggingPUPPIAK8
     123
     124  GenParticleFilter
    126125
    127126  TreeWriter
     
    12431242
    12441243
     1244#####################
     1245# Photon efficiency #
     1246#####################
     1247
     1248module Efficiency PhotonEfficiency {
     1249
     1250  ## input particles
     1251  set InputArray PhotonIsolation/photons
     1252  ## output particles
     1253  set OutputArray photons
     1254  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
     1255  # efficiency formula for photons
     1256  set EfficiencyFormula {                      (pt <= 10.0) * (0.00) + \
     1257                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.9635) + \
     1258         (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0)  * (0.9624) + \
     1259         (abs(eta) > 4.0)                                   * (0.00)}
     1260
     1261}
     1262
    12451263
    12461264#####################
     
    12621280
    12631281}
    1264 
    1265 
    1266 #####################
    1267 # Photon Id Loose   #
    1268 #####################
    1269 
    1270 module PhotonID PhotonLooseID {
    1271 
    1272   ## input particles
    1273   set InputPhotonArray PhotonIsolation/photons
    1274 
    1275   ## gen particles
    1276   set InputGenArray GenParticleFilter/filteredParticles
    1277 
    1278  
    1279   ## min pt to be considered, make sure this threshold is higher than threshold in particle filter
    1280   set PTMin  5.0
    1281 
    1282   ## this decides whether prompt or not (to be tuned, since FS and delphes have different isolation profiles)
    1283   set RelIsoMax 0.3
    1284 
    1285   ## output particles
    1286   set OutputArray photons
    1287 
    1288 
    1289   set PromptFormula {
    1290                           (pt <= 5.0) * (0.00) +
    1291                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 5.0 && pt <= 10.0) * (0.0454) +
    1292                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 5.0 && pt <= 10.0) * (0.0802) +
    1293                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 5.0 && pt <= 10.0) * (0.2095) +
    1294                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 5.0 && pt <= 10.0) * (0.0829) +
    1295                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 5.0 && pt <= 10.0) * (0.1510) +
    1296                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 10.0) * (0.0654) +
    1297                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 5.0 && pt <= 10.0) * (0.0127) +
    1298                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 5.0 && pt <= 10.0) * (0.0100) +
    1299                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 5.0 && pt <= 10.0) * (0.0249) +
    1300                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 10.0 && pt <= 20.0) * (0.6489) +
    1301                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 10.0 && pt <= 20.0) * (0.6649) +
    1302                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 10.0 && pt <= 20.0) * (0.6969) +
    1303                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 10.0 && pt <= 20.0) * (0.4604) +
    1304                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.5952) +
    1305                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.4665) +
    1306                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.2457) +
    1307                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 10.0 && pt <= 20.0) * (0.1481) +
    1308                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.2373) +
    1309                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 20.0 && pt <= 30.0) * (0.8522) +
    1310                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 20.0 && pt <= 30.0) * (0.8158) +
    1311                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 20.0 && pt <= 30.0) * (0.7961) +
    1312                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 20.0 && pt <= 30.0) * (0.5824) +
    1313                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 30.0) * (0.7721) +
    1314                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 30.0) * (0.7198) +
    1315                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 30.0) * (0.6192) +
    1316                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 20.0 && pt <= 30.0) * (0.6148) +
    1317                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30.0) * (0.5670) +
    1318                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.8954) +
    1319                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 30.0 && pt <= 40.0) * (0.8646) +
    1320                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 30.0 && pt <= 40.0) * (0.8257) +
    1321                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 30.0 && pt <= 40.0) * (0.6548) +
    1322                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.8366) +
    1323                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.8247) +
    1324                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 30.0 && pt <= 40.0) * (0.8556) +
    1325                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 30.0 && pt <= 40.0) * (0.8302) +
    1326                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40.0) * (0.6896) +
    1327                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 40.0 && pt <= 60.0) * (0.9275) +
    1328                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 40.0 && pt <= 60.0) * (0.9133) +
    1329                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 40.0 && pt <= 60.0) * (0.8604) +
    1330                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 40.0 && pt <= 60.0) * (0.7076) +
    1331                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 60.0) * (0.9058) +
    1332                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 60.0) * (0.9142) +
    1333                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 40.0 && pt <= 60.0) * (0.9426) +
    1334                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 40.0 && pt <= 60.0) * (0.9151) +
    1335                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 60.0) * (0.7297) +
    1336                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 60.0 && pt <= 90.0) * (0.9427) +
    1337                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 60.0 && pt <= 90.0) * (0.9333) +
    1338                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 60.0 && pt <= 90.0) * (0.8679) +
    1339                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 60.0 && pt <= 90.0) * (0.7325) +
    1340                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 60.0 && pt <= 90.0) * (0.9460) +
    1341                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 60.0 && pt <= 90.0) * (0.9486) +
    1342                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 60.0 && pt <= 90.0) * (0.9682) +
    1343                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 60.0 && pt <= 90.0) * (0.9407) +
    1344                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 90.0) * (0.7476) +
    1345                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 90.0 && pt <= 120.0) * (0.9471) +
    1346                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 90.0 && pt <= 120.0) * (0.9381) +
    1347                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 90.0 && pt <= 120.0) * (0.8795) +
    1348                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 90.0 && pt <= 120.0) * (0.7483) +
    1349                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 90.0 && pt <= 120.0) * (0.9584) +
    1350                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 90.0 && pt <= 120.0) * (0.9522) +
    1351                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 90.0 && pt <= 120.0) * (0.9728) +
    1352                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 90.0 && pt <= 120.0) * (0.9626) +
    1353                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 120.0) * (0.7983) +
    1354                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 120.0 && pt <= 160.0) * (0.9507) +
    1355                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 120.0 && pt <= 160.0) * (0.9413) +
    1356                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 120.0 && pt <= 160.0) * (0.8737) +
    1357                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 120.0 && pt <= 160.0) * (0.7335) +
    1358                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 120.0 && pt <= 160.0) * (0.9538) +
    1359                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 120.0 && pt <= 160.0) * (0.9435) +
    1360                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 120.0 && pt <= 160.0) * (0.9874) +
    1361                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 120.0 && pt <= 160.0) * (0.9812) +
    1362                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 160.0) * (0.8104) +
    1363                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 160.0 && pt <= 240.0) * (0.9490) +
    1364                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 160.0 && pt <= 240.0) * (0.9431) +
    1365                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 160.0 && pt <= 240.0) * (0.8814) +
    1366                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 160.0 && pt <= 240.0) * (0.7073) +
    1367                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 160.0 && pt <= 240.0) * (0.9371) +
    1368                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 160.0 && pt <= 240.0) * (0.9497) +
    1369                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 160.0 && pt <= 240.0) * (0.9930) +
    1370                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 160.0 && pt <= 240.0) * (0.9715) +
    1371                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 240.0) * (0.8121) +
    1372                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 240.0 && pt <= 300.0) * (0.9392) +
    1373                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 240.0 && pt <= 300.0) * (0.9413) +
    1374                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 240.0 && pt <= 300.0) * (0.8588) +
    1375                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 240.0 && pt <= 300.0) * (0.6232) +
    1376                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 240.0 && pt <= 300.0) * (0.8993) +
    1377                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 240.0 && pt <= 300.0) * (0.9485) +
    1378                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 240.0 && pt <= 300.0) * (0.9854) +
    1379                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 240.0 && pt <= 300.0) * (1.0000) +
    1380                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 240.0 && pt <= 300.0) * (0.8108) +
    1381                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 300.0 && pt <= 500.0) * (0.9227) +
    1382                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 300.0 && pt <= 500.0) * (0.9213) +
    1383                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 300.0 && pt <= 500.0) * (0.8613) +
    1384                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 300.0 && pt <= 500.0) * (0.5265) +
    1385                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 300.0 && pt <= 500.0) * (0.8159) +
    1386                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 300.0 && pt <= 500.0) * (0.9031) +
    1387                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 300.0 && pt <= 500.0) * (0.9791) +
    1388                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 300.0 && pt <= 500.0) * (1.0000) +
    1389                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 300.0 && pt <= 500.0) * (1.0000) +
    1390                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 500.0 && pt <= 14000.0) * (0.8381) +
    1391                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 500.0 && pt <= 14000.0) * (0.8670) +
    1392                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 500.0 && pt <= 14000.0) * (0.8603) +
    1393                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 500.0 && pt <= 14000.0) * (0.1091) +
    1394                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 500.0 && pt <= 14000.0) * (0.2821) +
    1395                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 500.0 && pt <= 14000.0) * (0.6818) +
    1396                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 500.0 && pt <= 14000.0) * (1.0000) +
    1397                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1398                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1399                           (abs(eta) > 3.0) * (0.00)
    1400   }
    1401 
    1402   set NonPromptFormula {
    1403                           (pt <= 5.0) * (0.00) +
    1404                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 5.0 && pt <= 10.0) * (0.0056) +
    1405                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 5.0 && pt <= 10.0) * (0.0331) +
    1406                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 5.0 && pt <= 10.0) * (0.1038) +
    1407                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 5.0 && pt <= 10.0) * (0.0287) +
    1408                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 5.0 && pt <= 10.0) * (0.0263) +
    1409                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 10.0) * (0.0081) +
    1410                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 5.0 && pt <= 10.0) * (0.0048) +
    1411                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 5.0 && pt <= 10.0) * (0.0040) +
    1412                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 5.0 && pt <= 10.0) * (0.0096) +
    1413                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 10.0 && pt <= 20.0) * (0.0858) +
    1414                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 10.0 && pt <= 20.0) * (0.1253) +
    1415                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 10.0 && pt <= 20.0) * (0.2731) +
    1416                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 10.0 && pt <= 20.0) * (0.1178) +
    1417                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.0613) +
    1418                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.0265) +
    1419                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.0212) +
    1420                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 10.0 && pt <= 20.0) * (0.0160) +
    1421                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.0530) +
    1422                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 20.0 && pt <= 30.0) * (0.1755) +
    1423                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 20.0 && pt <= 30.0) * (0.1583) +
    1424                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 20.0 && pt <= 30.0) * (0.2300) +
    1425                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 20.0 && pt <= 30.0) * (0.1228) +
    1426                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 30.0) * (0.0747) +
    1427                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 30.0) * (0.0331) +
    1428                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 30.0) * (0.0221) +
    1429                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 20.0 && pt <= 30.0) * (0.0237) +
    1430                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30.0) * (0.0876) +
    1431                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.1902) +
    1432                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 30.0 && pt <= 40.0) * (0.1746) +
    1433                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 30.0 && pt <= 40.0) * (0.2052) +
    1434                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 30.0 && pt <= 40.0) * (0.1182) +
    1435                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.0854) +
    1436                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.0465) +
    1437                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 30.0 && pt <= 40.0) * (0.0360) +
    1438                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 30.0 && pt <= 40.0) * (0.0275) +
    1439                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40.0) * (0.0803) +
    1440                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 40.0 && pt <= 60.0) * (0.1431) +
    1441                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 40.0 && pt <= 60.0) * (0.1403) +
    1442                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 40.0 && pt <= 60.0) * (0.1394) +
    1443                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 40.0 && pt <= 60.0) * (0.0842) +
    1444                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 60.0) * (0.0731) +
    1445                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 60.0) * (0.0446) +
    1446                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 40.0 && pt <= 60.0) * (0.0348) +
    1447                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 40.0 && pt <= 60.0) * (0.0206) +
    1448                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 60.0) * (0.0624) +
    1449                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 60.0 && pt <= 90.0) * (0.0623) +
    1450                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 60.0 && pt <= 90.0) * (0.0645) +
    1451                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 60.0 && pt <= 90.0) * (0.0760) +
    1452                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 60.0 && pt <= 90.0) * (0.0475) +
    1453                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 60.0 && pt <= 90.0) * (0.0374) +
    1454                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 60.0 && pt <= 90.0) * (0.0245) +
    1455                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 60.0 && pt <= 90.0) * (0.0178) +
    1456                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 60.0 && pt <= 90.0) * (0.0131) +
    1457                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 90.0) * (0.0478) +
    1458                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 90.0 && pt <= 120.0) * (0.0107) +
    1459                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 90.0 && pt <= 120.0) * (0.0147) +
    1460                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 90.0 && pt <= 120.0) * (0.0498) +
    1461                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 90.0 && pt <= 120.0) * (0.0233) +
    1462                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 90.0 && pt <= 120.0) * (0.0150) +
    1463                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 90.0 && pt <= 120.0) * (0.0100) +
    1464                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 90.0 && pt <= 120.0) * (0.0097) +
    1465                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 90.0 && pt <= 120.0) * (0.0092) +
    1466                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 120.0) * (0.0519) +
    1467                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 120.0 && pt <= 160.0) * (0.0054) +
    1468                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 120.0 && pt <= 160.0) * (0.0070) +
    1469                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 120.0 && pt <= 160.0) * (0.0386) +
    1470                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 120.0 && pt <= 160.0) * (0.0101) +
    1471                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 120.0 && pt <= 160.0) * (0.0087) +
    1472                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 120.0 && pt <= 160.0) * (0.0055) +
    1473                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 120.0 && pt <= 160.0) * (0.0115) +
    1474                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 120.0 && pt <= 160.0) * (0.0052) +
    1475                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 160.0) * (0.0606) +
    1476                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 160.0 && pt <= 240.0) * (0.0036) +
    1477                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 160.0 && pt <= 240.0) * (0.0057) +
    1478                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 160.0 && pt <= 240.0) * (0.0416) +
    1479                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 160.0 && pt <= 240.0) * (0.0050) +
    1480                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 160.0 && pt <= 240.0) * (0.0034) +
    1481                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 160.0 && pt <= 240.0) * (0.0031) +
    1482                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 160.0 && pt <= 240.0) * (0.0129) +
    1483                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 160.0 && pt <= 240.0) * (0.0136) +
    1484                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 240.0) * (0.0883) +
    1485                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 240.0 && pt <= 300.0) * (0.0028) +
    1486                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 240.0 && pt <= 300.0) * (0.0051) +
    1487                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 240.0 && pt <= 300.0) * (0.0526) +
    1488                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 240.0 && pt <= 300.0) * (0.0006) +
    1489                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 240.0 && pt <= 300.0) * (0.0026) +
    1490                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 240.0 && pt <= 300.0) * (0.0036) +
    1491                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 240.0 && pt <= 300.0) * (0.0189) +
    1492                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 240.0 && pt <= 300.0) * (0.0346) +
    1493                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 240.0 && pt <= 300.0) * (0.1151) +
    1494                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 300.0 && pt <= 500.0) * (0.0029) +
    1495                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 300.0 && pt <= 500.0) * (0.0071) +
    1496                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 300.0 && pt <= 500.0) * (0.0774) +
    1497                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 300.0 && pt <= 500.0) * (0.0017) +
    1498                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 300.0 && pt <= 500.0) * (0.0010) +
    1499                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 300.0 && pt <= 500.0) * (0.0051) +
    1500                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 300.0 && pt <= 500.0) * (0.0252) +
    1501                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 300.0 && pt <= 500.0) * (0.0758) +
    1502                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 300.0 && pt <= 500.0) * (0.1282) +
    1503                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 500.0 && pt <= 14000.0) * (0.0084) +
    1504                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 500.0 && pt <= 14000.0) * (0.0215) +
    1505                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 500.0 && pt <= 14000.0) * (0.1074) +
    1506                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1507                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1508                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 500.0 && pt <= 14000.0) * (0.0068) +
    1509                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 500.0 && pt <= 14000.0) * (0.0250) +
    1510                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1511                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1512                           (abs(eta) > 3.0) * (0.00)
    1513   }
    1514  
    1515   set FakeFormula {
    1516                           (pt <= 5.0) * (0.00) +
    1517                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 5.0 && pt <= 10.0) * (0.0092) +
    1518                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 5.0 && pt <= 10.0) * (0.0865) +
    1519                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 5.0 && pt <= 10.0) * (0.2503) +
    1520                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 5.0 && pt <= 10.0) * (0.3277) +
    1521                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 5.0 && pt <= 10.0) * (0.2397) +
    1522                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 10.0) * (0.0999) +
    1523                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 5.0 && pt <= 10.0) * (0.0589) +
    1524                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 5.0 && pt <= 10.0) * (0.0374) +
    1525                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 5.0 && pt <= 10.0) * (0.0256) +
    1526                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 10.0 && pt <= 20.0) * (0.2002) +
    1527                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 10.0 && pt <= 20.0) * (0.2972) +
    1528                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 10.0 && pt <= 20.0) * (0.4924) +
    1529                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 10.0 && pt <= 20.0) * (0.3284) +
    1530                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.2417) +
    1531                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.1044) +
    1532                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.0586) +
    1533                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 10.0 && pt <= 20.0) * (0.0362) +
    1534                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.0667) +
    1535                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 20.0 && pt <= 30.0) * (0.2450) +
    1536                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 20.0 && pt <= 30.0) * (0.2372) +
    1537                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 20.0 && pt <= 30.0) * (0.2893) +
    1538                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 20.0 && pt <= 30.0) * (0.1780) +
    1539                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 30.0) * (0.1115) +
    1540                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 30.0) * (0.0480) +
    1541                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 30.0) * (0.0287) +
    1542                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 20.0 && pt <= 30.0) * (0.0294) +
    1543                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30.0) * (0.0843) +
    1544                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.1607) +
    1545                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 30.0 && pt <= 40.0) * (0.1826) +
    1546                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 30.0 && pt <= 40.0) * (0.1645) +
    1547                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 30.0 && pt <= 40.0) * (0.1225) +
    1548                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.0608) +
    1549                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.0323) +
    1550                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 30.0 && pt <= 40.0) * (0.0265) +
    1551                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 30.0 && pt <= 40.0) * (0.0263) +
    1552                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40.0) * (0.0692) +
    1553                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 40.0 && pt <= 60.0) * (0.1367) +
    1554                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 40.0 && pt <= 60.0) * (0.1694) +
    1555                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 40.0 && pt <= 60.0) * (0.1327) +
    1556                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 40.0 && pt <= 60.0) * (0.0719) +
    1557                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 60.0) * (0.0349) +
    1558                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 60.0) * (0.0232) +
    1559                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 40.0 && pt <= 60.0) * (0.0232) +
    1560                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 40.0 && pt <= 60.0) * (0.0159) +
    1561                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 60.0) * (0.0459) +
    1562                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 60.0 && pt <= 90.0) * (0.0604) +
    1563                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 60.0 && pt <= 90.0) * (0.1698) +
    1564                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 60.0 && pt <= 90.0) * (0.0678) +
    1565                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 60.0 && pt <= 90.0) * (0.0447) +
    1566                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 60.0 && pt <= 90.0) * (0.0107) +
    1567                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 60.0 && pt <= 90.0) * (0.0095) +
    1568                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 60.0 && pt <= 90.0) * (0.0074) +
    1569                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 60.0 && pt <= 90.0) * (0.0062) +
    1570                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 90.0) * (0.0266) +
    1571                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 90.0 && pt <= 120.0) * (0.0123) +
    1572                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 90.0 && pt <= 120.0) * (0.0275) +
    1573                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 90.0 && pt <= 120.0) * (0.0526) +
    1574                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 90.0 && pt <= 120.0) * (0.0180) +
    1575                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 90.0 && pt <= 120.0) * (0.0061) +
    1576                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 90.0 && pt <= 120.0) * (0.0066) +
    1577                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 90.0 && pt <= 120.0) * (0.0177) +
    1578                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 90.0 && pt <= 120.0) * (0.0023) +
    1579                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 120.0) * (0.0036) +
    1580                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 120.0 && pt <= 160.0) * (0.0101) +
    1581                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 120.0 && pt <= 160.0) * (0.0101) +
    1582                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 120.0 && pt <= 160.0) * (0.0629) +
    1583                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 120.0 && pt <= 160.0) * (0.0294) +
    1584                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 120.0 && pt <= 160.0) * (0.0037) +
    1585                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 120.0 && pt <= 160.0) * (0.0000) +
    1586                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 120.0 && pt <= 160.0) * (0.0014) +
    1587                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 120.0 && pt <= 160.0) * (0.0000) +
    1588                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 160.0) * (0.0301) +
    1589                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 160.0 && pt <= 240.0) * (0.0390) +
    1590                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 160.0 && pt <= 240.0) * (0.0181) +
    1591                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 160.0 && pt <= 240.0) * (0.0784) +
    1592                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 160.0 && pt <= 240.0) * (0.0209) +
    1593                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 160.0 && pt <= 240.0) * (0.0097) +
    1594                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 160.0 && pt <= 240.0) * (0.0017) +
    1595                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 160.0 && pt <= 240.0) * (0.0000) +
    1596                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 160.0 && pt <= 240.0) * (0.0000) +
    1597                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 240.0) * (0.0123) +
    1598                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1599                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 240.0 && pt <= 300.0) * (0.0256) +
    1600                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 240.0 && pt <= 300.0) * (0.0385) +
    1601                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1602                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1603                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1604                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1605                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1606                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1607                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 300.0 && pt <= 500.0) * (0.0093) +
    1608                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 300.0 && pt <= 500.0) * (0.0137) +
    1609                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 300.0 && pt <= 500.0) * (0.0556) +
    1610                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1611                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 300.0 && pt <= 500.0) * (0.0032) +
    1612                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1613                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 300.0 && pt <= 500.0) * (0.0222) +
    1614                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1615                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1616                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 500.0 && pt <= 14000.0) * (0.0435) +
    1617                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 500.0 && pt <= 14000.0) * (0.0667) +
    1618                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1619                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1620                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1621                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1622                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1623                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 500.0 && pt <= 14000.0) * (1.0000) +
    1624                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1625                           (abs(eta) > 3.0) * (0.00)
    1626   }
    1627 
    1628 
    1629 
    1630 
    1631 }
    1632 
    1633 
    1634 
    1635 
    1636 
    1637 #####################
    1638 # Photon Id Tight   #
    1639 #####################
    1640 
    1641 module PhotonID PhotonTightID {
    1642 
    1643   ## input particles
    1644   set InputPhotonArray PhotonIsolation/photons
    1645 
    1646   ## gen particles
    1647   set InputGenArray GenParticleFilter/filteredParticles
    1648 
    1649  
    1650   ## min pt to be considered, make sure this threshold is higher than threshold in particle filter
    1651   set PTMin  5.0
    1652 
    1653   ## this decides whether prompt or not (to be tuned, since FS and delphes have different isolation profiles)
    1654   set RelIsoMax 0.3
    1655 
    1656   ## output particles
    1657   set OutputArray photons
    1658 
    1659   set PromptFormula {
    1660                           (pt <= 5.0) * (0.00) +
    1661                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 5.0 && pt <= 10.0) * (0.0283) +
    1662                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 5.0 && pt <= 10.0) * (0.0488) +
    1663                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 5.0 && pt <= 10.0) * (0.1689) +
    1664                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 5.0 && pt <= 10.0) * (0.0461) +
    1665                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 5.0 && pt <= 10.0) * (0.0512) +
    1666                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 10.0) * (0.0220) +
    1667                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 5.0 && pt <= 10.0) * (0.0042) +
    1668                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 5.0 && pt <= 10.0) * (0.0027) +
    1669                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 5.0 && pt <= 10.0) * (0.0069) +
    1670                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 10.0 && pt <= 20.0) * (0.4877) +
    1671                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 10.0 && pt <= 20.0) * (0.5161) +
    1672                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 10.0 && pt <= 20.0) * (0.6377) +
    1673                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 10.0 && pt <= 20.0) * (0.2684) +
    1674                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.3159) +
    1675                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.2393) +
    1676                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.0865) +
    1677                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 10.0 && pt <= 20.0) * (0.0378) +
    1678                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.1272) +
    1679                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 20.0 && pt <= 30.0) * (0.7255) +
    1680                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 20.0 && pt <= 30.0) * (0.6751) +
    1681                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 20.0 && pt <= 30.0) * (0.7367) +
    1682                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 20.0 && pt <= 30.0) * (0.3868) +
    1683                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 30.0) * (0.5364) +
    1684                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 30.0) * (0.4758) +
    1685                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 30.0) * (0.3500) +
    1686                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 20.0 && pt <= 30.0) * (0.3435) +
    1687                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30.0) * (0.4122) +
    1688                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.7528) +
    1689                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 30.0 && pt <= 40.0) * (0.7179) +
    1690                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 30.0 && pt <= 40.0) * (0.7497) +
    1691                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 30.0 && pt <= 40.0) * (0.4773) +
    1692                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.6226) +
    1693                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.6053) +
    1694                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 30.0 && pt <= 40.0) * (0.6562) +
    1695                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 30.0 && pt <= 40.0) * (0.6261) +
    1696                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40.0) * (0.5612) +
    1697                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 40.0 && pt <= 60.0) * (0.8039) +
    1698                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 40.0 && pt <= 60.0) * (0.7947) +
    1699                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 40.0 && pt <= 60.0) * (0.7843) +
    1700                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 40.0 && pt <= 60.0) * (0.5655) +
    1701                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 60.0) * (0.7533) +
    1702                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 60.0) * (0.7679) +
    1703                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 40.0 && pt <= 60.0) * (0.8353) +
    1704                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 40.0 && pt <= 60.0) * (0.7880) +
    1705                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 60.0) * (0.6282) +
    1706                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 60.0 && pt <= 90.0) * (0.8470) +
    1707                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 60.0 && pt <= 90.0) * (0.8446) +
    1708                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 60.0 && pt <= 90.0) * (0.7862) +
    1709                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 60.0 && pt <= 90.0) * (0.6286) +
    1710                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 60.0 && pt <= 90.0) * (0.8542) +
    1711                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 60.0 && pt <= 90.0) * (0.8579) +
    1712                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 60.0 && pt <= 90.0) * (0.8929) +
    1713                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 60.0 && pt <= 90.0) * (0.8416) +
    1714                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 90.0) * (0.6528) +
    1715                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 90.0 && pt <= 120.0) * (0.8707) +
    1716                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 90.0 && pt <= 120.0) * (0.8687) +
    1717                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 90.0 && pt <= 120.0) * (0.8057) +
    1718                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 90.0 && pt <= 120.0) * (0.6636) +
    1719                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 90.0 && pt <= 120.0) * (0.8818) +
    1720                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 90.0 && pt <= 120.0) * (0.8773) +
    1721                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 90.0 && pt <= 120.0) * (0.8986) +
    1722                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 90.0 && pt <= 120.0) * (0.8924) +
    1723                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 120.0) * (0.7387) +
    1724                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 120.0 && pt <= 160.0) * (0.8886) +
    1725                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 120.0 && pt <= 160.0) * (0.8819) +
    1726                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 120.0 && pt <= 160.0) * (0.8153) +
    1727                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 120.0 && pt <= 160.0) * (0.6434) +
    1728                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 120.0 && pt <= 160.0) * (0.8760) +
    1729                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 120.0 && pt <= 160.0) * (0.8637) +
    1730                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 120.0 && pt <= 160.0) * (0.9447) +
    1731                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 120.0 && pt <= 160.0) * (0.9438) +
    1732                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 160.0) * (0.7750) +
    1733                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 160.0 && pt <= 240.0) * (0.8929) +
    1734                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 160.0 && pt <= 240.0) * (0.8895) +
    1735                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 160.0 && pt <= 240.0) * (0.8281) +
    1736                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 160.0 && pt <= 240.0) * (0.6062) +
    1737                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 160.0 && pt <= 240.0) * (0.8341) +
    1738                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 160.0 && pt <= 240.0) * (0.8703) +
    1739                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 160.0 && pt <= 240.0) * (0.9746) +
    1740                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 160.0 && pt <= 240.0) * (0.9620) +
    1741                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 240.0) * (0.7808) +
    1742                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 240.0 && pt <= 300.0) * (0.8866) +
    1743                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 240.0 && pt <= 300.0) * (0.8882) +
    1744                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 240.0 && pt <= 300.0) * (0.8217) +
    1745                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 240.0 && pt <= 300.0) * (0.4747) +
    1746                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 240.0 && pt <= 300.0) * (0.7404) +
    1747                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 240.0 && pt <= 300.0) * (0.8682) +
    1748                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 240.0 && pt <= 300.0) * (0.9766) +
    1749                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 240.0 && pt <= 300.0) * (1.0000) +
    1750                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 240.0 && pt <= 300.0) * (0.7393) +
    1751                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 300.0 && pt <= 500.0) * (0.8654) +
    1752                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 300.0 && pt <= 500.0) * (0.8735) +
    1753                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 300.0 && pt <= 500.0) * (0.8295) +
    1754                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 300.0 && pt <= 500.0) * (0.3144) +
    1755                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 300.0 && pt <= 500.0) * (0.5566) +
    1756                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 300.0 && pt <= 500.0) * (0.7813) +
    1757                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 300.0 && pt <= 500.0) * (0.9749) +
    1758                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 300.0 && pt <= 500.0) * (0.9375) +
    1759                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 300.0 && pt <= 500.0) * (1.0000) +
    1760                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 500.0 && pt <= 14000.0) * (0.7800) +
    1761                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 500.0 && pt <= 14000.0) * (0.8373) +
    1762                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 500.0 && pt <= 14000.0) * (0.8382) +
    1763                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 500.0 && pt <= 14000.0) * (0.0273) +
    1764                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 500.0 && pt <= 14000.0) * (0.0769) +
    1765                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 500.0 && pt <= 14000.0) * (0.6818) +
    1766                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 500.0 && pt <= 14000.0) * (1.0000) +
    1767                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1768                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1769                           (abs(eta) > 3.0) * (0.00)
    1770   }
    1771  
    1772   set NonPromptFormula {
    1773                           (pt <= 5.0) * (0.00) +
    1774                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 5.0 && pt <= 10.0) * (0.0020) +
    1775                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 5.0 && pt <= 10.0) * (0.0124) +
    1776                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 5.0 && pt <= 10.0) * (0.0649) +
    1777                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 5.0 && pt <= 10.0) * (0.0134) +
    1778                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 5.0 && pt <= 10.0) * (0.0077) +
    1779                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 10.0) * (0.0045) +
    1780                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 5.0 && pt <= 10.0) * (0.0045) +
    1781                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 5.0 && pt <= 10.0) * (0.0038) +
    1782                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 5.0 && pt <= 10.0) * (0.0043) +
    1783                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 10.0 && pt <= 20.0) * (0.0351) +
    1784                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 10.0 && pt <= 20.0) * (0.0516) +
    1785                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 10.0 && pt <= 20.0) * (0.1754) +
    1786                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 10.0 && pt <= 20.0) * (0.0561) +
    1787                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.0185) +
    1788                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.0142) +
    1789                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.0174) +
    1790                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 10.0 && pt <= 20.0) * (0.0122) +
    1791                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.0201) +
    1792                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 20.0 && pt <= 30.0) * (0.0875) +
    1793                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 20.0 && pt <= 30.0) * (0.0760) +
    1794                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 20.0 && pt <= 30.0) * (0.1383) +
    1795                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 20.0 && pt <= 30.0) * (0.0572) +
    1796                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 30.0) * (0.0236) +
    1797                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 30.0) * (0.0132) +
    1798                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 30.0) * (0.0110) +
    1799                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 20.0 && pt <= 30.0) * (0.0100) +
    1800                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30.0) * (0.0425) +
    1801                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.1038) +
    1802                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 30.0 && pt <= 40.0) * (0.0935) +
    1803                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 30.0 && pt <= 40.0) * (0.1231) +
    1804                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 30.0 && pt <= 40.0) * (0.0539) +
    1805                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.0309) +
    1806                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.0188) +
    1807                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 30.0 && pt <= 40.0) * (0.0159) +
    1808                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 30.0 && pt <= 40.0) * (0.0108) +
    1809                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40.0) * (0.0419) +
    1810                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 40.0 && pt <= 60.0) * (0.0810) +
    1811                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 40.0 && pt <= 60.0) * (0.0791) +
    1812                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 40.0 && pt <= 60.0) * (0.0806) +
    1813                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 40.0 && pt <= 60.0) * (0.0399) +
    1814                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 60.0) * (0.0317) +
    1815                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 60.0) * (0.0225) +
    1816                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 40.0 && pt <= 60.0) * (0.0172) +
    1817                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 40.0 && pt <= 60.0) * (0.0099) +
    1818                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 60.0) * (0.0305) +
    1819                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 60.0 && pt <= 90.0) * (0.0337) +
    1820                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 60.0 && pt <= 90.0) * (0.0354) +
    1821                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 60.0 && pt <= 90.0) * (0.0369) +
    1822                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 60.0 && pt <= 90.0) * (0.0225) +
    1823                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 60.0 && pt <= 90.0) * (0.0202) +
    1824                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 60.0 && pt <= 90.0) * (0.0132) +
    1825                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 60.0 && pt <= 90.0) * (0.0095) +
    1826                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 60.0 && pt <= 90.0) * (0.0057) +
    1827                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 90.0) * (0.0196) +
    1828                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 90.0 && pt <= 120.0) * (0.0056) +
    1829                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 90.0 && pt <= 120.0) * (0.0072) +
    1830                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 90.0 && pt <= 120.0) * (0.0222) +
    1831                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 90.0 && pt <= 120.0) * (0.0113) +
    1832                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 90.0 && pt <= 120.0) * (0.0087) +
    1833                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 90.0 && pt <= 120.0) * (0.0055) +
    1834                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 90.0 && pt <= 120.0) * (0.0047) +
    1835                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 90.0 && pt <= 120.0) * (0.0058) +
    1836                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 120.0) * (0.0251) +
    1837                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 120.0 && pt <= 160.0) * (0.0026) +
    1838                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 120.0 && pt <= 160.0) * (0.0034) +
    1839                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 120.0 && pt <= 160.0) * (0.0138) +
    1840                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 120.0 && pt <= 160.0) * (0.0036) +
    1841                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 120.0 && pt <= 160.0) * (0.0049) +
    1842                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 120.0 && pt <= 160.0) * (0.0025) +
    1843                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 120.0 && pt <= 160.0) * (0.0049) +
    1844                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 120.0 && pt <= 160.0) * (0.0016) +
    1845                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 160.0) * (0.0234) +
    1846                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 160.0 && pt <= 240.0) * (0.0016) +
    1847                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 160.0 && pt <= 240.0) * (0.0018) +
    1848                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 160.0 && pt <= 240.0) * (0.0152) +
    1849                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 160.0 && pt <= 240.0) * (0.0018) +
    1850                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 160.0 && pt <= 240.0) * (0.0015) +
    1851                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 160.0 && pt <= 240.0) * (0.0011) +
    1852                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 160.0 && pt <= 240.0) * (0.0051) +
    1853                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 160.0 && pt <= 240.0) * (0.0030) +
    1854                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 240.0) * (0.0463) +
    1855                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 240.0 && pt <= 300.0) * (0.0009) +
    1856                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 240.0 && pt <= 300.0) * (0.0013) +
    1857                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 240.0 && pt <= 300.0) * (0.0115) +
    1858                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1859                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 240.0 && pt <= 300.0) * (0.0012) +
    1860                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 240.0 && pt <= 300.0) * (0.0019) +
    1861                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 240.0 && pt <= 300.0) * (0.0087) +
    1862                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 240.0 && pt <= 300.0) * (0.0094) +
    1863                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 240.0 && pt <= 300.0) * (0.0712) +
    1864                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 300.0 && pt <= 500.0) * (0.0007) +
    1865                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 300.0 && pt <= 500.0) * (0.0011) +
    1866                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 300.0 && pt <= 500.0) * (0.0173) +
    1867                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1868                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 300.0 && pt <= 500.0) * (0.0003) +
    1869                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 300.0 && pt <= 500.0) * (0.0014) +
    1870                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 300.0 && pt <= 500.0) * (0.0101) +
    1871                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 300.0 && pt <= 500.0) * (0.0152) +
    1872                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 300.0 && pt <= 500.0) * (0.1026) +
    1873                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 500.0 && pt <= 14000.0) * (0.0010) +
    1874                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 500.0 && pt <= 14000.0) * (0.0026) +
    1875                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 500.0 && pt <= 14000.0) * (0.0261) +
    1876                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1877                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1878                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 500.0 && pt <= 14000.0) * (0.0027) +
    1879                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1880                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1881                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1882                           (abs(eta) > 3.0) * (0.00)
    1883   }
    1884  
    1885   set FakeFormula {
    1886                           (pt <= 5.0) * (0.00) +
    1887                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 5.0 && pt <= 10.0) * (0.0028) +
    1888                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 5.0 && pt <= 10.0) * (0.0295) +
    1889                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 5.0 && pt <= 10.0) * (0.1131) +
    1890                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 5.0 && pt <= 10.0) * (0.2226) +
    1891                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 5.0 && pt <= 10.0) * (0.1545) +
    1892                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 10.0) * (0.0555) +
    1893                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 5.0 && pt <= 10.0) * (0.0253) +
    1894                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 5.0 && pt <= 10.0) * (0.0255) +
    1895                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 5.0 && pt <= 10.0) * (0.0129) +
    1896                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 10.0 && pt <= 20.0) * (0.0919) +
    1897                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 10.0 && pt <= 20.0) * (0.1518) +
    1898                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 10.0 && pt <= 20.0) * (0.3449) +
    1899                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 10.0 && pt <= 20.0) * (0.2604) +
    1900                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.1814) +
    1901                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.0603) +
    1902                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.0252) +
    1903                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 10.0 && pt <= 20.0) * (0.0205) +
    1904                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.0271) +
    1905                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 20.0 && pt <= 30.0) * (0.1401) +
    1906                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 20.0 && pt <= 30.0) * (0.1331) +
    1907                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 20.0 && pt <= 30.0) * (0.1863) +
    1908                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 20.0 && pt <= 30.0) * (0.1275) +
    1909                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 30.0) * (0.0762) +
    1910                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 30.0) * (0.0241) +
    1911                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 30.0) * (0.0138) +
    1912                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 20.0 && pt <= 30.0) * (0.0131) +
    1913                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30.0) * (0.0444) +
    1914                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 30.0 && pt <= 40.0) * (0.0849) +
    1915                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 30.0 && pt <= 40.0) * (0.0940) +
    1916                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 30.0 && pt <= 40.0) * (0.0961) +
    1917                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 30.0 && pt <= 40.0) * (0.0754) +
    1918                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 30.0 && pt <= 40.0) * (0.0348) +
    1919                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 30.0 && pt <= 40.0) * (0.0173) +
    1920                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 30.0 && pt <= 40.0) * (0.0176) +
    1921                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 30.0 && pt <= 40.0) * (0.0124) +
    1922                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40.0) * (0.0386) +
    1923                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 40.0 && pt <= 60.0) * (0.0783) +
    1924                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 40.0 && pt <= 60.0) * (0.0763) +
    1925                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 40.0 && pt <= 60.0) * (0.0774) +
    1926                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 40.0 && pt <= 60.0) * (0.0461) +
    1927                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 40.0 && pt <= 60.0) * (0.0232) +
    1928                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 40.0 && pt <= 60.0) * (0.0159) +
    1929                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 40.0 && pt <= 60.0) * (0.0165) +
    1930                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 40.0 && pt <= 60.0) * (0.0121) +
    1931                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 60.0) * (0.0273) +
    1932                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 60.0 && pt <= 90.0) * (0.0302) +
    1933                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 60.0 && pt <= 90.0) * (0.0827) +
    1934                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 60.0 && pt <= 90.0) * (0.0339) +
    1935                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 60.0 && pt <= 90.0) * (0.0198) +
    1936                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 60.0 && pt <= 90.0) * (0.0050) +
    1937                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 60.0 && pt <= 90.0) * (0.0053) +
    1938                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 60.0 && pt <= 90.0) * (0.0040) +
    1939                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 60.0 && pt <= 90.0) * (0.0042) +
    1940                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 90.0) * (0.0189) +
    1941                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 90.0 && pt <= 120.0) * (0.0084) +
    1942                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 90.0 && pt <= 120.0) * (0.0142) +
    1943                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 90.0 && pt <= 120.0) * (0.0283) +
    1944                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 90.0 && pt <= 120.0) * (0.0090) +
    1945                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 90.0 && pt <= 120.0) * (0.0024) +
    1946                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 90.0 && pt <= 120.0) * (0.0062) +
    1947                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 90.0 && pt <= 120.0) * (0.0103) +
    1948                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 90.0 && pt <= 120.0) * (0.0000) +
    1949                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 120.0) * (0.0036) +
    1950                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 120.0 && pt <= 160.0) * (0.0056) +
    1951                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 120.0 && pt <= 160.0) * (0.0058) +
    1952                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 120.0 && pt <= 160.0) * (0.0343) +
    1953                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 120.0 && pt <= 160.0) * (0.0110) +
    1954                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 120.0 && pt <= 160.0) * (0.0028) +
    1955                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 120.0 && pt <= 160.0) * (0.0000) +
    1956                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 120.0 && pt <= 160.0) * (0.0014) +
    1957                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 120.0 && pt <= 160.0) * (0.0000) +
    1958                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 160.0) * (0.0150) +
    1959                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 160.0 && pt <= 240.0) * (0.0228) +
    1960                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 160.0 && pt <= 240.0) * (0.0023) +
    1961                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 160.0 && pt <= 240.0) * (0.0196) +
    1962                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 160.0 && pt <= 240.0) * (0.0000) +
    1963                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 160.0 && pt <= 240.0) * (0.0075) +
    1964                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 160.0 && pt <= 240.0) * (0.0009) +
    1965                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 160.0 && pt <= 240.0) * (0.0000) +
    1966                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 160.0 && pt <= 240.0) * (0.0000) +
    1967                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 240.0) * (0.0000) +
    1968                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1969                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1970                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1971                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1972                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1973                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1974                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1975                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1976                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 240.0 && pt <= 300.0) * (0.0000) +
    1977                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1978                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1979                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1980                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1981                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1982                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1983                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 300.0 && pt <= 500.0) * (0.0222) +
    1984                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1985                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 300.0 && pt <= 500.0) * (0.0000) +
    1986                           (abs(eta) > 0.0 && abs(eta) <= 0.8) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1987                           (abs(eta) > 0.8 && abs(eta) <= 1.4) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1988                           (abs(eta) > 1.4 && abs(eta) <= 1.48) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1989                           (abs(eta) > 1.48 && abs(eta) <= 1.6) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1990                           (abs(eta) > 1.6 && abs(eta) <= 2.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1991                           (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1992                           (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1993                           (abs(eta) > 2.8 && abs(eta) <= 2.9) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1994                           (abs(eta) > 2.9 && abs(eta) <= 3.0) * (pt > 500.0 && pt <= 14000.0) * (0.0000) +
    1995                           (abs(eta) > 3.0) * (0.00)
    1996   }
    1997  
    1998 
    1999 }
    2000 
    2001 
    2002 
    2003 
    2004 
    20051282
    20061283######################
     
    40873364#  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
    40883365
    4089   add Branch PhotonLooseID/photons PhotonLoose Photon
    4090   add Branch PhotonTightID/photons PhotonTight Photon
    4091 
     3366  add Branch PhotonEfficiency/photons Photon Photon
    40923367  add Branch ElectronEfficiency/electrons Electron Electron
    40933368  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    40943369  add Branch MuonTightIdEfficiency/muons MuonTight Muon
    40953370
    4096 #  add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
     3371  add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
    40973372  add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
    40983373  add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
  • classes/DelphesClasses.h

    rfbad4c7 r792092a  
    275275  Float_t SumPt; // isolation variable
    276276
    277   Int_t Status; // 1: prompt -- 2: non prompt -- 3: fake
    278 
    279277  static CompBase *fgCompare; //!
    280278  const CompBase *GetCompare() const { return fgCompare; }
     
    553551
    554552  Int_t IsConstituent;
     553
    555554  Int_t IsFromConversion;
    556555
  • modules/ModulesLinkDef.h

    rfbad4c7 r792092a  
    5757#include "modules/TaggingParticlesSkimmer.h"
    5858#include "modules/PileUpJetID.h"
    59 #include "modules/PhotonID.h"
    6059#include "modules/ConstituentFilter.h"
    6160#include "modules/StatusPidFilter.h"
     
    110109#pragma link C++ class TaggingParticlesSkimmer+;
    111110#pragma link C++ class PileUpJetID+;
    112 #pragma link C++ class PhotonID+;
    113111#pragma link C++ class ConstituentFilter+;
    114112#pragma link C++ class StatusPidFilter+;
  • modules/TreeWriter.cc

    rfbad4c7 r792092a  
    463463    const TLorentzVector &position = candidate->Position;
    464464
     465
    465466    pt = momentum.Pt();
    466467    cosTheta = TMath::Abs(momentum.CosTheta());
     
    475476    entry->PT = pt;
    476477    entry->E = momentum.E();
     478
    477479    entry->T = position.T()*1.0E-3/c_light;
    478480
     
    487489
    488490    entry->EhadOverEem = candidate->Eem > 0.0 ? candidate->Ehad/candidate->Eem : 999.9;
    489 
    490     // 1: prompt -- 2: non prompt -- 3: fake
    491     entry->Status = candidate->Status;
    492491
    493492    FillParticles(candidate, &entry->Particles);
Note: See TracChangeset for help on using the changeset viewer.