Fork me on GitHub

Changeset 4689481 in git


Ignore:
Timestamp:
May 25, 2018, 8:28:31 PM (7 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
f37555a
Parents:
c3547ed (diff), fbad4c7 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.
Message:

Merge branch 'photonId'

Files:
2 added
5 edited

Legend:

Unmodified
Added
Removed
  • Makefile

    rc3547ed r4689481  
    379379        modules/TaggingParticlesSkimmer.h \
    380380        modules/PileUpJetID.h \
     381        modules/PhotonID.h \
    381382        modules/ConstituentFilter.h \
    382383        modules/StatusPidFilter.h \
     
    825826        classes/DelphesFactory.h \
    826827        classes/DelphesCylindricalFormula.h \
     828        external/ExRootAnalysis/ExRootResult.h \
     829        external/ExRootAnalysis/ExRootFilter.h \
     830        external/ExRootAnalysis/ExRootClassifier.h
     831tmp/modules/PhotonID.$(ObjSuf): \
     832        modules/PhotonID.$(SrcSuf) \
     833        modules/PhotonID.h \
     834        classes/DelphesClasses.h \
     835        classes/DelphesFactory.h \
     836        classes/DelphesFormula.h \
    827837        external/ExRootAnalysis/ExRootResult.h \
    828838        external/ExRootAnalysis/ExRootFilter.h \
     
    10771087        tmp/modules/PdgCodeFilter.$(ObjSuf) \
    10781088        tmp/modules/PhotonConversions.$(ObjSuf) \
     1089        tmp/modules/PhotonID.$(ObjSuf) \
    10791090        tmp/modules/PileUpJetID.$(ObjSuf) \
    10801091        tmp/modules/PileUpMerger.$(ObjSuf) \
     
    19411952        @touch $@
    19421953
     1954modules/PhotonID.h: \
     1955        classes/DelphesModule.h
     1956        @touch $@
     1957
    19431958external/fastjet/tools/Pruner.hh: \
    19441959        external/fastjet/ClusterSequence.hh \
  • cards/CMS_PhaseII/CMS_PhaseII_200PU_v03.tcl

    rc3547ed r4689481  
    22#  Main authors: Michele Selvaggi (CERN)
    33#
    4 #  Released on: Jun 26 - 2017
     4#  Released on: May 2018
    55#
    6 #  Version: v02
     6#  Version: v04
    77#
    88#
     
    5454  EFlowFilterCHS
    5555
     56  GenParticleFilter
    5657  PhotonFilter
    57 
    5858  PhotonCloner
     59
    5960  PhotonIsolation
    6061  PhotonIsolationCHS
    61   PhotonEfficiency
    6262  PhotonEfficiencyCHS
    6363
     64  PhotonLooseID
     65  PhotonTightID
     66 
    6467  ElectronCloner
    6568  ElectronIsolation
     
    121124  TauTaggingPUPPI
    122125  TauTaggingPUPPIAK8
    123 
    124   GenParticleFilter
    125126
    126127  TreeWriter
     
    12421243
    12431244
    1244 #####################
    1245 # Photon efficiency #
    1246 #####################
    1247 
    1248 module 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 
    12631245
    12641246#####################
     
    12801262
    12811263}
     1264
     1265
     1266#####################
     1267# Photon Id Loose   #
     1268#####################
     1269
     1270module 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
     1641module 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
    12822005
    12832006######################
     
    33644087#  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
    33654088
    3366   add Branch PhotonEfficiency/photons Photon Photon
     4089  add Branch PhotonLooseID/photons PhotonLoose Photon
     4090  add Branch PhotonTightID/photons PhotonTight Photon
     4091
    33674092  add Branch ElectronEfficiency/electrons Electron Electron
    33684093  add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
    33694094  add Branch MuonTightIdEfficiency/muons MuonTight Muon
    33704095
    3371   add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
     4096#  add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
    33724097  add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
    33734098  add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
  • classes/DelphesClasses.h

    rc3547ed r4689481  
    275275  Float_t SumPt; // isolation variable
    276276
     277  Int_t Status; // 1: prompt -- 2: non prompt -- 3: fake
     278
    277279  static CompBase *fgCompare; //!
    278280  const CompBase *GetCompare() const { return fgCompare; }
     
    551553
    552554  Int_t IsConstituent;
    553 
    554555  Int_t IsFromConversion;
    555556
  • modules/ModulesLinkDef.h

    rc3547ed r4689481  
    5757#include "modules/TaggingParticlesSkimmer.h"
    5858#include "modules/PileUpJetID.h"
     59#include "modules/PhotonID.h"
    5960#include "modules/ConstituentFilter.h"
    6061#include "modules/StatusPidFilter.h"
     
    109110#pragma link C++ class TaggingParticlesSkimmer+;
    110111#pragma link C++ class PileUpJetID+;
     112#pragma link C++ class PhotonID+;
    111113#pragma link C++ class ConstituentFilter+;
    112114#pragma link C++ class StatusPidFilter+;
  • modules/TreeWriter.cc

    rc3547ed r4689481  
    463463    const TLorentzVector &position = candidate->Position;
    464464
    465 
    466465    pt = momentum.Pt();
    467466    cosTheta = TMath::Abs(momentum.CosTheta());
     
    476475    entry->PT = pt;
    477476    entry->E = momentum.E();
    478 
    479477    entry->T = position.T()*1.0E-3/c_light;
    480478
     
    489487
    490488    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;
    491492
    492493    FillParticles(candidate, &entry->Particles);
Note: See TracChangeset for help on using the changeset viewer.