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