PndGemFindClusters.h

Go to the documentation of this file.

//* $Id: */

// -------------------------------------------------------------------------
// -----                     PndGemFindClusters header file                -----
// -----                  Created 15/02/2009 by R. Karabowicz          -----
// -------------------------------------------------------------------------


#ifndef PNDGEMFINDCLUSTERS_H
#define PNDGEMFINDCLUSTERS_H 1


#include "FairTask.h"

#include "TStopwatch.h"

#include <map>
#include <set>


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

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

class PndGemFindClusters : public FairTask
{


 public :

  PndGemFindClusters();


  PndGemFindClusters(Int_t iVerbose);


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


  virtual ~PndGemFindClusters();


  virtual void Exec(Option_t* opt);



 private:

  PndGemDigiPar*    fDigiPar;     
  TClonesArray*     fDigis;       
  TClonesArray*     fClusters;    
  std::map <Double_t, Int_t> fTimeOrderedDigis;
  std::vector <DigiCluster>  fDigiClusters;

  Int_t             fTNofEvents;
  Int_t             fTNofDigis;
  Int_t             fTNofClusters;

  TStopwatch fTimer;

  virtual void SetParContainers();


  virtual InitStatus Init();


  virtual InitStatus ReInit();


  Int_t SortDigis();

  //  /** Find clusters  **/
  //void FindClusters();//not used

  Int_t CreateClusters();

  Int_t WriteClusters();

  // /** Compare digi to clusters **/
  //Bool_t CompareDigiToClusters(Int_t digiNumber);

  Bool_t CompareDigiToClustersDigis(Int_t digiNumber);


  //  /** Analyze clusters **/
  //void AnalyzeClusters();//not used


  //void ClearClusters();//not used

  //void ClearClusters2();

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

  void JoinTwoClusters(Int_t clus1, Int_t clus2);

  void AddDigiToCluster(Int_t digiNr, Int_t clusNr);

  virtual void Finish();


  ClassDef(PndGemFindClusters,1);

};

#endif