CbmRichRingFitterBase.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CbmRichRingFitterBase.h

#ifndef CBM_RICH_RING_FITTER_BASE
#define CBM_RICH_RING_FITTER_BASE

#include "CbmRichRingLight.h"
class CbmRichRingFitterBase
{
public:

   CbmRichRingFitterBase(){ }

  virtual ~CbmRichRingFitterBase() { }

  virtual void DoFit(
        CbmRichRingLight* ring) = 0;

protected:

  virtual void CalcChi2(
        CbmRichRingLight* ring)
  {
     int nofHits = ring->GetNofHits();
      if ( nofHits < 4 ) {
         ring->SetChi2(-1.);
         return;
      }

      float chi2 = 0.;
      float r = ring->GetRadius();
      float xc = ring->GetCenterX();
      float yc = ring->GetCenterY();

      for (int i = 0; i < nofHits; i++) {
         float xh = ring->GetHit(i).fX;
         float yh = ring->GetHit(i).fY;
         float d = r - sqrt((xc - xh)*(xc - xh) + (yc - yh)*(yc - yh));

         chi2 += d*d;
      }
      ring->SetChi2(chi2);
  }

  static const int MAX_NOF_HITS_IN_RING = 400; // maximum possible number of hits
};

#endif