DTSegmentUpdator.h

Go to the documentation of this file.

#ifndef DTSegment_DTSegmentUpdator_h
#define DTSegment_DTSegmentUpdator_h

/* C++ Headers */
#include <vector>

#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"

#include "FWCore/Framework/interface/ESHandle.h"

#include "Geometry/DTGeometry/interface/DTGeometry.h"

/* ====================================================================== */

/* Collaborating Class Declarations */
class DTSegmentCand;
class DTRecSegment2D;
class DTRecSegment4D;
class DTLinearFit;
class DTRecHitBaseAlgo;
class DTChamberRecSegment2D;

namespace edm{class EventSetup; class ParameterSet;}

/* Class DTSegmentUpdator Interface */
class DTSegmentUpdator{

  public:

    DTSegmentUpdator(const edm::ParameterSet& config) ;

    ~DTSegmentUpdator() ;

    /* Operations */ 

    bool fit(DTSegmentCand* seg) const;

    void fit(DTRecSegment2D* seg) const;

    void fit(DTRecSegment4D* seg) const;

    void update(DTRecSegment4D* seg, const bool calcT0 = false) const;

    void update(DTRecSegment2D* seg) const;

    void calculateT0corr(DTRecSegment2D* seg) const;
    void calculateT0corr(DTRecSegment4D* seg) const;

    void setES(const edm::EventSetup& setup);

  protected:

  private:
    DTLinearFit* theFitter; // the linear fitter
    DTRecHitBaseAlgo* theAlgo; // the algo for hit reconstruction
    edm::ESHandle<DTGeometry> theGeom; // the geometry

    void updateHits(DTRecSegment2D* seg,
                    GlobalPoint &gpos,
                    GlobalVector &gdir,
                    const int step=2) const;

   //rejects bad hits (due to deltas) for phi segment
   void rejectBadHits(DTChamberRecSegment2D* ) const;

    void fit(const std::vector<float>& x,
             const std::vector<float>& y, 
             const std::vector<float>& sigy,
             LocalPoint& pos,
             LocalVector& dir,
             AlgebraicSymMatrix& covMat,
             double& chi2) const;

    void Fit4Var(const std::vector<float>& xfit,
                 const std::vector<float>& yfit,
                 const std::vector<int>& lfit,
                 const std::vector<double>& tfit,
                 const int nptfit,
                 float& cminf, double& vminf,
                 double& chi2fit) const;

    bool vdrift_4parfit;
    double T0_hit_resolution;
    bool perform_delta_rejecting;
    bool debug;

};

#endif // DTSegment_DTSegmentUpdator_h