MdCornerClippingCorrecter.h

Go to the documentation of this file.

#ifndef _MdCornerClippingCorrecter_MdCornerClippingCorrecter_h_
#define _MdCornerClippingCorrecter_MdCornerClippingCorrecter_h_
#include <fwk/VModule.h>
#include <mevt/Module.h>
#include <mdet/Module.h>
#include <utl/Vector.h>

namespace MdCornerClippingCorrecterAG {

  class MdCornerClippingCorrecter : public fwk::VModule {
    public:
      MdCornerClippingCorrecter();
      virtual ~MdCornerClippingCorrecter();


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

    private:
      // parametrization pcc(theta, dphi) = b(theta) + m(theta)|sin(dphi)|
      // m(theta) = fPar0 + fPar1*sec(theta)
      // b(theta) = fPar2 + fPar3*sec(theta)
      double fPar0 = 0;
      double fPar1 = 0;
      double fPar2 = 0;
      double fPar3 = 0;
      bool fComputeCornerClippingProbability = true;

      //double GetSinDeltaPhi(const mdet::Module& mDet, const utl::Vector& axis, const utl::CoordinateSystemPtr& coreCS, const bool OutputVerbose);
      //double GetCornerClippingProbability(const mdet::Module& mDet, const utl::Vector& axis, const utl::CoordinateSystemPtr& coreCS, const bool OutputVerbose);
      //double GetNumberOfEstimatedMuonsCorrected(const mevt::Module& mEvt, const mdet::Module& mDet, const utl::Vector& axis, const utl::CoordinateSystemPtr& cs, const bool OutputVerbose);
      void FillCornerClippingProbability(const evt::Event& event, mevt::Module& module, const mdet::Module& mdetModule);

      REGISTER_MODULE("MdCornerClippingCorrecter", MdCornerClippingCorrecter);
  };

}

#endif