OpFlashAnaAlg.h

Go to the documentation of this file.

#ifndef OPFLASHANAALG_H
#define OPFLASHANAALG_H

#include <vector>
#include <memory>

#include "TTree.h"
#include "TH1.h"

#include "lardataobj/RecoBase/OpHit.h"
#include "lardataobj/RecoBase/OpFlash.h"

namespace opdet{

  class OpFlashAnaAlg{

  public:

    OpFlashAnaAlg(){ fMakeOpDetPEHist=false; }
    void SetOpFlashTree(TTree*, bool makeOpDetPEHist=true);
    void SetOpHitTree(TTree*);
    
    void FillOpFlashes(const std::vector<recob::OpFlash>&);
    void FillOpHits(const std::vector<recob::OpHit>&);
    
  private:

    struct FlashAnaStruct{
      double        FlashTime;
      double        FlashTimeWidth;
      double        FlashAbsTime;
      double        FlashY;
      double        FlashYWidth;
      double        FlashZ;
      double        FlashZWidth;
      double        FlashFastToTotal;
      double        FlashTotalPE;
      unsigned int  FlashFrame;
      int           FlashOnBeamTime;
      bool          FlashInBeamFrame;

      TH1D*         FlashOpDetPEHist;

      std::string   LeafList;
      FlashAnaStruct()
      : LeafList("time/D:timewidth/D:abstime/D:y/D:ywidth/D:z/D:zwidth/D:fasttototal/D:totalpe/D:onbeamtime/I:frame/i:inbeamframe/O")
      { FlashOpDetPEHist = new TH1D(); }
      ~FlashAnaStruct() { delete FlashOpDetPEHist; }
    };

    FlashAnaStruct fOpFlashAnaStruct;    
    bool fMakeOpDetPEHist;

    TTree* fOpFlashTree;
    void   FillOpFlashTree(const std::vector<recob::OpFlash>&);
    

    struct HitAnaStruct{
      double        HitPeakTime;
      double        HitPeakTimeAbs;
      double        HitWidth;
      double        HitArea;
      double        HitAmplitude;
      double        HitFastToTotal;
      double        HitPE;
      unsigned int  HitFrame;
      int           HitOpChannel;

      std::string   LeafList;
      HitAnaStruct()
      : LeafList("time/D:abstime/D:width/D:area/D:amplitude/D:fasttototal/D:pe/D:frame/I:opchannel/i")
      { }
    };
    HitAnaStruct fOpHitAnaStruct;

    TTree* fOpHitTree;
    void   FillOpHitTree(const std::vector<recob::OpHit>&);

  };
  
}

#endif