BSEmcOnlineBumpSplitter.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//  $Id:$
//
// Description:
//  Class BSEmcOnlineBumpSplitter
//      Implementation of OnlineBumpSplitter which splits clusters based on
//      local maxima, determined in a 0th order approximation, i.e. by looking
//      for the two most energetic digis and checking if they are neighbours.
//
// Environment:
//  Software developed for the PANDA experiment at GSI.
//
// Author List:
//      Marcel Tiemens
//------------------------------------------------------------------------
//#pragma once
#ifndef BSEMCONLINEBUMPSPLITTER_HH
#define BSEMCONLINEBUMPSPLITTER_HH

#include <map>
#include <memory>
#include <string>
#include <vector>

#include "Rtypes.h"
#include "RtypesCore.h"

#include "FairTask.h"

#include "PndContainerI.h"
#include "PndProcess.h"

#include "BSEmcClusteringData.h"
#include "BSEmcCrystalPositionPar.h"
#include "BSEmcDigi.h"
#include "BSEmcGeoNeighbouringRelationPar.h"
#include "BSEmcPositionAlgo.h"
#include "BSEmcSubCluster.h"

class BSEmcCluster;
class BSEmcCrystalPositionPar;
class BSEmcDigi;
class BSEmcGeoNeighbouringRelationPar;
class BSEmcSubCluster;
class PndParameterRegister;
class TBuffer;
class TClass;
class TMemberInspector;
struct BSEmcClusteringData;
template <class T>
class PndContainerI;
template <class T>
class PndOutputContainerI;

class BSEmcOnlineBumpSplitter : public PndProcess<BSEmcClusteringData> {
 public:
  BSEmcOnlineBumpSplitter();
  virtual ~BSEmcOnlineBumpSplitter();
  virtual void SetDetectorName(const std::string &t_detectorName) /*override*/;
  virtual void SetData(BSEmcClusteringData *t_data) /*override*/;
  virtual void SetupParameters(const PndParameterRegister *t_paramRegister) /*override*/;
  virtual void Process() /*override*/;
  virtual void PostProcess();
  std::vector<BSEmcSubCluster *> Subcluster(const std::vector<const BSEmcDigi *> &t_digiCont, const std::vector<BSEmcCluster *> &t_clusterCont);

  void SetMinimumSeedEnergy(Double_t minE) { fSeedEnergyThreshold = minE; };

  void SetPositionParName(const std::string &t_parName) { fPositionParName = t_parName; }
  void SetNeighbouringRelationParName(const std::string &t_parName) { fNeighbouringRelationParName = t_parName; }

 protected:
  BSEmcSubCluster *AddBump();
  void SingleBumpCluster(BSEmcCluster *, Int_t, const std::vector<const BSEmcDigi *> &t_digiCont);

 private:
  PndContainerI<BSEmcDigi> *fDigiArray{nullptr};
  PndOutputContainerI<BSEmcCluster> *fClusterArray{nullptr};
  PndOutputContainerI<BSEmcSubCluster> *fSubClusterArray{nullptr};
  std::vector<BSEmcSubCluster *> fSubClusterVec{};
  std::string fNeighbouringRelationParName{""};
  BSEmcGeoNeighbouringRelationPar *fNeighbouringRelationPar{nullptr};
  std::string fPositionParName{""};
  BSEmcCrystalPositionPar *fPositionPar{nullptr};
  std::unique_ptr<BSEmcPositionAlgo> fPositionProcess{new BSEmcPositionAlgo};
  Double_t fSeedEnergyThreshold{0.03}; 

  ClassDef(BSEmcOnlineBumpSplitter, 1);
};
#endif /*BSEMCONLINEBUMPSPLITTER_HH*/