JRunAnalyzer Class Reference

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

#include <JRunAnalyzer.hh>

Public Member Functions
	JRunAnalyzer ()
	Default constructor. More...
	JRunAnalyzer (string file, string detectorFile, JLimit_t nTimeslices, JLimit_t nSummaryslices, JLimit_t nEvents, bool pmtanalysis)
	Constructor. More...
	~JRunAnalyzer ()
	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 JRunAnalyzer.hh.

Constructor & Destructor Documentation

JRunAnalyzer::JRunAnalyzer ( )

inline

Default constructor.

Definition at line 44 of file JRunAnalyzer.hh.

{}

JRunAnalyzer::JRunAnalyzer ( 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	Option for analysis of PMT data

Definition at line 56 of file JRunAnalyzer.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);
  } 

JRunAnalyzer::~JRunAnalyzer ( )

inline

Destructor.

Definition at line 80 of file JRunAnalyzer.hh.

{}

Member Function Documentation

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

inline

Definition at line 88 of file JRunAnalyzer.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_Triggered_over_Snapshot_hits -> Fill((Double_t) event.size<JDAQTriggeredHit> () / event.size<JDAQSnapshotHit > ());
      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);
        }
      }

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

        JDAQTriggerMask trigger_mask(event.getTriggerMask(*hit));

        if(trigger_mask.hasTriggerBit(getTriggerBit<JTrigger3DMuon>())) counter_3dmuon++; 
         
        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);

          if(trigger_mask.hasTriggerBit(getTriggerBit<JTrigger3DMuon>())){
            histograms.h_trigger.h_Triggered_hits_3dmuon_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());
        }
      }

      if(counter_3dmuon != 0) histograms.h_trigger.h_Triggered_hits_3dmuon->Fill(counter_3dmuon);

      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 JRunAnalyzer::iterateSummarysliceTree ( JTreeScanner< JDAQSummaryslice > &  scanner )

inline

Definition at line 179 of file JRunAnalyzer.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, nHRVcount);
          }
        
          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;    
            const double factor = 1.0e-3;   // [kHz]
              
            for (int i = 0 ; i < NUMBER_OF_PMTS ; i++){

              rate += frame->getRate(i, factor);

              h2->Fill(i , frame->getRate(i));

              PMT_rate_quantiles[string][floor-1][i].put(frame->getRate(i, factor)); 

              histograms.h_summary.h_pmt_rate_distribution->Fill(frame->getRate(i, factor), frame->getWeight(i, factor));

            }

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

            double rate = 0;    
            const double factor = 1.0e-3;   // [kHz]
                  
            for (int i = 0 ; i < NUMBER_OF_PMTS ; i++){

              rate += frame->getRate(i, factor);

              histograms.h_summary.h_pmt_rate_distribution->Fill(frame->getRate(i, factor), frame->getWeight(i, factor));
            }

            histograms.h_summary.h_dom_rate_distribution->Fill(rate);

            DOM_rate_quantiles[string].resize(modulesPerString);
          
            DOM_rate_quantiles[string][floor-1].put(rate);
            
          }
          
        }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 JRunAnalyzer::iterateTimesliceTree ( JTreeScanner< T > &  scanner )

inline

Definition at line 300 of file JRunAnalyzer.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 JRunAnalyzer::readSummaryData ( )

inline

Definition at line 366 of file JRunAnalyzer.hh.

                         {

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

template<class T >

void JRunAnalyzer::readTimesliceData ( )

inline

Definition at line 383 of file JRunAnalyzer.hh.

                           {

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

void JRunAnalyzer::readEvents ( )

inline

Definition at line 400 of file JRunAnalyzer.hh.

                    {

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

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

JRA_Histograms JRunAnalyzer::getHistograms ( )

inline

Definition at line 415 of file JRunAnalyzer.hh.

                                {

    return histograms ;
  }

Member Data Documentation

string JRunAnalyzer::inputFile

private

Definition at line 29 of file JRunAnalyzer.hh.

JDetector JRunAnalyzer::detector

private

Definition at line 30 of file JRunAnalyzer.hh.

JModuleRouter* JRunAnalyzer::router

private

Definition at line 31 of file JRunAnalyzer.hh.

JRA_Histograms JRunAnalyzer::histograms

private

Definition at line 32 of file JRunAnalyzer.hh.

int JRunAnalyzer::modulesPerString

private

Definition at line 33 of file JRunAnalyzer.hh.

JLimit_t JRunAnalyzer::numberOfTimeslices

private

Definition at line 34 of file JRunAnalyzer.hh.

JLimit_t JRunAnalyzer::numberOfSummaryslices

private

Definition at line 35 of file JRunAnalyzer.hh.

JLimit_t JRunAnalyzer::numberOfEvents

private

Definition at line 36 of file JRunAnalyzer.hh.

bool JRunAnalyzer::pmt_analysis

private

Definition at line 37 of file JRunAnalyzer.hh.

---

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

• JRunAnalyzer.hh