pmtbaseline_ana.hh

Go to the documentation of this file.

#ifndef PMTBASELINE_ANA_HH
#define PMTBASELINE_ANA_HH

#include "ana_base.hh"
#include <map>
#include <TH1D.h>

class pmtbaseline_ana : public ana_base{

public:

  pmtbaseline_ana();

  virtual ~pmtbaseline_ana(){};

  virtual bool initialize();

  virtual bool analyze(storage_manager* storage);

  virtual bool finalize();

  void set_range_pulse_start_time(double v1, double v2)   {_cut_tstart.first=v1;   _cut_tstart.second=v2;   };

  void set_range_pulse_start_time_reco(double v1, double v2)   {_cut_tstart_reco.first=v1;   _cut_tstart_reco.second=v2;   };

  void set_range_pulse_end_time(double v1, double v2)     {_cut_tend.first=v1;     _cut_tend.second=v2;     };

  void set_range_pulse_amp(double v1, double v2)          {_cut_amp.first=v1;      _cut_amp.second=v2;      };

  void set_range_pulse_charge(double v1,double v2)        {_cut_charge.first=v1;   _cut_charge.second=v2;   };

  void set_range_ped_mean(double v1, double v2)           {_cut_pedbase.first=v1;  _cut_pedbase.second=v2;  };

  void set_range_ped_rms(double v1, double v2)            {_cut_pedrms.first=v1;   _cut_pedrms.second=v2;   };

  void set_range_event_id(PMT::word_t v1, PMT::word_t v2) {_cut_event_id.first=v1; _cut_event_id.second=v2; };

  void set_range_channel(PMT::ch_number_t v1, PMT::ch_number_t v2) {_cut_channels.first=v1; _cut_channels.second=v2;};

  void set_range_sum_charge(double v1, double v2)         {_cut_sum_charge.first=v1; _cut_sum_charge.second=v2;};

  void set_range_sum_peak(double v1, double v2)           {_cut_sum_peak.first=v1; _cut_sum_peak.second=v2;};

  void set_range_npulse(double v1, double v2)             {_cut_npulse.first=v1; _cut_npulse.second=v2;};

  //
  // Getters
  //

  const std::pair<double,double> range_pulse_start_time() const   {return _cut_tstart;};

  const std::pair<double,double> range_pulse_start_time_reco() const   {return _cut_tstart_reco;};

  const std::pair<double,double> range_pulse_end_time() const     {return _cut_tend;};

  const std::pair<double,double> range_pulse_amp() const          {return _cut_amp;};

  const std::pair<double,double> range_pulse_charge() const       {return _cut_charge;};

  const std::pair<double,double> range_ped_mean() const           {return _cut_pedbase;};

  const std::pair<double,double> range_ped_rms() const            {return _cut_pedrms;};

  const std::pair<PMT::word_t,PMT::word_t> range_event_id() const {return _cut_event_id;};

  const std::pair<PMT::word_t,PMT::word_t> range_channel() const  {return _cut_channels;};

  const std::pair<double,double> range_sum_charge() const         {return _cut_sum_charge;};

  const std::pair<double,double> range_sum_peak() const           {return _cut_sum_peak;};

  const std::pair<UInt_t,UInt_t> range_npulse() const         {return _cut_npulse;};

protected:

  void reset_cuts();

  std::map<PMT::ch_number_t,TH1D*> _hArray_PedMean, _hArray_PedRMS, _hArray_Charge, _hArray_Peak;

  std::pair<double,double> _cut_tstart, _cut_tstart_reco, _cut_tend, _cut_amp, _cut_charge, _cut_pedbase, _cut_pedrms;

  std::pair<double,double> _cut_sum_charge, _cut_sum_peak;

  std::pair<UInt_t,UInt_t> _cut_npulse;  

  std::pair<PMT::word_t,PMT::word_t> _cut_event_id;

  std::pair<PMT::ch_number_t,PMT::ch_number_t> _cut_channels;

};

#endif
 // end of doxygen group 

---