HoughTransNCPi0.h

Go to the documentation of this file.

//
// HoughTrans
//
// Mark's hough transform.
//
// Created:  M. Dorman -- May, 2005
//
// $Author: rhatcher $ 
//
// $Revision: 1.2 $
// 
// $Name:  $
//
// $Id: HoughTransNCPi0.h,v 1.2 2008/05/02 15:46:25 rhatcher Exp $
//



#ifndef houghtransNCPi0_h
#define houghtransNCPi0_h

#include "TH2F.h"
#include "CandNtupleSR/NtpSRStrip.h"

class HoughTransNCPi0 {

 public:
  ~HoughTransNCPi0();
  HoughTransNCPi0();
  HoughTransNCPi0(Int_t numSamps,Float_t lowGrad,Float_t highGrad);
  HoughTransNCPi0(Int_t numSamps,Float_t lowGrad,Float_t highGrad,Float_t lowIcept,Float_t highIcept);
  void SetVtxz(Float_t vertz);
  Float_t GetVtxz();
  void FillHough(NtpSRStrip* stp);
  void FillHough(Float_t zPos, Float_t tPos);
 // void FillHoughfixed(Float_t zPos, Float_t tPos);
  void FillHoughweight(Float_t zPos, Float_t tPos, Float_t energy);
  void DrawHough();
  Int_t EntriesHough();
  void ResetHough();
  void DeleteHough();
  Int_t GetPeakHeight();
  Int_t GetPeakGradBin();
  Int_t GetPeakIceptBin();
  Float_t GetPeakGradVal();
  Float_t GetPeakIceptVal();
  Float_t GetCos();
  Float_t GetGradientweight();
  Float_t GetTransverseVtxCoord();
  Float_t GetIntercept();
  Float_t GetRms(Float_t frac);

  // the arguement frac should take values from [0,1] and
  // represents the fraction of fullness a bin must have
  // relative to the peak to be included in the RMS 
  // calculation (I use 0.75)

  // before calling GetRms() it is worth checking that the
  // peak height in Hough space is not zero, in fact 'low'
  // values for the peak height mean that the hits that 
  // contribute to the peak (and hence the 'track' found)
  // may not actually be a track at all

  TH2F* houghSpace;
  Float_t vtxz;

 private:

  //TH2F houghSpace;
  //Float_t vtxz;

};

#endif

---