hDict[i]["h_mass_passed_1b"].Fill(tchain.mjj,bWeight(SFs,1))
hDict[i]["h_mass_passed_2b"].Fill(tchain.mjj,bWeight(SFs,2))
hDict[i]["h_mass_passed_le1b"].Fill(tchain.mjj,bWeight(SFs,1))
hDict[i]["h_mass_passed_le1b"].Fill(tchain.mjj,bWeight(SFs,2))
return hDict
if __name__ == '__main__':
rt.gROOT.SetBatch()
setTDRStyle.setTDRStyle()
###################################################################
parser = OptionParser(usage=usage)
parser.add_option('-f','--flavour',dest="flavour",type="string",default="none",
help="Name of the signal flavour")
parser.add_option('-m','--model',dest="model",type="string",default="qq",
help="Name of the signal model")
parser.add_option('-s','--su',dest='su',type = 'string',default='central',help='central/up/down')
(options,args) = parser.parse_args()
flavour = options.flavour
model = options.model
su = options.su
print "selected flavour:",flavour
import ROOT as rt
from rootTools import tdrstyle as setTDRStyle
rt.gROOT.SetBatch()
setTDRStyle.setTDRStyle()
from ROOT import *
import math
import sys,os
treename = 'deepntuplizer/tree'
VariableList = ['n_sv','nsv','sv_pt','sv_phi','sv_etarel','sv_phirel','sv_deltaR','sv_mass','sv_ntracks','sv_chi2','sv_ndf','sv_normchi2','sv_dxy','sv_dxyerr','sv_dxysig','sv_d3d','sv_d3derr','sv_costhetasvpv','sv_enratio','npv','rho','ntreInt','event_no','jet_no','gen_pt','Delta_gen_pt','isB','isGBB','isBB','isLeptonicB','isLeptonicB_C','isC','isGCC','isCC','isUD','isS','isG','isUndefined','genDecay','isPhysB','isPhysGBB','isPhysBB','isPhysLeptonicB','isPhysLeptonicB_C','isPhysC','isPhysGCC','isPhysCC','isPhysUD','isPhysS','isPhysG','isPhysUndefined','jet_pt','jet_corr_pt','jet_eta','jet_phi','jet_mass','jet_energy','jet_looseID','jet_qgI','QG_ptD','QG_axis2','QG_mult','y_multiplicity','y_charged_multiplicity','y_neutral_multiplicity','y_ptD','y_axis1','y_axis2','y_pt_dr_log','muons_numbers','electrons_number','muons_isLooseMuon','muons_isTightMuon','muons_isSoftMuon','muons_isHighPtMuon','muons_pt','muons_relEta','muons_relPhi','muons_energy','electrons_pt','electrons_relEta','electrons_relPhi','electrons_energy','gen_pt_Recluster','gen_pt_WithNu','Delta_gen_pt_recluster','Delta_gen_pt_WithNu','nCpfcan','Cpf_can_pt','Cpfca_phi','Cpfcan_phi','cpfcan_ptrel','cpfcan_erel','Cpfcan_phirel','Cpfcan_etarel','Cpfcan_deltaR','Cpfcan_puppiw','Cpfcan_dxy','Cpfcan_dxyerrinv','Cpfcan_dxyig','Cpfcan_dz','Cpfcan_VTX_ass','Cpfcan_fromPV','Cpfcan_drminsv','Cpfcan_vertex_rho','Cpfcan_vertex_phirel','Cpfcan_vertex_etarel','Cpfcan_chi2','Cpfcan_quality','nNpfcand','Npfcan_pt','Npfcan_eta','Npfcan_phi','Npfcan_ptrel','Npfcan_erel','Npfcan_puppiw','Npfcan_phirel','Npfcan_etarel','Npfcan_deltaR','Npfcan_isGamma']
VariableList=['n_sv','nsv','sv_pt','sv_phi','sv_etarel','sv_phirel','sv_deltaR','sv_mass','sv_ntracks','sv_chi2','sv_ndf','sv_normchi2','sv_dxy','sv_dxyerr','sv_dxysig','sv_d3d','sv_d3derr','sv_costhetasvpv','sv_enratio','jet_pt','jet_corr_pt','jet_eta','jet_phi','jet_mass','jet_energy']
VariableList = ['n_sv','sv_pt','sv_phi','sv_eta','sv_deltaR','sv_mass','sv_ntracks','sv_chi2','jet_corr_pt','jet_eta','jet_phi','jet_mass','jet_energy']
VariableList = ['n_sv', 'nsv', 'sv_pt', 'sv_eta', 'sv_phi', 'sv_etarel', 'sv_phirel', 'sv_deltaR', 'sv_mass', 'sv_ntracks', 'sv_chi2', 'sv_ndf', 'sv_normchi2', 'sv_dxy', 'sv_dxyerr', 'sv_dxysig', 'sv_d3d', 'sv_d3derr', 'sv_d3dsig', 'sv_costhetasvpv', 'sv_enratio', 'npv', 'rho', 'ntrueInt', 'event_no', 'jet_no', 'gen_pt', 'Delta_gen_pt', 'isB', 'isGBB', 'isBB', 'isLeptonicB', 'isLeptonicB_C', 'isC', 'isGCC', 'isCC', 'isUD', 'isS', 'isG', 'isUndefined', 'genDecay', 'isPhysB', 'isPhysGBB', 'isPhysBB', 'isPhysLeptonicB', 'isPhysLeptonicB_C', 'isPhysC', 'isPhysGCC', 'isPhysCC', 'isPhysUD', 'isPhysS', 'isPhysG', 'isPhysUndefined', 'jet_pt', 'jet_corr_pt', 'jet_eta', 'jet_phi', 'jet_mass', 'jet_energy', 'jet_looseId', 'jet_qgl', 'QG_ptD', 'QG_axis2', 'QG_mult', 'y_multiplicity', 'y_charged_multiplicity', 'y_neutral_multiplicity', 'y_ptD', 'y_axis1', 'y_axis2', 'y_pt_dr_log', 'muons_number', 'electrons_number', 'muons_isLooseMuon', 'muons_isTightMuon', 'muons_isSoftMuon', 'muons_isHighPtMuon', 'muons_pt', 'muons_relEta', 'muons_relPhi', 'muons_energy', 'electrons_pt', 'electrons_relEta', 'electrons_relPhi', 'electrons_energy', 'gen_pt_Recluster', 'gen_pt_WithNu', 'Delta_gen_pt_Recluster', 'Delta_gen_pt_WithNu', 'pfCombinedInclusiveSecondaryVertexV2BJetTags', 'pfCombinedMVAV2BJetTags', 'pfDeepCSVJetTags_probb', 'pfDeepCSVJetTags_probbb', 'pfDeepCSVJetTags_probc', 'pfDeepCSVJetTags_probudsg', 'pfJetBProbabilityBJetTags', 'pfJetProbabilityBJetTags', 'softPFElectronBJetTags', 'softPFMuonBJetTags', 'n_Cpfcand', 'nCpfcand', 'Cpfcan_pt', 'Cpfcan_eta', 'Cpfcan_phi', 'Cpfcan_ptrel', 'Cpfcan_erel', 'Cpfcan_phirel', 'Cpfcan_etarel', 'Cpfcan_deltaR', 'Cpfcan_puppiw', 'Cpfcan_dxy', 'Cpfcan_dxyerrinv', 'Cpfcan_dxysig', 'Cpfcan_dz', 'Cpfcan_VTX_ass', 'Cpfcan_fromPV', 'Cpfcan_drminsv', 'Cpfcan_vertex_rho', 'Cpfcan_vertex_phirel', 'Cpfcan_vertex_etarel', 'Cpfcan_BtagPf_trackMomentum', 'Cpfcan_BtagPf_trackEta', 'Cpfcan_BtagPf_trackEtaRel', 'Cpfcan_BtagPf_trackPtRel', 'Cpfcan_BtagPf_trackPPar', 'Cpfcan_BtagPf_trackDeltaR', 'Cpfcan_BtagPf_trackPtRatio', 'Cpfcan_BtagPf_trackPParRatio', 'Cpfcan_BtagPf_trackSip3dVal', 'Cpfcan_BtagPf_trackSip3dSig', 'Cpfcan_BtagPf_trackSip2dVal', 'Cpfcan_BtagPf_trackSip2dSig', 'Cpfcan_BtagPf_trackDecayLen', 'Cpfcan_BtagPf_trackJetDistVal', 'Cpfcan_BtagPf_trackJetDistSig', 'Cpfcan_isMu', 'Cpfcan_isEl', 'Cpfcan_lostInnerHits', 'Cpfcan_numberOfPixelHits', 'Cpfcan_chi2', 'Cpfcan_quality', 'n_Npfcand', 'nNpfcand', 'Npfcan_pt', 'Npfcan_eta', 'Npfcan_phi', 'Npfcan_ptrel', 'Npfcan_erel', 'Npfcan_puppiw', 'Npfcan_phirel', 'Npfcan_etarel', 'Npfcan_deltaR', 'Npfcan_isGamma', 'Npfcan_HadFrac', 'Npfcan_drminsv', 'trackJetPt', 'jetNTracks', 'TagVarCSV_jetNSecondaryVertices', 'TagVarCSV_trackSumJetEtRatio', 'TagVarCSV_trackSumJetDeltaR', 'TagVarCSV_vertexCategory', 'TagVarCSV_trackSip2dValAboveCharm', 'TagVarCSV_trackSip2dSigAboveCharm', 'TagVarCSV_trackSip3dValAboveCharm', 'TagVarCSV_trackSip3dSigAboveCharm', 'n_TagVarCSV_jetNSelectedTracks', 'TagVarCSV_jetNSelectedTracks', 'TagVarCSVTrk_trackPtRel', 'TagVarCSVTrk_trackDeltaR', 'TagVarCSVTrk_trackPtRatio', 'TagVarCSVTrk_trackSip3dSig', 'TagVarCSVTrk_trackSip2dSig', 'TagVarCSVTrk_trackDecayLenVal', 'TagVarCSVTrk_trackJetDistVal', 'n_TagVarCSV_jetNTracksEtaRel', 'TagVarCSV_jetNTracksEtaRel', 'TagVarCSV_trackEtaRel', 'trackPParRatio', 'trackSip2dVal', 'trackSip3dVal', 'trackMomentum', 'trackEta', 'trackPPar', 'n_StoredVertices', 'NStoredVertices', 'TagVarCSV_vertexMass', 'TagVarCSV_vertexNTracks', 'TagVarCSV_vertexEnergyRatio', 'TagVarCSV_vertexJetDeltaR', 'TagVarCSV_flightDistance2dVal', 'TagVarCSV_flightDistance2dSig', 'TagVarCSV_flightDistance3dVal', 'TagVarCSV_flightDistance3dSig']
Sample={}
#Sample['bstar'] = {
#'0500':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_500GeV.root',
#'1000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_1000GeV.root',
#'2000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_2000GeV.root',
#'3000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_3000GeV.root',
#'4000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_4000GeV.root',
#'5000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_5000GeV.root',
#'6000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_6000GeV.root',
#'7000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_7000GeV.root',
#'8000':'/eos/cms/store/group/phys_exotica/dijet/Dijet13TeV/TylerW/BTV/bstar/bstar_8000GeV.root',
def setStyle(self):
print "Drawer::setStyle"
setTDRStyle.setTDRStyle()
gStyle.SetTitleSize(0.06, "XYZ")
gStyle.SetLabelFont(42, "XYZ")
gStyle.SetLabelSize(0.05, "XYZ")
gStyle.SetCanvasBorderMode(0)
gStyle.SetFrameBorderMode(0)
gStyle.SetCanvasColor(0)
gStyle.SetPadTickX(1)
gStyle.SetPadTickY(1)
gStyle.SetPadLeftMargin(0.15)
gStyle.SetPadRightMargin(0.05)
gStyle.SetPadTopMargin(0.05)
gStyle.SetPadBottomMargin(0.15)
TGaxis.SetMaxDigits(3)
gROOT.ForceStyle()
gROOT.Reset()
tdrstyle.setTDRStyle()
gROOT.SetStyle('tdrStyle')
gROOT.ForceStyle()
#change the CMS_lumi variables (see CMS_lumi.py)
CMS_lumi.writeExtraText = 1
CMS_lumi.extraText = "Simulation Preliminary"
CMS_lumi.lumi_sqrtS = "13 TeV" # used with iPeriod = 0, e.g. for simulation-only plots (default is an empty string)
iPos = 11
if( iPos==0 ): CMS_lumi.relPosX = 0.12
H_ref = 800;
W_ref = 800;
W = W_ref
H = H_ref