(Anti)Neutrino Directionality Study

Calculating the average antineutrino direction in the PROSPECT detector of antineturinos produced in the HFIR core. Paper in progress.

Summary

The analysis code first reads analyzed PROSPECT data passed through the analysis framework PROSPECT2x and creates histograms to track the average displacement between the prompt and delayed events. The x and y positions are approximated as the midpoint of the segment in which the events occurred. The inoperative segments requires a more complex selection of events, which will be described in the paper. This also leads to a more complex error calculation which is outlined in the code. After the events are counted, the background subtraction is carried to get the IBD signal events which are then put through the error calculations and angle reconstructions.

Requirements

ROOT 5+
C++ compliler (I'm using g++ in the instructions)
Analyzed and calibrated PROSPECT data

How To Use (Debian)

# Clone this repository
git clone https://github.com/shayshunk/NeutrinoDirectionality

# Go into the repository
cd NeutrinoDirectionality

# Install dependencies (quick way below but not recommended)
# Go here to find the CERN recommended way: https://root.cern/install/
sudo snap install root-framework

# Edit the header NeutrinoDirectionality.h to point to your data files
# Line 26-29 contain my local paths but yours will be different

# Run the code
# I'm using g++ because it's straightforward but you can configure another compiler
# You can also run in macro mode but it'll be slower
g++ NeutrinoDirectionality.cc -o Directionality `root-config --cflags --glibs`
./Directionality

# Make the plots
# Just use macro mode, it's fast enough that there's no time lost
# Compiling changes the plot aspect ratio for some reason
root NeutrinoDirectionalityPlots.cc

The above is an explanation on how to run the two main programs that calculate the average antineutrino angles and then plot them. The other piece of code AveragePromptLocation.cc attempts to calculate the average prompt location in the detector by using a simulation with no missing segments. It requires a bit more tweaking to run and is less refined.

Customization Options

A few command line options are available for customizing the output of the code. The follwing variables in NeutrinoDirectionality.cc let you define how verbose you want the output of the values calculated to be.

DETECTOR_VERBOSITY 
LIVETIME_VERBOSITY 
IBDCOUNT_VERBOSITY 
MEAN_VERBOSITY 
ANGLES_STATISTICS
COVARIANCE_VERBOSITY

Each option is accessed by adding a flag after the executable while running in the terminal. Example: ./Directionality -M -C -S. The options are are as follows:

-D will set DETECTOR_VERBOSITY to true. It will print the detector configurations for each period of PROSPECT data. These are available on the PROSPECT docdb.
-L sets LIVETIME_VERBOSITY to true. It will print the total livetime for all Reactor On and Reactor Off events.
-I sets IBD_COUNT VERBOSITY to true. It will print the total and effective IBDs in each direction for each dataset. Effective IBDs are calculated through Poisson statistics.
-M sets MEAN_VERBOSITY to true. It will print the p components and respective errors that are calculated prior to the addition of systematics. These components are then used to calculate the tangents of the angles. More details in the paper.
-S sets ANGLE_STATISTICS to true. It will print only the statistical errors for all angles.
-C sets COVARIANCE_VERBOSITY to true. It will print out the semi-major axes of the 1-sigma ellipse that the plots show. Note that these are not the errors on the angles individually but due to covariance in the angle calculations, the 68% C.L. ellipse combining both angles is a different shape and size than just using the errors on the angles as the semi-major axes.

The other option is contained in Formatting.h. I added a few quick functions that return a certain formatting (bold/underline) or color for more aesthetically pleasing output. These only work on Linux terminals. If working on another platform or the output simply looks jumbled or unpleasant, turn off the special formatting on line 4 by setting it to 0.

#define COLORFUL_FORMATTING 1

Sample Output

This is the final output of NeutrinoDirectionality.cc with the current dataset and angle calculation methods. All verbosity flags above are set to 0 except for COLORFUL_FORMATTING.

This the final plot as of now, though it is to subject to change as the paper develops.

shayshunk / NeutrinoDirectionality