TIcDetector.cxx

Go to the documentation of this file.

//  TIcDetector.cxx.
#include <iostream>
#include <sstream>

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

#include "TIcDetector.h"
#include "TSubEvent.h"
#include "THit.h"

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

TIcDetector::TIcDetector(const Char_t *name, UInt_t channels)
  : THiDetector(name, channels)
{
  stringstream plate;
  fPlateNames.assign(channels, name);

  for (unsigned int i = 0; i < channels; i++) {
    plate << i;
    fPlateNames[i].push_back('_');
    fPlateNames[i] += plate.str().data();
    plate.str(string());
  }
}

TIcDetector::~TIcDetector()
{
}

void TIcDetector::Print(const Char_t *option) const
{
  cout << "IC detector:" << endl;
  THiDetector::Print(option);
}

void TIcDetector::SetSignals(const TSubEvent &subEvent)
{
  const vector<TDataItem<UShort_t> > &adcs = subEvent.GetAdcs();

  for (UInt_t i = 0; i < adcs.size(); i++)
    if (fCalib.GetAdc().TestChannel(adcs[i].GetChannel())) {
      fAdcs.push_back(&adcs[i]);
      fHit = true;
    }
}

void TIcDetector::FindHits(vector<THit> *hits) const
{
  THit p;

  for (UInt_t i = 0; i < fEnergies.size(); i++) {
    p.SetDetector(fPlateNames[fEnergies[i]->GetChannel()].data());
    p.SetEnergy(fEnergies[i]->GetValue());
    hits->push_back(p);
  }
}

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

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

Bool_t TIcDetector::PlotEvT(TH2F *hist, Float_t time,
                            const TCutG *cut, Bool_t option) const
{
  Bool_t status = false;
  
  for (UInt_t i = 0; i < fEnergies.size() &&
         fEnergies[i]->GetChannel() == 0; i++) {
    Float_t energy = fEnergies[i]->GetValue();
    if (cut && cut->IsInside(time, energy))
      status = true;
    if (option || status)
      hist->Fill(time, energy);
  }
  return status;
}

Bool_t TIcDetector::PlotAnodeSum(TH1F *hist, Int_t lower,
                                 Int_t upper, Bool_t option) const
{
  Bool_t status = false;
  Float_t sum = 0;

  for (UInt_t i = 0; i < fEnergies.size(); i++)
    sum += fEnergies[i]->GetValue();
  if (lower != upper && sum >= lower && sum <= upper)
    status = true;
  if (option || status)
    hist->Fill(sum);
  return status;
}

Bool_t TIcDetector::PlotEdE(TH2F *hist, UShort_t echan, UShort_t dechan,
                            const TCutG *cut, Bool_t option) const
{
  Bool_t status = false;
  Float_t e  = 0;
  Float_t de = 0;

  for (UInt_t i = 0; i < fEnergies.size(); i++)
    if (fEnergies[i]->GetChannel() == echan)
      e = fEnergies[i]->GetValue();
    else if (fEnergies[i]->GetChannel() == dechan)
      de = fEnergies[i]->GetValue();
  if (cut && cut->IsInside(e, de))
    status = true;
  if ((option || status) && e && de)
    hist->Fill(e, de);
  return status;
}

---