MD2ADST.h

Go to the documentation of this file.

#ifndef _MD2ADST_h_
#define _MD2ADST_h_

#include <utl/CoordinateSystem.h>
#include <evt/ShowerMRecData.h>
#include <mdet/MDetector.h>
#include <mevt/MEvent.h>
#include <adst/RecEvent.h>


class RecEvent;
namespace evt {
  class Event;
}

namespace otoa {

  class Config;


  /*
    \author Martin Maur
    \date 14 Apr 2012
    \author Diego Ravignani
    \date 31 Aug 2015
  */

  class MD2ADST {
  public:
    MD2ADST(const otoa::Config& config);

    void Convert(const evt::Event& inEvent, RecEvent& outEvent) const;

  private:
    const otoa::Config& fConfig;

    // Make master detector data available across the board
    const det::Detector* fDetector;
    const mdet::MDetector* fMdDetector;
    utl::CoordinateSystemPtr fReferenceCS;

    inline const otoa::Config& GetConfig() const { return fConfig; }

    void FillMEvent(const evt::Event& e, MDEvent& me) const;
    void FillRecShower(const evt::ShowerMRecData& sr, MdRecShower& mr) const;
    MdRecCounter MakeCounter(const mevt::Counter& c, const utl::CoordinateSystemPtr& cs) const;
    // passing the mdet::Module to lookup master data
    MdRecModule MakeModule(const mevt::Module& em, const mdet::Module& dm, const utl::CoordinateSystemPtr& cs) const;
    // no option but to pass mdet::Module again, the needed master data is here (but should be at the channel level!)
    MdRecChannel MakeChannel(const mevt::Channel& c, const mdet::Module& m) const;

    // MakeSimCounter still needs to be reorganised in the same way done for reconstructed data
    // A coordinate system is provided to find the counter coordinates relative to the shower plane
    MdSimCounter MakeSimCounter(const mevt::Counter& c, const utl::CoordinateSystemPtr& showerPlaneCS) const;

    utl::CoordinateSystemPtr GetShowerPlaneCS(const evt::ShowerMRecData& sr) const;

    void AddParticle(const utl::Particle& p, const mdet::Counter& c, const int moduleId, const int scintId, MdSimScintillator& ms) const;

  };

}


#endif