TECALTestbeamModule.cxx

Go to the documentation of this file.

#include "TDatum.hxx"
#include "TPrincipal.h"
#include "TRealDatum.hxx"

#include "TECALTestbeamModule.hxx"


#include "../cvstags/TECALTestbeamModule.cxx"

ND::TECALTestbeamModule::TECALTestbeamModule(const char* name,
                                             const char* title) {
  SetNameTitle(name, title);
  // Enable this module by default:
  fIsEnabled = kTRUE;
  fDescription = "ECALTestbeam recon output";
  fCVSTagName = CVSTAG;
  fCVSID = CVSID;
}

ND::TECALTestbeamModule::~TECALTestbeamModule() {}

Bool_t ND::TECALTestbeamModule::ProcessFirstEvent(ND::TND280Event& event) {
  return true;
}

void ND::TECALTestbeamModule::InitializeModule() {}

void ND::TECALTestbeamModule::InitializeBranches() {
  fOutputTree->Branch("Cerenkov1Lo", &Cerenkov1Lo, "Cerenkov1Lo[23]/I");
  fOutputTree->Branch("Cerenkov1Hi", &Cerenkov1Hi, "Cerenkov1Hi[23]/I");
  fOutputTree->Branch("Cerenkov2Lo", &Cerenkov2Lo, "Cerenkov2Lo[23]/I");
  fOutputTree->Branch("Cerenkov2Hi", &Cerenkov2Hi, "Cerenkov2Hi[23]/I");
  fOutputTree->Branch("TOF", &TOF, "TOF[23]/I");
  fOutputTree->Branch("TriggerWord", &TriggerWord, "TriggerWord/I");

  fOutputTree->Branch("TestBeamPID", &PIDResult, "PIDResult/I");
  fOutputTree->Branch("TestBeamMomentum", &Momentum, "Momentum/D");
  fOutputTree->Branch("TestBeamAngle", &Angle, "Angle/I");
}

bool ND::TECALTestbeamModule::FillTree(ND::TND280Event& event) {
  ND::THandle<ND::TIntegerDatum> Ch1Lo =
      event.Get<ND::TIntegerDatum>("Cerenkov1Lo");
  ND::THandle<ND::TIntegerDatum> Ch2Lo =
      event.Get<ND::TIntegerDatum>("Cerenkov2Lo");
  ND::THandle<ND::TIntegerDatum> Ch1Hi =
      event.Get<ND::TIntegerDatum>("Cerenkov1Hi");
  ND::THandle<ND::TIntegerDatum> Ch2Hi =
      event.Get<ND::TIntegerDatum>("Cerenkov2Hi");
  ND::THandle<ND::TIntegerDatum> TOFLo =
      event.Get<ND::TIntegerDatum>("TOFHighGain");
  ND::THandle<ND::TIntegerDatum> TOFHi =
      event.Get<ND::TIntegerDatum>("TOFLowGain");
  ND::THandle<ND::TIntegerDatum> TWord =
      event.Get<ND::TIntegerDatum>("TriggerWord");

  for (int i = 0; i < 23; ++i) {
     Cerenkov1Lo[i] = Ch1Lo->at(i);
     Cerenkov1Hi[i] = Ch2Lo->at(i);
     Cerenkov2Lo[i] = Ch1Hi->at(i);
     Cerenkov2Hi[i] = Ch2Hi->at(i);
     TOF[i] = TOFLo->at(i);
  }

  TriggerWord = TWord->GetValue();

  ND::THandle<ND::TIntegerDatum> pidResultDatum =
      event.Get<ND::TIntegerDatum>("testBeamPID");
  ND::THandle<ND::TIntegerDatum> angleDatum =
      event.Get<ND::TIntegerDatum>("testBeamAngle");
  ND::THandle<ND::TRealDatum> momentumDatum =
      event.Get<ND::TRealDatum>("testBeamMomentum");

  PIDResult = pidResultDatum->GetValue();
  Momentum = momentumDatum->GetValue();
  Angle = angleDatum->GetValue();

  return true;
}