led_pulse_analyzer Class Reference
[Analysis, Analysis]

This is a class to visualize reconstructed pulse information together with the raw waveform on a histogram. More...

#include <led_pulse_analyzer.hh>

Inheritance diagram for led_pulse_analyzer:


Public Member Functions
	led_pulse_analyzer ()
	Default constructor.
virtual	~led_pulse_analyzer ()
	Default destructor.
virtual bool	initialize ()
	Initialization method to be called before the analysis event loop.
virtual bool	analyze (storage_manager *storage)
	Analyze a data event-by-event.
virtual bool	finalize ()
	Finalize method to be called after all events processed.
void	clear_event ()
	Clear event.
void	set_pulse_type (DATA_STRUCT::DATA_TYPE type)
void	set_ref_channel (PMT::ch_number_t ch)
void	set_frame_size (PMT::word_t frame_size)
void	set_time_cut (PMT::word_t time_cut)
void	set_output_file (TFile *fout)
	A setter for analysis output file poitner.
void	set_verbosity (MSG::Level level)
	Setter for the verbosity level.
MSG::Level	get_verbosity () const
	Getter for the verbosity level.
const std::string	class_name () const
	Getter for the class name.
Protected Attributes
TFile *	_fout
	Analysis output file pointer.
char	_buf [200]
	char buffer for message manipulation
bool	_verbosity [MSG::MSG_TYPE_MAX]
	holder for enabled message levels
MSG::Level	_verbosity_level
	holder for specified verbosity level
std::string	_name
	class name holder
Private Attributes
PMT::word_t	_event_id
	event id holder
PMT::word_t	_ch_frame_id
PMT::word_t	_ch_sample_id
PMT::word_t	_disc_id
TTree *	_ch_tree
int	_ch
int	_nsample
double	_wf_charge
double	_wf_peak
double	_wf_ts
double	_wf_tm
double	_wf_te
double	_ped_mean
double	_ped_rms
led_pulse_selector	_selector

Detailed Description

This is a class to visualize reconstructed pulse information together with the raw waveform on a histogram.

This helps to debug/understand pulse reconstruction algorithm output. It is meant to be used interactively and not in a batch mode as it does not store any objects in an output. It allows a user to specify a range for various pulse-wise and event-wise parameters so that one can filter out irrelevant pulses from visualization routine (which can take time). See Analysis/mac/led_pulse_analyzer.py for an example usage.

Definition at line 45 of file led_pulse_analyzer.hh.

Constructor & Destructor Documentation

led_pulse_analyzer::led_pulse_analyzer ( )

Default constructor.

Definition at line 6 of file led_pulse_analyzer.cc.

References _ch_tree, and decoder_base::_name.

00006                                        : _selector() { 
00007   _name="led_pulse_analyzer"; 
00008   _ch_tree = 0;
00009 }

virtual led_pulse_analyzer::~led_pulse_analyzer ( ) [inline, virtual]

Default destructor.

Definition at line 53 of file led_pulse_analyzer.hh.

00053 {};

Member Function Documentation

bool led_pulse_analyzer::analyze ( storage_manager * storage ) [virtual]

Analyze a data event-by-event.

Implements ana_base.

Definition at line 35 of file led_pulse_analyzer.cc.

References _ch, _ch_frame_id, _ch_sample_id, _ch_tree, _disc_id, _ped_mean, _ped_rms, _selector, _wf_charge, _wf_peak, _wf_te, _wf_tm, _wf_ts, led_pulse_selector::analyze(), MSG::ERROR, storage_manager::get_data(), pulse_selector_base::pulse_type(), led_pulse_selector::select_pulse(), and Message::send().

00035                                                          {
00036   
00037   pulse_collection *pulses = (pulse_collection*)(storage->get_data(_selector.pulse_type()));
00038 
00039   if(pulses->size()==0) {
00040     Message::send(MSG::ERROR,__FUNCTION__,"No Pulse Found!");
00041     return false;
00042   }
00043 
00044   if(!_selector.analyze(storage)) return false;
00045 
00046   // Fill data for this event
00047   for(pulse_collection::const_iterator iter(pulses->begin());
00048       iter!=pulses->end();
00049       ++iter){
00050 
00051     if( !( _selector.select_pulse(&(*iter))) ) continue;
00052       
00053     _ch           = (*iter).channel_number();
00054     _ch_frame_id  = (*iter).frame_id();
00055     _ch_sample_id = (*iter).timeslice();
00056     _disc_id      = (int)((*iter).disc_id());
00057     _wf_charge    = (*iter).charge();
00058     _wf_peak      = (*iter).pulse_peak();
00059     _wf_ts        = (*iter).start_time();
00060     _wf_tm        = (*iter).max_time();
00061     _wf_te        = (*iter).end_time();
00062     _ped_mean     = (*iter).ped_mean();
00063     _ped_rms      = (*iter).ped_rms();
00064     _ch_tree->Fill();
00065     
00066   }
00067   
00068   return true;
00069 }

const std::string decoder_base::class_name ( ) const [inline, inherited]

Getter for the class name.

Definition at line 46 of file decoder_base.hh.

References decoder_base::_name.

00046 {return _name;};

void led_pulse_analyzer::clear_event ( )

Clear event.

bool led_pulse_analyzer::finalize ( ) [virtual]

Finalize method to be called after all events processed.

Implements ana_base.

Definition at line 71 of file led_pulse_analyzer.cc.

References _ch_tree, and ana_base::_fout.

00071                                   {
00072 
00073   // This function is called at the end of event loop.
00074   // Do all variable finalization you wish to do here.
00075   // If you need, you can store your ROOT class instance in the output
00076   // file. You have an access to the output file through "_fout" pointer.
00077   //amplitude->Write();
00078   _fout->cd();
00079   _ch_tree->Write();
00080   return true;
00081 }

MSG::Level decoder_base::get_verbosity ( ) const [inline, inherited]

Getter for the verbosity level.

Definition at line 43 of file decoder_base.hh.

References decoder_base::_verbosity_level.

00043 {return _verbosity_level;};

bool led_pulse_analyzer::initialize ( ) [virtual]

Initialization method to be called before the analysis event loop.

Implements ana_base.

Definition at line 11 of file led_pulse_analyzer.cc.

References _ch, _ch_frame_id, _ch_sample_id, _ch_tree, _disc_id, _event_id, _nsample, _ped_mean, _ped_rms, _selector, _wf_charge, _wf_peak, _wf_te, _wf_tm, _wf_ts, and led_pulse_selector::initialize().

00011                                     {
00012 
00013   if(!_ch_tree){
00014 
00015     _ch_tree = new TTree("ch_tree","");
00016     _ch_tree->Branch("event_id",&_event_id,"event_id/I");
00017     _ch_tree->Branch("frame_id",&_ch_frame_id,"frame_id/I");
00018     _ch_tree->Branch("sample_id",&_ch_sample_id,"sample_id/I");
00019     _ch_tree->Branch("disc_id",&_disc_id,"disc_id/I");
00020     _ch_tree->Branch("ch",&_ch,"ch/I");
00021     _ch_tree->Branch("charge",&_wf_charge,"charge/D");
00022     _ch_tree->Branch("peak",&_wf_peak,"peak/D");
00023     _ch_tree->Branch("t_start",&_wf_ts,"t_start/D");
00024     _ch_tree->Branch("t_max",&_wf_tm,"t_max/D");
00025     _ch_tree->Branch("t_end",&_wf_te,"t_end/D");
00026     _ch_tree->Branch("ped_mean",&_ped_mean,"ped_mean/D");
00027     _ch_tree->Branch("ped_rms",&_ped_rms,"ped_rms/D");
00028     _ch_tree->Branch("nsample",&_nsample,"nsample/I");
00029 
00030   }
00031 
00032   return _selector.initialize();
00033 }

void led_pulse_analyzer::set_frame_size ( PMT::word_t frame_size ) [inline]

Definition at line 71 of file led_pulse_analyzer.hh.

References _selector, and led_pulse_selector::set_frame_size().

00071 {_selector.set_frame_size(frame_size);};

void ana_base::set_output_file ( TFile * fout ) [inline, inherited]

A setter for analysis output file poitner.

Definition at line 46 of file ana_base.hh.

References ana_base::_fout.

00046 {_fout=fout;};

void led_pulse_analyzer::set_pulse_type ( DATA_STRUCT::DATA_TYPE type ) [inline]

Definition at line 67 of file led_pulse_analyzer.hh.

References _selector, and pulse_selector_base::set_pulse_type().

00067 {_selector.set_pulse_type(type);};

void led_pulse_analyzer::set_ref_channel ( PMT::ch_number_t ch ) [inline]

Definition at line 69 of file led_pulse_analyzer.hh.

References _selector, and led_pulse_selector::set_ref_channel().

00069 {_selector.set_ref_channel(ch);};

void led_pulse_analyzer::set_time_cut ( PMT::word_t time_cut ) [inline]

Definition at line 73 of file led_pulse_analyzer.hh.

References _selector, and led_pulse_selector::set_time_cut().

00073 {_selector.set_time_cut(time_cut);};

void decoder_base::set_verbosity ( MSG::Level level ) [inherited]

Setter for the verbosity level.

Reimplemented in ana_processor.

Definition at line 11 of file decoder_base.cc.

References decoder_base::_verbosity, decoder_base::_verbosity_level, MSG::DEBUG, MSG::ERROR, MSG::INFO, MSG::MSG_TYPE_MAX, MSG::NORMAL, and MSG::WARNING.

Referenced by algo_slow_readout_decoder::algo_slow_readout_decoder(), compare_content(), decoder_base::decoder_base(), decoder_manager::decoder_manager(), decoder_manager::initialize(), main(), xmit_event_search::run(), bin_word_search::run(), and ana_processor::set_verbosity().

00011                                               {
00012 
00013   _verbosity_level=level;
00014 
00015   for(size_t i=(size_t)(MSG::DEBUG); i<(size_t)(MSG::MSG_TYPE_MAX); ++i)
00016     _verbosity[i]=false;
00017 
00018   switch(level){
00019   case MSG::DEBUG:
00020     _verbosity[MSG::DEBUG]=true;
00021   case MSG::INFO:
00022     _verbosity[MSG::INFO]=true;
00023   case MSG::NORMAL:
00024     _verbosity[MSG::NORMAL]=true;
00025   case MSG::WARNING:
00026     _verbosity[MSG::WARNING]=true;
00027   case MSG::ERROR:
00028     _verbosity[MSG::ERROR]=true;
00029   case MSG::MSG_TYPE_MAX:
00030     break;
00031   }
00032 
00033 }