JRunAnalyzer Class Reference

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

#include <JRunAnalyzer.hh>

Classes
struct	value_type

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 ()
void	writeToFile (TFile &f)

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

Detailed Description

Class dedicated to the analysis of KM3NeT runs.

Definition at line 29 of file JRunAnalyzer.hh.

Constructor & Destructor Documentation

JRunAnalyzer::JRunAnalyzer ( )

inline

Default constructor.

Definition at line 45 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 57 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);
    highest_floor  = getRangeOfFloors(detector).getUpperLimit();
  } 

JRunAnalyzer::~JRunAnalyzer ( )

inline

Destructor.

Definition at line 81 of file JRunAnalyzer.hh.

{}

Member Function Documentation

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

inline

Definition at line 89 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;
      double tmin = numeric_limits<double>::max();
      double tmax = 0;
      
      for (JDAQEvent::const_iterator<JDAQTriggeredHit> hit = event.begin<JDAQTriggeredHit>(); hit != event.end<JDAQTriggeredHit>(); ++hit) {

        if(hit->getT() < tmin) tmin = hit->getT();
        if(hit->getT() > tmax) tmax = hit->getT();
        
        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());
        }
      }

      histograms.h_trigger.h_event_duration->Fill(tmax - tmin);

      histograms.h_trigger.h_pmt_distribution_triggered_hits->Scale(1. / (Double_t) event.size<JDAQTriggeredHit> ());
     
      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 197 of file JRunAnalyzer.hh.

                                                                             {

    map<int, value_type>               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 ();

          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, (1./distance(scanner.begin(), scanner.end()))*100);   
          }
        
          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, factor));

              PMT_rate_quantiles[string][floor][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(highest_floor + 1);
          
            DOM_rate_quantiles[string][floor].put(rate);
            
          } 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(highest_floor + 1);
          
            DOM_rate_quantiles[string][floor].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 <= highest_floor; j++) {
      
        if (i->second[j].getCount() > 0) histograms.h_summary.h_rate_summary->Fill(i->first, j, i->second[j].getMean());
      }
    }

    if(pmt_analysis == true){

      for(map<int, value_type>::const_iterator i = PMT_rate_quantiles.begin(); i != PMT_rate_quantiles.end(); ++i) {

        for (int j = 1; j <= highest_floor; 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, 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 314 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 380 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 397 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 414 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 429 of file JRunAnalyzer.hh.

                                {

    return histograms ;
  }

void JRunAnalyzer::writeToFile ( TFile &  f )

inline

Definition at line 438 of file JRunAnalyzer.hh.

                             {

    f.mkdir("Detector");
    f   .cd("Detector");

    if (histograms.h_summary.h_fifo_per_dom)          histograms.h_summary.h_fifo_per_dom          -> Write();
    if (histograms.h_summary.h_daq_status_per_dom)    histograms.h_summary.h_daq_status_per_dom    -> Write();
    if (histograms.h_summary.h_hrv_per_dom)           histograms.h_summary.h_hrv_per_dom           -> Write();
    if (histograms.h_summary.h_rate_summary)          histograms.h_summary.h_rate_summary          -> Write();
    if (histograms.h_summary.h_pmt_rate_distribution) histograms.h_summary.h_pmt_rate_distribution -> Write();
    if (histograms.h_summary.h_dom_rate_distribution) histograms.h_summary.h_dom_rate_distribution -> Write();

    histograms.Write_histogram_table_to_file(f , MAKE_STRING("Detector"), histograms.h_timeslice.h_dom_mean_rates);

    for (typename vector < JManager < string , TH2D >* >::const_iterator i = histograms.h_timeslice.m_pmt_tot_distributions.begin();
         i != histograms.h_timeslice.m_pmt_tot_distributions.end() ; ++i){

      if ((*i)){

        for (typename JManager < string , TH2D >::const_iterator j = (*i) -> begin() ; j != (*i) -> end() ; ++j){

          j->second->Scale(1., "width");
        }
      }
    }

    histograms.Write_manager_table_in_key_dir (f , histograms.h_timeslice.m_pmt_tot_distributions);
    histograms.Write_manager_table_in_key_dir (f , histograms.h_timeslice.m_pmt_rate_distributions);

    if (histograms.h_summary.m_summary_rate_distribution) histograms.Write_manager_in_key_dir (f , histograms.h_summary.m_summary_rate_distribution);

    if (histograms.h_summary.m_mean_summary_rate)              histograms.Write_manager_in_key_dir( f , histograms.h_summary.m_mean_summary_rate);
    if (histograms.h_summary.m_mean_summary_rate_distribution) histograms.Write_manager_in_key_dir( f , histograms.h_summary.m_mean_summary_rate_distribution);

    histograms.h_timeslice.Fill_mean_ToT_histograms();

    histograms.Write_manager_table_in_key_dir (f , histograms.h_timeslice.m_mean_ToT);
    histograms.Write_manager_table_in_key_dir (f , histograms.h_timeslice.m_mean_ToT_distribution);

    f.mkdir ( MAKE_STRING ("JDAQEvent").c_str());
    f.cd    ("JDAQEvent");

    if (histograms.h_trigger.h_Trigger_bit_event) histograms.h_trigger.h_Trigger_bit_event -> Write();
    if (histograms.h_trigger.h_Trigger_bit_hit)   histograms.h_trigger.h_Trigger_bit_hit   -> Write();
    if (histograms.h_trigger.h_Triggered_hits)        histograms.h_trigger.h_Triggered_hits        -> Write();
    if (histograms.h_trigger.h_Triggered_hits_3dmuon) histograms.h_trigger.h_Triggered_hits_3dmuon -> Write();
    if (histograms.h_trigger.h_Triggered_hits_3dmuon_per_module) histograms.h_trigger.h_Triggered_hits_3dmuon_per_module -> Write();
    if (histograms.h_trigger.h_Snapshot_hits)         histograms.h_trigger.h_Snapshot_hits         -> Write();
    if (histograms.h_trigger.h_Triggered_over_Snapshot_hits)       histograms.h_trigger.h_Triggered_over_Snapshot_hits  -> Write();
    if (histograms.h_trigger.h_event_duration)  histograms.h_trigger.h_event_duration -> Write();
    if (histograms.h_trigger.h_Number_of_overlays)  histograms.h_trigger.h_Number_of_overlays -> Write();
    if (histograms.h_trigger.h_pmt_distribution_triggered_hits) {histograms.h_trigger.h_pmt_distribution_triggered_hits -> Write();}
    if (histograms.h_trigger.h_pmt_distribution_snapshot_hits)  {histograms.h_trigger.h_pmt_distribution_snapshot_hits  -> Write();}
    if (histograms.h_trigger.h_tot_distribution_triggered_hits) { histograms.h_trigger.h_tot_distribution_triggered_hits->Scale(1, "width") ; histograms.h_trigger.h_tot_distribution_triggered_hits-> Write();}
    if (histograms.h_trigger.h_tot_distribution_snapshot_hits)  {histograms.h_trigger.h_tot_distribution_snapshot_hits->Scale(1, "width") ; histograms.h_trigger.h_tot_distribution_snapshot_hits -> Write();}
    if (histograms.h_trigger.h_Triggered_hits_per_module) { histograms.h_trigger.h_Triggered_hits_per_module-> Write();}
    if (histograms.h_trigger.h_Snapshot_hits_per_module)  { histograms.h_trigger.h_Snapshot_hits_per_module -> Write();}
    if (histograms.h_trigger.m_Snapshot_hits_per_pmt) histograms.Write_manager_in_key_dir(f , histograms.h_trigger.m_Snapshot_hits_per_pmt);
  }

Member Data Documentation

string JRunAnalyzer::inputFile

private

Definition at line 31 of file JRunAnalyzer.hh.

JModuleRouter* JRunAnalyzer::router

private

Definition at line 32 of file JRunAnalyzer.hh.

JRA_Histograms JRunAnalyzer::histograms

private

Definition at line 33 of file JRunAnalyzer.hh.

int JRunAnalyzer::highest_floor

private

Definition at line 34 of file JRunAnalyzer.hh.

JLimit_t JRunAnalyzer::numberOfTimeslices

private

Definition at line 35 of file JRunAnalyzer.hh.

JLimit_t JRunAnalyzer::numberOfSummaryslices

private

Definition at line 36 of file JRunAnalyzer.hh.

JLimit_t JRunAnalyzer::numberOfEvents

private

Definition at line 37 of file JRunAnalyzer.hh.

bool JRunAnalyzer::pmt_analysis

private

Definition at line 38 of file JRunAnalyzer.hh.

---

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

• JRunAnalyzer.hh