Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-11-30 01:55:03

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 )