PndEmcMergePreclusters.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

//-----------------------------------------------------------------------
// File and Version Information:
// $Id: $
// 
//---------------------------------------------------------------------
//#pragma once
#ifndef PNDEMCMERGEPRECLUSTERS_H
#define PNDEMCMERGEPRECLUSTERS_H

#include "FairTask.h"
#include <vector>       

#include "FairTSBufferFunctional.h"
#include "TStopwatch.h"
#include "TH1.h"

class TClonesArray;
class TObjectArray;
class PndEmcDigi;
class PndEmcCluster;
class PndEmcGeoPar;
class PndEmcRecoPar;
class BinaryFunctor;


class PndEmcMergePreclusters : public FairTask
{
public:

  // Constructors

  PndEmcMergePreclusters(Int_t verbose=0, Bool_t storeclusters=kTRUE);

  // Destructor

  virtual ~PndEmcMergePreclusters( );

  virtual InitStatus Init();

  virtual void Exec(Option_t* opt);

  virtual void FinishTask();

  void SetStorageOfData(Bool_t val);    

  void StoreClusterBaseDigis(Bool_t val = kTRUE) { fStoreClusterBase = kTRUE; }
 
  void SetTimebunchCutTime(Double_t time) { fTimebunchCutTime = time; fAutoTime=kFALSE; } 
  void SetPreclusterTimeThreshold(Double_t dt) { fClusterActiveTime = dt; }
  void SetAutoDetermineTimecuts(Bool_t autotime) {fAutoTime=!autotime; }

  void SetClusterMinimumEnergy(Double_t minE) { fClusterEnergyCut = minE; }
  void EnableRemovalOfLowEnergyClusters(Bool_t enable) {fRemoveLowEclus = enable;} 
  void SetPositionMethod(Int_t method) { fPosMethod = method; } 
  void SetNeighbourMethod(Int_t nbmethod) { fNbMethod = nbmethod; } 
protected:
  virtual void FinishClusters();
  void FinishCluster(PndEmcCluster* tmpcluster);
  virtual void RemoveLowEnergyClusters();
  
private:
    TClonesArray* fDigiArray;
    TClonesArray* fPreclusterArray;

    TClonesArray* fClusterArray;    

    PndEmcGeoPar*     fGeoPar;       
    PndEmcRecoPar*    fRecoPar;      
    Double_t fClusterEnergyCut;
    
    Double_t fTimebunchCutTime; 
    Double_t fClusterActiveTime; 
    BinaryFunctor* fClusterFunctor;

    std::vector<Double_t> fClusterPosParam;

    Int_t fPosMethod;
    Int_t fNbMethod;
    Int_t evtCounter;
    Int_t precCounter;
    Int_t fNrOfPres;
    Int_t nMrgProg;
    Int_t nTotClusters;
    Int_t fRemovedClusters;
    Int_t nTotDigis;
    Int_t wrongConnection;
    Double_t CNtotRtime;
    Double_t CNtotCtime;

    TH1D* hDx;
    TH1D* hDy;
    TH1D* hDz;
    TH1D* hTimeDifference;

    Bool_t fStoreClusters;
    Bool_t fStoreClusterBase;   //restore digis in case of a timebased run
    Bool_t fRemoveLowEclus;
    Bool_t fAutoTime; // toggle manual overwriting of timebunch cut time

    virtual void SetParContainers();

    TStopwatch fTimer;

    ClassDef(PndEmcMergePreclusters, 2)
};
#endif // PNDEMCMERGEPRECLUSTERS_HH