HarpoKalmanFitting.h

Go to the documentation of this file.

 #ifndef HARPOKALMANFITTING_H
#define HARPOKALMANFITTING_H

//
// HarpoKalmanFitting
//
// Track fitting with Kalman filter
//

#include "TObject.h"
#include "HarpoAnalyse.h"
#include "HarpoEvent.h"
#include "HarpoRecoEvent.h"
#include "HarpoKalman.h"
#include "TProfile.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TCanvas.h"
#include "TPolyLine.h"
#include "TMatrixD.h"
#include "TGraphErrors.h"

#define NTRACK 100


class HarpoKalmanFitting : public HarpoKalman {

 public:
  
  virtual ~HarpoKalmanFitting() {}

  HarpoKalmanFitting() {}

  void   Init();
  void   Save(char *mode=NULL);
  void   process(); 
  void   print(); 

  TH1F* hNcl;
  /* TGraph* GetGraph(Int_t i, Int_t plane = 0) {  */
  /*   if(plane<0 || plane>1) return 0; */
  /*   if(fId[i][plane]<0) return 0; */
  /*   if(i>=0 && i<fNtr[plane]) return fGraphs[plane][fId[i][plane]]; */
  /*      //if(i>=0 && i<fNtr[plane]) return fGraphs[plane][i]; */
  /*   else return 0; */
  /* } */
  Int_t GetNtr(Int_t plane){if(plane<0 || plane>1) return 0; return fNtr[plane];}
  Int_t GetId(Int_t i,Int_t plane) {return fId[i][plane];}

 private:

  Int_t      InitPlane(TClonesArray* clArray, Int_t plane);
  TMatrixD    GetTrackEnd(TClonesArray* clArray, Int_t plane, Int_t iTr);
  void     FindTrack(TClonesArray* clArray, Int_t icl0, Int_t icl1, TMatrixD Corig, Int_t plane, Int_t &color);
  Int_t    FitTrack(TClonesArray* clArray, Int_t plane, HarpoRecoKalmanTracks* track);
  //  void     SpliceTracks(Int_t plane);
  Double_t GetQtrack(Int_t itr, Int_t plane);
  Double_t GetQtracks(Int_t plane);


  Int_t fId[NTRACK][2];
  Double_t fQcommon[NTRACK][NTRACK];


  TH2F*   fHistQQtestVsDist;
  TH2F*   fHistQQmin;

  /* TMatrixD H; */
  /* TMatrixD Ht; */
  /* TMatrixD Q; */
  /* TMatrixD G; */
  /* TMatrixD V; */
  /* TMatrixD I; */
  /* TMatrixD Corig; */
  
  Int_t fNclMin; 
  Int_t fNclMin2; 


  ClassDef(HarpoKalmanFitting,1)  //Run AnalysePrint
};

#endif