PndEmcMakeCluster.h

Go to the documentation of this file.

//-----------------------------------------------------------------------
// File and Version Information:
// $Id: $
// Software developed for the PANDA Detector at GSI.
// Author List:
//  Jan Zhong
//  Marcel Tiemens
//  Áron Kripkó
//---------------------------------------------------------------------
//#pragma once
#ifndef PNDEMCMAKECLUSTER_H
#define PNDEMCMAKECLUSTER_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 PndEmcDigiPar;
class PndEmcRecoPar;
class BinaryFunctor;

class PndEmcMakeCluster : public PndPersistencyTask {
 public:
  // Constructors
  PndEmcMakeCluster(Int_t verbose = 0, Bool_t storeclusters = kTRUE);
  // Destructor
  virtual ~PndEmcMakeCluster();

  virtual InitStatus Init();
  virtual void Exec(Option_t *opt);

  virtual void ExecuteTasks(Option_t *option)
  {
    PndPersistencyTask::ExecuteTasks(option);
    FinishClusters();
  }

  void SetStorageOfData(Bool_t val); 

  void StoreClusterBaseDigis(Bool_t val = kTRUE) { fStoreClusterBase = val; }

  void SetClusterActiveTime(Double_t time) { fClusterActiveTime = time; } 

  void SetClusterMinimumEnergy(Double_t minE) { fClusterEnergyCut = minE; }
  void EnableRemovalOfLowEnergyClusters(Bool_t enable) { fRemoveLowEclus = enable; }
  void EnableBetterNeutralReconstruction(Bool_t enable) { fMerge = enable; }

 protected:
  virtual void SetParContainers();
  virtual void FinishClusters();
  virtual void RemoveLowEnergyClusters();

 private:
  void FinishCluster(PndEmcCluster *tmpcluster);

  void cleansortmclist(std::vector<Int_t> &newlist, TClonesArray *mcTrackArray);
  // don't allow copying (-Weffc++)

  TClonesArray *fDigiArray;

  TClonesArray *fHitArray;
  TClonesArray *fMCTrackArray;

  TClonesArray *fClusterArray; 

  PndEmcGeoPar *fGeoPar;   
  PndEmcDigiPar *fDigiPar; 
  PndEmcRecoPar *fRecoPar; 
  Int_t fVerbose;

  Double_t fDigiEnergyTresholdBarrel;   
  Double_t fDigiEnergyTresholdFWD;      
  Double_t fDigiEnergyTresholdBWD;      
  Double_t fDigiEnergyTresholdShashlyk; 
  Double_t fClusterEnergyCut;           
  Double_t fMaxECut;

  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;

  std::vector<Double_t> fClusterPosParam;

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

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

  TStopwatch fTimer;

  ClassDef(PndEmcMakeCluster, 3)
};
#endif // PNDEMCMAKECLUSTER_HH