CEvent/StationRecData.h

Go to the documentation of this file.

#ifndef _cevt_StationRecData_h_
#define _cevt_StationRecData_h_

#include <utl/TimeStamp.h>
#include <utl/ShadowPtr.h>
//#include <utl/ShadowPtr_fwd.h>
#include <utl/TimeDistribution.h>
#include <utl/MultiTimeDistribution.h>

#include <cevt/StationConstants.h>


static const char CvsId_cevt_StationRecData[] =
  "$Id$";


namespace cevt {

  class Station;

  class StationRecData {
  public:
    typedef StationConstants::SignalComponent SignalComponent;

    void SetTotalHits(const unsigned int nhits, const SignalComponent source=StationConstants::eTotal) 
    {fTotalHits[source] = nhits;} 

    void SetTotalHitsPileUpCorrected(const double nhits, const SignalComponent source=StationConstants::eTotal) 
    {fTotalHitsPileUpCorrected[source] = nhits;} 

    unsigned int GetTotalHits(const SignalComponent source=StationConstants::eTotal) const
    { return fTotalHits[source]; } 

    double GetTotalHitsPileUpCorrected(const SignalComponent source=StationConstants::eTotal) const
    { return  fTotalHitsPileUpCorrected[source]; } 


    // Keep these two methods for backward compatibility
    void SetTotalMuonHitsPileUpCorrected(const double nhits) 
    { fTotalMuonHitsPileUpCorrected = nhits; }

    double GetTotalMuonHitsPileUpCorrected() const 
    { return fTotalMuonHitsPileUpCorrected; }


    // Set's and Get's for fiducial area stuff

    void SetFiducialArea(const double & fiducialArea) 
    { fFiducialArea = fiducialArea; }

    void SetHasFiducialAreaInfo(const bool & flag) 
    { fHasFiducialAreaInfo = flag; }

    bool HasFiducialAreaInfo(void) const 
    { return fHasFiducialAreaInfo; }

    double GetFiducialArea(void) const 
    { return fFiducialArea; }

    void SetNumberPadsInFiducialArea(const unsigned int nPads) 
    {fNumberPadsInFiducialArea = nPads;}

    void SetNumberPadsOutFiducialArea(const unsigned int nPads) 
    {fNumberPadsOutFiducialArea = nPads;}

    unsigned int GetNumberPadsInFiducialArea(void) const 
    {return fNumberPadsInFiducialArea;}

    unsigned int GetNumberPadsOutFiducialArea(void) const 
    {return fNumberPadsOutFiducialArea;}


    void SetTotalHitsInFiducialArea(const unsigned int nhits, const SignalComponent source=StationConstants::eTotal) 
    {fTotalHitsInFiducialArea[source] = nhits;}

    void SetTotalHitsOutFiducialArea(const unsigned int nhits, const SignalComponent source=StationConstants::eTotal) 
    {fTotalHitsOutFiducialArea[source] = nhits;}

    void SetTotalHitsInFiducialAreaPileUpCorrected(const double nhits, const SignalComponent source=StationConstants::eTotal) 
    {fTotalHitsInFiducialAreaPileUpCorrected[source] = nhits;}

    void SetTotalHitsOutFiducialAreaPileUpCorrected(const double nhits, const SignalComponent source=StationConstants::eTotal) 
    {fTotalHitsOutFiducialAreaPileUpCorrected[source] = nhits;}


    unsigned int GetTotalHitsInFiducialArea(const SignalComponent source=StationConstants::eTotal) const
    { return fTotalHitsInFiducialArea[source]; }

    unsigned int GetTotalHitsOutFiducialArea(const SignalComponent source=StationConstants::eTotal) const
    { return fTotalHitsOutFiducialArea[source]; }

    double GetTotalHitsInFiducialAreaPileUpCorrected(const SignalComponent source=StationConstants::eTotal) const 
    { return fTotalHitsInFiducialAreaPileUpCorrected[source]; }

    double GetTotalHitsOutFiducialAreaPileUpCorrected(const SignalComponent source=StationConstants::eTotal) const 
    { return fTotalHitsOutFiducialAreaPileUpCorrected[source]; }


    // Set's and Get's for mass overburden stuff
    void SetHasSlantMassInfo(const bool & flag) 
    { fHasSlantMassInfo = flag; }
    void SetAverageSlantMass(const double & mass) 
    {fAverageSlantMass = mass;}
    void SetAverageSlantMassInFiducialArea(const double & mass) 
    {fAverageSlantMassInFiducialArea = mass;}
    void SetAverageSlantMassOutFiducialArea(const double & mass) 
    {fAverageSlantMassOutFiducialArea = mass;}

    bool HasSlantMassInfo(void) const 
    { return fHasSlantMassInfo; }
    double GetAverageSlantMass(void) const 
    { return fAverageSlantMass; }
    double GetAverageSlantMassInFiducialArea(void) const
    { return fAverageSlantMassInFiducialArea; } 
    double GetAverageSlantMassOutFiducialArea(void) const
    { return fAverageSlantMassOutFiducialArea; }


  private:
    StationRecData();
    ~StationRecData() {;}

    unsigned int fTotalHits[StationConstants::eLastSource];
    double fTotalHitsPileUpCorrected[StationConstants::eLastSource];
    double fTotalMuonHitsPileUpCorrected;

    bool fHasFiducialAreaInfo;
    double fFiducialArea;
    unsigned int fNumberPadsInFiducialArea;
    unsigned int fNumberPadsOutFiducialArea;

    unsigned int fTotalHitsInFiducialArea[StationConstants::eLastSource];
    unsigned int fTotalHitsOutFiducialArea[StationConstants::eLastSource];
    double fTotalHitsInFiducialAreaPileUpCorrected[StationConstants::eLastSource];
    double fTotalHitsOutFiducialAreaPileUpCorrected[StationConstants::eLastSource];

    bool fHasSlantMassInfo;
    double fAverageSlantMass;
    double fAverageSlantMassInFiducialArea;
    double fAverageSlantMassOutFiducialArea;

  public:
    // Time distribution of hits in fiducial area

    utl::TimeDistributionI&  
    GetHitsInFiducialArea(const SignalComponent source = StationConstants::eTotal)
    { return fHitsInFiducialArea->GetTimeDistribution(source); }
    const utl::TimeDistributionI& 
    GetHitsInFiducialArea(const SignalComponent source = StationConstants::eTotal) const
    { return fHitsInFiducialArea->GetTimeDistribution(source); }

    void MakeHitsInFiducialArea(const SignalComponent source = StationConstants::eTotal);

    bool HasHitsInFiducialArea(const SignalComponent source = StationConstants::eTotal) const
    { return fHitsInFiducialArea->HasLabel(source); }

    const utl::MultiTimeDistributionI& GetMultiHitsInFiducialArea()
    { return *fHitsInFiducialArea; }

    void ClearHitsInFiducialArea()
    { fHitsInFiducialArea->Clear(); }

    typedef utl::MultiTimeDistributionI::Iterator HitsInFiducialAreaIterator;
    typedef utl::MultiTimeDistributionI::ConstIterator ConstHitsInFiducialAreaIterator;

    HitsInFiducialAreaIterator HitsInFiducialAreaBegin() { return fHitsInFiducialArea->Begin(); }
    HitsInFiducialAreaIterator HitsInFiducialAreaEnd()   { return fHitsInFiducialArea->End(); }
    ConstHitsInFiducialAreaIterator HitsInFiducialAreaBegin() const { return fHitsInFiducialArea->Begin(); }
    ConstHitsInFiducialAreaIterator HitsInFiducialAreaEnd() const { return fHitsInFiducialArea->End(); }

    // Time distribution of hits outside fiducial area

    utl::TimeDistributionI&  
    GetHitsOutFiducialArea(const SignalComponent source = StationConstants::eTotal)
    { return fHitsOutFiducialArea->GetTimeDistribution(source); }
    const utl::TimeDistributionI& 
    GetHitsOutFiducialArea(const SignalComponent source = StationConstants::eTotal) const
    { return fHitsOutFiducialArea->GetTimeDistribution(source); }

    void MakeHitsOutFiducialArea(const SignalComponent source = StationConstants::eTotal);

    bool HasHitsOutFiducialArea(const SignalComponent source = StationConstants::eTotal) const
    { return fHitsOutFiducialArea->HasLabel(source); }

    const utl::MultiTimeDistributionI& GetMultiHitsOutFiducialArea()
    { return *fHitsOutFiducialArea; }
    
    void ClearHitsOutFiducialArea()
    { fHitsOutFiducialArea->Clear(); }

    typedef utl::MultiTimeDistributionI::Iterator HitsOutFiducialAreaIterator;
    typedef utl::MultiTimeDistributionI::ConstIterator ConstHitsOutFiducialAreaIterator;

    HitsOutFiducialAreaIterator HitsOutFiducialAreaBegin() { return fHitsOutFiducialArea->Begin(); }
    HitsOutFiducialAreaIterator HitsOutFiducialAreaEnd()   { return fHitsOutFiducialArea->End(); }
    ConstHitsOutFiducialAreaIterator HitsOutFiducialAreaBegin() const { return fHitsOutFiducialArea->Begin(); }
    ConstHitsOutFiducialAreaIterator HitsOutFiducialAreaEnd() const { return fHitsOutFiducialArea->End(); }

  private:

    utl::InitializedShadowPtr<utl::MultiTimeDistributionI> fHitsInFiducialArea;
    utl::InitializedShadowPtr<utl::MultiTimeDistributionI> fHitsOutFiducialArea;

    friend class cevt::Station;
    friend class utl::ShadowPtr<StationRecData>;

  }; 

} // end of namespace cevt

#endif

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// compile-command: "make -C .. -k"
// End: