RunAnalyzer Class Reference

Class dedicated to the analysis of KM3NeT runs. More...

#include <RunAnalyzer.hh>

Public Member Functions
	RunAnalyzer ()
	Default constructor. More...
	RunAnalyzer (string file, string detectorFile, JLimit_t nTimeslices, JLimit_t nSummaryslices, JLimit_t nEvents, bool pmtanalysis)
	Constructor. More...
	~RunAnalyzer ()
	Destructor. More...
void	iterateEventTree (JTreeScanner< JDAQEvent > &scanner)
void	iterateSummarysliceTree (JTreeScanner< JDAQSummaryslice > &scanner)
template<class T >
void	iterateTimesliceTree (JTreeScanner< T > &scanner)
void	readSummaryData ()
template<class T >
void	readTimesliceData ()
void	readEvents ()
JRA_Histograms	getHistograms ()

Private Attributes
string	inputFile
JDetector	detector
JModuleRouter *	router
JRA_Histograms	histograms
int	modulesPerString
JLimit_t	numberOfTimeslices
JLimit_t	numberOfSummaryslices
JLimit_t	numberOfEvents
bool	pmt_analysis

Detailed Description

Class dedicated to the analysis of KM3NeT runs.

Author

rgruiz, adomi

Definition at line 27 of file RunAnalyzer.hh.

Constructor & Destructor Documentation

RunAnalyzer::RunAnalyzer ( )

inline

Default constructor.

Definition at line 44 of file RunAnalyzer.hh.

{}

RunAnalyzer::RunAnalyzer ( string  file, string  detectorFile, JLimit_t  nTimeslices, JLimit_t  nSummaryslices, JLimit_t  nEvents, bool  pmtanalysis  )

inline

Constructor.

Parameters

file	Path to the file to be analyzed
detectorFile	path to a detector file
nTimeslices	Number of timeslices to be read
nSummaryslices	Number of summary slices to be read
nEvents	Number of frames to be read
pmtanalysis	PMT analysis option

Definition at line 56 of file RunAnalyzer.hh.

                                                                                                                                       :
    inputFile                       (file),
    router                          (NULL),
    numberOfTimeslices       (nTimeslices),
    numberOfSummaryslices (nSummaryslices),
    numberOfEvents               (nEvents),
    pmt_analysis             (pmtanalysis)
  {
    try {
      load (detectorFile, detector);
    }
    catch(const JException & error) {
      cerr << "FATAL ERROR. Could not open detector file '" << detectorFile << "'." << endl;
      exit(1);
    }
    
    histograms       = JRA_Histograms(detector);
    router           = new JModuleRouter(detector);
    modulesPerString = JDETECTOR::getNumberOfModules (detector) / getNumberOfStrings (detector);
  } 

RunAnalyzer::~RunAnalyzer ( )

inline

Destructor.

Definition at line 80 of file RunAnalyzer.hh.

{}

Member Function Documentation

void RunAnalyzer::iterateEventTree ( JTreeScanner< JDAQEvent > &  scanner )

inline

Definition at line 88 of file RunAnalyzer.hh.

                                                               {
  
    while (scanner.hasNext()) {

      JDAQEvent event = *(scanner.next());
  
      histograms.h_trigger.h_Snapshot_hits      -> Fill((Double_t) event.size<JDAQSnapshotHit > ());
      histograms.h_trigger.h_Triggered_hits     -> Fill((Double_t) event.size<JDAQTriggeredHit> ());
      histograms.h_trigger.h_Number_of_overlays -> Fill(event.getOverlays());

      for (unsigned int i = 0; i != NUMBER_OF_TRIGGER_BITS; ++i) {

        if (event.hasTriggerBit(i)) {

          histograms.h_trigger.h_Trigger_bit_event -> Fill((Double_t) i);
        }
      }

      for (JDAQEvent::const_iterator<JDAQTriggeredHit> hit = event.begin<JDAQTriggeredHit>(); hit != event.end<JDAQTriggeredHit>(); ++hit) {

        if (router->hasModule(hit->getModuleID())) {

          const JModule& module = router->getModule (hit->getModuleID());
        
          histograms.h_trigger.h_pmt_distribution_triggered_hits->Fill(hit->getPMT());
          histograms.h_trigger.h_tot_distribution_triggered_hits->Fill(hit->getToT());

          int String = module.getString();
          int Floor  = module.getFloor();

          histograms.h_trigger.h_Triggered_hits_per_module -> Fill(String , Floor);
        
          for (unsigned int i = 0; i != NUMBER_OF_TRIGGER_BITS; ++i) {
          
            if (hit -> hasTriggerBit(i)) {

              histograms.h_trigger.h_Trigger_bit_hit->Fill((Double_t) i);
            }
          }
        }else{
          FATAL("JModuleRouter trying to access non existing identifier: "<< hit->getModuleID());
        }
      }
      
      for (JDAQEvent::const_iterator<JDAQSnapshotHit> hit = event.begin<JDAQSnapshotHit>() ; hit != event.end<JDAQSnapshotHit>() ; ++hit){

        if (router->hasModule(hit->getModuleID())) {

          const JModule& module = router->getModule (hit->getModuleID());
        
          int String = module.getString();
          int Floor  = module.getFloor();
          int pmt    = hit->  getPMT();

          if(pmt_analysis == true){
            
            (*histograms.h_trigger.m_Snapshot_hits_per_pmt)[MAKE_STRING("Detector/DU" + to_string(String))] -> Fill(pmt, Floor);

          }

          histograms.h_trigger.h_pmt_distribution_snapshot_hits                                           -> Fill(hit->getPMT());       
          histograms.h_trigger.h_tot_distribution_snapshot_hits                                           -> Fill(hit->getToT());
          histograms.h_trigger.h_Snapshot_hits_per_module                                                 -> Fill(String, Floor);

          
          }else{
            FATAL("JModuleRouter trying to access non existing identifier: "<< hit->getModuleID());
        }
      }
    }
  }

void RunAnalyzer::iterateSummarysliceTree ( JTreeScanner< JDAQSummaryslice > &  scanner )

inline

Definition at line 166 of file RunAnalyzer.hh.

                                                                             {
    
    std::map <int,vector<vector <JQuantile>>> PMT_rate_quantiles;
    std::map <int,vector<JQuantile>>          DOM_rate_quantiles;

    while (scanner.hasNext()){

      JDAQSummaryslice slice = *(scanner.next());

      for (JDAQSummaryslice::const_iterator frame = slice.begin() ; frame != slice.end() ; ++ frame) {

        if (router->hasModule(frame->getModuleID())) {

          const JModule& module = router->getModule (frame->getModuleID());

          int string = module.getString();
          int floor  = module.getFloor ();
        
          PMT_rate_quantiles[string].resize(modulesPerString , vector <JQuantile>(NUMBER_OF_PMTS)); 

          JDAQFrameStatus status = frame -> getDAQFrameStatus();
          int nFIFOcount         = status.countFIFOStatus();
          int nHRVcount          = status.countHighRateVeto();

          if (nHRVcount > 0) {
            histograms.h_summary.h_hrv_per_dom->Fill(string , floor);
          }
        
          if (status.testDAQStatus() == false) {
            histograms.h_summary.h_daq_status_per_dom->Fill(string , floor);      
          }
        
          histograms.h_summary.h_fifo_per_dom->Fill(string , floor , nFIFOcount); 

          if(pmt_analysis == true){
          
            TH2D*  h2   = (*histograms.h_summary.m_summary_rate_distribution)[MAKE_STRING("Detector/DU" + to_string(string) + "/F" + to_string(floor))];

            double rate = 0;    

            for (int i = 0 ; i < NUMBER_OF_PMTS ; i++){

              rate += frame->getRate(i);

              h2->Fill(i , frame->getRate(i) * 1e-3);

              PMT_rate_quantiles[string][floor-1][i].put(frame->getRate(i)*1e-3); 

            }

            DOM_rate_quantiles[string].resize(modulesPerString);
          
            DOM_rate_quantiles[string][floor-1].put(rate * 1e-3);
            
          } else {

            double rate = 0;    
                  
            for (int i = 0 ; i < NUMBER_OF_PMTS ; i++){
              rate += frame->getRate(i);
              histograms.h_summary.h_pmt_rate_distribution->Fill(frame->getRate(i)*1e-3);
            }
            
            histograms.h_summary.h_dom_rate_distribution->Fill(rate*1e-3);

            DOM_rate_quantiles[string].resize(modulesPerString);
          
            DOM_rate_quantiles[string][floor-1].put(rate * 1e-3);
            
          }
          
        }else{
          FATAL("JModuleRouter trying to access non existing identifier: "<< frame->getModuleID());
        }      
      }
    }
    
  
    for (std::map<int , vector<JQuantile> >::const_iterator i = DOM_rate_quantiles.begin() ; i!= DOM_rate_quantiles.end() ; ++i) {
      
      for (int j=1 ; j < modulesPerString + 1 ; j++) {
      
        if (i->second[j-1].getCount() > 0) histograms.h_summary.h_rate_summary->Fill(i->first, j, i->second[j-1].getMean());
      }
    }

    
    if(pmt_analysis == true){
    
      for (std::map<int , vector < vector<JQuantile> > >::const_iterator i = PMT_rate_quantiles.begin() ; i!= PMT_rate_quantiles.end() ; ++i) {
    
        for (int j=0 ; j < modulesPerString ; j++){
      
          for (int k=0 ; k < NUMBER_OF_PMTS ; k++){

            if (i -> second[j][k].getCount() > 0){

              (*histograms.h_summary.m_mean_summary_rate)             [MAKE_STRING("Detector/DU" + to_string(i->first))]->Fill(k, j+1, i->second[j][k].getMean());
              (*histograms.h_summary.m_mean_summary_rate_distribution)[MAKE_STRING("Detector/DU" + to_string(i->first))]->Fill(        i->second[j][k].getMean());
            }
          }
        }
      }
    }
  }

template<class T >

void RunAnalyzer::iterateTimesliceTree ( JTreeScanner< T > &  scanner )

inline

Definition at line 279 of file RunAnalyzer.hh.

                                                           {
    
    int ts_type = JIndexOf<JDAQTimesliceTypes_t , T>::value;

    std::map < int , std::map <int , JQuantile> > DOM_rate_quantiles;

    double inverseFrameTimeSec = 1. / (1.0e-9 * getFrameTime());

    while (scanner.hasNext()){
      
      T slice = *(scanner.next());

      for(auto & frame : slice) {
        if (router->hasModule(frame.getModuleID())) {

          const JModule& module = router->getModule (frame.getModuleID());
          double  rate   = frame.numberOfHits * inverseFrameTimeSec;
          int     string = module.getString();
          int     floor  = module.getFloor ();

          DOM_rate_quantiles[string][floor].put(rate);

          vector <int>  pmt_hit_count (NUMBER_OF_PMTS , 0) ;

          if(pmt_analysis == true){
          
            TH2D* h2 = (*histograms.h_timeslice.m_pmt_tot_distributions[ts_type])[MAKE_STRING("Detector/DU" + to_string(string) + "/F" + to_string(floor) + "/" + to_string(frame.getModuleID()))];
            
            for (JDAQSuperFrame::const_iterator hit = frame.begin() ; hit != frame.end() ; ++hit){
          
            h2 -> Fill(hit->getPMT() , hit->getToT()) ;
            
            pmt_hit_count[hit->getPMT()]++;
            
            }
            
            for (int pmt = 0 ; pmt != NUMBER_OF_PMTS ; ++pmt) {
            
              (*histograms.h_timeslice.m_pmt_rate_distributions[ts_type])[MAKE_STRING("Detector/DU" + to_string(string) + "/F" + to_string(floor))] -> Fill(pmt , 1e-3 * pmt_hit_count[pmt] * inverseFrameTimeSec);
                  
            }

          }
            
          }else{
            FATAL("JModuleRouter trying to access non existing identifier: "<< frame.getModuleID());
          }      
        }      
    }
      
    
    for (std::map< int , std::map<int,JQuantile>>::const_iterator i = DOM_rate_quantiles.begin() ; i != DOM_rate_quantiles.end() ; i++){

      for (std::map<int , JQuantile>::const_iterator j = i->second.begin() ; j != i->second.end() ; j++){

        if (j->second.getCount() > 0) histograms.h_timeslice.h_dom_mean_rates[ts_type] -> Fill (i->first , j->first , j->second.getMean() ) ;
      }
    }
  }

void RunAnalyzer::readSummaryData ( )

inline

Definition at line 345 of file RunAnalyzer.hh.

                         {

    JTreeScanner <JDAQSummaryslice> scanner(inputFile, numberOfSummaryslices);
    if (scanner.hasNext()) {
      histograms.initialize_summary_histograms();
      scanner.rewind();
      iterateSummarysliceTree(scanner);
    }
  }

template<class T >

void RunAnalyzer::readTimesliceData ( )

inline

Definition at line 362 of file RunAnalyzer.hh.

                           {

    JTreeScanner <T> scanner(inputFile, numberOfTimeslices);
    if(scanner.hasNext()) {
      
      histograms.initialize_timeslice_histograms <T>();
      scanner.rewind();
      iterateTimesliceTree(scanner);
    }
  }

void RunAnalyzer::readEvents ( )

inline

Definition at line 379 of file RunAnalyzer.hh.

                    {

    JTreeScanner <JDAQEvent> scanner(inputFile, numberOfEvents);
        
    if(scanner.hasNext()) {

      histograms.initialize_trigger_histograms();
      scanner.rewind();
      iterateEventTree(scanner);
    }
  }

JRA_Histograms RunAnalyzer::getHistograms ( )

inline

Definition at line 394 of file RunAnalyzer.hh.

                                {

    return histograms ;
  }

Member Data Documentation

string RunAnalyzer::inputFile

private

Definition at line 29 of file RunAnalyzer.hh.

JDetector RunAnalyzer::detector

private

Definition at line 30 of file RunAnalyzer.hh.

JModuleRouter* RunAnalyzer::router

private

Definition at line 31 of file RunAnalyzer.hh.

JRA_Histograms RunAnalyzer::histograms

private

Definition at line 32 of file RunAnalyzer.hh.

int RunAnalyzer::modulesPerString

private

Definition at line 33 of file RunAnalyzer.hh.

JLimit_t RunAnalyzer::numberOfTimeslices

private

Definition at line 34 of file RunAnalyzer.hh.

JLimit_t RunAnalyzer::numberOfSummaryslices

private

Definition at line 35 of file RunAnalyzer.hh.

JLimit_t RunAnalyzer::numberOfEvents

private

Definition at line 36 of file RunAnalyzer.hh.

bool RunAnalyzer::pmt_analysis

private

Definition at line 37 of file RunAnalyzer.hh.

---

The documentation for this class was generated from the following file:

• RunAnalyzer.hh