/FLOWER

soFtware for LOW-Energy Reconstruction

Primary LanguageC++

FLOWER: soFtware for LOW-Energy Reconstruction

An energy reconstruction script for low-energy events that works with (hk-)BONSAI and other low-energy reconstruction tools - it just requires a positional vertex.

Compiling

Make sure the following are sourced

  • WCSim ($WCSIMDIR set)

Then

export FLOWERDIR=/path/to/FLOWER
make

Running the WCSim/BONSAI script bundled here

Make sure the following are sourced

  • WCSim ($WCSIMDIR set)
  • BONSAI ($BONSAIDIR set)
  • FLOWER ($FLOWERDIR set)

Then

export PATH=$FLOWERDIR/rootflower:$PATH
rootflower -b -q flower_with_bonsai.C+'("/path/to/wcsim/file.root",1,"SuperK")'

Running in other code

  • Construct a class member, giving it a detectorname (it sets default values and determines the exact energy correction factors) and detector geometry
WCSimFlower(const char * detectorname, WCSimRootGeom * geom, int verbose);
  • Override default values of parameters with these methods
void SetDarkRate(double darkrate);
void SetNPMTs   (int npmts);
void SetNWorkingPMTs(int nworkingpmts);
void SetNeighbourDistance(double neighbourdistance);
void SetShortDuration(double shortduration);
void SetLongDuration(double longduration);
void SetTopBottomDistance(double hi, double lo);
  • Every event/trigger, send it a list of hit times, hit tube IDs, and reconstructed vertex position (x,y,z)
double GetEnergy(std::vector<float> times, std::vector<int> tubeIds, float * vertex);