QwSIS3320_Accumulator.h

Go to the documentation of this file.

/**
 * \class QwSIS3320_Accumulator QwSIS3320_Accumulator.h
 *
 * \brief SIS3320 sampling ADC accumulator
 *
 * \author  W. Deconinck
 * \date  2009-09-04 18:06:23
 * \ingroup QwCompton
 *
 * The QwSIS3320_Accumulator should allow convenient access to the accumulator
 * data collected with the SIS3320 for the Compton photon detector.  This class
 * implements its own sum, difference, and ratio methods inherited from the
 * general VQwDataElement.
 *
 */

#ifndef __QwSIS3320_Accumulator__
#define __QwSIS3320_Accumulator__

// System headers
#include <iostream>
#include <vector>
#include <numeric>

// ROOT headers
#include <TTree.h>

// Boost math library for random number generation
#include <boost/random.hpp>

// Qweak headers
#include "VQwDataElement.h"

class QwSIS3320_Accumulator: public VQwDataElement {

  public:

    QwSIS3320_Accumulator(TString name = "") {
      SetElementName(name);
      SetNumberOfDataWords(3);
      fMinValue = 0; fMaxValue = 0;
      fMinTime = 0; fMaxTime = 0;
    };
    virtual ~QwSIS3320_Accumulator() { };

    Int_t GetMinValue() const { return fMinValue; };
    Int_t GetMaxValue() const { return fMaxValue; };
    Int_t GetMinTime() const { return fMinTime; };
    Int_t GetMaxTime() const { return fMaxTime; };

    Double_t GetNumberOfSamples() const { return fNumberOfSamples; };
    Double_t GetAccumulatorSum() const { return fAccumulatorSum; };
    Double_t GetAccumulatorAvg() const {
      if (fAccumulatorSum > 0)
        return fAccumulatorSum / fNumberOfSamples;
      else return 0.0;
    };

    void  ClearEventData() { fAccumulatorSum = 0; fNumberOfSamples = 0; };
    Int_t ProcessEvBuffer(UInt_t* buffer, UInt_t num_words_left, UInt_t subelement = 0);
    void  ProcessEvent() { };

    const QwSIS3320_Accumulator operator/ (const Double_t &value) const;
    const QwSIS3320_Accumulator operator* (const Double_t &value) const;
    const QwSIS3320_Accumulator operator+ (const Double_t &value) const;
    const QwSIS3320_Accumulator operator- (const Double_t &value) const;
    QwSIS3320_Accumulator& operator/= (const Double_t &value);
    QwSIS3320_Accumulator& operator*= (const Double_t &value);
    QwSIS3320_Accumulator& operator+= (const Double_t &value);
    QwSIS3320_Accumulator& operator-= (const Double_t &value);
    const QwSIS3320_Accumulator operator+ (const QwSIS3320_Accumulator &value) const;
    const QwSIS3320_Accumulator operator- (const QwSIS3320_Accumulator &value) const;
    QwSIS3320_Accumulator& operator=  (const QwSIS3320_Accumulator &value);
    QwSIS3320_Accumulator& operator+= (const QwSIS3320_Accumulator &value);
    QwSIS3320_Accumulator& operator-= (const QwSIS3320_Accumulator &value);
    void Sum(const QwSIS3320_Accumulator &value1, const QwSIS3320_Accumulator &value2);
    void Difference(const QwSIS3320_Accumulator &value1, const QwSIS3320_Accumulator &value2);
    void Ratio(const QwSIS3320_Accumulator &numer, const QwSIS3320_Accumulator &denom);

    void  ConstructHistograms(TDirectory *folder, TString &prefix);
    void  FillHistograms();

    void  ConstructBranchAndVector(TTree *tree, TString &prefix, std::vector<Double_t> &values);
    void  FillTreeVector(std::vector<Double_t> &values) const;

    // Output stream operator<< for an accumulator
    friend std::ostream& operator<< (std::ostream& stream, const QwSIS3320_Accumulator& a);

    // Define operator+=, operator-=, etc, for accumulator asymmetries

  protected:

    Double_t fNumberOfSamples;  ///< Number of accumulated samples
    Double_t fAccumulatorSum; ///< Accumulator sum
    Double_t fAccumulatorAvg; ///< Accumulator average

  private:

    int fMinValue, fMaxValue; //! Value-based accumulator limits
    int fMinTime, fMaxTime; //! Time-based accumulator limits

    // Ntuple array indices
    size_t fTreeArrayIndex; //! Index of this data element in tree
    size_t fTreeArrayNumEntries; //! Number of entries from this data element

    static const unsigned int INDEX_NUM;
    static const unsigned int INDEX_SUM;

};

// Output stream operator<< for the accumulators
inline std::ostream& operator<< (std::ostream& stream, const QwSIS3320_Accumulator& a) {
  stream << a.GetAccumulatorSum() << " (" << a.GetNumberOfSamples() << ")";
  return stream;
}

#endif // __QwSIS3320_Accumulator__