Cobraa Coincident-background reactor antineutrino analysis is an analysis for reactor discovery via the simulation and evaluation of coincident events. It can also optionally use reconstruction of coincident events.
It's a toolchain which handles full-detector simulation through to sensitivity analysis for the current NEO detector design options: 16m/22m with Gd-water or Gd-WbLS.
Adapted from the WATCHMAN Watchmakers package (Author: Marc Bergevin, LLNL), it constitutes a significant overhaul of the original Watchmakers code. Retains the convenience of the directory organisation and macro/job production but incorporates streamlining and a full analysis update.
To be used alongside the FRED (BONSAI) reconstruction or CoRe (BONSAI) pair reconstruction.
Analysis performs full evaluation of coincidences for both signal and background and optimises sensitivity to a reactor as a function of cuts in up to 7 dimensions:
- Prompt energy threshold
- Delayed energy threshold
- Fiducial volume
- Time between triggers dT
- Distance between triggers dR (FRED only)
- Fit quality (BONSAI timing goodness)
- Maximum prompt energy
To install:
$git clone https://github.com/ait-watchman/cobraa
$cd cobraa
$./configure
$source env_wm.sh
Basic use:
-
Create macros and simulation jobs for a small-ish number of events (16m/Gd-water/20% photocoverage by default)
cobraa -m -j -e 2500 --lightSimWater
-
Run jobs (locally by default)
source job/job*.sh
-
Run FRED or CoRe reconstruction (this stage is not incorporated into Cobraa)
-
Merge ntuple root files from reconstruction
cobraa -M
-
Map coincidences after cuts
cobraa --coincidences [--core]
-
Calculate rates and optimise signal significance
cobraa --sensitivity [--core]
See Wiki for more details