Circe.h

Go to the documentation of this file.

#ifndef TRACKFIT_CIRCE_H
#define TRACKFIT_CIRCE_H
#include <vector>
#include "Math/Minimizer.h"
namespace recobase { 
  class CellHit; 
  class Prong; 
}

namespace trk
{
  class Measurement {
  public:
    std::vector<const recobase::CellHit*> fHits;  
    std::vector<double>                   fW;     
  };
  
  class Circe : public ROOT::Math::IMultiGenFunction {
  public:
    static const unsigned int kNmax = 20;
  public:
    Circe() { }

    void   SeedProngs(const Measurement* meas, 
                      unsigned int       nprong,
                      double*            seed) const;
    double Fit(int nprog, const Measurement* m);
    void   MakeProngs(std::vector<recobase::Prong>& prong);

    // Complete the IMultiGenFunction interface
    unsigned int                   NDim()                  const;
    double                         DoEval(const double* p) const;
    ROOT::Math::IMultiGenFunction* Clone()                 const;
    
  private:
    const Measurement* fM;            
    unsigned int       fN;            
    double             fVtx[3];       
    double             fTheta[kNmax]; 
    double             fPhi[kNmax];   
    std::vector<int>   fProngAssn;    
  };
}

#endif

---