BSEmcFwEndcapBlindDeconvProcess.h

Go to the documentation of this file.

#ifndef BSEMCFWENDCAPBLINDDECOFEATUREEXTRACTIONPROCESS_HH
#define BSEMCFWENDCAPBLINDDECOFEATUREEXTRACTIONPROCESS_HH

#include <algorithm>
#include <memory>
#include <string>
#include <vector>

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

#include "PndInputContainer.h"
#include "PndOutputContainerI.h"
#include "PndProcess.h"

#include "BSEmcDigi.h"
#include "BSEmcDigitizationData.h"
#include "BSEmcMultiWaveform.h"

#include "BSEmcPSAFPGABlindDeconv.h"
#include "BSEmcPSAFPGABaselineFinder.h"
#include "BSEmcPSAFPGAPeakfinder.h"
#include "BSEmcPSAFPGAPileUpDetection.h"
#include "BSEmcPSAFPGAPeakfinderTime.h"

class PndParameterRegister;
class BSEmcAbsPSA;
class BSEmcWaveform;
class TBuffer;
class TClass;
class TMemberInspector;
struct BSEmcDigitizationData;
template <class T>
class PndOutputContainerI;

class BSEmcFwEndcapBlindDeconvProcess : public PndProcess<BSEmcDigitizationData> {
 public:
  BSEmcFwEndcapBlindDeconvProcess();
  virtual ~BSEmcFwEndcapBlindDeconvProcess();

  virtual void SetDetectorName(const std::string &t_detectorName) /*override*/;
  virtual void SetData(BSEmcDigitizationData *t_data) /*override*/;
  virtual void SetupParameters(const PndParameterRegister *t_parameterRegister) /*override*/;
  virtual void Process() /*override*/;

  std::vector<BSEmcDigi *> ExtractDigis(const std::vector<const BSEmcMultiWaveform *> &t_waveforms) const;
  void SetWaveformBranchName(const TString &t_branchname) { fWaveformBranchName = t_branchname; }

  Bool_t IsVptt(const Int_t t_detId) const;
  void SetVpttIds(const std::vector<Int_t> &t_vpttIds)
  {
    fVpttIds = t_vpttIds;
    std::sort(fVpttIds.begin(), fVpttIds.end());
  }

 protected:
  TString fHitBranchName{""};
  TString fWaveformBranchName{""};

  PndOutputContainerI<BSEmcDigi> *fDigiArray{nullptr};
  PndInputContainer<BSEmcMultiWaveform> fWaveformArray{};

  Double_t fEnergyDigiThreshold{0};

  BSEmcPSAFPGABlindDeconv fBlindDeconvoluterAPD{};
  BSEmcPSAFPGABaselineFinder fBaselineFinderAPD{};
  BSEmcPSAFPGAPeakfinder fPeakFinderAPD{};
  BSEmcPSAFPGAPeakfinderTime fPeakTimeFinderAPD{};
  BSEmcPSAFPGAPileUpDetection fPileUpDetectorAPD{};

  BSEmcPSAFPGABlindDeconv fBlindDeconvoluterVPTT{};
  BSEmcPSAFPGABaselineFinder fBaselineFinderVPTT{};
  BSEmcPSAFPGAPeakfinder fPeakFinderVPTT{};
  BSEmcPSAFPGAPeakfinderTime fPeakTimeFinderVPTT{};
  BSEmcPSAFPGAPileUpDetection fPileUpDetectorVPTT{};

  std::vector<Int_t> fVpttIds{};

  ClassDef(BSEmcFwEndcapBlindDeconvProcess, 1);
};

#endif /*BSEMCFWENDCAPBLINDDECOFEATUREEXTRACTIONPROCESS_HH*/