CSCSegAlgoShowering.h

Go to the documentation of this file.

#ifndef CSCSegment_CSCSegAlgoShowering_h
#define CSCSegment_CSCSegAlgoShowering_h

#include <DataFormats/CSCRecHit/interface/CSCRecHit2D.h>
#include <DataFormats/CSCRecHit/interface/CSCSegment.h>
#include <Geometry/CSCGeometry/interface/CSCChamber.h>
#include <FWCore/ParameterSet/interface/ParameterSet.h>

#include <vector>

class CSCSegAlgoShowering {

 public:

  typedef std::vector<const CSCRecHit2D*> ChamberHitContainer;

  explicit CSCSegAlgoShowering(const edm::ParameterSet& ps);

  virtual ~CSCSegAlgoShowering();

  CSCSegment showerSeg( const CSCChamber* aChamber, const ChamberHitContainer& rechits );


 private:

  bool isHitNearSegment(const CSCRecHit2D* h) const;
  bool addHit(const CSCRecHit2D* hit, int layer);
  void updateParameters(void);
  bool hasHitOnLayer(int layer) const;
  void compareProtoSegment(const CSCRecHit2D* h, int layer);
  CLHEP::HepMatrix derivativeMatrix(void) const;
  AlgebraicSymMatrix weightMatrix(void) const;
  AlgebraicSymMatrix calculateError(void) const;
  void flipErrors(AlgebraicSymMatrix&) const;

  void pruneFromResidual();

  // Member variables
  const std::string myName; 
  const CSCChamber* theChamber;

  ChamberHitContainer protoSegment;
  float       protoSlope_u;
  float       protoSlope_v;
  LocalPoint  protoIntercept;       
  double      protoChi2;
  LocalVector protoDirection;

  // input from .cfi file
  bool   debug;
  int    minHitsPerSegment;
  double dRPhiFineMax;
  double dPhiFineMax;
  float  tanPhiMax;
  float  tanThetaMax;
  float  chi2Max;
  float  maxRatioResidual;
//  float  maxDR;
  float  maxDTheta;
  float  maxDPhi;

};
#endif