TGrDetector.cxx

Go to the documentation of this file.

//  TGrDetector.cxx.
#include <iostream>

#include <TCutG.h>
#include <TH2F.h>

#include "TGrDetector.h"
#include "TCalibration.h"
#include "TSubEvent.h"

using std::cout;
using std::endl;

TGrDetector::TGrDetector(const Char_t *name, UInt_t channels)
  : TDetector(name, channels), fCalib(channels)
{
  fAdcs.reserve(channels);
  fLeTdcs.reserve(channels);
  fCfTdcs.reserve(channels);
}

TGrDetector::TGrDetector(const TGrDetector &rhs) : TDetector(rhs)
{
  fAdcs   = rhs.fAdcs;
  fLeTdcs = rhs.fLeTdcs;
  fCfTdcs = rhs.fCfTdcs;
  fCalib  = rhs.fCalib;
}

TGrDetector::~TGrDetector()
{
}

TGrDetector & TGrDetector::operator=(const TGrDetector &rhs)
{
  if (&rhs != this) {
    TDetector::operator=(rhs);
    fAdcs   = rhs.fAdcs;
    fLeTdcs = rhs.fLeTdcs;
    fCfTdcs = rhs.fCfTdcs;
    fCalib  = rhs.fCalib;
  }
  return *this;
}

void TGrDetector::Clear(const Char_t *option)
{
  TDetector::Clear();
  fAdcs.clear();
  fLeTdcs.clear();
  fCfTdcs.clear();
}

void TGrDetector::Print(const Char_t *option) const
{
  TDetector::Print(option);
  fCalib.Print(option);
  if (option[0] == 'a')
    PrintData();
}

void TGrDetector::PrintData() const
{
  cout << "TGrDetector data:" << endl;
  for (vector<const TDataItem<UShort_t> *>::const_iterator it =
         fAdcs.begin(), end = fAdcs.end(); it < end; it++)
    (*it)->Print("adc:    ");
  for (vector<const TDataItem<UShort_t> *>::const_iterator it =
         fLeTdcs.begin(), end = fLeTdcs.end(); it < end; it++)
    (*it)->Print("leTdc:  ");
  for (vector<const TDataItem<UShort_t> *>::const_iterator it =
         fCfTdcs.begin(), end = fCfTdcs.end(); it < end; it++)
    (*it)->Print("cfTdc:  ");
  for (vector<const TDataItem<Float_t> *>::const_iterator it = 
         fEnergies.begin(), end = fEnergies.end(); it < end; it++)
    (*it)->Print("energy: ");
  for (vector<const TDataItem<Float_t> *>::const_iterator it =
         fTimes.begin(), end = fTimes.end(); it < end; it++)
    (*it)->Print("time:   ");
}

void TGrDetector::SetSignals(const TSubEvent &subEvent)
{
}

void TGrDetector::Calibrate()
{
  fCalib.GetAdc().Calibrate(fAdcs, fEnergies);
  fCalib.GetLeTdc().Calibrate(fLeTdcs, fTimes);
  fCalib.GetCfTdc().Calibrate(fCfTdcs, fTimes);
}

Bool_t
TGrDetector::PlotEnergies(TH2F *hist, const TCutG *cut, Bool_t option) const
{
  Bool_t status = false;

  for (UInt_t i = 0; i < fEnergies.size(); i++) {
    Float_t val = fEnergies[i]->GetValue();
    UShort_t chan = fCalib.GetAdc().GetHistogramMap(fEnergies[i]->GetChannel());
    if (cut && cut->IsInside(val, chan))
      status = true;
    if (option || status)
      hist->Fill(val, chan);
  }
  return status;
}

Bool_t TGrDetector::PlotChannels(TH1F *hist, Int_t lower,
                                 Int_t upper, Bool_t option) const
{
  Bool_t status = false;

  for (UInt_t i = 0; i < fEnergies.size(); i++) {
    UShort_t channel =
      fCalib.GetAdc().GetHistogramMap(fEnergies[i]->GetChannel());
    if (lower != upper && channel >= lower && channel <= upper)
      status = true;
    if (option || status)
      hist->Fill(channel);
  }
  return status;
}

---