MuonProductionDepthFinder.h

Go to the documentation of this file.

#ifndef _MuonProductionDepthFinderGL_MuonProductionDepthFinder_h_
#define _MuonProductionDepthFinderGL_MuonProductionDepthFinder_h_

#include <fwk/VModule.h>
#include <evt/Event.h>
#include <utl/Vector.h>
#include <utl/Point.h>
#include <utl/AxialVector.h>
#include <utl/TimeStamp.h>

#include <sevt/Station.h>
#include <sevt/StationConstants.h>

#include <atm/ProfileResult.h>

namespace evt {
  class Event;
}

namespace sevt {
  class SEvent;
};


namespace MuonProductionDepthFinderGL {

  class MuonProductionDepthFinder : public fwk::VModule {

  public:
    MuonProductionDepthFinder() { }
    virtual ~MuonProductionDepthFinder() { }

    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event& event);
    fwk::VModule::ResultFlag Finish();

  private:
    void InitCoordinateSys(const evt::Event& event,
                           double dX, double dY,
                           double dTheta, double dPhi);
    void SetMuonProductionDepthHist(evt::Event& event);
    bool IsContained(const sevt::SEvent& sevent) const;

    atm::ProfileResult fSlantDepthProfile;
    std::vector<double> fHLAY;
    std::vector<double> fAATM;
    std::vector<double> fBATM;
    std::vector<double> fCATM;
    double fSlantDepth;

    utl::Point fCore;
    utl::TimeStamp fCoreTime;

    utl::CoordinateSystemPtr fLocalCS;
    utl::CoordinateSystemPtr fShowerCS;
    double fTheta;
    double fPhi;

    double fSignalThreshold;
    bool fUseConstRadiusCut;
    double fRadiusCut;
    utl::TabulatedFunction fTblFuncRadiusCut;
    double fTimeShift;
    double fSignalMin;

    double fXmuMin;
    double fXmuMax;
    double fNXmuBins;

    bool fUsePionPathCorrection;

    REGISTER_MODULE("MuonProductionDepthFinderGL", MuonProductionDepthFinder);

  };

}


#endif

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