import ROOT import os import math import numpy as np import copy from array import array from Workflow_Handler import Workflow_Handler #Supress the opening of many Canvas's ROOT.gROOT.SetBatch(True) isBDT_with_Wmass = False # If true, pT(pi) and ET(gamma) in the BDT are normalized to Wmass myWF = Workflow_Handler("Signal","Data",isBDT_with_Wmass) ttbar_sig_calib_file = ROOT.TFile.Open("ttbar_signal_ratio.root") ttbar_sig_calib = ttbar_sig_calib_file.Get("ttbar_signal_ratio") ##Bools for rundom SF variation random_mu_SF = False #------if True, muon scale factors are sampled from a Gaussian random_ele_SF = False #------if True, electron scale factors are sampled from a Gaussian random_ph_SF = False #------if True, photon scale factors are sampled from a Gaussian ############################################################################ # # #-------------------------- Integrated luminosity -------------------------# # # ############################################################################ #Normalize to this luminsity, in fb-1 #luminosity_norm = 36.46
import ROOT import math import numpy as np ##normalize to the 2016 lumi luminosity_norm = 5. from Workflow_Handler import Workflow_Handler myWF = Workflow_Handler("Signal_H800_A300") def is_Event_selected(jet_btag,jet_pt): """Save events according to some basic selection criteria""" btag_cut = jet_btag[0] > 0.97 and jet_btag[1] > 0.97 and jet_btag[2] > 0.97 and jet_btag[3] > 0.97 #btag_cut = jet_btag[0] > 0.89 and jet_btag[1] > 0.89 and jet_btag[2] > 0.97 and jet_btag[3] > 0.97 #pt_cut = jet_pt[0] > 165. and jet_pt[1] > 130. and jet_pt[2] > 130. and jet_pt[3] > 110. #pt_cut = jet_pt[3] > 50. return btag_cut #and pt_cut ##Here starts the program Norm_Map = myWF.get_normalizations_map() steps_cut1 = 20 cut1_init = 70. cut1_stepsize = 10. steps_cut2 = 20 cut2_init = 70. cut2_stepsize = 10. steps_cut3 = 20
import ROOT import os from Workflow_Handler import Workflow_Handler myWF = Workflow_Handler("Signal_H800_A300") if not os.path.exists("plots"): os.makedirs("plots") ##Do all the scaling to this luminosity, in fb-1 luminosity_norm = 36. signal_magnify = 100. ##These are the histograms to be merged list_histos = [ "h_jet1pt", "h_jet2pt", "h_jet3pt", "h_jet4pt", "h_jet1Btag", "h_jet2Btag", "h_jet3Btag", "h_jet4Btag", "h_delta_Phi_pair", "h_delta_Eta_pair", "h_Events", "h_jet1eta", "h_jet2eta" ] ##Here starts the program Norm_Map = myWF.get_normalizations_map() ##Get the handlers for all the histos h_QCD = dict() h_QCD[list_histos[0]] = ROOT.TH1F(list_histos[0], "p_{t} of 1st jet", 200, 0., 500.) h_QCD[list_histos[1]] = ROOT.TH1F(list_histos[1], "p_{t} of 2nd jet", 200, 0., 500.) h_QCD[list_histos[2]] = ROOT.TH1F(list_histos[2], "p_{t} of 3rd jet", 200, 0., 500.)
import ROOT import math import numpy as np ##normalize to the 2016 lumi luminosity_norm = 5. from Workflow_Handler import Workflow_Handler myWF = Workflow_Handler("Signal_H800_A300") def is_Event_selected(jet_btag, jet_pt): """Save events according to some basic selection criteria""" btag_cut = jet_btag[0] > 0.97 and jet_btag[1] > 0.97 and jet_btag[ 2] > 0.97 and jet_btag[3] > 0.97 #btag_cut = jet_btag[0] > 0.89 and jet_btag[1] > 0.89 and jet_btag[2] > 0.97 and jet_btag[3] > 0.97 #pt_cut = jet_pt[0] > 165. and jet_pt[1] > 130. and jet_pt[2] > 130. and jet_pt[3] > 110. #pt_cut = jet_pt[3] > 50. return btag_cut #and pt_cut ##Here starts the program Norm_Map = myWF.get_normalizations_map() steps_cut1 = 20 cut1_init = 70. cut1_stepsize = 10. steps_cut2 = 20
import ROOT import math import numpy as np import sys #sys.path.insert(0, '/afs/cern.ch/user/r/rselvati/wpigamma/CMSSW_8_0_28/src/StandardModel/WPiGamma/test/scripts') from scripts.create_normalization_table import get_xsec_fromsample ##normalize to the 2016 lumi luminosity_norm = 36.46 from Workflow_Handler import Workflow_Handler myWF = Workflow_Handler("Signal", "Data", isMedium=True) #-----Some selection bools------# is_pi_gamma = False is_ele_gamma = True is_mu_pT = False is_ele_pT = False is_deltaphi_mu_pi = False is_deltaphi_ele_pi = False is_piIso_03 = False is_piIso_05 = False def is_Event_selected(Wmass, nBjets, lepton_iso, ele_gamma_InvMass): #,lep_pT):#,pi_pT,gamma_eT): """Save events according to some basic selection criteria""" bjet_cut = nBjets >= 2. mass_cut_down = Wmass >= 50.
import ROOT from Workflow_Handler import Workflow_Handler myWF = Workflow_Handler("Signal_H500_A200") ##Do all the scaling to this luminosity, in fb-1 luminosity_norm = 10. signal_magnify = 100. ##These are the histograms to be merged list_histos = ["h_jet1pt", "h_jet2pt", "h_jet3pt", "h_jet4pt", "h_jet1Btag", "h_jet2Btag", "h_jet3Btag", "h_jet4Btag", "h_delta_Phi_pair", "h_delta_Eta_pair", "h_Events", "h_jet1eta", "h_jet2eta"] ##Here starts the program Norm_Map = myWF.get_normalizations_map() ##Get the handlers for all the histos h_QCD = dict() h_QCD[list_histos[0]] = ROOT.TH1F(list_histos[0], "p_{t} of 1st jet", 200, 0., 500.) h_QCD[list_histos[1]] = ROOT.TH1F(list_histos[1], "p_{t} of 2nd jet", 200, 0., 500.) h_QCD[list_histos[2]] = ROOT.TH1F(list_histos[2], "p_{t} of 3rd jet", 200, 0., 500.) h_QCD[list_histos[3]] = ROOT.TH1F(list_histos[3], "p_{t} of 4th jet", 200, 0., 500.) h_QCD[list_histos[4]] = ROOT.TH1F(list_histos[4], "Btag value of 1st jet", 50, 0.89, 1.) h_QCD[list_histos[5]] = ROOT.TH1F(list_histos[5], "Btag value of 2nd jet", 50, 0.89, 1.) h_QCD[list_histos[6]] = ROOT.TH1F(list_histos[6], "Btag value of 3rd jet", 50, 0.89, 1.) h_QCD[list_histos[7]] = ROOT.TH1F(list_histos[7], "Btag value of 4th jet", 50, 0.89, 1.) h_QCD[list_histos[8]] = ROOT.TH1F(list_histos[8], "#Delta_{#phi} between the two jet pairs", 50, 0., 6.28) h_QCD[list_histos[9]] = ROOT.TH1F(list_histos[9], "#Delta_{#eta} between the two jet pairs", 50, -10., 10.) h_QCD[list_histos[10]] = ROOT.TH1F(list_histos[10], "Events after pre-selection steps", 6, 0., 6.) h_QCD[list_histos[11]] = ROOT.TH1F(list_histos[11], "#eta of 1st jet", 50, -10., 10.) h_QCD[list_histos[12]] = ROOT.TH1F(list_histos[12], "#eta of 2nd jet", 50, -10., 10.)
'runningEra_option', help= 'Type <<0>> for 2016, <<1>> for 2017, <<2>> for 2018, <<3>> for combination 2016+2017' ) args = p.parse_args() # Switch from muon to electron channel, and from 2016 to 2017 if args.isBDT_option == "preselection": isBDT = False if args.isBDT_option == "BDT": isBDT = True runningEra = int(args.runningEra_option) #---------------------------------# isBDT_with_Wmass = False # If true, pT(pi) and ET(gamma) in the BDT are normalized to Wmass myWF = Workflow_Handler("Signal", "Data", isBDT_with_Wmass, runningEra) #Bools for rundom SF variation random_mu_SF = False #------if True, muon scale factors are sampled from a Gaussian random_ele_SF = False #------if True, electron scale factors are sampled from a Gaussian random_ph_SF = False #------if True, photon scale factors are sampled from a Gaussian ############################################################################ # # #-------------------------- Integrated luminosity -------------------------# # # ############################################################################ #Normalize to this luminsity, in fb-1 luminosity_norm_2016 = 35.86