OpFlashAnaAlg.cxx

Go to the documentation of this file.

#include "larana/OpticalDetector/OpFlashAnaAlg.h"

#include "TTree.h"

void opdet::OpFlashAnaAlg::SetOpFlashTree(TTree* tree, bool makeOpDetPEHist)
{
  fOpFlashTree = tree;
  fMakeOpDetPEHist = makeOpDetPEHist;

  fOpFlashTree->Branch(
    fOpFlashTree->GetName(), &fOpFlashAnaStruct, fOpFlashAnaStruct.LeafList.c_str());
  if (fMakeOpDetPEHist) {
    fOpFlashAnaStruct.FlashOpDetPEHist->SetName("hOpDetPEs");
    fOpFlashAnaStruct.FlashOpDetPEHist->SetTitle("PEs per Optical Detector; OpDet Number; PE");
    fOpFlashTree->Branch("flash_OpDetPEHist", "TH1D", &(fOpFlashAnaStruct.FlashOpDetPEHist));
  }
}

void opdet::OpFlashAnaAlg::SetOpHitTree(TTree* tree)
{
  fOpHitTree = tree;
  fOpHitTree->Branch(fOpHitTree->GetName(), &fOpHitAnaStruct, fOpHitAnaStruct.LeafList.c_str());
}

void opdet::OpFlashAnaAlg::FillOpFlashes(const std::vector<recob::OpFlash>& flashVector)
{
  if (fOpFlashTree) FillOpFlashTree(flashVector);
}

void opdet::OpFlashAnaAlg::FillOpHits(const std::vector<recob::OpHit>& hitVector)
{
  if (fOpHitTree) FillOpHitTree(hitVector);
}

void opdet::OpFlashAnaAlg::FillOpFlashTree(const std::vector<recob::OpFlash>& flashVector)
{
  for (auto const& flash : flashVector) {

    fOpFlashAnaStruct.FlashTime = flash.Time();
    fOpFlashAnaStruct.FlashTimeWidth = flash.TimeWidth();
    fOpFlashAnaStruct.FlashAbsTime = flash.AbsTime();
    fOpFlashAnaStruct.FlashFrame = flash.Frame();
    fOpFlashAnaStruct.FlashY = flash.YCenter();
    fOpFlashAnaStruct.FlashYWidth = flash.YWidth();
    fOpFlashAnaStruct.FlashZ = flash.ZCenter();
    fOpFlashAnaStruct.FlashZWidth = flash.ZWidth();
    fOpFlashAnaStruct.FlashInBeamFrame = flash.InBeamFrame();
    fOpFlashAnaStruct.FlashOnBeamTime = flash.OnBeamTime();
    fOpFlashAnaStruct.FlashFastToTotal = flash.FastToTotal();
    fOpFlashAnaStruct.FlashTotalPE = flash.TotalPE();

    if (fMakeOpDetPEHist) {
      int n_opdets = flash.WireCenters().size();
      fOpFlashAnaStruct.FlashOpDetPEHist->SetBins(n_opdets, -0.5, (float)n_opdets - 0.5);
      for (int i = 0; i < n_opdets; i++)
        fOpFlashAnaStruct.FlashOpDetPEHist->SetBinContent(i + 1, flash.PE(i));
    }

    fOpFlashTree->Fill();
  }
}

void opdet::OpFlashAnaAlg::FillOpHitTree(const std::vector<recob::OpHit>& hitVector)
{
  for (auto const& hit : hitVector) {

    fOpHitAnaStruct.HitPeakTime = hit.PeakTime();
    fOpHitAnaStruct.HitPeakTimeAbs = hit.PeakTimeAbs();
    fOpHitAnaStruct.HitWidth = hit.Width();
    fOpHitAnaStruct.HitArea = hit.Area();
    fOpHitAnaStruct.HitAmplitude = hit.Amplitude();
    fOpHitAnaStruct.HitFastToTotal = hit.FastToTotal();
    fOpHitAnaStruct.HitPE = hit.PE();
    fOpHitAnaStruct.HitFrame = hit.Frame();
    fOpHitAnaStruct.HitOpChannel = hit.OpChannel();

    fOpHitTree->Fill();
  }
}