EcalTBHodoscopeRecInfoAlgo.h

Go to the documentation of this file.

#ifndef RecoTBCalo_EcalTBHodoscopeReconstructor_EcalTBHodoscopeRecInfoAlgo_HH
#define RecoTBCalo_EcalTBHodoscopeReconstructor_EcalTBHodoscopeRecInfoAlgo_HH

#include "TBDataFormats/EcalTBObjects/interface/EcalTBHodoscopeRawInfo.h"
#include "TBDataFormats/EcalTBObjects/interface/EcalTBHodoscopeRecInfo.h"
#include "Geometry/EcalTestBeam/interface/EcalTBHodoscopeGeometry.h"

#include <vector>
#include <cmath>


class EcalTBHodoscopeRecInfoAlgo {

 public:
  EcalTBHodoscopeRecInfoAlgo();

  explicit EcalTBHodoscopeRecInfoAlgo(int fitMethod, const std::vector<double>& planeShift, const std::vector<double>& zPosition);

  ~EcalTBHodoscopeRecInfoAlgo() 
    {
      if (myGeometry_)
    delete myGeometry_;
    };

  EcalTBHodoscopeRecInfo reconstruct(const EcalTBHodoscopeRawInfo& hodoscopeRawInfo) const;

 private:
  class BeamTrack {
  public:
    float x;
    float xS;
    float xQ;

    bool operator<(BeamTrack &b2)
      {
    return ( fabs(xS) < fabs(b2.xS) );
      }
    
    BeamTrack(float x0,float xs,float xq):x(x0),xS(xs),xQ(xq) {}

    ~BeamTrack() {}
  private:
    BeamTrack() {}
  };

  //Methods taken from h4ana. They can change in a future version

  void clusterPos(float &x, float &xQuality,const int& ipl, const int& xclus, const int& wclus) const ;

  void fitHodo(float &x, float &xQuality,
           const int& ipl, const int& nclus, const std::vector<int>& xclus, const std::vector<int>& wclus) const ;

  void fitLine(float &x, float &xSlope, float &xQuality,
           const int& ipl1, const int& nclus1, const std::vector<int>& xclus1, const std::vector<int>& wclus1,
           const int& ipl2, const int& nclus2, const std::vector<int>& xclus2, const std::vector<int>& wclus2) const;

  int fitMethod_;

  std::vector<double> planeShift_;
  std::vector<double> zPosition_;

  //for the moment mantaining it here
  EcalTBHodoscopeGeometry* myGeometry_;

};

#endif