Detector.hh File Reference

#include "JDetector/JDetector.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JTools/JQuantile.hh"
#include <tuple>

Go to the source code of this file.

Functions
JDetector	loadDetector (string detectorFile)
	Loads the content of a detector file in a JDetector object. More...
vector< double >	getDetFile_t0_differences (int strNr, JDetector detector, int ref_pmt, int tgt_pmt)
	Loops over the floors in the string chosen by the user. More...
vector< tuple< double, double, double > >	get_t0_offsets_vs_depth (int strNr, const JDetector &detector)
double	get_PMT_distance (JPMT src_pmt, JPMT tgt_pmt)
	Calculate the distance between two JPMTs in the detector. More...
void	Print_Detector (JDetector &detector, int strNr)
	Prints the t0 of all the PMTs in a DU. More...

Function Documentation

loadDetector()

JDetector loadDetector ( string  detectorFile )

inline

Loads the content of a detector file in a JDetector object.

Author

rgruiz

Parameters

detectorFile

The name of a .detx

Returns

a JDetector object with the information of the file name

Definition at line 28 of file Detector.hh.

                                                  {
 
  JDetector detector;
 
  try {
 
    load(detectorFile, detector);
 
  }
 
  catch(const JException& error) {
 
    cerr << "FATAL ERROR. Could not open detector file '" << detectorFile << "'." << endl ;
 
    exit(1) ;
 
  }
 
  return detector ;
 
}

getDetFile_t0_differences()

vector<double> getDetFile_t0_differences ( int  strNr, JDetector  detector, int  ref_pmt, int  tgt_pmt  )

inline

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.

Parameters

detector	The name of a .detx
strNr	String number
ref_pmt	number of reference pmt
tgt_pmt	number of target pmt

Returns

a JDetector object with the information of the file name

Definition at line 66 of file Detector.hh.

                                                                                                           {
 
  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 ;
 
}

get_t0_offsets_vs_depth()

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

inline

Definition at line 83 of file Detector.hh.

                                                                                                                 {
 
  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 ;
 
}

get_PMT_distance()

double get_PMT_distance ( JPMT  src_pmt, JPMT  tgt_pmt  )

inline

Calculate the distance between two JPMTs in the detector.

Parameters

src_pmt	One pmt
tgt_pmt	Other pmt

Returns

distance between both JPMTs

Definition at line 133 of file Detector.hh.

                                                           {
    
  const JPosition3D pmt1_pos = src_pmt.getPosition() ;
  
  const JPosition3D pmt2_pos = tgt_pmt.getPosition() ;
  
  const JVector3D dist(pmt1_pos - pmt2_pos) ;
  
  return dist.getLength() ;
    
}

Print_Detector()

void Print_Detector ( JDetector &  detector, int  strNr  )

inline

Prints the t0 of all the PMTs in a DU.

Parameters

detector	A JDetector
strNr	The DU number

Definition at line 152 of file Detector.hh.

                                                           {
 
  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;
 
    } 
 
  }
 
}