Reproducing the nue First Analysis results¶
Set up the appropriate release¶
- Setup the first analysis branch (TODO: simplify when next tag is available and s/PRIVATE/PUBLIC/ below)
setup_nova -r R15-11-03-ana.a -b maxopt newrel -t R15-11-03-ana.a mytestrel cd mytestrel srt_setup -a addpkg_svn -b CAFAna R15-11-03-ana-br make
Ignore the warning about sam_web_client
. If you find yourself needing to use a test release or to make a commit to this branch please consult NOvA Snapshot and Release Policy
Making the ND Data/MC plot included in the nue paper¶
- Run
plot_paper_nd_datamc.C
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/plot_paper_nd_datamc.C
Making the FD Data/MC plot included in the nue paper¶
- Run
plot_fd_datamc.C
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/plot_fd_datamc.C false
Assumes the extrapolation file already exists. Points to the FirstAna/nue_first_ana_sens_newposc3.root file currently.
Note: Running this with option "true" will produce a large array of data/MC comparison plots but assumes a local file
produced first. In order to make this input file:
- Run
make_nue_fd_datamc.C
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_nue_fd_datamc.C false
This takes a very long time.
Producing extrapolation file and event counts¶
- Run
make_nue_firstana_extrap.C
andmake_nue_firstana_exposure.C
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_nue_firstana_extrap.C cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_nue_firstana_exposure.C true true true
- Wait many hours (TODO about how long?)
- Or, here's one I prepared earlier:
/nova/ana/nu_e_ana/FirstAna/nue_first_ana_sens_newposc3.root # MC /nova/ana/nu_e_ana/FirstAna/nue_first_ana_sens_data_v2.root # data
The old files with the wrong detector positions used in the originally presented results are available under oldpos/
Making contours/slices¶
- Run
contours.C
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/contours.C nue_first_ana_sens_newposc3.root nue_first_ana_sens_data_v2.root lid cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/contours.C nue_first_ana_sens_newposc3.root nue_first_ana_sens_data_v2.root lem
Applying Feldman-Cousins corrections¶
- Use a batch system to run 1000 instances of this
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_surfs.C nue_first_ana_sens_newposc3.root lid 10000 0.5 true true cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_surfs.C nue_first_ana_sens_newposc3.root lem 10000 0.5 true true
- and this
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_slices.C nue_first_ana_sens_all_newposc3.root lid 10000 0.5 true true cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_slices.C nue_first_ana_sens_all_newposc3.root lem 10000 0.5 true true
hadd
the results into 8 "fccol" files. Official versions in
/nova/ana/nu_e_ana/FirstAna/fc/expts/
- Convert experiments into surfaces
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_surfs.C nue_first_ana_sens_newposc3.root lid 0 0.5 true true cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_surfs.C nue_first_ana_sens_newposc3.root lem 0 0.5 true true cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_slices.C nue_first_ana_sens_newposc3.root lid 0 0.5 true true cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_slices.C nue_first_ana_sens_newposc3.root lem 0 0.5 true true
- Official versions in
/nova/ana/nu_e_ana/FirstAna/fc/surfs/
- Make corrected contours and slices
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/contours.C nue_first_ana_sens_newposc3.root nue_first_ana_sens_data_v2.root lid 0.5 fcsurf_nh_lid_0.5_syst_osc.root fcsurf_ih_lid_0.5_syst_osc.root fcslice_nh_lid_0.5_syst_osc.root fcslice_ih_lid_0.5_syst_osc.root cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/contours.C nue_first_ana_sens_newposc3.root nue_first_ana_sens_data_v2.root lem 0.5 fcsurf_nh_lem_0.5_syst_osc.root fcsurf_ih_lem_0.5_syst_osc.root fcslice_nh_lem_0.5_syst_osc.root fcslice_ih_lem_0.5_syst_osc.root
You can control whether these macros make the paper-style plots or not by editing the gPaper
variable at the top.
Compute significance to reject theta13=0¶
- Run make_fc_zero_rej_hybrid.C
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_zero_rej_hybrid.C nue_first_ana_sens_newposc3.root lid 6 100000 0.5 true true cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_fc_zero_rej_hybrid.C nue_first_ana_sens_newposc3.root lem 11 100000 0.5 true true
Making the text file for theorists¶
- Run
make_theorist_text_file.C
cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_theorist_text_file.C nue_first_ana_sens_newposc3.root LID > table_lid.txt cafe -bq $SRT_PRIVATE_CONTEXT/CAFAna/nue/FirstAna/make_theorist_text_file.C nue_first_ana_sens_newposc3.root LEM > table_lem.txt