Fork me on GitHub

Changeset 2f9fb43 in git


Ignore:
Timestamp:
Dec 16, 2015, 10:49:11 AM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
7538969
Parents:
fdfe022
Message:

validated card for upgrade PhaseII studies

File:
1 edited

Legend:

Unmodified
Added
Removed
  • cards/CMS_PhaseII_140PU_conf4.tcl

    rfdfe022 r2f9fb43  
    1 #######################################
     1######################################
    22# Order of execution of various modules
    33#######################################
     
    2121  Calorimeter
    2222 
     23  ElectronFilter
     24 
    2325  TrackPileUpSubtractor
    2426   
     
    2628  EFlowMergerAllTracks
    2729  EFlowMerger
    28   RunPUPPI
    2930 
    3031  NeutrinoFilter
    3132 
    32   RhoGridFastJet
    33   RhoPuppiGridFastJet
    34 
     33  Rho
     34 
    3535  GenJetFinder
    3636  RecoJetFinder
    37   PuppiJetFinder
     37 
     38  JetPileUpSubtractor
     39 
     40  JetFlavorAssociation
    3841   
    39   RecoJetPileUpSubtractor
    40   PuppiJetPileUpSubtractor
    41    
    42   RecoJetFlavorAssociation
    43   PuppiJetFlavorAssociation
    44    
    45   RecoJetPileUpID
    46   PuppiJetPileUpID
    47  
    4842  RecoJetBTaggingLoose
    4943  RecoJetBTaggingMedium
    5044  RecoJetBTaggingTight
    5145
    52   PuppiJetBTaggingLoose
    53   PuppiJetBTaggingMedium
    54   PuppiJetBTaggingTight
    55 
    5646  GenMissingET
    5747  RecoMissingET
    58   JetPileUpIDMissingET
    59   PuppiMissingET
    60  
     48 
    6149  GenScalarHT
    6250  RecoScalarHT
    63   JetPileUpIDScalarHT
    64   PuppiScalarHT
    65    
     51   
    6652  PhotonEfficiency
    6753  PhotonIsolation
     
    9480
    9581  # pre-generated minbias input file
    96   set PileUpFile MinBias.pileup
     82  set PileUpFile /home/fynu/mvidal/scratch/Delphes_CMS/MinBias100K_14TeV.pileup
    9783
    9884  # Get rid of beam spot from http://red-gridftp11.unl.edu/Snowmass/MinBias100K_14TeV.pileup ...
     
    323309  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
    324310
    325   set ECalEnergyMin 0.5
    326   set HCalEnergyMin 0.5
    327 
    328   set ECalEnergySignificanceMin 2.0
    329   set HCalEnergySignificanceMin 2.0
     311  set ECalEnergyMin 0.0
     312  set HCalEnergyMin 0.0
     313
     314  set ECalEnergySignificanceMin 0.0
     315  set HCalEnergySignificanceMin 0.0
    330316
    331317  set SmearTowerCenter false
     
    402388}
    403389
     390
     391#################
     392# Electron filter
     393#################
     394
     395module PdgCodeFilter ElectronFilter {
     396  set InputArray Calorimeter/eflowTracks
     397  set OutputArray electrons
     398  set Invert true
     399  add PdgCode {11}
     400  add PdgCode {-11}
     401}
     402
     403
    404404##########################
    405405# Track pile-up subtractor
     
    409409# add InputArray InputArray OutputArray
    410410  add InputArray Calorimeter/eflowTracks eflowTracks
    411   add InputArray ElectronEnergySmearing/electrons electrons
     411  add InputArray ElectronFilter/electrons electrons
    412412  add InputArray MuonMomentumSmearing/muons muons
    413413
     
    511511}
    512512
     513#############
     514# Rho pile-up
     515#############
     516
     517module FastJetFinder Rho {
     518#  set InputArray Calorimeter/towers
     519  set InputArray EFlowMerger/eflow
     520
     521  set ComputeRho true
     522  set RhoOutputArray rho
     523
     524  # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
     525  set AreaAlgorithm 5
     526
     527  # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     528  set JetAlgorithm 4
     529  set ParameterR 0.4
     530  set GhostEtaMax 5.0
     531  set RhoEtaMax 5.0
     532
     533  add RhoEtaRange -5.0 -4.0
     534  add RhoEtaRange -4.0 -2.5
     535  add RhoEtaRange -2.5 2.5
     536  add RhoEtaRange 2.5 4.0
     537  add RhoEtaRange 4.0 5.0
     538
     539  set JetPTMin 0.0
     540}
     541
    513542
    514543#################################
     
    576605
    577606
    578 ################
    579 ## Jet finder ##
    580 ################
     607
     608############
     609# Jet finder
     610############
    581611
    582612module FastJetFinder RecoJetFinder {
    583   set InputArray   EFlowMerger/eflow
    584   set OutputArray  jets
     613#  set InputArray Calorimeter/towers
     614  set InputArray EFlowMerger/eflow
     615
     616  set OutputArray jets
     617
    585618  # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
    586   set AreaAlgorithm 1
    587   # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    588   set JetAlgorithm  6
    589   set ParameterR    0.4
    590   set JetPTMin      10.0
    591 }
    592 
    593 
    594 
    595 #####################
    596 ## Make puppi jets ##
    597 #####################
    598 
    599 module FastJetFinder PuppiJetFinder {
    600   set InputArray RunPUPPI/PuppiParticles
    601   set OutputArray jets
    602   # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
    603   set AreaAlgorithm 1
     619  set AreaAlgorithm 5
     620
    604621  # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
    605622  set JetAlgorithm 6
    606   set ParameterR   0.4
    607   set JetPTMin     10.0
    608 }
    609 
     623  set ParameterR 0.4
     624
     625  set JetPTMin 5.0
     626}
    610627
    611628
     
    614631###########################
    615632
    616 module JetPileUpSubtractor RecoJetPileUpSubtractor {
     633module JetPileUpSubtractor JetPileUpSubtractor {
    617634  set JetInputArray RecoJetFinder/jets
    618   set RhoInputArray RhoGridFastJet/rho
     635  set RhoInputArray Rho/rho
    619636
    620637  set OutputArray jets
    621638
    622   set JetPTMin 20.0
    623 }
    624 
    625 
    626 ###########################
    627 # Jet Pile-Up Subtraction
    628 ###########################
    629 
    630 module JetPileUpSubtractor PuppiJetPileUpSubtractor {
    631   set JetInputArray PuppiJetFinder/jets
    632   set RhoInputArray RhoPuppiGridFastJet/rho
    633 
    634   set OutputArray jets
    635 
    636   set JetPTMin 20.0
    637 }
    638 
     639  set JetPTMin 5.0
     640}
    639641
    640642
     
    643645########################
    644646
    645 module JetFlavorAssociation RecoJetFlavorAssociation {
     647module JetFlavorAssociation JetFlavorAssociation {
    646648 
    647649  set PartonInputArray Delphes/partons
    648650  set ParticleInputArray Delphes/allParticles
    649651  set ParticleLHEFInputArray Delphes/allParticlesLHEF
    650   set JetInputArray RecoJetPileUpSubtractor/jets
     652  set JetInputArray JetPileUpSubtractor/jets
    651653 
    652654  set DeltaR 0.5
     
    657659
    658660
    659 
    660 ########################
    661 # Jet Flavor Association
    662 ########################
    663 
    664 module JetFlavorAssociation PuppiJetFlavorAssociation {
    665  
    666   set PartonInputArray Delphes/partons
    667   set ParticleInputArray Delphes/allParticles
    668   set ParticleLHEFInputArray Delphes/allParticlesLHEF
    669   set JetInputArray PuppiJetPileUpSubtractor/jets
    670  
    671   set DeltaR 0.5
    672   set PartonPTMin 1.0
    673   set PartonEtaMax 2.5
    674 
    675 }
    676 
    677 
    678 ###########################
    679 # Jet Pile-Up ID
    680 ###########################
    681 
    682 module PileUpJetID RecoJetPileUpID {
    683  
    684   set JetInputArray RecoJetPileUpSubtractor/jets
    685  
    686   # Using constituents does not make sense with Charged hadron subtraction                                                                                                           
    687   # In 0 mode, dR cut used instead                                                                                                                                                   
    688   set UseConstituents 0
    689 
    690   set TrackInputArray TrackMerger/tracks
    691   set NeutralInputArray NeutralTowerMerger/eflowTowers
    692  
    693   set OutputArray jets
    694   set NeutralsInPassingJets eflowTowers
    695  
    696   set ParameterR 0.4
    697 
    698   set MeanSqDeltaRMaxBarrel 0.07
    699   set BetaMinBarrel 0.13
    700   set MeanSqDeltaRMaxEndcap 0.07
    701   set BetaMinEndcap 0.15
    702   set MeanSqDeltaRMaxForward 0.01
    703  
    704   set JetPTMin 10.0
    705 
    706 }
    707 
    708 ###########################
    709 # Jet Pile-Up ID
    710 ###########################
    711 
    712 module PileUpJetID PuppiJetPileUpID {
    713  
    714   set JetInputArray PuppiJetPileUpSubtractor/jets
    715  
    716   # Using constituents does not make sense with Charged hadron subtraction                                                                                                           
    717   # In 0 mode, dR cut used instead                                                                                                                                                   
    718   set UseConstituents 0
    719 
    720   set TrackInputArray RunPUPPI/puppiTracks
    721   set NeutralInputArray RunPUPPI/puppiNeutrals
    722  
    723   set OutputArray jets
    724   set NeutralsInPassingJets eflowTowers
    725  
    726   set ParameterR 0.4
    727 
    728   set MeanSqDeltaRMaxBarrel 0.07
    729   set BetaMinBarrel 0.13
    730   set MeanSqDeltaRMaxEndcap 0.07
    731   set BetaMinEndcap 0.15
    732   set MeanSqDeltaRMaxForward 0.01
    733  
    734   set JetPTMin 10.0
    735 
    736 }
    737661#############
    738662# b-tagging #
     
    740664module BTagging RecoJetBTaggingLoose {
    741665
    742   set JetInputArray RecoJetPileUpSubtractor/jets
     666  set JetInputArray JetPileUpSubtractor/jets
    743667
    744668  set BitNumber 0
     
    920844module BTagging RecoJetBTaggingMedium {
    921845
    922   set JetInputArray RecoJetPileUpSubtractor/jets
     846  set JetInputArray JetPileUpSubtractor/jets
    923847
    924848  set BitNumber 1
     
    11001024module BTagging RecoJetBTaggingTight {
    11011025
    1102   set JetInputArray RecoJetPileUpSubtractor/jets
     1026  set JetInputArray JetPileUpSubtractor/jets
    11031027
    11041028  set BitNumber 2
     
    12811205
    12821206
    1283 #############
    1284 # b-tagging #
    1285 #############
    1286 module BTagging PuppiJetBTaggingLoose {
    1287 
    1288   set JetInputArray PuppiJetPileUpSubtractor/jets
    1289 
    1290   set BitNumber 0
    1291 
    1292   add EfficiencyFormula {0}      {(pt <= 20.0) * (0.000) + \
    1293                                   (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.096) + \
    1294                                   (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.112) + \
    1295                                   (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0812) + \
    1296                                   (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.0908) + \
    1297                                   (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0792) + \
    1298                                   (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.0858) + \
    1299                                   (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.0917) + \
    1300                                   (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.084) + \
    1301                                   (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.0906) + \
    1302                                   (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.0989) + \
    1303                                   (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1022) + \
    1304                                   (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1035) + \
    1305                                   (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1095) + \
    1306                                   (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1201) + \
    1307                                   (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1348) + \
    1308                                   (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1482) + \
    1309                                   (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1629) + \
    1310                                   (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.1775) + \
    1311                                   (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.2002) + \
    1312                                   (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.1977) + \
    1313                                   (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2084) + \
    1314                                   (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2195) + \
    1315                                   (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2424) + \
    1316                                   (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2909) + \
    1317                                   (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3457) + \
    1318                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1319                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.074765) + \
    1320                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.100053) + \
    1321                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.071492) + \
    1322                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.084796) + \
    1323                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.076927) + \
    1324                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.08424) + \
    1325                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.093118) + \
    1326                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.084629) + \
    1327                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.092977) + \
    1328                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.10206) + \
    1329                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.102344) + \
    1330                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.098435) + \
    1331                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.105507) + \
    1332                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.112841) + \
    1333                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.126329) + \
    1334                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.140759) + \
    1335                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.153193) + \
    1336                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.107869) + \
    1337                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.119527) + \
    1338                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.08688) + \
    1339                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.089324) + \
    1340                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.097172) + \
    1341                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.000) + \
    1342                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.000) + \
    1343                                   (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
    1344                                   (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.02) +\
    1345                                   (abs(eta) > 4.0) * (0.000) \
    1346   }
    1347 
    1348   add EfficiencyFormula {4}          { (pt <= 20.0) * (0.000) + \
    1349                                        (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.387) + \
    1350                                        (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.448) + \
    1351                                        (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.408) + \
    1352                                        (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.427) + \
    1353                                        (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.408) + \
    1354                                        (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.425) + \
    1355                                        (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.426) + \
    1356                                        (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.4) + \
    1357                                        (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.415) + \
    1358                                        (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.416) + \
    1359                                        (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.405) + \
    1360                                        (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.387) + \
    1361                                        (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.39) + \
    1362                                        (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.389) + \
    1363                                        (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.389) + \
    1364                                        (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.381) + \
    1365                                        (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.381) + \
    1366                                        (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.367) + \
    1367                                        (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.369) + \
    1368                                        (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.326) + \
    1369                                        (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.335) + \
    1370                                        (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.326) + \
    1371                                        (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.341) + \
    1372                                        (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.403) + \
    1373                                        (abs(eta) <= 1.8) * (pt > 2000.0) * (0.47) + \
    1374                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1375                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.2497) + \
    1376                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.31891) + \
    1377                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.273) + \
    1378                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.28445) + \
    1379                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.28036) + \
    1380                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.28453) + \
    1381                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.29495) + \
    1382                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.27024) + \
    1383                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.28912) + \
    1384                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.29048) + \
    1385                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.27507) + \
    1386                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.23327) + \
    1387                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.2493) + \
    1388                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.2416) + \
    1389                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.26652) + \
    1390                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.24852) + \
    1391                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.26927) + \
    1392                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.19302) + \
    1393                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.19433) + \
    1394                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.17523) + \
    1395                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.14981) + \
    1396                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.16666) + \
    1397                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1398                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1399                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
    1400                                        (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.29*tanh(pt*0.00997 - 0.143)) + \
    1401                                        (abs(eta) > 4.0) * (0.000)
    1402   }
    1403 
    1404    # efficiency formula for b-jets
    1405   add EfficiencyFormula {5}          { (pt <= 20.0) * (0.000) + \
    1406                                        (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.75) + \
    1407                                        (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.827) + \
    1408                                        (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
    1409                                        (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.853) + \
    1410                                        (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.855) + \
    1411                                        (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.862) + \
    1412                                        (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.868) + \
    1413                                        (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.865) + \
    1414                                        (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.863) + \
    1415                                        (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.857) + \
    1416                                        (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.851) + \
    1417                                        (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.838) + \
    1418                                        (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.831) + \
    1419                                        (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.817) + \
    1420                                        (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.796) + \
    1421                                        (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.772) + \
    1422                                        (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.76) + \
    1423                                        (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.743) + \
    1424                                        (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.703) + \
    1425                                        (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.638) + \
    1426                                        (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.605) + \
    1427                                        (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.572) + \
    1428                                        (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.541) + \
    1429                                        (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.567) + \
    1430                                        (abs(eta) <= 1.8) * (pt > 2000.0) * (0.603) + \
    1431                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1432                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.6063) + \
    1433                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.7188) + \
    1434                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.72) + \
    1435                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.7365) + \
    1436                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.7462) + \
    1437                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.7454) + \
    1438                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.7415) + \
    1439                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.727) + \
    1440                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.7112) + \
    1441                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.7112) + \
    1442                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.6754) + \
    1443                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.6359) + \
    1444                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.6348) + \
    1445                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.6115) + \
    1446                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.5585) + \
    1447                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.5608) + \
    1448                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.5208) + \
    1449                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.456) + \
    1450                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.4524) + \
    1451                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.388) + \
    1452                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.3928) + \
    1453                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.3823) + \
    1454                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1455                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1456                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
    1457                                        (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
    1458                                        (abs(eta) > 4.0) * (0.000)
    1459   }
    1460 }
    1461 
    1462 
    1463 #############
    1464 # b-tagging #
    1465 #############
    1466 module BTagging PuppiJetBTaggingMedium {
    1467 
    1468   set JetInputArray PuppiJetPileUpSubtractor/jets
    1469 
    1470   set BitNumber 1
    1471  
    1472   add EfficiencyFormula {0}       { (pt <= 20.0) * (0.000) + \
    1473                                     (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00654) + \
    1474                                     (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00921) + \
    1475                                         (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.00573) + \
    1476                                         (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00694) + \
    1477                                         (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0062) + \
    1478                                         (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00708) + \
    1479                                         (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00779) + \
    1480                                         (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00693) + \
    1481                                         (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00777) + \
    1482                                         (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00862) + \
    1483                                         (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01038) + \
    1484                                         (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01189) + \
    1485                                         (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01279) + \
    1486                                         (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.01452) + \
    1487                                         (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01696) + \
    1488                                         (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01958) + \
    1489                                         (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.02253) + \
    1490                                         (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.01787) + \
    1491                                         (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.02154) + \
    1492                                         (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.01839) + \
    1493                                         (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.01987) + \
    1494                                         (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.02351) + \
    1495                                         (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.02937) + \
    1496                                         (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.04001) + \
    1497                                         (abs(eta) <= 1.8) * (pt > 2000.0) * (0.0542) + \
    1498                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1499                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.004236) + \
    1500                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.006653) + \
    1501                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.005512) + \
    1502                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.00754) + \
    1503                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.005813) + \
    1504                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.006439) + \
    1505                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.008063) + \
    1506                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.00647) + \
    1507                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.007583) + \
    1508                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.008543) + \
    1509                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.01034) + \
    1510                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.011253) + \
    1511                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.012945) + \
    1512                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.014474) + \
    1513                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.017361) + \
    1514                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.020912) + \
    1515                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.023139) + \
    1516                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.010756) + \
    1517                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.012569) + \
    1518                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
    1519                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.006428) + \
    1520                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00887) + \
    1521                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1522                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1523                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
    1524                                         (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
    1525                                         (abs(eta) > 4.0) * (0.000) \
    1526 }
    1527 
    1528  # efficiency formula for c-jets (misidentification rate)
    1529   add EfficiencyFormula {4}           { (pt <= 20.0) * (0.000) + \
    1530                                         (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.1102) + \
    1531                                         (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.1344) + \
    1532                                         (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.1025) + \
    1533                                         (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.1025) + \
    1534                                         (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.1063) + \
    1535                                         (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.1087) + \
    1536                                         (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.1124) + \
    1537                                         (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.109) + \
    1538                                         (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.111) + \
    1539                                         (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.1091) + \
    1540                                         (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1087) + \
    1541                                         (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1091) + \
    1542                                         (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1107) + \
    1543                                         (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1061) + \
    1544                                         (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1017) + \
    1545                                         (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1) + \
    1546                                         (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.0966) + \
    1547                                         (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0697) + \
    1548                                         (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0679) + \
    1549                                         (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0503) + \
    1550                                         (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0514) + \
    1551                                         (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0481) + \
    1552                                         (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0667) + \
    1553                                         (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0861) + \
    1554                                         (abs(eta) <= 1.8) * (pt > 2000.0) * (0.092) + \
    1555                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1556                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.03331) + \
    1557                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.04361) + \
    1558                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.03863) + \
    1559                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.04287) + \
    1560                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.04431) + \
    1561                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.04452) + \
    1562                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.04339) + \
    1563                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0436) + \
    1564                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0456) + \
    1565                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.05138) + \
    1566                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.04794) + \
    1567                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.04004) + \
    1568                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.04713) + \
    1569                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.04515) + \
    1570                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.05314) + \
    1571                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.05143) + \
    1572                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.05936) + \
    1573                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.02357) + \
    1574                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.03222) + \
    1575                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
    1576                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.02621) + \
    1577                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.01709) + \
    1578                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1579                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1580                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
    1581                                         (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.1898*tanh(pt*0.00997 - 0.143)) + \
    1582                                         (abs(eta) > 4.0) * (0.000)
    1583   }
    1584 
    1585  # efficiency formula for b-jets
    1586   add EfficiencyFormula {5}           { (pt <= 20.0) * (0.000) + \
    1587                                         (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.536) + \
    1588                                         (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.6439) + \
    1589                                         (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.6504) + \
    1590                                         (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.6716) + \
    1591                                         (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.6841) + \
    1592                                         (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.6896) + \
    1593                                         (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.6916) + \
    1594                                         (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.6882) + \
    1595                                         (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.6838) + \
    1596                                         (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.6715) + \
    1597                                         (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.6554) + \
    1598                                         (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.6366) + \
    1599                                         (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.6192) + \
    1600                                         (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.595) + \
    1601                                         (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.5551) + \
    1602                                         (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.5138) + \
    1603                                         (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.4884) + \
    1604                                         (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.4009) + \
    1605                                         (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.3459) + \
    1606                                         (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.2523) + \
    1607                                         (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2404) + \
    1608                                         (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2198) + \
    1609                                         (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2263) + \
    1610                                         (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2614) + \
    1611                                         (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3194) + \
    1612                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1613                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.3254) + \
    1614                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.4339) + \
    1615                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.4499) + \
    1616                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.4716) + \
    1617                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.4766) + \
    1618                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.4788) + \
    1619                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.4863) + \
    1620                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.4891) + \
    1621                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.462) + \
    1622                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.4583) + \
    1623                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.4247) + \
    1624                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.3775) + \
    1625                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.3734) + \
    1626                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.3348) + \
    1627                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.2939) + \
    1628                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.285) + \
    1629                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.2421) + \
    1630                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.1565) + \
    1631                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.1522) + \
    1632                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.1231) + \
    1633                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.1607) + \
    1634                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.1323) + \
    1635                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1636                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1637                                         (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
    1638                                             (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
    1639                                         (abs(eta) > 4.0) * (0.000)
    1640  }
    1641 }
    1642 
    1643 #############
    1644 # b-tagging #
    1645 #############
    1646 module BTagging PuppiJetBTaggingTight {
    1647 
    1648   set JetInputArray PuppiJetPileUpSubtractor/jets
    1649 
    1650   set BitNumber 2
    1651 
    1652   add EfficiencyFormula {0}          { (pt <= 20.0) * (0.000) + \
    1653                                        (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.000164) + \
    1654                                        (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.000235) + \
    1655                                        (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.000266) + \
    1656                                        (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.000329) + \
    1657                                        (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.000309) + \
    1658                                        (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.000309) + \
    1659                                        (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.000546) + \
    1660                                        (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.000499) + \
    1661                                        (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.000642) + \
    1662                                        (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.000742) + \
    1663                                        (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.000928) + \
    1664                                        (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.001323) + \
    1665                                        (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.001392) + \
    1666                                        (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00154) + \
    1667                                        (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.002094) + \
    1668                                        (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.002427) + \
    1669                                        (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.002927) + \
    1670                                        (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.001854) + \
    1671                                        (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002355) + \
    1672                                        (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002297) + \
    1673                                        (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.002433) + \
    1674                                        (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002706) + \
    1675                                        (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.003602) + \
    1676                                        (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.004987) + \
    1677                                        (abs(eta) <= 1.8) * (pt > 2000.0) * (0.007414) + \
    1678                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1679                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.000573) + \
    1680                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.000574) + \
    1681                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.000521) + \
    1682                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.000786) + \
    1683                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.000539) + \
    1684                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.000673) + \
    1685                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.000934) + \
    1686                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.000781) + \
    1687                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.000949) + \
    1688                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.000977) + \
    1689                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.001168) + \
    1690                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.000879) + \
    1691                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.000812) + \
    1692                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.001215) + \
    1693                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.001679) + \
    1694                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.001893) + \
    1695                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.002723) + \
    1696                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.003555) + \
    1697                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.003881) + \
    1698                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
    1699                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.005563) + \
    1700                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.007611) + \
    1701                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1702                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1703                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
    1704                                        (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
    1705                                        (abs(eta) > 4.0) * (0.000) \
    1706    }
    1707 
    1708  # efficiency formula for c-jets (misidentification rate)
    1709   add EfficiencyFormula {4}     { (pt <= 20.0) * (0.000) + \
    1710                                   (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00531) + \
    1711                                        (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00567) + \
    1712                                        (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0064) + \
    1713                                        (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00673) + \
    1714                                        (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.00766) + \
    1715                                        (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00729) + \
    1716                                        (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00674) + \
    1717                                        (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00824) + \
    1718                                        (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00888) + \
    1719                                        (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00919) + \
    1720                                        (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01021) + \
    1721                                        (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01041) + \
    1722                                        (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01027) + \
    1723                                        (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00999) + \
    1724                                        (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01047) + \
    1725                                        (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01014) + \
    1726                                        (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.01021) + \
    1727                                        (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.00601) + \
    1728                                        (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0054) + \
    1729                                        (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.00487) + \
    1730                                        (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.00519) + \
    1731                                        (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.00469) + \
    1732                                        (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.00651) + \
    1733                                        (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.01299) + \
    1734                                        (abs(eta) <= 1.8) * (pt > 2000.0) * (0.00897) + \
    1735                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1736                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.01014) + \
    1737                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.01288) + \
    1738                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.01392) + \
    1739                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.01533) + \
    1740                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.01508) + \
    1741                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.01579) + \
    1742                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.01106) + \
    1743                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.01346) + \
    1744                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.01315) + \
    1745                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.01156) + \
    1746                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0082) + \
    1747                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.00439) + \
    1748                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.00744) + \
    1749                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.00685) + \
    1750                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.00755) + \
    1751                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.00706) + \
    1752                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.00428) + \
    1753                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.01031) + \
    1754                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.00976) + \
    1755                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
    1756                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.00374) + \
    1757                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00854) + \
    1758                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1759                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1760                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
    1761                                        (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
    1762                                        (abs(eta) > 4.0) * (0.000)
    1763 }
    1764 
    1765   # efficiency formula for b-jets
    1766   add EfficiencyFormula {5}          { (pt <= 20.0) * (0.000) + \
    1767                                        (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.2426) + \
    1768                                        (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.327) + \
    1769                                        (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.3559) + \
    1770                                        (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.3704) + \
    1771                                        (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.3824) + \
    1772                                        (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.3844) + \
    1773                                        (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.3848) + \
    1774                                        (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.3862) + \
    1775                                        (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.3778) + \
    1776                                        (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.3622) + \
    1777                                        (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.3299) + \
    1778                                        (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.2889) + \
    1779                                        (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.2815) + \
    1780                                        (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.253) + \
    1781                                        (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.221) + \
    1782                                        (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1963) + \
    1783                                        (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1739) + \
    1784                                        (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0992) + \
    1785                                        (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0788) + \
    1786                                        (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0581) + \
    1787                                        (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0534) + \
    1788                                        (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0521) + \
    1789                                        (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0626) + \
    1790                                        (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0826) + \
    1791                                        (abs(eta) <= 1.8) * (pt > 2000.0) * (0.1022) + \
    1792                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
    1793                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.1562) + \
    1794                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.2499) + \
    1795                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.2956) + \
    1796                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.3128) + \
    1797                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.3147) + \
    1798                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.3222) + \
    1799                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.304) + \
    1800                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.3051) + \
    1801                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.2657) + \
    1802                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.2578) + \
    1803                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.2087) + \
    1804                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.1634) + \
    1805                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.1651) + \
    1806                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.1353) + \
    1807                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.109) + \
    1808                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0799) + \
    1809                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0699) + \
    1810                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.054) + \
    1811                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0718) + \
    1812                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.0746) + \
    1813                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.0803) + \
    1814                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.0882) + \
    1815                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
    1816                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
    1817                                        (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
    1818                                        (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
    1819                                        (abs(eta) > 4.0) * (0.000)
    1820   }
    1821 
    1822 }
    1823 
    1824 
    18251207
    18261208#####################
     
    18331215  set ParticleInputArray Delphes/allParticles
    18341216  set PartonInputArray Delphes/partons
    1835   set JetInputArray RecoJetPileUpSubtractor/jets
     1217  set JetInputArray JetPileUpSubtractor/jets
    18361218 
    18371219  set DeltaR 0.4
     
    18661248}
    18671249
    1868 module Merger JetPileUpIDMissingET {
    1869   add InputArray TrackPileUpSubtractor/eflowTracks
    1870   add InputArray RecoJetPileUpID/eflowTowers
    1871   set MomentumOutputArray momentum
    1872 }
    1873 
    1874 module Merger PuppiMissingET {
    1875   add InputArray RunPUPPI/PuppiParticles
    1876   set MomentumOutputArray momentum
    1877 }
    18781250
    18791251
     
    18931265 set EnergyOutputArray energy
    18941266}
    1895 
    1896 module Merger JetPileUpIDScalarHT {
    1897   add InputArray TrackPileUpSubtractor/eflowTracks
    1898   add InputArray RecoJetPileUpID/eflowTowers
    1899   set MomentumOutputArray energy
    1900 }
    1901 
    1902 module Merger PuppiScalarHT {
    1903  # add InputArray InputArray
    1904  add InputArray RunPUPPI/PuppiParticles
    1905  set EnergyOutputArray energy
    1906 }
    1907 
    19081267
    19091268#####################
     
    19391298 
    19401299  # select a rho for the isolation
    1941   set RhoInputArray RhoGridFastJet/rho
     1300  set RhoInputArray Rho/rho
    19421301 
    19431302  # output array
     
    19481307 
    19491308  # minimum pT 
    1950   set PTMin     0.5
     1309  set PTMin     1.0
    19511310 
    19521311  # iso ratio to cut
     
    20101369  set IsolationInputArray EFlowMergerAllTracks/eflow
    20111370 
    2012   set RhoInputArray RhoGridFastJet/rho
     1371  set RhoInputArray Rho/rho
    20131372 
    20141373  set OutputArray electrons
    20151374 
    20161375  set DeltaRMax 0.3
    2017   set PTMin 0.5
     1376  set PTMin 1.0
    20181377  set PTRatioMax 9999.
    20191378
    20201379}
     1380
     1381
    20211382
    20221383###################
     
    20381399                          (abs(eta) <= 4.00) * (pt >  50. && pt <= 70.)  * (0.98) + \
    20391400                          (abs(eta) <= 4.00) * (pt > 70.0 )  * (1.00) + \
    2040                          (abs(eta) > 4.00)  * (0.00)
     1401                          (abs(eta) > 4.00)  * (0.00)
    20411402 }
    20421403
    20431404}
     1405
     1406
    20441407
    20451408##################
     
    20531416  set IsolationInputArray EFlowMergerAllTracks/eflow
    20541417 
    2055   set RhoInputArray RhoGridFastJet/rho
     1418  set RhoInputArray Rho/rho
    20561419 
    20571420  set OutputArray muons
    20581421 
    20591422  set DeltaRMax 0.3
    2060   set PTMin 0.5
     1423  set PTMin 1.0
    20611424  set PTRatioMax 9999.
    20621425
     
    20991462  #add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
    21001463 
    2101   add Branch RhoGridFastJet/rho RhoGridFastJet Rho
    2102   add Branch RhoPuppiGridFastJet/rho RhoPuppiGridFastJet Rho
    2103  
     1464  add Branch Rho/rho Rho Rho
     1465   
    21041466  ## Jets 
    2105   #add Branch GenJetFinder/jets GenJet Jet
    2106   #add Branch RecoJetFinder/jets RecoJet Jet
    2107   #add Branch PuppiJetFinder/jets PuppiJet Jet
    2108  
    2109  
    2110   add Branch RecoJetPileUpID/jets RecoJetPileUpID Jet
    2111   add Branch PuppiJetPileUpID/jets PuppiJetPileUpID Jet
     1467  add Branch GenJetFinder/jets GenJet Jet
     1468  add Branch JetPileUpSubtractor/jets Jet Jet
    21121469 
    21131470  ## MET
    21141471  add Branch GenMissingET/momentum GenMissingET MissingET
    2115   add Branch RecoMissingET/momentum RecoMissingET MissingET
    2116   add Branch JetPileUpIDMissingET/momentum JetPileUpIDMissingET MissingET
    2117   add Branch PuppiMissingET/momentum PuppiMissingET MissingET
    2118 
     1472  add Branch RecoMissingET/momentum MissingET MissingET
     1473 
    21191474  ## HT
    21201475  add Branch GenScalarHT/energy GenHT ScalarHT
    21211476  add Branch RecoScalarHT/energy RecoHT ScalarHT
    2122   add Branch JetPileUpIDScalarHT/energy JetPileUpIDHT ScalarHT
    2123   add Branch PuppiScalarHT/energy PuppiHT ScalarHT
    2124 
     1477 
    21251478  ## photons and leptons
    21261479  add Branch ElectronIsolation/electrons Electron Electron
     
    21281481  add Branch MuonIsolation/muons Muon Muon
    21291482
    2130 }
     1483
     1484
     1485}
Note: See TracChangeset for help on using the changeset viewer.