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

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 ;

}

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 ;

}

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 ;

}

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() ;
    
}

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;

    } 

  }
 
}