Project

General

Profile

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 and make_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