MTDCPEBase.h

Go to the documentation of this file.

#ifndef RecoLocalFastTime_FTLClusterizer_MTDCPEBase_H
#define RecoLocalFastTime_FTLClusterizer_MTDCPEBase_H 1

//-----------------------------------------------------------------------------
// \class        MTDCPEBase
//-----------------------------------------------------------------------------

#include <utility>
#include <memory>
#include <vector>
#include "TMath.h"

#include "MTDClusterParameterEstimator.h"

#include "Geometry/MTDGeometryBuilder/interface/MTDGeometry.h"
#include "Geometry/MTDGeometryBuilder/interface/MTDGeomDetUnit.h"
#include "Geometry/MTDGeometryBuilder/interface/ProxyMTDTopology.h"
#include "Geometry/MTDGeometryBuilder/interface/RectangularMTDTopology.h"
#include "Geometry/CommonDetUnit/interface/GeomDetType.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementPoint.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementError.h"
#include "DataFormats/GeometrySurface/interface/GloballyPositioned.h"


#include <unordered_map>

#include <iostream>


class MTDCPEBase : public MTDClusterParameterEstimator
{
public:
   struct DetParam
   {
      DetParam() {}
      const MTDGeomDetUnit * theDet;
      const ProxyMTDTopology * theTopol;
      const RectangularMTDTopology * theRecTopol;

      GeomDetType::SubDetector thePart;      
      Local3DPoint theOrigin;
      float theThickness;
      float thePitchX;
      float thePitchY;
   };
   
   struct ClusterParam
   {
      ClusterParam(const FTLCluster & cl) : theCluster(&cl) {}

      virtual ~ClusterParam() = default;

      const FTLCluster * theCluster;
   };
   
public:
   MTDCPEBase(edm::ParameterSet const& conf, const MTDGeometry& geom);  
   
   
   inline ReturnType getParameters(const FTLCluster & cl,
                                   const GeomDetUnit    & det ) const override
   {
      
      DetParam const & dp = detParam(det);
      ClusterParam cp(cl);
      setTheClu( dp, cp );
      auto tuple = std::make_tuple(
                   localPosition(dp, cp),
                   localError(dp, cp),
                   clusterTime(dp, cp),
                   clusterTimeError(dp, cp)
                   );
      return tuple;
   }
   
   //--------------------------------------------------------------------------
   // In principle we could use the track too to get an angle if needed
   //--------------------------------------------------------------------------
   inline ReturnType getParameters(const FTLCluster & cl,
                                   const GeomDetUnit    & det,
                                   const LocalTrajectoryParameters & ltp ) const override
   {
     return getParameters(cl,det);
   }
   
   
   
private:
   //--------------------------------------------------------------------------
   // This is where the action happens.
   //--------------------------------------------------------------------------
   virtual LocalPoint localPosition(DetParam const & dp, ClusterParam & cp) const;
   virtual LocalError localError   (DetParam const & dp, ClusterParam & cp) const;
   virtual TimeValue  clusterTime(DetParam const & dp, ClusterParam & cp) const;
   virtual TimeValueError  clusterTimeError(DetParam const & dp, ClusterParam & cp) const;
   
   void fillDetParams();
   
protected:   
   //---------------------------------------------------------------------------
   //  Data members
   //---------------------------------------------------------------------------
   
   //--- Global quantities   
   const MTDGeometry & geom_;          // geometry
   
protected:
   
   void  setTheClu( DetParam const & dp, ClusterParam & cp ) const ;   
   
   //---------------------------------------------------------------------------
   //  Cluster-level services.
   //---------------------------------------------------------------------------
   
   DetParam const & detParam(const GeomDetUnit & det) const;
   
   using DetParams=std::vector<DetParam>;
   DetParams m_DetParams;
   
};

#endif