JRA_Histograms Class Reference

#include <JRunHistograms.hh>

Public Member Functions
	JRA_Histograms ()
	JRA_Histograms (const JDetector &detector)
void	initialize_summary_histograms ()
template<class T >
void	initialize_timeslice_histograms ()
void	initialize_trigger_histograms ()
template<class T >
void	Write_histogram_table_to_file (TFile &f, string dirname, vector< vector< T * > > table)
template<class T >
void	Write_histogram_table_to_file (TFile &f, string dirname, vector< T * > table)
template<class T , class V >
void	Write_manager_to_file (TFile &f, string dirname, JManager< T, V > *table)
template<class T , class V >
void	Replace_wildcard_in_name (JManager< T, V > *manager, char wc= '%')
template<class T , class V >
void	Write_manager_in_key_dir (TFile &f, JManager< T, V > *manager)
template<class T , class V >
void	Write_manager_table_in_key_dir (TFile &f, vector< JManager< T, V > * > table)
void	Write_to_file (TFile &f)

Public Attributes
SummaryHistograms	h_summary
TimesliceHistograms	h_timeslice
TriggerHistograms	h_trigger
std::set< int >	du_ids
int	modules_per_string

Detailed Description

Definition at line 740 of file JRunHistograms.hh.

Constructor & Destructor Documentation

JRA_Histograms::JRA_Histograms ( )

inline

Definition at line 750 of file JRunHistograms.hh.

{}

JRA_Histograms::JRA_Histograms ( const JDetector &  detector )

inline

Definition at line 752 of file JRunHistograms.hh.

                                            {

    du_ids             = getStringIDs(detector);
    modules_per_string = JDETECTOR::getNumberOfModules (detector) / getNumberOfStrings (detector);
    h_summary          = SummaryHistograms();
    h_timeslice        = TimesliceHistograms();
    h_trigger          = TriggerHistograms ();
  }

Member Function Documentation

void JRA_Histograms::initialize_summary_histograms ( )

inline

Definition at line 766 of file JRunHistograms.hh.

                                       {
    
    h_summary.initialize(du_ids , modules_per_string);

  }

template<class T >

void JRA_Histograms::initialize_timeslice_histograms ( )

inline

Definition at line 778 of file JRunHistograms.hh.

                                         {
    
    const int       index = JIndexOf<JDAQTimesliceTypes_t , T>::value ;
    const string   prefix = "KM3NETDAQ::JDAQ" ;
    string        ts_name = T::Class_Name();
    string::size_type pos = ts_name.find(prefix);

    if (pos != string::npos) ts_name.replace(ts_name.find(prefix) , prefix.length() , "");
      
    h_timeslice.initialize(du_ids , modules_per_string , index , ts_name);
  }

void JRA_Histograms::initialize_trigger_histograms ( )

inline

Definition at line 793 of file JRunHistograms.hh.

                                       {
    
    h_trigger.initialize (du_ids , modules_per_string);
  }

template<class T >

void JRA_Histograms::Write_histogram_table_to_file ( TFile &  f, string  dirname, vector< vector< T * > >  table  )

inline

Definition at line 806 of file JRunHistograms.hh.

                                                                                                 {

    if(f.GetDirectory(dirname.c_str()) == 0) f.mkdir (dirname.c_str());

    f.cd (dirname.c_str());

    for (int i=0 ; i < (int)table.size(); i++){

      for (int j=0 ; j< (int)table[i].size(); j++){

        if (table[i][j]) table [i][j] -> Write();
      }
    }
  }

template<class T >

void JRA_Histograms::Write_histogram_table_to_file ( TFile &  f, string  dirname, vector< T * >  table  )

inline

Definition at line 829 of file JRunHistograms.hh.

                                                                                      {

    if(f.GetDirectory(dirname.c_str()) == 0) f.mkdir (dirname.c_str());

    f.cd (dirname.c_str());

    for (int i=0 ; i < (int)table.size(); i++){

      if (table[i]) table[i] -> Write();
    }
  }

template<class T , class V >

void JRA_Histograms::Write_manager_to_file ( TFile &  f, string  dirname, JManager< T, V > *  table  )

inline

Definition at line 849 of file JRunHistograms.hh.

                                                                                    {

    if(f.GetDirectory(dirname.c_str()) == 0) f.mkdir (dirname.c_str());

    f.cd (dirname.c_str());
    
    for (typename JManager < T , V >::const_iterator i = table -> begin() ; i != table -> end() ; ++i){
      
      i -> second -> Write();
    }
  }

template<class T , class V >

void JRA_Histograms::Replace_wildcard_in_name ( JManager< T, V > *  manager, char  wc = '%'  )

inline

Definition at line 867 of file JRunHistograms.hh.

                                                                            {
    
    for (typename JManager < T , V >::const_iterator i = manager -> begin() ; i != manager -> end() ; ++i){

      if (i -> second -> GetTitle()){

        std::string     buffer = i -> second -> GetTitle();
        string::size_type ipos = buffer.find(wc);
        
        if (ipos!=std::string::npos){
          
          ostringstream os;
          
          os << i -> first ;
          
          buffer.replace(ipos, 1, os.str());
          
          i -> second -> SetTitle(buffer.c_str());
        }
      }
    }
  }

template<class T , class V >

void JRA_Histograms::Write_manager_in_key_dir ( TFile &  f, JManager< T, V > *  manager  )

inline

Definition at line 896 of file JRunHistograms.hh.

                                                                      {
    
    for (typename JManager < T , V >::const_iterator i = manager -> begin() ; i != manager -> end() ; ++i){
      
      std::string fullpath = MAKE_STRING(i->second->GetName());

      int          pos = fullpath.rfind ('/');
      std::string name = fullpath.substr (pos + 1);
      std::string path = fullpath.substr (0 , pos);

      if (f.GetDirectory(path.c_str()) == 0) f.mkdir (path.c_str());

      f.cd(path.c_str());

      i -> second -> SetName(name.c_str());
      i -> second -> Write();
    }
  }

template<class T , class V >

void JRA_Histograms::Write_manager_table_in_key_dir ( TFile &  f, vector< JManager< T, V > * >  table  )

inline

Definition at line 921 of file JRunHistograms.hh.

                                                                                        {

    for (typename vector < JManager < T , V >* >::const_iterator i = table.begin() ; i != table.end() ; ++i){

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

          std::string fullpath = MAKE_STRING(j-> second -> GetName());

          int          pos = fullpath.rfind  ('/');
          std::string name = fullpath.substr (pos + 1);
          std::string path = fullpath.substr (0 , pos);

          if (f.GetDirectory(path.c_str()) == 0) f.mkdir (path.c_str());

          f.cd(path.c_str());

          j -> second -> SetName(name.c_str());
          j -> second -> Write();
        }
      }
    }
  }

void JRA_Histograms::Write_to_file ( TFile &  f )

inline

Definition at line 950 of file JRunHistograms.hh.

                               {

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

    if (h_summary.h_fifo_per_dom)          h_summary.h_fifo_per_dom          -> Write();
    if (h_summary.h_daq_status_per_dom)    h_summary.h_daq_status_per_dom    -> Write();
    if (h_summary.h_hrv_per_dom)           h_summary.h_hrv_per_dom           -> Write();
    if (h_summary.h_rate_summary)          h_summary.h_rate_summary          -> Write();
    if (h_summary.h_pmt_rate_distribution) h_summary.h_pmt_rate_distribution -> Write();
    if (h_summary.h_dom_rate_distribution) h_summary.h_dom_rate_distribution -> Write();
  
    Write_histogram_table_to_file(f , MAKE_STRING("Detector/DOM_mean_rates")          , h_timeslice.h_dom_mean_rates);

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

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

          for(int k=1 ; k < j -> second -> GetNbinsX() ; ++k){

            for(int t=1 ; t < j -> second -> GetNbinsY() ; ++t){

              double width = j -> second -> GetYaxis() -> GetBinWidth(t);
              
              j -> second -> SetBinContent(k,t,j -> second -> GetBinContent(k,t) / width);
              j -> second -> SetBinError(k,t,j -> second -> GetBinError(k,t) / width);
            }
          } 
        }      
      }
    }

    Write_manager_table_in_key_dir (f , h_timeslice.m_pmt_tot_distributions);
    Write_manager_table_in_key_dir (f , h_timeslice.m_pmt_rate_distributions);
    if (h_summary.m_summary_rate_distribution) Write_manager_in_key_dir (f , h_summary.m_summary_rate_distribution);

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

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

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

    if (h_trigger.h_Trigger_bit_event) h_trigger.h_Trigger_bit_event -> Write(); 
    if (h_trigger.h_Trigger_bit_hit)   h_trigger.h_Trigger_bit_hit   -> Write(); 
    if (h_trigger.h_Trigger_time)      h_trigger.h_Trigger_time      -> Write();

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

    if (h_trigger.h_Triggered_hits)      h_trigger.h_Triggered_hits     -> Write(); 
    if (h_trigger.h_Snapshot_hits)       h_trigger.h_Snapshot_hits      -> Write();
    if (h_trigger.h_Number_of_overlays)  h_trigger.h_Number_of_overlays -> Write();

    f.mkdir ( MAKE_STRING ("JDAQEvent/hits_pmt_distributions").c_str());
    f.cd    ("JDAQEvent/hits_pmt_distributions");
    
    if (h_trigger.h_pmt_distribution_triggered_hits) {h_trigger.h_pmt_distribution_triggered_hits -> Write();}
    if (h_trigger.h_pmt_distribution_snapshot_hits)  {h_trigger.h_pmt_distribution_snapshot_hits  -> Write();}

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

    if (h_trigger.h_tot_distribution_triggered_hits) {normalizeBW1D (h_trigger.h_tot_distribution_triggered_hits) ; h_trigger.h_tot_distribution_triggered_hits-> Write();}
    if (h_trigger.h_tot_distribution_snapshot_hits)  {normalizeBW1D (h_trigger.h_tot_distribution_snapshot_hits)  ; h_trigger.h_tot_distribution_snapshot_hits -> Write();}

    f.cd ("JDAQEvent");
    
    if (h_trigger.h_Triggered_hits_per_module) { h_trigger.h_Triggered_hits_per_module-> Write();}
    if (h_trigger.h_Snapshot_hits_per_module)  { h_trigger.h_Snapshot_hits_per_module -> Write();}
    
    if (h_trigger.m_Snapshot_hits_per_pmt) Write_manager_in_key_dir(f , h_trigger.m_Snapshot_hits_per_pmt);
  }

Member Data Documentation

SummaryHistograms JRA_Histograms::h_summary

Definition at line 744 of file JRunHistograms.hh.

TimesliceHistograms JRA_Histograms::h_timeslice

Definition at line 745 of file JRunHistograms.hh.

TriggerHistograms JRA_Histograms::h_trigger

Definition at line 746 of file JRunHistograms.hh.

std::set<int> JRA_Histograms::du_ids

Definition at line 747 of file JRunHistograms.hh.

int JRA_Histograms::modules_per_string

Definition at line 748 of file JRunHistograms.hh.

---

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

• JRunHistograms.hh