MdMuonCounter.h

Go to the documentation of this file.

#ifndef _MdMuonCounterAG_MdMuonCounter_h_
#define _MdMuonCounterAG_MdMuonCounter_h_

#include <utl/config.h> // before all

#include <fwk/VModule.h>
#include <mevt/MEvent.h>
#include <mevt/Counter.h>
#include <mdet/MDetector.h>
#include <evt/ShowerRecData.h>

#include <utl/Trace-fwd.h>
#include <utl/MessageLoggerConfig.h>
#include <utl/UnitsConfig.h>

#include <string>
#include <iostream>
#include <memory>

#include <boost/utility.hpp>


namespace MdMuonCounterAG {

  class VMuonCounterFunctor;


  class MdMuonCounter : public boost::noncopyable, public fwk::VModule {

  public:
    enum Strategy {
      eGap,             //< Require a gap in 1s.
      eGapRelax,        //< Require at least a number of gaps in between the 1s.
      eAmountInWindow,  //< Only check quantity of 1s in window.
      eConsecutive,     //< Check consecutiveness of 1s.
      eAmountGlobal     //< Count ones and divide by constant.
    };
    static const char* const kStrategyTags[];

    MdMuonCounter() { }
    ~MdMuonCounter() { }
    VModule::ResultFlag Init();
    VModule::ResultFlag Run(evt::Event& e);
    VModule::ResultFlag Finish();

  private:
    utl::MessageLoggerConfig fLog;

    unsigned int fWindowSize;

    unsigned int fNOnes;

    unsigned int fNGaps;

    double fNOnesPerPatternMatch;

    unsigned int fNMinOnes;

    bool fConvertTo80Mhz;

    Strategy fStrategy;

    std::unique_ptr<VMuonCounterFunctor> fCounter;

    utl::UnitsConfig fUnits;

    void DowngradeTrace(utl::TraceB& trace, unsigned int dataReduction = 4);

    unsigned int FillChannelRecData(mevt::Channel& ch);

    void FillModuleRecData(mevt::Module& m, const mdet::Module&) /*const*/;

    std::vector<unsigned int> GetPatternMatchBins(mevt::Module& m) /*const*/;

    void FillChannelsOn(mevt::Module& m, const std::vector<unsigned int>&, const mdet::Module&) const;

    double EstimateNumberOfMuons(const utl::TraceUI&, const size_t segmentation) const;

    utl::TimeInterval GetTraceOffset(const mevt::Module& m, const int counterId) const;

    utl::TimeStamp GetSdTraceStartTime(const evt::Event&, const int counterId) const;

    double ComputeSignalT50(mevt::Counter&, const utl::Vector& rAxis) const;

    double GetMedian(const std::vector<double> v) const;

    double GetTimeCorrection(const utl::Vector& dposition, const utl::Vector& axis) const;

    bool fSetTimeStamps;

    bool fSyncCountingHisto;

    typedef std::map<int, double> TModuleDelay;
    TModuleDelay fModuleDelay;

    REGISTER_MODULE("MdMuonCounterAG", MdMuonCounter);

  };

}


#endif