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