HarpoMatching.h

Go to the documentation of this file.

#ifndef HARPOMATCHING
#define HARPOMATCHING

//
// NewHarpoMatching
//
// Generic class for HARPO Detectors 
//
#include <cstring>
#include "TMatrixD.h"
#include "TObject.h"
#include "TH2D.h"
#include "TH2F.h"
#include "TH1F.h"
#include "TVirtualPad.h"
#include "TGraph2D.h"
#include "TGNumberEntry.h"
#include "HarpoConfig.h"
#include "HarpoAnalyse.h"

/* #define facXY 1.075 */
/* #define facTronc 0.6 */
/* #define facRec 0.8 */
/* #define chi2TroncNcut 0.5 */
//#define NADC 511
//#define NCHAN 288
//#define NCLUSTERMAX 3000
//#define CUTCLUSTER 500;
#define NTESTMAT 50
#define NTRACK 100


class HarpoMatching : public HarpoAnalyse
{
  
 public:

  HarpoMatching(){}

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



  //fontion
  /* void ResetMatching(); */
  void Matching(TH1F* hX, TH1F* hY, Double_t &coverageX, Double_t &coverageY, Double_t &chi2);

  Int_t InitEventNew();
  /* void MatchingNew(Int_t i1, Int_t i2); */
  /* void AddGoodMatch(Int_t i1, Int_t i2, Int_t trX, Int_t trY, Double_t chi2); */



  TH1* GetHist(Int_t i, Int_t match = -1);
  void DisplayAnalysis(TRootEmbeddedCanvas* ecTab,TGListBox* info);
  void DrawMatching(HarpoRecoMatching* match, TVirtualPad* c);
  void ConfigFrame(TGMainFrame* fMain, Int_t id);
  void SetConfig();

  //  Int_t GetNgoodmatch(){return ngoodmatch;}

  /* Int_t QtTx[NTRACK]; */
  /* Int_t QtTy[NTRACK]; */
  /* Int_t iXtestmatch[ NTESTMAT],iYtestmatch[NTESTMAT];//, indgoodm[NTRACK]; */
  /* Double_t  chi2TroncN[50],/\* DecTr[NTRACK],*\/ coveX[ NTESTMAT],coveY[ NTESTMAT]; */
  /* Double_t QtTrX[NADC][NTRACK], QtTrY[NADC][NTRACK]; */
  /* Int_t iMatchX[NTRACK], iMatchY[NTRACK]; */
  /* Int_t iMatchTrX[NTESTMAT], iMatchTrY[NTESTMAT]; */
  /* Int_t mTrX, mTrY; */
  /* Int_t nbinNZTrX[NTRACK]; */
  /* Int_t nbinNZTrY[NTRACK]; */
  /* Int_t Qtrack_x[NTRACK],Qtrack_y[NTRACK]; */
  /* Double_t  Qtotal_x[NTRACK][NCLUSTERMAX],Qtotal_y[NTRACK][NCLUSTERMAX]; */
  // Double_t ThetaX,ThetaY;

  /* TGraph2D* GetGraph3D(Int_t i){if(i<ngoodmatch) return fGraph3D[i]; return 0;} */


 private:

  void DoMatchingNew();
  void DoMatching();

  void MakeCompoundTracks(Int_t iTr, ULong64_t mask, Int_t plane);
  Int_t GetNcom(ULong64_t iTr, ULong64_t mask, Int_t plane);
  TH1F* GetCompoundHist(ULong64_t mask, Int_t plane);
  Int_t CombineTracks(ULong64_t mask, Int_t plane);

  Double_t FacXY;
  Double_t FacTronc;
  Double_t FacRec;
  Double_t Chi2TroncNcut;
  Int_t fNcomMin;
  TGNumberEntry* fChooseFacXY;
  TGNumberEntry* fChooseFacTronc;
  TGNumberEntry* fChooseFacRec;
  TGNumberEntry* fChooseChi2TroncNcut;
  TGNumberEntry* fChooseNcomMin;

  TH2F* hQxyVsNcom;
  TH1F* hCoverageX;
  TH1F* hCoverageY;
  TH1F* hChi2;
  TH2F* hQxyVsNcomEvent;
  TH1F* hCoverageXEvent;
  TH1F* hCoverageYEvent;
  TH1F* hChi2Event;
  TH1F* hQvsTX[NTRACK];
  TH1F* hQvsTY[NTRACK];

  static const Int_t kNmaxCompoundTracks = NTRACK;
  ULong64_t fCompoundTrack[2][kNmaxCompoundTracks];
  Int_t fNcompoundTracks[2];
  Int_t fTrMax[2];
  ULong64_t fMask[NTRACK][2];

  /* TGraph2D* fGraph3D[NTESTMAT]; */
  //Int_t fMatchId[NTRACK][NTRACK];

  ClassDef(HarpoMatching,1)  //Event structure
    };

#endif