HarpoKalmanTracking.h

Go to the documentation of this file.

 #ifndef HARPOKALMANTRACKING_H
#define HARPOKALMANTRACKING_H

//
// HarpoKalmanTracking
//
// Dummy analise which print all run events.
//

#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 HarpoKalmanTracking : public HarpoKalman {

 public:
  
  virtual ~HarpoKalmanTracking() {}

  HarpoKalmanTracking() {}

  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];}

  void SetCanvas(TVirtualPad* c, Int_t plane){fCanvas[plane] = c;}
  void SetHist(TH2F* h, Int_t plane){fHist[plane] = h;}

 private:

  Int_t      InitPlane(TClonesArray* clArray, Int_t plane);
  Int_t    GetMapEdges(TClonesArray* clArray, Int_t plane, Int_t &iMin, Int_t &iMax, Int_t &jMin, Int_t &jMax);
  void     FindTrack(TClonesArray* clArray, Int_t icl0, Int_t icl1, TMatrixD Corig, Int_t plane, Int_t &color);
  Int_t    AddTrack(TClonesArray* clArray, Int_t plane);
  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];
  TGraph* fGraphs[2][NTRACK];

  Double_t fQtot;
  Double_t fQused;


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

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

#endif