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;
