import sys, os from ROOT import gROOT gROOT.ProcessLine(".x " + os.environ['MYSW_DIR'] + "/Utils/rootlogon.C") from ROOT import Main maker = Main.Maker() # BNBComisc maker.SetInputFile( "/Users/deltutto/RealWork/CCInclusiveEventSelection/Files/Input/ubxsec_output_mc_bnbcosmic_mcc8.9_500k.root" ) maker.SetOutputFile( "/Users/deltutto/RealWork/CCInclusiveEventSelection/Files/Output/ubxsecana_output_bnbcosmic_mcc8.9_500k.root" ) maker.SetEntries(50000) maker.SetInitialEntry(0) maker.SetBeamSpillStart(3.2) maker.SetBeamSpillEnd(4.8) maker.SetFlashShift(0.) maker.SetGainCalibration(198) maker.SetCalculatePOT(True) maker.SetIsData(False) maker.SetExtraWeight(1.028) # Flux correction maker.SetTargetFluxSystematic("total")
import sys, os, math from ROOT import gROOT gROOT.ProcessLine(".x " + os.environ['MYSW_DIR'] + "/Utils/rootlogon.C") from ROOT import Main #gROOT.SetBatch(kTRUE) analyser = Main.CC1uNPAna() #analyser.SetBNBCosmicFile ("/uboone/data/users/jiangl/ubxsec_static/v06_26_01_22_Feb/ubxsecana_output_mc_bnbcosmic_ubcodev06_26_01_22_detsyst_CV.root"); analyser.SetBNBCosmicFile( "/uboone/data/users/kirby/ubxsec_static/v06_26_01_22_Apr/detsyst/ubxsecana_output_mc_bnbcosmic_ubcodev06_26_01_22.root" ) #analyser.SetInTimeCosmicFile ("/Users/deltutto/CCInclusiveFiles/Output/ubxsecana_output_bnbcosmic_mcc8.7_test6.root") # Just a placeholder analyser.SetDirtFile( "/uboone/data/users/jiangl/ubxsec_static/v06_26_01_22_Apr19/ubxsecana_output_mc_bnbdirt_ubcodev06_26_01_22.root" ) analyser.SetBNBONFile( "/uboone/data/users/kirby/ubxsec_static/v06_26_01_22_Apr/ubxsecana_output_data_onbeam_ubcodev06_26_01_22.root" ) analyser.SetEXTBNBFile( "/uboone/data/users/kirby/ubxsec_static/v06_26_01_22_Apr/ubxsecana_output_data_offbeam_ubcodev06_26_01_22.root" ) analyser.SetBNBPOT(1.582e+20) analyser.SetBNBONTriggers(31865413.0) analyser.SetEXTBNBTriggers(71687142.0) analyser.SetPrefix("CV")
import sys, os, math from ROOT import gROOT gROOT.ProcessLine(".x " + os.environ['MYSW_DIR'] + "/Utils/rootlogon.C") from ROOT import Main analyser = Main.Analyse() analyser.SetBNBCosmicFile( "/Users/deltutto/CCInclusiveFiles/Output/ubxsecana_output_mc_bnbcosmic_ubcodev06_26_01_22__v4.root" ) # tune 1 analyser.SetInTimeCosmicFile( "/Users/deltutto/CCInclusiveFiles/Output/ubxsecana_output_bnbcosmic_mcc8.7_test6.root" ) # Just a placeholder analyser.SetBNBONFile( "/Users/deltutto/CCInclusiveFiles/Output/ubxsecana_output_data_overlay_ubcodev06_26_01_22__v4.root" ) # fake data! Overlay analyser.SetEXTBNBFile( "/Users/deltutto/CCInclusiveFiles/Output/ubxsecana_output_data_extbnb_run1_ubcodev06_26_01_22__v4.root" ) analyser.SetBNBPOT(9.1371E20) #88601 1.65319e+20 #9.6181e+20 analyser.SetBNBONTriggers( 213802895.68115577) #36749130.19319 #213802895.68115577 analyser.SetEXTBNBTriggers(72299264) # 35388924 bnbon triggers corrspond to 1.592e+20 POT # and for bnbext we have 72299264, bnbext stays the same # for bnbon we have to increase the triggers to 35388924/1.592e+20*9.6181e+20 = 213802895.68115577
import sys, os from ROOT import gROOT gROOT.ProcessLine(".x " + os.environ['MYSW_DIR'] + "/Utils/rootlogon.C") from ROOT import Main maker = Main.Simplified_Gen() calc_syst_list = ["PerfectReconstruction_neutronsDetected"] #calc_syst_list = ["PerfectReconstruction_neutronsUndetected"] #calc_syst_list = ["OurSetup_neutronsDetected"] #calc_syst_list = ["OurSetup_neutronsUndetected"] #calc_syst_list = ["DUNE_FastMC"] for systname in calc_syst_list: #maker.SetInputFile("/pnfs/genie/persistent/users/jiangl/GENIE_sample/v2_12_0_alt/gntp.30.gst.2120alt.root"); #maker.SetInputFile("/pnfs/genie/persistent/users/jiangl/GENIE_sample/v2_12_0_alt_anumu/gntp.all60.gst.anumu.root"); #numu #maker.SetInputFile("/pnfs/genie/persistent/users/jiangl/GENIE_sample/v2_12_0_alt_nue/gntp.all60.gst.nue.root"); maker.SetInputFile("/pnfs/genie/persistent/users/jiangl/GENIE_sample/v2_12_0_alt_anue/gntp.all60.gst.anue.root"); maker.SetOutputFile("/genie/app/users/jiangl/NuPRISM_output/output_simplifiedgen_anue.root"); maker.SetNeutrinoFlavor(1); # the options are 1(numu), 2(anumu), 3(nue, 4(anue) maker.SetInteractionType(4); # the options are 1 (QE), 2 (RES), 3 (DIS) 4 (MEC) 5(COH) 0 ALL;