Detector.hh

Go to the documentation of this file.

#ifndef __DETECTOR__
#define __DETECTOR__


#include "JDetector/JDetector.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JTools/JQuantile.hh"

#include <tuple>

using namespace JDETECTOR ;
using namespace std;


/**
 * \author rgruiz
 */



/**
 * Loads the content of a detector file in a JDetector object
 *
 * \param detectorFile The name of a .detx
 * \return a JDetector object with the information of the file name
 */
inline JDetector loadDetector(string detectorFile){

  JDetector detector;

  try {

    load(detectorFile, detector);

  }

  catch(const JException& error) {

    cerr << "FATAL ERROR. Could not open detector file '" << detectorFile << "'." << endl ;

    exit(1) ;

  }

  return detector ;

}



/**
 * Loops over the floors in the string chosen by the user. At each 
 * iteration "i" this function calculates:
 * dt = t0_tgt[i+1] - t0_ref[i]
 * where t0_tgt is the t0 of the target PMT and t0_ref 
 * is the t0 of the reference PMT (the closest to the nanobeacon). 
 * The result is stored in a vector of dimension nFloors - 1 and this vector is returned.
 * 
 * \param detector The name of a .detx
 * \param strNr String number 
 * \param ref_pmt number of reference pmt
 * \param tgt_pmt number of target pmt
 * \return a JDetector object with the information of the file name
 */
inline vector<double> getDetFile_t0_differences(int strNr , JDetector detector , int ref_pmt , int tgt_pmt){

  int nFloors = getNumberOfFloors(detector) ;

  vector<double> T0_diff (nFloors - 1) ;

  for (int i=1 ; i< nFloors ; i++) {

    T0_diff[i-1]  = detector.getModule(JModuleLocation(strNr , i+1 )).getPMT(tgt_pmt).getT0() - detector.getModule(JModuleLocation(strNr , i )).getPMT(ref_pmt).getT0() ;

  }

  return T0_diff ;

}


inline vector < tuple <double , double , double> > get_t0_offsets_vs_depth(int strNr , const JDetector& detector){

  int nFloors = getNumberOfFloors(detector) ;

  vector< tuple <double , double , double> > t0_vs_depth (nFloors) ;

  JTOOLS::JQuantile Q ("") ;

  for (int i = 1 ; i<nFloors+1 ; i++){

    JModule module = detector.getModule(JModuleLocation(strNr , i)) ;

    double depth = module.getPosition().getZ() ;

    JTOOLS::JQuantile q ("") ;

    for (int j = 0 ; j < getNumberOfPMTs(module) ; j++){

      JPMT pm = module.getPMT(j) ;

      q.put(pm.getT0()) ;

    }

    tuple < double , double , double > t (depth , q.getMean() , q.getSTDev()) ;
    
    t0_vs_depth[i-1] = t ;

    Q.put(q.getMean()) ;

  }
  
  for (auto & t: t0_vs_depth){
    
    get<1>(t) -= Q.getMean() ; 

  }

  return t0_vs_depth ;

}


/**
 * Calculate the distance between two JPMTs in the detector
 *
 * \param src_pmt One pmt
 * \param tgt_pmt Other pmt
 * \return distance between both JPMTs
 */  
inline double get_PMT_distance(JPMT src_pmt , JPMT tgt_pmt){
    
  const JPosition3D pmt1_pos = src_pmt.getPosition() ;
  
  const JPosition3D pmt2_pos = tgt_pmt.getPosition() ;
  
  const JVector3D dist(pmt1_pos - pmt2_pos) ;
  
  return dist.getLength() ;
    
}


/**
 * Prints the t0 of all the PMTs in a DU
 *
 * \param detector A JDetector
 * \param strNr The DU number
 */
inline void Print_Detector(JDetector& detector , int strNr){

  int nFloors = 18 ;

  int nPMTs = 31 ;

  JModuleRouter moduleRouter(detector) ;

  for(int i=1 ; i < nFloors ; i++){

    JModule dom = detector.getModule(JModuleLocation(strNr , i+1)) ;

    for(int j=0 ; j < nPMTs ; j++){

      cout <<dom.getPMT(j).getT0() << endl;

    } 

  }
 
}

#endif