TMcpDetector.cxx

Go to the documentation of this file.

//  TMcpDetector.cxx.
#include <iostream>

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

#include "TMcpDetector.h"
#include "TSubEvent.h"
#include "THit.h"

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

TMcpDetector::TMcpDetector(const Char_t *name, UInt_t channels)
  : THiDetector(name, channels)
{
}

TMcpDetector::~TMcpDetector()
{
}

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

void TMcpDetector::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 TMcpDetector::FindHits(vector<THit> *hits) const
{
  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);
  }
}

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

  if (fEnergies.size() == GetChannels()) {

    Float_t Ia = fEnergies[0]->GetValue();
    Float_t Ib = fEnergies[1]->GetValue();
    Float_t Ic = fEnergies[2]->GetValue();
    Float_t Id = fEnergies[3]->GetValue();
    Float_t L = 50;

    // from Mike's thesis, page 25
    Float_t x = L/2 * ((Ib + Ic) - (Ia + Id)) / (Ia + Ib + Ic + Id);
    Float_t y = L/2 * ((Ia + Ib) - (Ic + Id)) / (Ia + Ib + Ic + Id);
    cout << Ia << " " << Ib << " " << Ic << " " << Id << endl;
    cout << x << " " << y << endl;
    if (cut && cut->IsInside(x, y))
      status = true;
    if (option || status)
      hist->Fill(x, y);
  }
  return status;
}

---