def main(): '''Main function which can run the hadronic and the semi-leptonic decay. The study is selected at the beginning with the parameters. Nothing should be changed in the function''' # parameters to choose the study ecm = 162.6 m_w = 79.385 width = '2' reconstruction = "direct" reco_level = 'reco' # NO NEED reco_type = '' FSI = 'BE' # Analysis histdict = gethists() filename = "/afs/cern.ch/user/s/slablack/fcc/heppy/ee_ZZ_qqqqBkg_162_50K_ee_kt_SM.root" print("Reading from file {}".format(filename)) rootfile = TFile(filename) tree = rootfile.Get("events") # output file outdir = "/afs/cern.ch/user/s/slablack/fcc/output_wmass_test/output_ee_ZZ_qqqqBkg_162_diWmasscut_Bkg_50K_ee_kt_SM_{}" #.format(loop) # outdir = "./output/hadronic/FSI_study/shift/woTreatment/W_mass_{}_M{}_width{}".format(ecm, str(m_w)[:2], width) # Function call w_mass_hadronic(tree, ecm, m_w, histdict, reconstruction, reco_level, reco_type) histogram.save_root_file(histdict, outdir)
def main(): '''Main function which can run the hadronic and the semi-leptonic decay. The study is selected at the beginning with the parameters. Nothing should be changed in the function''' #input and output file from runWmass.sh fileinput6 = sys.argv[1] fileoutput6 = sys.argv[2] # parameters to choose the study ecm = 162.6 m_w = 79.385 width = '2' reconstruction = "direct" reco_level = 'reco' # NO NEED reco_type = '' FSI = 'BE' # Analysis histdict = gethists() filename = fileinput6 print("Reading from file {}".format(filename)) rootfile = TFile(filename) tree = rootfile.Get("events") # output file outdir = fileoutput6 print("Reading from file {} output file {}".format(filename,outdir)) # outdir = "./output/hadronic/FSI_study/shift/woTreatment/W_mass_{}_M{}_width{}".format(ecm, str(m_w)[:2], width) # Function call w_mass_hadronic(tree, ecm, m_w, histdict, reconstruction, reco_level, reco_type) histogram.save_root_file(histdict, outdir)
def main(): '''Main function . The study is selected at the beginning with the parameters. Nothing should be changed in functions ''' # parameters to choose the study nevts = "50K" ecm = 365 m_w = 80.385 # Analysis filename = "../heppy_output/semi-leptonic/ee_WW_qqll_{}_{}_M{}_FSR_nonFSR.root".format( ecm, nevts, str(m_w)[:2]) rootfile = TFile(filename) tree = rootfile.Get("events") # output file outdir = "./output/plots/semi-leptonic/ISR_FSR_study/{}GeV_FSR_nonFSR".format( ecm) # histogram.create_output_dir(outdir) # Function call histdict = gethists() mis_energy(tree, histdict) histogram.save_root_file(histdict, outdir)
def main(): ''' Main function, define paramters and process the correlation matrix definition''' # parameters to choose the study ecm = 365 #GeV m_w = 80.385 #GeV width = 2 reco_type = '_pcut' # input filename = "../heppy_output/hadronic/ref_templates/ee_WW_qqqq_{}_M{}_width{}.root".format(ecm, str(m_w)[:2], width) rootfile = TFile(filename) tree = rootfile.Get("events") # Output outdir = "./output/pythia_2018/hadronic/correlation_matrix/corre_matrix_{}_M{}_width{}{}".format(ecm, str(m_w)[:2], width, reco_type) histdict = gethists() coefficients(tree, histdict, reco_type) histogram.save_root_file(histdict, outdir) print_corr_coef_into_file(histdict, outdir) h_list = ["h_uncert_phi"] for h_ in h_list: fit(h_, outdir)
def main(): '''Main function which can run the hadronic and the semi-leptonic decay. The study is selected at the beginning with the parameters. Nothing should be changed in the function''' # parameters to choose the study nevts = "50K" ecm = 365 m_w = 80.385 width = '3' reco_level = "reco" # Analysis histdict = gethists() filename = "../heppy_output/pythia_2018/semi-leptonic/onlyMuons/width_variation/ee_WW_qqll_{}_{}_M{}_width{}.root".format( ecm, nevts, str(m_w)[:2], width) print("Reading from file {}".format(filename)) rootfile = TFile(filename) tree = rootfile.Get("events") # output file outdir = "./output/pythia_2018/semi-leptonic/width_variation/W_mass_{}_{}_M{}_width{}_direct".format( ecm, nevts, str(m_w)[:2], width) w_mass_semi_leptonic(ecm, tree, histdict, reco_level) histogram.save_root_file(histdict, outdir)
def main(): # input indir = '../heppy_output/reconnectionStudy/' infile = indir + 'ee_WW_qqqq_300_50K_M80_genLevel_ee_kt_woReconnection_woISR_woFSR_woRemnant_perfectPairing.root' rootfile = TFile(infile) tree = rootfile.Get("events") # output outdir = "./output/testHadronicChannel/genLevelStudy/reconnection/angleStudy" histdict = gethists() angle_study(tree, histdict) histogram.save_root_file(histdict, outdir)