PndGemFindClustersTB.h

Go to the documentation of this file.

//* $Id: */

// -------------------------------------------------------------------------
// -----                     PndGemFindClustersTB header file                -----
// -----                  Created 15/06/2013 by R. Karabowicz          -----
// -------------------------------------------------------------------------

#ifndef PNDGEMFINDCLUSTERSTB_H
#define PNDGEMFINDCLUSTERSTB_H 1

#include "FairTask.h"
#include "FairTSBufferFunctional.h"

#include "TStopwatch.h"

#include <map>
#include <set>

class TClonesArray;
class PndGemDigiPar;
class PndGemSensor;
class PndGemStation;

struct DigiClusterTB {
  Int_t detId;
  std::vector<Int_t> digiNr;
  std::vector<Int_t> chanNr;
  std::vector<Double_t> sigADC;
  std::vector<Double_t> sigTDC;
  Double_t cluTDC;
  Double_t cluADC;
  Double_t cluPos;
  Int_t cluPMn;
  Int_t cluPMx;
  Double_t cluMPs;
  Double_t cluMVl;
};

class PndGemFindClustersTB : public FairTask {

 public:
  PndGemFindClustersTB();

  PndGemFindClustersTB(Int_t iVerbose);

  PndGemFindClustersTB(const char *name, Int_t iVerbose = 0);

  virtual ~PndGemFindClustersTB();

  virtual void Exec(Option_t *opt);

  void SetInBranchName(TString ibn) { fInBranchName = ibn.Data(); }

 private:
  PndGemDigiPar *fDigiPar; 
  TClonesArray *fDigis;    
  TClonesArray *fClusters; 
  Int_t fTNofEvents;
  Int_t fTNofDigis;
  Int_t fTNofClusters;

  BinaryFunctor *fFunctor; 

  Double_t fPrepTime;
  Double_t fCreateTime;
  Double_t fSortTime;
  Double_t fCheckTime;
  Double_t fAnaTime;
  Double_t fWriteTime;
  Double_t fAllTime;

  TString fInBranchName;

  std::vector<DigiClusterTB> fDigiClusters;

  TStopwatch fTimer;

  virtual void SetParContainers();

  virtual InitStatus Init();

  virtual InitStatus ReInit();

  Int_t CreateClusters();

  Int_t WriteClusters();

  void AnalyzeClusters();

  void SortClusters();
  void SortCluster(Int_t clus);
  void PrintClusters();
  void PrintCluster(Int_t clus);

  void CheckClusters();
  void CheckCluster(Int_t clus);

  void JoinTwoClusters(Int_t clus1, Int_t clus2);

  void AddDigiToCluster(Int_t digiNr, Int_t clusNr);

  Bool_t CompareDigiToClusters(Int_t digiNumber);

  virtual void Finish();

  ClassDef(PndGemFindClustersTB, 1);
};

#endif