Wiki¶

Simulation of nucleon decay events at DUNE FD¶

To simulate nucleon decay events a GENIE-based generation has been implemented. A total of 68 single-nucleon exclusive decay channels listed in the 2016 update of the PDG have been included, the list includes 2-, 3- and 5-body decays.

The following is a quick guide on how to generated nucleon decay events using DUNE software platform. Before you starting you most have your computing privileges and accounts, including an account on the dune gpvm nodes (e.g. dunegpvm01.fnal.gov) set. The following assumes you are working on one of the dune gpvm nodes.

1. First do
   

   source /grid/fermiapp/products/dune/setup_dune.sh

2. Then, setup dunetpc (chose your version you want to setup, to check the list of versions available you can do "ups list -aK+ dunetpc"
   

    setup dunetpc v06_82_00 -q "prof:e15"

   
   This will automatically setup larsoft as well, dunetpc is particular implementation of LArsoft. You should then be able to run any fhicl file. (For more information related to dunetpc please visit dunetpc wiki page https://cdcvs.fnal.gov/redmine/projects/dunetpc/wiki/_Tutorial_)
3. Then, you can generate nucleon decay events using the corresponding fhicl file under /dunetpc/fcl/dunefd/gen/ndk/prodndk_dune10kt_1x2x6.fcl by doing
   

   lar -c prodndk_dune10kt_1x2x6.fcl -n 10 

   
   This will generated 10 events of proton decay p --> K + nu_bar using DUNE FD geometry 1x2x6 geometry configuration, the decay mode by default is 19 following the PDG convention. To generated a different decay mode you will need to add these lines to your fhicl file
   

   physics.producers.generator.DecayMode: 1
   physics.producers.generator.DecayedNucleon: 2212
   

   
   This would generated p-->pi_zero+e^{+}. Full decay mode list can be found here. (Nucleon decay mode list)
4. Once you have finished with the generation of events you will need to following the standard geant4 simulation, detector simulation and reconstruction using the corresponding fhicl files. You can find them at

/dunetpc/fcl/dunefd/g4/standard_g4_dune10kt_1x2x6.fcl
/dunetpc/fcl/dunefd/detsim/standard_detsim_dune10kt_1x2x6.fcl
/dunetpc/fcl/dunefd/reco/standard_reco_dune10kt_1x2x6.fcl

dunendk Analysis module¶

To analyze the reconstructed LArSoft files an analysis module has been developed, the module will save reconstructed objects, tracks, showers, PIDA values, etc and truth information into a ROOT Tree for further analysis
To use the analysis module a quick guide is described below

1. Follow the two first steps from the previous section
2. Now, you should make a new directory where you intend to put your code (Do not use your home directory)
   

   mkdir mydunendkAna
   cd mydunendkAna
   

3. Then, inside the new directory, do
   

   mrb newDev
   

4. Do the following to use the local release
   

   source localProducts_XXXX/setup
   

5. Now, pull down dunendk
   

   cd srcs
   mrb g dunendk
   

   
   IMPORTANT NOTE: dunendk repository depends on dunetpc please check the dunetpc version needed to pull down & compile dunendk, you can find that info at /ups/product_deps if a different version was setup an error message will appear
6. To the build directory in your development area,
   

   cd mydunendkAna/build*
   

   
   and then, if this is your first time building with these local packages, do
   

   mrbsetenv
   

   
   Once properly set up, you can build/install
   

   mrb i -j4
   

   
   After everything is built and installed, do
   

   mrb i -j4
   

   
   Now you can run the fcl file to do analysis using NDKAna_module.cc
   

   cd srcs/dunendk/fcl/
   lar -c NDKAnalysis.fcl myfile_reco.root