The scripts and tools to complete the salt phase signal extraction are contained on the KEG data area in /data/data003/dgrant/analysis . This directory contains the version 3.0 release of the tsigex code and the macros/rootfiles created for signal extraction (in sub-directories drg and parameterizer).
Whether one is working with MC files or MCPROD processed neutrino datafiles it is often the case that there are many hundreds of files containing the events one is interested in. To simplify, we combine the events that meet a specific set of cuts from each of these files into a single ROOT file for the data processing. There are a series of macros to perform this combining, located in /data/data003/dgrant/analysis/parametizer/scripts/raw_data, with the name combine_xx.C . These macros are designed to accept a list of files to combine (using a txt file) and user defined cuts. The resultant ROOT files are stored in /data/data003/dgrant/analysis/parametizer/roots/raw_gen. We note that the file /data/data003/dgrant/analysis/parametizer/jeanne_runlist.txt lists all neutrino candidate events for the 1st salt phase open dataset.
Once one has a ROOT file containing a set of events that are interesting for study, it is useful to next histogram these events. In the case where the events are neutrino MC events of a specific flavour, we wish to fit the histograms to specific parameterized functions (as listed in chapter 5 of my PhD thesis). Macros to perform this histogramming and fits are found in /data/data003/dgrant/analysis/parametizer/scripts/fits.
The results of the parameterized fit functions are used to generate a final set of PDFs histograms for use in the signal extraction code. These histograms are generated and added to a single root file using the macros in /data/data003/dgrant/analysis/parametizer/scripts/hist_gen. The resultant root file may be found in /data/data003/dgrant/analysis/parametizer/roots/hist_prod.
Each of the MC neutrino event datasets were adjusted according to the determined detector resolution functions. These datasets produced new event histograms which were fit with the parametric functions to give new systematic PDFs that in turn were used in the signal extraction. The scripts to perform this fit are located in /data/data003/dgrant/analysis/parametizer/scripts/systematics, and the resultant root files are found in /data/data003/dgrant/dgrant/analysis/parametizer/roots/systematics.
The neutrino ROOT datafile in /data/data003/dgrant/analysis/drg/roots called neutrino_data_dec19.root. This directory also contains the other ROOT files used in the signal extraction process (the titles of these tiles makes them fairly self-explainatory).
We note that the signal extraction is completed using c++ code which is compiled against the tsigex library. This code is stored in /data/data003/dgrant/analysis/drg/code. The code may be compiled with the makefiles stored in /data/data003/dgrant/analysis/drg/makefiles, and the resultant executables are written to /data/data003/dgrant/analysis/drg/executables. When performing signal extraction tsigex uses a card file, found in /data/data003/dgrant/analysis/drg/datfiles. We note that these card files contain the official cuts used in the sig-ex process and the PDFs.
The /data/data003/dgrant/analysis/drg/scripts directory contains the macros used to compare the data histograms to the sum of the extracted CC, NC, ES and External neutrons with the appropriate PDFs. We also note that this directory contains the CORSET and CORGEN programs, re-written from the fortran version into a c++ macro. These macros are used with the scripts in /data/data003/dgrant/analysis/drg/systematics_param to compute the systematic uncertainties due to parameterization of the histograms to form PDFs. This process is described in detail in chapter 6 of the thesis.