EEmcSectorFit.h

#ifndef __EEmcSectorFit_h__
#define __EEmcSectorFit_h__

#include "TMinuit.h"
#include <vector>
class TH1F;
class EEmcSectorFit : public TMinuit
{
public:
  EEmcSectorFit(Int_t maxGammas=10);

  virtual ~EEmcSectorFit();

  void print() const; 

  void SetHistograms(TH1F *u, TH1F *v){ mSMD[0]=u; mSMD[1]=v; }  

  Double_t FitFunc( Double_t x, Int_t plane ) const;

  virtual Int_t Eval(Int_t np,Double_t* gr,Double_t& x2,Double_t* p,Int_t flg); 

  Double_t Residual( Int_t x, Int_t plane) const;
  Double_t Residual( Int_t x, Int_t plane, Int_t dx, Int_t side=0) const; 

  Int_t MaxStrip(Int_t plane) const;

  void AddCandidate( Double_t yield, Double_t sigma, Double_t u, Double_t v );

  void Draw(Option_t *opts);

  void InitParameters();

  void TryPermutations();
  
  void Clear(Option_t *opts="");

  Int_t numberOfCandidates(){ return (Int_t)yield.size(); }

  void  GetCandidate( Int_t i, Double_t &nmips, Double_t &width, Double_t &u, Double_t &v ){ nmips=yield[i]; width=sigma[i]; u=umean[i]; v=vmean[i]; }
  void  GetLastCandidate( Double_t &nmips, Double_t &width, Double_t &u, Double_t &v ){ nmips=yield.back(); width=sigma.back(); u=umean.back(); v=vmean.back(); }
 
  TH1F *histo(Int_t plane) { return mSMD[plane]; }
  Double_t chi2() const { return mChi2; }
  Int_t ndf() const { return mNDF; }  

  void AddFits(TH1F *u, TH1F *v);

  Bool_t doPermutations;

protected:

  TH1F *mSMD[2];

  std::vector<Double_t> yield;
  std::vector<Double_t> sigma;
  std::vector<Double_t> umean;
  std::vector<Double_t> vmean;

  Double_t mChi2;
  Int_t mNDF; 

  ClassDef(EEmcSectorFit,1);

};

#endif