PndEmcMakeClusterOnline.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 PNDEMCMAKECLUSTERONLINE_H
#define PNDEMCMAKECLUSTERONLINE_H

#include <PndPersistencyTask.h>
#include <vector>       

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

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

class PndEmcMakeClusterOnline : public PndPersistencyTask
{
public:

  // Constructors

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

  // Destructor

  virtual ~PndEmcMakeClusterOnline( );

  virtual InitStatus Init();

  virtual void Exec(Option_t* opt);

  virtual void FinishTask();

  void SetStorageOfData(Bool_t val);    

  void StoreClusterBaseDigis(Bool_t val = kTRUE) { fStoreClusterBase = val; }
 
  void SetTimebunchCutTime(Double_t time) { fTimebunchCutTime = time; } 
  void SetClusterActiveTime(Double_t dt) { fClusterActiveTime = dt; }

  void SetClusterMinimumEnergy(Double_t minE) { fClusterEnergyCut = minE; }
  void EnableRemovalOfLowEnergyClusters(Bool_t enable) {fRemoveLowEclus = enable;} 
  void SetAutoDetermineTimecuts(Bool_t autodet) {fAutoDetermine = autodet;}

protected:
  virtual void FinishClusters();
  void FinishCluster(PndEmcCluster* tmpcluster);
  virtual void RemoveLowEnergyClusters();
  
private:
    TClonesArray* fDigiArray;

    TClonesArray* fClusterArray;    

    PndEmcGeoPar*     fGeoPar;       
    PndEmcRecoPar*    fRecoPar;      
    Double_t fDigiEnergyTresholdBarrel; 
    Double_t fDigiEnergyTresholdFWD; 
    Double_t fDigiEnergyTresholdBWD; 
    Double_t fDigiEnergyTresholdShashlyk;
    Double_t fClusterEnergyCut; 
    Double_t fTimebunchCutTime; // Defines how long clusters are kept open in timebased reconstruction
    Double_t fClusterActiveTime; 
    BinaryFunctor* fDigiFunctor;

    Int_t fNrOfEvents;
    Int_t fNrOfDigis;
    Int_t nOnlProg;
    Int_t digiCounter;
    Int_t evtCounter;
    Double_t OtotRtime;
    Double_t OtotCtime;

    std::vector<Double_t> fClusterPosParam;

    Bool_t fStoreClusters;
    Bool_t fStoreClusterBase;   //re-store digis in case of a timebased run
    Bool_t fRemoveLowEclus;
    Bool_t fAutoDetermine;

    TH1I* hClusMultiplicity;
    TH1I* hEventsPerCluster;
    TH1D* hTimeDifference;

    virtual void SetParContainers();

    TStopwatch fTimer;

    ClassDef(PndEmcMakeClusterOnline, 1)
};
#endif // PNDEMCMAKECLUSTERONLINE_HH