TDetector.cxx

Go to the documentation of this file.

//  TDetector.cxx.
#include <iostream>
#include <iomanip>

#include "TDetector.h"
#include "THit.h"

using std::cout;
using std::cerr;
using std::endl;
using std::setw;

TDetector::TDetector(const Char_t *name, UInt_t channels)
  : fHit(false), fChannels(channels), fName(name)
{
  fEnergies.reserve(channels);
  fTimes.reserve(channels);
}

TDetector::TDetector(const TDetector &rhs)
{
  vector<const TDataItem<Float_t> *>::const_iterator it;
  vector<const TDataItem<Float_t> *>::const_iterator end;
  for (it = rhs.fEnergies.begin(), end = rhs.fEnergies.end(); it < end; it++)
    fEnergies.push_back(*it);
  for (it = rhs.fTimes.begin(), end = rhs.fTimes.end(); it < end; it++)
    fTimes.push_back(*it);
  fCoordinates = rhs.fCoordinates;
  fHit         = rhs.fHit;
  fChannels    = rhs.fChannels;
  fName        = rhs.fName;
}

TDetector::~TDetector()
{
  Clear();
}

TDetector & TDetector::operator=(const TDetector &rhs)
{
  if (&rhs != this) {
    vector<const TDataItem<Float_t> *>::const_iterator it;
    vector<const TDataItem<Float_t> *>::const_iterator end;
    Clear();
    for (it = rhs.fEnergies.begin(), end = rhs.fEnergies.end(); it < end; it++)
      fEnergies.push_back(*it);
    for (it = rhs.fTimes.begin(), end = rhs.fTimes.end(); it < end; it++)
      fTimes.push_back(*it);
    fCoordinates = rhs.fCoordinates;
    fHit         = rhs.fHit;
    fChannels    = rhs.fChannels;
    fName        = rhs.fName;
  }
  return *this;
}

void TDetector::Clear(const Char_t *option)
{
  vector<const TDataItem<Float_t> *>::iterator it;
  vector<const TDataItem<Float_t> *>::iterator end;
  for (it = fEnergies.begin(), end = fEnergies.end(); it < end; it++)
    delete *it;
  fEnergies.clear();
  for (it = fTimes.begin(), end = fTimes.end(); it < end; it++)
    delete *it;
  fTimes.clear();
  fHit = false;
}

void TDetector::Print(const Char_t *option) const
{
  if (! fChannels)
    cout << "Print: WARNING - detector channels set to 0" << endl;
  if (fName.empty())
    cout << "Print: WARNING - detector name not set"      << endl;
  else {
    cout << "Detector " << fName << " with " << fChannels;
    fChannels != 1 ? cout << " channels" << endl : cout << " channel" << endl;
  }
  fCoordinates.Print();
}

void TDetector::PrintData() const
{
}

void TDetector::Calibrate()
{
}

Bool_t
TDetector::PlotEnergies(TH2F *hist, const TCutG *cut, Bool_t option) const
{
  return false;
}

Bool_t TDetector::PlotTimes(TH2F *hist, const TCutG *cut, Bool_t option) const
{
  return false;
}

Bool_t TDetector::PlotEvT(TH2F *hist, Float_t time,
                          const TCutG *cut, Bool_t option) const
{
  return false;
}

void TDetector::FindHits(vector<THit> *hits) const
{
}

Float_t TDetector::GetMaxEnergy() const
{
  vector<const TDataItem<Float_t> *>::const_iterator max =
    max_element(fEnergies.begin(), fEnergies.end(), GreaterValue<Float_t>());
  return (max != fEnergies.end()) ? (*max)->GetValue() : 0;
}

---