e4nu Collaboration - (e,e'1pi) with Brittany

[jtenavidal]

@BrittanyKCohen96

Here is the command I use to get the comparisons to your data:

MCLOC=/Users/juliatenavidal/Desktop/Postdoc/e4nu/FinalPionProductionAnalysis/single_pion/simulation/
DATALOC=/Users/juliatenavidal/Desktop/Postdoc/e4nu/FinalPionProductionAnalysis/single_pion/data/
OUTLOC=/Users/juliatenavidal/Desktop/Postdoc/e4nu/FinalPionProductionAnalysis/single_pion/plots/
OBS=RecoEvPion


./plot_e4nuanalysis --mc_location $MCLOC --output_location $OUTLOC --output_name clas6analysis_1pim_1GeV --input_mc_files e4nuanalysis_1pim_GEM21_11a_Rarita_LFG_Q2_01_1GeV_eCarbon_NoRad --model_names GEM21_11a_Rarita --observable_list $OBS --data_name "CLAS6 1#pi^{-} 1GeV" --analysis_id "1pim" --store_root true --data_location $DATALOC --input_data_file e4nuanalysis_1pimanalysis_e_on_1000060120_1161MeV_4MaxBkgMult_clas6data --data_name "CLAS6 1#pi^{-} 1GeV" --analysis_id "1pim" --store_root true --nofsi_file e4nuanalysis_1pim_GEM21_11a_Rarita_LFG_Q2_01_1GeV_eCarbon_NoRad_NoFSI --add-systematics BkgClosure:1,Photon:1 --observable_list_y RecoQ2 --observable_y_cuts 0.1,0.2,0.4,0.6,0.8

[jtenavidal]

Multiple instances of SetTreeBranch have been removed. Now, to add a new branch, one must only change the PlottingUtils.cxx and PlottingUtils.h

Update code for results!

[jtenavidal]

Mott scaling fixed. With the new code, by default you will not scale by Q4

[jtenavidal]

Next steps:
1 - edit https://github.com/e4nu/e4nuanalysiscode/blob/develop/PionSep2025/ConfFiles/mc_conf/run_clas6mc_1pimanalysis_eC12_1GeV.sh
https://github.com/e4nu/e4nuanalysiscode/blob/develop/PionSep2025/ConfFiles/mc_conf/run_clas6mc_1pimanalysis_eC12_2GeV.sh
https://github.com/e4nu/e4nuanalysiscode/blob/develop/PionSep2025/ConfFiles/mc_conf/run_clas6mc_1pimanalysis_eC12_4GeV.sh
and
https://github.com/e4nu/e4nuanalysiscode/blob/develop/PionSep2025/ConfFiles/mc_conf/run_clas6mc_1pipanalysis_eC12_1GeV.sh
https://github.com/e4nu/e4nuanalysiscode/blob/develop/PionSep2025/ConfFiles/mc_conf/run_clas6mc_1pipanalysis_eC12_2GeV.sh
https://github.com/e4nu/e4nuanalysiscode/blob/develop/PionSep2025/ConfFiles/mc_conf/run_clas6mc_1pipanalysis_eC12_4GeV.sh

Change Line 49 line to compute the true reconstructed instead (ComputeTrueAccCorr->ComputeTrueRecoAccCorr):
./e4nuanalysis --conf-file ${conf_file} --root-file ${path_mc_files}${RadiatedFiles[0]} --output-file ${path_output}${OutputFilesRadiated[0]} --analysis-type ComputeTrueRecoAccCorr --xsec-file ${path_xsec}${XSecFiles[0]} --rad-corr true

2 - Re-run Radiative part only (no need to re-run all of it). You can either copy paste the configuration and the line above or comment out the lines you don't want to run.
3 - Run the 2D cross section for Q2 vs W. W ranges 1, 1.45, 1.6 and above
4 - Run the 2D cross section for Q2 vs W. W ranges 1, 1.2 ,1.45, 1.6 and above
For the 2D plots, we want the slices and the coloured 2D plots.

5 - re-run closure test step to add all observables of interest. For instance, edit this file:
https://github.com/e4nu/e4nuanalysiscode/blob/develop/PionSep2025/ConfFiles/mc_conf/clas6mc_1pimanalysis_eC12_1GeV.txt
Change ObservableList RecoW,ECal,pim_mom to -> ObservableList RecoW,ECal,pim_mom,RecoQ2,pip_mom,pim_theta, ...
You run it in line 45 from all run_clas6mc_ files.

---

For the plots of section 4.2.5 in the technote, which show the per-sector variation, these are automatically generated in the new version of the plotting scripts. The plots are stored in XSecPerSector and have the following naming scheme
(output name given in the commandline) + "dataonly_stack_eventrate"+observable_name+"_persector.root"

I merged it. Now you should have access. I had to change the binning for some of my observables, the rest is untouched.

[jtenavidal]

It would be good to run the predictions with the Geant 4 FSI model.

You can add it in the configuration file as follows:
https://github.com/e4nu/e4nuanalysiscode/blob/master/ConfFiles/mc_conf/run_clas6mc_1p1pimanalysis_eC12_2GeV.sh

declare -a InputFiles=(    
    "Rarita/Carbon/2GeV/"
    "Rarita/Carbon/NoFsi/2GeV/"
    "Dipole/Carbon/2GeV/"
-->    "Geant4/Carbon/2GeV/"
)

declare -a OutputFiles=(
    "e4nuanalysis_1p1pim_GEM21_11a_Rarita_LFG_Q2_04_2GeV_eCarbon_NoRad"
    "e4nuanalysis_1p1pim_GEM21_11a_Rarita_LFG_Q2_04_2GeV_eCarbon_NoRad_NoFSI"
    "e4nuanalysis_1p1pim_GEM21_11a_Dipole_LFG_Q2_04_2GeV_eCarbon_NoRad"
 -->   "e4nuanalysis_1p1pim_GEM21_11d_Rarita_LFG_Q2_04_2GeV_eCarbon_NoRad"
)

declare -a XSecFiles=(
    "GEM21_11a_Rarita_Q2_04_C.gst.root"
    "GEM21_11a_Rarita_Q2_04_C.gst.root"
    "GEM21_11a_Dipole_Q2_04_C.gst.root"
 -->   "GEM21_11d_Rarita_Q2_04_C.gst.root"
)

2 - To run the closure test for the background, you want to use this script:
https://github.com/e4nu/e4nuanalysiscode/blob/master/src/plotting_apps/plot_BkgCheck.cxx

To run it:
./plot_bkg --observable ECal --input-file YourFile_closuretest.root

The file named as bakground_debug_${OBS}_syst.root you will need to add the systeamtic error to your data using the --bkg-systematics option.

For instance my command looks like this:

./plot_bkg --input-file $MCLOC/e4nuanalysis_1p1pim_GEM21_11a_Rarita_LFG_Q2_01_1GeV_eCarbon_NoRad_closuretest.root --observable "$OBS"

mv "bakground_debug_${OBS}syst.root" "$MCLOC/clas6analysis_1p1pim_1GeV_bakground_debug${OBS}_syst.root"

./plot_e4nuanalysis --mc_location $MCLOC --output_location $OUTLOC --output_name clas6analysis_1p1pim_1GeV
--input_mc_files e4nuanalysis_1p1pim_GEM21_11a_Rarita_LFG_Q2_01_1GeV_eCarbon_NoRad,e4nuanalysis_1p1pim_GEM21_11d_Rarita_LFG_Q2_01_1GeV_eCarbon_NoRad,e4nuanalysis_1p1pim_GEM21_11a_Dipole_LFG_Q2_01_1GeV_eCarbon_NoRad
--model_names GEM21_11a_Rarita,GEM21_11d_Rarita,Dipole
--observable_list $OBS --data_name "CLAS6 1p1#pi^{-} 1GeV"
--analysis_id "1p1pim" --store_root true --data_location $DATALOC --input_data_file e4nuanalysis_1p1pimanalysis_e_on_1000060120_1161MeV_5MaxBkgMult_clas6data
--data_name "CLAS6 1p1#pi^{-} 1GeV" --analysis_id "1p1pim" --store_root true --nofsi_file e4nuanalysis_1p1pim_GEM21_11a_Rarita_LFG_Q2_01_1GeV_eCarbon_NoRad_NoFSI
--input_rad_files e4nuanalysis_1p1pim_GEM21_11a_Rarita_LFG_Q2_01_1GeV_eCarbon_Rad \n
--bkg-systematics $MCLOC/clas6analysis_1p1pim_1GeV_bakground_debug_${OBS}_syst.root

It also plots in a separate file which will contain the comparison between the true background and the estimated.

[BrittanyKCohen96]

1GeV pi-:

AngleElvsHad:

Electron momentum:

Electron angle:

Pi- momentum:

Pi- angle:

RecoW:

RecoQ2:

[BrittanyKCohen96]

2GeV pi-:

AngleElvsHad:

Electron momentum:

Electron angle:

Pi- momentum:

Pi- angle:

RecoW:

RecoQ2:

[BrittanyKCohen96]

4GeV pi-:

AngleElvsHad:

Electron momentum:

Electron angle:

Pi- momentum:

Pi- angle:

RecoW:

RecoQ2:

[BrittanyKCohen96]

2D pi- plots:

1GeV:

2GeV:

4GeV:

[BrittanyKCohen96]

1GeV pi+ results:

AngleElvsHad:

Electron momentum:

Electron angle:

Pi+ momentum:

Pi+ angle:

RecoW:

RecoQ2:

[BrittanyKCohen96]

2GeV pi+ 1D results:

AngleElvsHad:

Electron momentum:

Electron angle:

Pi+ momentum:

Pi+ angle:

RecoW:

RecoQ2:

[BrittanyKCohen96]

4GeV pi+ 1D results:

AngleElvsHad:

Electron momentum:

Electron angle:

Pi+ momentum:

Pi+ angle:

RecoW:

RecoQ2:

[BrittanyKCohen96]

2D plots for pi+:

1GeV:

2GeV:

4GeV: