HarpoVertexing.h

Go to the documentation of this file.

#ifndef HARPOVERTEXING_H
#define HARPOVERTEXING_H

//
// HarpoVertexing
//
// Dummy analysis which creates basic test histograms
//

#include "TObject.h"
#include "HarpoAnalyse.h"
#include "HarpoEvent.h"
#include "TProfile.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TGraph2D.h"
#include "TGNumberEntry.h"
#include "TGraphErrors.h"

class HarpoVertexing : public HarpoAnalyse {

 public:
  
  HarpoVertexing() {}

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


  Double_t GetFraction(Int_t ndet, Double_t x, Double_t z, Double_t &phi0, Double_t &phi1);
  void GetConnex(Int_t i, Int_t j, Int_t n, Int_t n0, HarpoDccMap* m, TH2F* h);
  //  TH1F* GetCircle(TH2F* hMap, Double_t x, Double_t z, Double_t phi0, Double_t radius);

  //  TGraph* ReduceData(Int_t ndet, Int_t rebin, Double_t &z0, Double_t &x0);
  /* void GetThrust(Int_t nPoints, Double_t* z, Double_t* x, Double_t &tz, Double_t &tx); */
  /* void GetThrust2(Int_t nPoints, Double_t* z, Double_t* x, Double_t* q, Double_t &zM, Double_t &xM, Double_t &tz, Double_t &tx); */
  //  void GetMoments(Int_t plane, Double_t &xM, Double_t &zM, Double_t &xxM, Double_t &zzM, Double_t &xzM, Double_t &xxxM, Double_t &zzzM, Double_t &xxzM, Double_t &xzzM);
  void DisplayAnalysis(TRootEmbeddedCanvas* ecTab, TGListBox* infobox);
  void ConfigFrame(TGMainFrame* fMain, Int_t id);
  void SetConfig();

  void FindVertexes();
  void RefineVertexPosition();
  void FillVertexAngles();

  TArrayD* GetPeaks2(TH1F* h);
  TArrayD* GetPeaks(Int_t n, Double_t *x, Double_t *w, Int_t plane, Int_t v);
  void     GetPeak(Int_t n, Double_t *x, Double_t *w, Double_t &muOut, Double_t &sigOut, Int_t plane, Int_t v);

  void FitVertex(TH2F* h, HarpoRecoVertex* vertex);
  TGraph* GetGraph(TH2F* h, Int_t n, Double_t x0, Double_t z0, Double_t px, Double_t pz);
  TArrayD* FitGraph(TGraph* g, const Int_t nFit);



 private:


  TGNumberEntry* fChooseRadiusX;
  TGNumberEntry* fChooseRadiusZ;
  TGNumberEntry* fChooseMinFrac;
  TGNumberEntry* fChooseSeparation;
  TGNumberEntry* fChooseQmin;

  Double_t fRadiusX;
  Double_t fRadiusZ;
  Double_t fMinFrac;
  Double_t fSeparation;
  Double_t fQmin;

  TH1F* fHist[2][10][3];
  TH1F* fHist2[2][10][3];

  TH2F* fHistProj2Draw[20];
  TH2F* fHistProj2D[20];
  TH1F* fHistMean[20];
  TH2F* fHistMap[20];
  TH1F* fHistDiff;

  Int_t fNbins;
  Int_t fNbinsTheta;

  Bool_t fDisplay;
  ClassDef(HarpoVertexing,1)  //Run AnalysePrint
};

#endif