Project

General

Profile

Muon Monitors g4numi » History » Version 1

Version 1/5 - Next » - Current version
Tyler Rehak, 02/15/2019 09:28 AM
rough outline


Muon Monitors g4numi

Get current version of g4numi-minerva

Generate g4numi_ files with whatever interesting parameter changes

Those are in /persistent/test

nu2mubatch translates g4numi_ files to muon_ files, outputs to /data

hadd muon_ files into 4 files to /persistent/nu2mu

/pnfs/minerva/persistent/users/trehak/flux/test/nu2mu

hadd muon_me000z200i_merged_0_0012.root muon_me000z200i_[0-4]_0012.root; hadd muon_me000z200i_merged_1_0012.root muon_me000z200i_[5-9]_0012.root; hadd muon_me000z200i_merged_2_0012.root muon_me000z200i_1[0-4]_0012.root; hadd muon_me000z200i_merged_3_0012.root muon_me000z200i_1[5-9]_0012.root;

g4numinu2mu translates muon_ files to hadmm_ files in /persistent

source setup_beamsim.sh

/minerva/app/users/trehak/g4numinu2mu/g4numi

python ProcessG4NuMInu2mu.py --infile /pnfs/minerva/persistent/users/trehak/flux/test/nu2mu/muon_me000z200i_merged_0_0012.root --infilename muon_me000z200i_merged_0_0012.root --run_number 120; python ProcessG4NuMInu2mu.py --infile /pnfs/minerva/persistent/users/trehak/flux/test/nu2mu/muon_me000z200i_merged_1_0012.root --infilename muon_me000z200i_merged_1_0012.root --run_number 121; python ProcessG4NuMInu2mu.py --infile /pnfs/minerva/persistent/users/trehak/flux/test/nu2mu/muon_me000z200i_merged_2_0012.root --infilename muon_me000z200i_merged_2_0012.root --run_number 122; python ProcessG4NuMInu2mu.py --infile /pnfs/minerva/persistent/users/trehak/flux/test/nu2mu/muon_me000z200i_merged_3_0012.root --infilename muon_me000z200i_merged_3_0012.root --run_number 123;

combine hadmm for each run using hadd to /data

Now plot the hadmm
source: hadmmntuples in /data/hadmmNtuple
saved plots in Documents/nu2muplots

hadmm->Draw("hadmmdata.mmpz0","hadmmdata.mmpz0>80000")
hadmm
>Draw("hadmmdata.mmpz1","hadmmdata.mmpz1>80000")
hadmm
>Draw("hadmmdata.mmpz2","hadmmdata.mmpz2>80000")
hadmm
>Draw("hadmmdata.mmxpos0","hadmmdata.mmxpos0>50000")
hadmm
>Draw("hadmmdata.mmxpos1","hadmmdata.mmxpos1>50000")
hadmm
>Draw("hadmmdata.mmxpos2","hadmmdata.mmxpos2>-50000")