TBgoDetector.cxx

Go to the documentation of this file.

//  TBgoDetector.cxx.
#include <iostream>

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

#include "TBgoDetector.h"
#include "TSubEvent.h"
#include "THit.h"

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

TBgoDetector::TBgoDetector(const Char_t *name, UInt_t channels)
  : TGrDetector(name, channels)
{
}

TBgoDetector::~TBgoDetector()
{
}

void TBgoDetector::Print(const Char_t *option) const
{
  cout << "BGO detector:" << endl;
  TGrDetector::Print(option);
}

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

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

void TBgoDetector::FindHits(vector<THit> *hits) const
{
  //
  // we only rebuild 1 hit here
  //
  THit p;

  if (fEnergies.size() && fEnergies[0]->GetValue() > 0) {
    p.SetEnergy(fEnergies[0]->GetValue());
    if (fTimes.size())
      p.SetTime(fTimes[0]->GetValue());
    p.SetDetector(GetName());
    p.SetPosition(fCoordinates.X(), fCoordinates.Y(), fCoordinates.Z());
    hits->push_back(p);
  }
}

void TBgoDetector::SetXYZ(Int_t x, Int_t y, Int_t z)
{
  fCoordinates.SetXYZ(x, y, z);
}

Bool_t TBgoDetector::PlotShiftedEnergies(TH2F *hist, Float_t valueShift,
                                         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());
    val += valueShift;
    if (cut && cut->IsInside(val, chan))
      status = true;
    if (option || status)
      hist->Fill(val, chan);
  }
  return status;
}

---