###################################### # Order of execution of various modules ####################################### set ExecutionPath { PileUpMerger ParticlePropagator ChargedHadronTrackingEfficiency ElectronTrackingEfficiency MuonTrackingEfficiency ChargedHadronMomentumSmearing ElectronEnergySmearing MuonMomentumSmearing TrackMerger Calorimeter ElectronFilter TrackPileUpSubtractor NeutralTowerMerger EFlowMergerAllTracks EFlowMerger NeutrinoFilter Rho GenJetFinder RecoJetFinder JetPileUpSubtractor JetFlavorAssociation RecoJetBTaggingLoose RecoJetBTaggingMedium RecoJetBTaggingTight GenMissingET RecoMissingET GenScalarHT RecoScalarHT PhotonEfficiency PhotonIsolation ElectronEfficiency ElectronIsolation MuonEfficiency MuonIsolation TauTagging StatusPidFilter TreeWriter } ############### # PileUp Merger ############### module PileUpMerger PileUpMerger { set InputArray Delphes/stableParticles set ParticleOutputArray stableParticles set VertexOutputArray vertices # pre-generated minbias input file set PileUpFile /home/fynu/mvidal/scratch/Delphes_CMS/MinBias100K_14TeV.pileup # Get rid of beam spot from http://red-gridftp11.unl.edu/Snowmass/MinBias100K_14TeV.pileup ... set InputBSX 2.44 set InputBSY 3.39 # replace it with beam spot from CMSSW files set OutputBSX 0.24 set OutputBSY 0.39 # average expected pile up set MeanPileUp 140 # maximum spread in the beam direction in m set ZVertexSpread 0.15 # maximum spread in time in s set TVertexSpread 1.5E-09 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))} set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)} } ##################################################################################### # Propagate particles in cylinder and divide charged particles in different classes # ##################################################################################### module ParticlePropagator ParticlePropagator { ## take particles after beam spot smearing set InputArray PileUpMerger/stableParticles ## produce independent output collection: all particles, only charged hadrons, electrons and muons set OutputArray stableParticles set ChargedHadronOutputArray chargedHadrons set ElectronOutputArray electrons set MuonOutputArray muons ## radius of the magnetic field coverage, in m set Radius 1.29 ## half-length of the magnetic field coverage, in m set HalfLength 3.00 ## magnetic field set Bz 3.8 } #################################### # Charged hadron tracking efficiency #################################### module Efficiency ChargedHadronTrackingEfficiency { ## particles after propagation set InputArray ParticlePropagator/chargedHadrons set OutputArray chargedHadrons # tracking efficiency formula for charged hadrons set EfficiencyFormula { (pt <= 0.2) * (0.00) + \ (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \ (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \ (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \ (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \ (abs(eta) > 4.0) * (0.00) } } ##################################### # Electron tracking efficiency - ID #################################### module Efficiency ElectronTrackingEfficiency { set InputArray ParticlePropagator/electrons set OutputArray electrons # tracking efficiency formula for electrons set EfficiencyFormula { (pt <= 0.2) * (0.00) + \ (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \ (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \ (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \ (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \ (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \ (abs(eta) > 4.0) * (0.00) } } ########################## # Muon tracking efficiency ########################## module Efficiency MuonTrackingEfficiency { set InputArray ParticlePropagator/muons set OutputArray muons # tracking efficiency formula for muons set EfficiencyFormula { (pt <= 0.2) * (0.00) + \ (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \ (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \ (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \ (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \ (abs(eta) > 4.0) * (0.00) } } ######################################## # Momentum resolution for charged tracks ######################################## module MomentumSmearing ChargedHadronMomentumSmearing { ## hadrons after having applied the tracking efficiency set InputArray ChargedHadronTrackingEfficiency/chargedHadrons set OutputArray chargedHadrons # resolution formula for charged hadrons set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \ (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \ (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \ (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80) } } ################################# # Energy resolution for electrons ################################# module EnergySmearing ElectronEnergySmearing { set InputArray ElectronTrackingEfficiency/electrons set OutputArray electrons # set ResolutionFormula {resolution formula as a function of eta and energy} set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.02) + \ (abs(eta) <= 2.5) * (energy > 2.5e1) * (energy*0.028) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (energy*0.032) + \ (abs(eta) > 3.0 && abs(eta) <= 5.0) * (energy*0.07) } } ############################### # Momentum resolution for muons ############################### module MomentumSmearing MuonMomentumSmearing { set InputArray MuonTrackingEfficiency/muons set OutputArray muons # resolution formula for muons set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.012) + \ (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \ (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.012) + \ (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.030) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.011) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.011) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.022) + \ (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.030) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \ (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \ (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \ (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80) } } ############## # Track merger ############## module Merger TrackMerger { # add InputArray InputArray add InputArray ChargedHadronMomentumSmearing/chargedHadrons add InputArray ElectronEnergySmearing/electrons add InputArray MuonMomentumSmearing/muons set OutputArray tracks } ############################################################## # Calorimeter : emulate calorimiter answer making caloTowers # ############################################################## module Calorimeter Calorimeter { ## particle from the propagation without any efficiency or smearing (for neutrals) set ParticleInputArray ParticlePropagator/stableParticles ## track after smearing and efficiency: used for charged particles set TrackInputArray TrackMerger/tracks ## output collections set TowerOutputArray towers set PhotonOutputArray photons set EFlowTrackOutputArray eflowTracks set EFlowPhotonOutputArray eflowPhotons set EFlowNeutralHadronOutputArray eflowNeutralHadrons set ECalEnergyMin 0.0 set HCalEnergyMin 0.0 set ECalEnergySignificanceMin 0.0 set HCalEnergySignificanceMin 0.0 set SmearTowerCenter false set pi [expr {acos(-1)}] # lists of the edges of each tower in eta and phi # each list starts with the lower edge of the first tower # the list ends with the higher edged of the last tower # 5 degrees towers set PhiBins {} for {set i -70} {$i <= 70} {incr i} { add PhiBins [expr {$i * $pi/70.0}] } foreach eta {-4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 -1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4} { add EtaPhiBins $eta $PhiBins } # 10 degrees towers set PhiBins {} for {set i -11} {$i <= 11} {incr i} { add PhiBins [expr {$i * $pi/11.0}] } foreach eta {-4.35 -4.175 -4 4.175 4.35 4.525} { add EtaPhiBins $eta $PhiBins } # 20 degrees towers set PhiBins {} for {set i -9} {$i <= 9} {incr i} { add PhiBins [expr {$i * $pi/9.0}] } foreach eta {-5 -4.7 -4.525 4.7 5} { add EtaPhiBins $eta $PhiBins } # default energy fractions {abs(PDG code)} {Fecal Fhcal} add EnergyFraction {0} {0.0 1.0} # energy fractions for e, gamma and pi0 add EnergyFraction {11} {1.0 0.0} add EnergyFraction {22} {1.0 0.0} add EnergyFraction {111} {1.0 0.0} # energy fractions for muon, neutrinos and neutralinos add EnergyFraction {12} {0.0 0.0} add EnergyFraction {13} {0.0 0.0} add EnergyFraction {14} {0.0 0.0} add EnergyFraction {16} {0.0 0.0} add EnergyFraction {1000022} {0.0 0.0} add EnergyFraction {1000023} {0.0 0.0} add EnergyFraction {1000025} {0.0 0.0} add EnergyFraction {1000035} {0.0 0.0} add EnergyFraction {1000045} {0.0 0.0} # energy fractions for K0short and Lambda add EnergyFraction {310} {0.3 0.7} add EnergyFraction {3122} {0.3 0.7} # set ECalResolutionFormula {resolution formula as a function of eta and energy} set ECalResolutionFormula {(abs(eta) <= 1.5) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \ (abs(eta) > 1.5 && abs(eta) <= 4.0) * sqrt(energy^2*0.008^2 + energy*0.092^2 + 0.088^2) + \ (abs(eta) > 4.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.97^2)} # set HCalResolutionFormula {resolution formula as a function of eta and energy} set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \ (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.050^2 + energy*0.706^2) + \ (abs(eta) > 3.0 && abs(eta) <= 4.9) * sqrt(energy^2*0.05^2 + energy*1.00^2)} } ################# # Electron filter ################# module PdgCodeFilter ElectronFilter { set InputArray Calorimeter/eflowTracks set OutputArray electrons set Invert true add PdgCode {11} add PdgCode {-11} } ########################## # Track pile-up subtractor ########################## module TrackPileUpSubtractor TrackPileUpSubtractor { # add InputArray InputArray OutputArray add InputArray Calorimeter/eflowTracks eflowTracks add InputArray ElectronFilter/electrons electrons add InputArray MuonMomentumSmearing/muons muons set VertexInputArray PileUpMerger/vertices # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution # Z vertex resolution in m set ZVertexResolution 0.0001 } #################### # Neutral tower merger #################### module Merger NeutralTowerMerger { # add InputArray InputArray add InputArray Calorimeter/eflowPhotons add InputArray Calorimeter/eflowNeutralHadrons set OutputArray eflowTowers } ################################## # Energy flow merger (all tracks) ################################## module Merger EFlowMergerAllTracks { # add InputArray InputArray add InputArray TrackMerger/tracks add InputArray Calorimeter/eflowPhotons add InputArray Calorimeter/eflowNeutralHadrons set OutputArray eflow } #################### # Energy flow merger #################### module Merger EFlowMerger { # add InputArray InputArray add InputArray TrackPileUpSubtractor/eflowTracks add InputArray Calorimeter/eflowPhotons add InputArray Calorimeter/eflowNeutralHadrons set OutputArray eflow } ########################### ### Run the puppi code ### ########################### module RunPUPPI RunPUPPI { ## input information set TrackInputArray TrackMerger/tracks set NeutralInputArray NeutralTowerMerger/eflowTowers set PVInputArray PileUpMerger/vertices set MinPuppiWeight 0.05 set UseExp false ## define puppi algorithm parameters (more than one for the same eta region is possible) add EtaMinBin 0. 2.5 2.5 3.0 3.0 add EtaMaxBin 2.5 3.0 3.0 10.0 10.0 add PtMinBin 0. 0.5 0.5 0.5 0.5 add ConeSizeBin 0.25 0.25 0.25 0.25 0.25 add RMSPtMinBin 0.1 0.5 0.5 0.5 0.5 add RMSScaleFactorBin 1.0 1.0 1.0 1.0 1.0 add NeutralMinEBin 0.2 1.0 1.0 1.5 1.5 add NeutralPtSlope 0.02 0.02 0.02 0.02 0.02 add ApplyCHS true true true true true add UseCharged true false false false false add ApplyLowPUCorr true true true true true add MetricId 5 5 1 5 1 ## output name set OutputArray PuppiParticles set OutputArrayTracks puppiTracks set OutputArrayNeutrals puppiNeutrals } ##################### # Neutrino Filter ##################### module PdgCodeFilter NeutrinoFilter { set InputArray Delphes/stableParticles set OutputArray filteredParticles set PTMin 0.0 add PdgCode {12} add PdgCode {14} add PdgCode {16} add PdgCode {-12} add PdgCode {-14} add PdgCode {-16} } ############# # Rho pile-up ############# module FastJetFinder Rho { # set InputArray Calorimeter/towers set InputArray EFlowMerger/eflow set ComputeRho true set RhoOutputArray rho # 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 set AreaAlgorithm 5 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt set JetAlgorithm 4 set ParameterR 0.4 set GhostEtaMax 5.0 set RhoEtaMax 5.0 add RhoEtaRange -5.0 -4.0 add RhoEtaRange -4.0 -2.5 add RhoEtaRange -2.5 2.5 add RhoEtaRange 2.5 4.0 add RhoEtaRange 4.0 5.0 set JetPTMin 0.0 } ################################# # Rho pile-up with std particles ################################# module FastJetGridMedianEstimator RhoGridFastJet { set InputArray EFlowMerger/eflow set RhoOutputArray rho # add GridRange rapmin rapmax drap dphi # rapmin - the minimum rapidity extent of the grid # rapmax - the maximum rapidity extent of the grid # drap - the grid spacing in rapidity # dphi - the grid spacing in azimuth add GridRange -5.0 -4.0 1.0 1.0 add GridRange -4.0 -2.5 1.0 1.0 add GridRange -2.5 2.5 1.0 1.0 add GridRange 2.5 4.0 1.0 1.0 add GridRange 4.0 5.0 1.0 1.0 } ################################### # Rho pile-up with PUPPI particles ################################### module FastJetGridMedianEstimator RhoPuppiGridFastJet { set InputArray RunPUPPI/PuppiParticles set RhoOutputArray rho # add GridRange rapmin rapmax drap dphi # rapmin - the minimum rapidity extent of the grid # rapmax - the maximum rapidity extent of the grid # drap - the grid spacing in rapidity # dphi - the grid spacing in azimuth add GridRange -5.0 -4.0 1.0 1.0 add GridRange -4.0 -2.5 1.0 1.0 add GridRange -2.5 2.5 1.0 1.0 add GridRange 2.5 4.0 1.0 1.0 add GridRange 4.0 5.0 1.0 1.0 } ##################### # MC truth jet finder ##################### module FastJetFinder GenJetFinder { set InputArray NeutrinoFilter/filteredParticles set OutputArray jets # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt set JetAlgorithm 6 set ParameterR 0.4 set JetPTMin 15.0 } ############ # Jet finder ############ module FastJetFinder RecoJetFinder { # set InputArray Calorimeter/towers set InputArray EFlowMerger/eflow set OutputArray jets # 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 set AreaAlgorithm 5 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt set JetAlgorithm 6 set ParameterR 0.4 set JetPTMin 5.0 } ########################### # Jet Pile-Up Subtraction ########################### module JetPileUpSubtractor JetPileUpSubtractor { set JetInputArray RecoJetFinder/jets set RhoInputArray Rho/rho set OutputArray jets set JetPTMin 5.0 } ######################## # Jet Flavor Association ######################## module JetFlavorAssociation JetFlavorAssociation { set PartonInputArray Delphes/partons set ParticleInputArray Delphes/allParticles set ParticleLHEFInputArray Delphes/allParticlesLHEF set JetInputArray JetPileUpSubtractor/jets set DeltaR 0.5 set PartonPTMin 1.0 set PartonEtaMax 2.5 } ############# # b-tagging # ############# module BTagging RecoJetBTaggingLoose { set JetInputArray JetPileUpSubtractor/jets set BitNumber 0 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.096) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.112) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0812) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.0908) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0792) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.0858) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.0917) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.084) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.0906) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.0989) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1022) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1035) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1095) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1201) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1348) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1482) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1629) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.1775) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.2002) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.1977) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2084) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2195) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2424) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2909) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3457) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.074765) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.100053) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.071492) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.084796) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.076927) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.08424) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.093118) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.084629) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.092977) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.10206) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.102344) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.098435) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.105507) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.112841) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.126329) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.140759) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.153193) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.107869) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.119527) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.08688) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.089324) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.097172) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.02) +\ (abs(eta) > 4.0) * (0.000) \ } add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.387) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.448) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.408) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.427) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.408) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.425) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.426) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.4) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.415) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.416) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.405) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.387) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.39) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.389) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.389) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.381) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.381) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.367) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.369) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.326) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.335) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.326) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.341) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.403) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.47) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.2497) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.31891) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.273) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.28445) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.28036) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.28453) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.29495) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.27024) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.28912) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.29048) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.27507) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.23327) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.2493) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.2416) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.26652) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.24852) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.26927) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.19302) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.19433) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.17523) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.14981) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.16666) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.29*tanh(pt*0.00997 - 0.143)) + \ (abs(eta) > 4.0) * (0.000) } # efficiency formula for b-jets add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.75) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.827) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.853) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.855) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.862) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.868) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.865) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.863) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.857) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.851) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.838) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.831) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.817) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.796) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.772) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.76) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.743) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.703) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.638) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.605) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.572) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.541) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.567) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.603) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.6063) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.7188) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.72) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.7365) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.7462) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.7454) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.7415) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.727) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.7112) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.7112) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.6754) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.6359) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.6348) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.6115) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.5585) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.5608) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.5208) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.456) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.4524) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.388) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.3928) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.3823) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \ (abs(eta) > 4.0) * (0.000) } } ############# # b-tagging # ############# module BTagging RecoJetBTaggingMedium { set JetInputArray JetPileUpSubtractor/jets set BitNumber 1 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00654) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00921) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.00573) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00694) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0062) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00708) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00779) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00693) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00777) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00862) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01038) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01189) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01279) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.01452) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01696) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01958) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.02253) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.01787) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.02154) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.01839) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.01987) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.02351) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.02937) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.04001) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.0542) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.004236) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.006653) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.005512) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.00754) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.005813) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.006439) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.008063) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.00647) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.007583) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.008543) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.01034) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.011253) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.012945) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.014474) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.017361) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.020912) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.023139) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.010756) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.012569) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.006428) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00887) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \ (abs(eta) > 4.0) * (0.000) \ } # efficiency formula for c-jets (misidentification rate) add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.1102) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.1344) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.1025) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.1025) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.1063) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.1087) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.1124) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.109) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.111) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.1091) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1087) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1091) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1107) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1061) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1017) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.0966) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0697) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0679) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0503) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0514) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0481) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0667) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0861) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.092) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.03331) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.04361) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.03863) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.04287) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.04431) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.04452) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.04339) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0436) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0456) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.05138) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.04794) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.04004) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.04713) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.04515) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.05314) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.05143) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.05936) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.02357) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.03222) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.02621) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.01709) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.1898*tanh(pt*0.00997 - 0.143)) + \ (abs(eta) > 4.0) * (0.000) } # efficiency formula for b-jets add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.536) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.6439) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.6504) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.6716) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.6841) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.6896) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.6916) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.6882) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.6838) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.6715) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.6554) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.6366) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.6192) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.595) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.5551) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.5138) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.4884) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.4009) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.3459) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.2523) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2404) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2198) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2263) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2614) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3194) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.3254) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.4339) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.4499) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.4716) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.4766) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.4788) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.4863) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.4891) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.462) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.4583) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.4247) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.3775) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.3734) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.3348) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.2939) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.285) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.2421) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.1565) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.1522) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.1231) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.1607) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.1323) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \ (abs(eta) > 4.0) * (0.000) } } ############# # b-tagging # ############# module BTagging RecoJetBTaggingTight { set JetInputArray JetPileUpSubtractor/jets set BitNumber 2 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.000164) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.000235) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.000266) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.000329) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.000309) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.000309) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.000546) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.000499) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.000642) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.000742) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.000928) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.001323) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.001392) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00154) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.002094) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.002427) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.002927) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.001854) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002355) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002297) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.002433) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002706) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.003602) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.004987) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.007414) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.000573) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.000574) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.000521) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.000786) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.000539) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.000673) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.000934) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.000781) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.000949) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.000977) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.001168) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.000879) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.000812) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.001215) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.001679) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.001893) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.002723) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.003555) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.003881) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.005563) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.007611) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \ (abs(eta) > 4.0) * (0.000) \ } # efficiency formula for c-jets (misidentification rate) add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00531) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00567) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0064) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00673) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.00766) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00729) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00674) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00824) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00888) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00919) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01021) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01041) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01027) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00999) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01047) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01014) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.01021) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.00601) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0054) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.00487) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.00519) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.00469) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.00651) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.01299) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.00897) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.01014) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.01288) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.01392) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.01533) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.01508) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.01579) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.01106) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.01346) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.01315) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.01156) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0082) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.00439) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.00744) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.00685) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.00755) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.00706) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.00428) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.01031) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.00976) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.00374) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00854) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \ (abs(eta) > 4.0) * (0.000) } # efficiency formula for b-jets add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \ (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.2426) + \ (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.327) + \ (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.3559) + \ (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.3704) + \ (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.3824) + \ (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.3844) + \ (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.3848) + \ (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.3862) + \ (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.3778) + \ (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.3622) + \ (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.3299) + \ (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.2889) + \ (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.2815) + \ (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.253) + \ (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.221) + \ (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1963) + \ (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1739) + \ (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0992) + \ (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0788) + \ (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0581) + \ (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0534) + \ (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0521) + \ (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0626) + \ (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0826) + \ (abs(eta) <= 1.8) * (pt > 2000.0) * (0.1022) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.1562) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.2499) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.2956) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.3128) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.3147) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.3222) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.304) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.3051) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.2657) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.2578) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.2087) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.1634) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.1651) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.1353) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.109) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0799) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0699) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.054) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0718) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.0746) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.0803) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.0882) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \ (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \ (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \ (abs(eta) > 4.0) * (0.000) } } ##################### ## tau tagging ## ##################### module TauTagging TauTagging { set ParticleInputArray Delphes/allParticles set PartonInputArray Delphes/partons set JetInputArray JetPileUpSubtractor/jets set DeltaR 0.4 set TauPTMin 1.0 set TauEtaMax 4.0 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt} # default efficiency formula (misidentification rate) add EfficiencyFormula {0} {(abs(eta)<1.8)*0.006+(abs(eta)>1.8)*0.015} # efficiency formula for tau-jets add EfficiencyFormula {15} {0.65} } ##################### ## MissingET ## ##################### ### -sum of all particles after filtering neutrinos module Merger GenMissingET { add InputArray NeutrinoFilter/filteredParticles set MomentumOutputArray momentum } module Merger RecoMissingET { add InputArray EFlowMerger/eflow set MomentumOutputArray momentum } ############## ## Scalar HT # ############## module Merger GenScalarHT { # add InputArray InputArray add InputArray NeutrinoFilter/filteredParticles set EnergyOutputArray energy } module Merger RecoScalarHT { # add InputArray InputArray add InputArray EFlowMerger/eflow set EnergyOutputArray energy } ##################### # Photon efficiency # ##################### module Efficiency PhotonEfficiency { ## input particles set InputArray Calorimeter/eflowPhotons ## output particles set OutputArray photons # set EfficiencyFormula {efficiency formula as a function of eta and pt} # efficiency formula for photons set EfficiencyFormula { (pt <= 10.0) * (0.00) + \ (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \ (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \ (abs(eta) > 4.0) * (0.00)} } #################### # Photon isolation # #################### module Isolation PhotonIsolation { # particle for which calculate the isolation set CandidateInputArray PhotonEfficiency/photons # isolation collection set IsolationInputArray EFlowMergerAllTracks/eflow # select a rho for the isolation set RhoInputArray Rho/rho # output array set OutputArray photons # isolation cone set DeltaRMax 0.3 # minimum pT set PTMin 1.0 # iso ratio to cut set PTRatioMax 9999. } ####################### # Electron efficiency # ####################### module Efficiency ElectronEfficiency { set InputArray TrackPileUpSubtractor/electrons set OutputArray electrons # set EfficiencyFormula {efficiency formula as a function of eta and pt} # efficiency formula for electrons set EfficiencyFormula { (pt <= 4.0) * (0.00) + \ (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \ (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \ (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \ (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \ (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \ (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \ (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \ (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \ (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \ (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \ (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \ (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \ (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \ (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \ (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \ (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \ (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \ (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \ (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \ (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \ (abs(eta) > 4.0) * (0.00) } } ###################### # Electron isolation # ###################### module Isolation ElectronIsolation { set CandidateInputArray ElectronEfficiency/electrons # isolation collection set IsolationInputArray EFlowMergerAllTracks/eflow set RhoInputArray Rho/rho set OutputArray electrons set DeltaRMax 0.3 set PTMin 1.0 set PTRatioMax 9999. } ################### # Muon efficiency # ################### module Efficiency MuonEfficiency { set InputArray TrackPileUpSubtractor/muons set OutputArray muons # set EfficiencyFormula {efficiency as a function of eta and pt} # efficiency formula for muons set EfficiencyFormula { (pt <= 2.0) * (0.00) + \ (abs(eta) <= 4.00) * (pt > 2.0 && pt <= 3.0) * (0.51) + \ (abs(eta) <= 4.00) * (pt > 3.0 && pt <= 4.0) * (0.85) + \ (abs(eta) <= 4.00) * (pt > 4.0 && pt <= 11.0) * (0.93) + \ (abs(eta) <= 4.00) * (pt > 11. && pt <= 50.) * (0.96) + \ (abs(eta) <= 4.00) * (pt > 50. && pt <= 70.) * (0.98) + \ (abs(eta) <= 4.00) * (pt > 70.0 ) * (1.00) + \ (abs(eta) > 4.00) * (0.00) } } ################## # Muon isolation # ################## module Isolation MuonIsolation { set CandidateInputArray MuonEfficiency/muons # isolation collection set IsolationInputArray EFlowMergerAllTracks/eflow set RhoInputArray Rho/rho set OutputArray muons set DeltaRMax 0.3 set PTMin 1.0 set PTRatioMax 9999. } ############################################################################################################### # StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 # ############################################################################################################### module StatusPidFilter StatusPidFilter { ## take the particles from Pythia8 not adding pile-up set InputArray Delphes/allParticles set OutputArray filteredParticles set PTMin 5.0 } ################## # ROOT tree writer ################## module TreeWriter TreeWriter { ## branch notation :