beamtime/cern-oct11/go4/RICH/CbmRichRingFitterBase.h (r4864/r3566)

Go to the documentation of this file.

#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