QwPMT_Channel.h

Go to the documentation of this file.

/*!
 * \file   QwPMT_Channel.h
 * \brief  PMT channel data element for tracking subsystem
 * \author P. M. King
 * \date   2007-05-08
 */
 
#pragma once
 
#include <vector>
#include "TTree.h"
 
//jpan: Mersenne Twistor: A 623-diminsionally equidistributed
//uniform pseudorandom number generator
#include "TRandom3.h"
 
#include "VQwDataElement.h"
 
 
/**
 * \class QwPMT_Channel
 * \ingroup QwTracking
 * \brief PMT channel data element for photomultiplier tube readout
 *
 * Handles data from photomultiplier tube channels including raw values,
 * calibration, and basic data element operations for PMT-based detectors
 * in the tracking system.
 */
class QwPMT_Channel: public VQwDataElement {
  /******************************************************************
   *  Class: QwPMT_Channel
   *
   ******************************************************************/
 public:
  /// Default constructor
  QwPMT_Channel() {
    // Prepare the random number generator.
    gRandom->SetSeed();
  };
  /// Copy constructor
  QwPMT_Channel(const QwPMT_Channel& source)
  : VQwDataElement(source)
  { };
  /// Constructor with name
  QwPMT_Channel(TString name) {
    InitializeChannel(name);
  };
  /// Virtual destructor
  ~QwPMT_Channel() override { };
 
  void  InitializeChannel(TString name){
    SetElementName(name);
    ClearEventData();
  };
 
  void  ClearEventData() override;
  void  RandomizeEventData(int helicity, int SlotNum, int ChanNum);
  void  EncodeEventData(std::vector<UInt_t> &TrigBuffer);
  Int_t ProcessEvBuffer(UInt_t* buffer, UInt_t num_words_left, UInt_t subelement=0) override{return 0;};
 
  void     SetValue(Double_t data) { fValue = data; };
  Double_t GetValue() const        { return fValue; };
  void SetSubbankID(const Int_t bank_index)         { fCrate = bank_index; };
  void SetModule(const Int_t slot_num)              { fModule = slot_num; };
  Int_t           GetSubbankID()     const { return fCrate; };
  Int_t           GetModule()        const { return fModule; };
 
  void  ProcessEvent();
 
  QwPMT_Channel& operator=  (const QwPMT_Channel &value);
  QwPMT_Channel& operator+= (const QwPMT_Channel &value);
  QwPMT_Channel& operator-= (const QwPMT_Channel &value);
  void Sum(const QwPMT_Channel &value1, const QwPMT_Channel &value2);
  void Difference(const QwPMT_Channel &value1, const QwPMT_Channel &value2);
 
  void  ConstructHistograms(TDirectory *folder, TString &prefix) override;
  void  FillHistograms() override;
 
  void  ConstructBranchAndVector(TTree *tree, TString &prefix, QwRootTreeBranchVector &values);
  void  FillTreeVector(QwRootTreeBranchVector &values) const;
 
  void PrintValue() const override;
  void PrintInfo() const override;
  void PrintErrorCounters() const override {};
 
 
 protected:
 
 
 private:
  static const Bool_t kDEBUG;
 
  /*  ADC Calibration                     */
  static const Double_t kPMT_VoltsPerBit;
 
  /*  Channel information data members    */
 
 
  /*  Ntuple array indices */
  size_t fTreeArrayIndex;
  size_t fTreeArrayNumEntries;
 
  /*  Event data members */
  Double_t fValue;
 
  Int_t fCrate;                      ///< ROC number
  Int_t fModule;                     ///< slot number
};