HarpoTrackingPh.h

Go to the documentation of this file.

 #ifndef HARPOTRACKINGPH_H
#define HARPOTRACKINGPH_H

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

#include "TObject.h"
#include "HarpoAnalyse.h"
#include "HarpoEvent.h"
#include "HarpoRecoEvent.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 HarpoTrackingPh : public HarpoAnalyse {

 public:
  
  virtual ~HarpoTrackingPh() {}

  HarpoTrackingPh() {}

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

  void SetQmin(Double_t val){fQmin = val;}
  void SetWidthMinT(Double_t val){fWidthMinT = val;}
  void SetWidthMinC(Double_t val){fWidthMinC = val;}
  void SetWidthMax(Double_t val){fWidthMax = val;}


 private:

  TArrayI* FindNext(TMatrixD X, TMatrixD C, Double_t angle, Int_t Ntr, TClonesArray* clArray, TArrayI* arr, Int_t ncl, Int_t plane, Int_t color, Int_t fill,Int_t smooth, Int_t q);
  TArrayI* FindNextClosest(TMatrixD X, TMatrixD C, Double_t angle, Int_t Ntr, TClonesArray* clArray, TArrayI* arr, Int_t ncl, Int_t plane, Int_t color, Int_t fill,Int_t smooth, Int_t q);
  //  TArrayI* FindNextClosest(TMatrixD X, TMatrixD C, Int_t Ncl, TClonesArray* clArray, TArrayI* arr, Int_t plane, Int_t color, Int_t fill,Int_t smooth);
  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(TArrayI* arr, TClonesArray* clArray, Int_t plane);
  Int_t    AddTrack(TClonesArray* clArray, Int_t plane);
  //  Int_t AddTrackOld(TArrayI* arr, TClonesArray* clArray, Int_t plane);
  //  void AddAllTracks(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);

  Double_t GetResolution(HarpoRecoClusters* cl);

  Int_t Ccl[NCHAN][20];
  Int_t NCcl[NCHAN];
  Int_t Tcl[NADC][20];
  Int_t NTcl[NADC];
  Int_t fNtr[2];

  Int_t fNclTrack[NTRACK];
  TGraph* fGraphs[2][NTRACK];


  Double_t fStartPointX[NTRACK][10];
  Double_t fStartPointZ[NTRACK][10];
  //  Double_t fStartAngle[NTRACK][10];
  Double_t fStartDirX[NTRACK][10];
  Double_t fStartDirZ[NTRACK][10];
  Double_t fEndPointX[NTRACK][10];
  Double_t fEndPointZ[NTRACK][10];
  Double_t fEndDirX[NTRACK][10];
  Double_t fEndDirZ[NTRACK][10];
  //  Double_t fEndAngle[NTRACK][10];
  Int_t fStartIndex[NTRACK];
  Int_t fEndIndex[NTRACK];
  Int_t fId[NTRACK][2];
  Double_t fQcommon[NTRACK][NTRACK];

  //  Double_t fMapTmp[NCHAN][NADC];
  ULong64_t fMapTmp[NCHAN][NADC];
  Double_t fQtot;
  Double_t fQused;

  Int_t used[4000];
  Int_t reused[4000];
  Int_t fTrack[4000];
  TVirtualPad* fCanvas[2];
  TH2F* fHist[2];

  TH1F* hDistTracks;
  TH1F* hDistMinTracks;
  TH1F*  hThetaTracks;
  TH2F*  hThetaTracksNtr;
  TH2F*  hThetaDist[10];

  TH2F*   fHistQQtestVsDist;
  TH2F*   fHistQQmin;

  /* TMatrixD H; */
  /* TMatrixD Ht; */
  /* TMatrixD Q; */
  /* TMatrixD G; */
  /* TMatrixD V; */
  /* TMatrixD I; */
  /* TMatrixD Corig; */
  
  Double_t fQmin;
  Double_t fWidthMinT;
  Double_t fWidthMinC;
  Double_t fWidthMax;
  Int_t fNclMin; 
  Int_t fNclMin2; 
  Double_t fMaxSlopeType;
  Double_t  fResKalman;
  Double_t  fScatKalman;
  Double_t  fResFinder;
  Double_t fScatFinder;


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

#endif