RecoParticleFlow/python/defaults_cfi.py

0001 #response config
0002
0003 ptbins = [
0004     15, 21, 28, 37, 49, 64, 84, 114, 153, 196, 272, 330, 395, 468, 548, 686, 846, 1032, 1248, 1588, 2000, 2500, 3103, 3832, 4713, 5777, 7000
0005     #10,24,32,43,56,74,97,133,174,245,300,362,430,507,592,686,846,1032,1248,1588,2000,2500,3000,4000,6000
0006 ]
0007 etabins = [0.0, 0.5, 1.3, 2.1, 2.5, 3.0, 5.0]
0008
0009 def response_distribution_name(iptbin, ietabin):
0010     #convert 0.5 -> "05"
0011     eta_string = "{0:.1f}".format(etabins[ietabin+1]).replace(".", "")
0012     return "reso_dist_{0:.0f}_{1:.0f}_eta{2}".format(ptbins[iptbin], ptbins[iptbin+1], eta_string)
0013
0014 def genjet_distribution_name(ietabin):
0015     eta_string = "{0:.1f}".format(etabins[ietabin+1]).replace(".", "")
0016     return "genjet_pt_eta{0}".format(eta_string)
0017
0018 def recojet_distribution_name(ietabin):
0019     eta_string = "{0:.1f}".format(etabins[ietabin+1]).replace(".", "")
0020     return "recojet_pt_eta{0}".format(eta_string)
0021
0022
0023 jetResponseDir = 'ParticleFlow/JetResponse/'
0024 genjetDir = 'ParticleFlow/GenJets/'
0025
0026 #offset config
0027
0028 etaBinsOffset = [-5.191, -4.889, -4.716, -4.538, -4.363, -4.191, -4.013, -3.839, -3.664, -3.489, -3.314, -3.139, -2.964, -2.853, -2.65,
0029    -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,
0030    -0.879, -0.783, -0.696, -0.609, -0.522, -0.435, -0.348, -0.261, -0.174, -0.087, 0,
0031    0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609, 0.696, 0.783, 0.879, 0.957, 1.044, 1.131, 1.218, 1.305, 1.392, 1.479,
0032    1.566, 1.653, 1.74, 1.83, 1.93, 2.043, 2.172, 2.322, 2.5, 2.65, 2.853, 2.964, 3.139, 3.314, 3.489, 3.664, 3.839, 4.013,
0033    4.191, 4.363, 4.538, 4.716, 4.889, 5.191]
0034
0035 muLowOffset = 0
0036 muHighOffset = 100
0037
0038 npvLowOffset = 0
0039 npvHighOffset = 100
0040 pvzHighOffset = 100
0041
0042 eBinsOffset = 1000
0043 eLowOffset = 0
0044 eHighOffset = 1000
0045
0046 '''
0047 Defined in PFCandidate.cc
0048 ParticleTypes (use absolute value)
0049 211 h,       // charged hadron
0050 11  e,       // electron
0051 13  mu,      // muon
0052 22  gamma,   // photon
0053 130 h0,      // neutral hadron
0054 1   h_HF,        // HF tower identified as a hadron
0055 2   egamma_HF    // HF tower identified as an EM particle
0056 '''
0057
0058 #map cmssw reco::PFCandidate::pdgId() particle type to our type
0059 candidateDict = {211:"chm", 11:"lep", 13:"lep", 22:"ne", 130:"nh", 1:"hfh", 2:"hfe"}
0060 #our particle types (note chu = unmatched charged hadrons)
0061 candidateType = ["chm", "chu", "nh", "ne", "hfh", "hfe", "lep"]
0062
0063 # map reco::PFCandidate::pdgId() to particle type for PFCand
0064 pdgIDDict = {211: "chargedHadron", 11:"electron", 13:"muon", 22:"photon", 130:"neutralHadron", 1:"HF_hadron", 2: "HF_EM_particle"}
0065
0066 offsetPlotBaseName = 'p_offset_eta'
0067 offsetDir = 'ParticleFlow/Offset/'
0068 offsetR = 0.4
0069
0070 offsetVariableType = ["npv", "mu"]
0071
0072 def offset_name( var, ivar, itype ) :
0073     return "{0}_{1}{2}_{3}".format( offsetPlotBaseName, var, ivar, itype )