FitDigitizerAlgo.h

Go to the documentation of this file.

/*
 * FitDigitizerAlgo.h
 *
 *  Created on: 20 Feb 2020
 *      Author: Valerio Formato (based on code by Junjing Wang)
 */

#ifndef HERD_FITDIGITIZERALGO_H_
#define HERD_FITDIGITIZERALGO_H_

// HerdSoftware headers
#include "algorithms/digitization/FitDigitizationTools/AttenuationTool.h"
#include "algorithms/digitization/FitDigitizationTools/DistributionChannelTool.h"
#include "algorithms/digitization/FitDigitizationTools/DistributionFiberTool.h"
#include "algorithms/digitization/FitDigitizationTools/PhotonTool.h"
#include "algorithms/digitization/FitDigitizationTools/SiPMTool.h"
#include "dataobjects/FitChannelInfoColl.h"
#include "dataobjects/FitGeoParamsColl.h"
#include "dataobjects/FitHitsColl.h"
#include "dataobjects/FitParticleHitsColl.h"

// EventAnalysis headers
#include "algorithm/Algorithm.h"

using namespace EA;

namespace Herd {

// clang-format off
// clang-format on

class FitDigitizerAlgo : public Algorithm {
public:
  FitDigitizerAlgo(const std::string &name);

  bool Initialize() override;

  bool Process() override;

private:
  StorePtr _evStore, _globStore;

  FitHitsColl _fitParticleHitsColl;
  observer_ptr<FitGeoParamsColl> _origFitGeoParamsColl;
  std::shared_ptr<FitGeoParamsColl> _digFitGeoParamsColl;
  std::shared_ptr<FitChannelInfoColl> _digFitChannelInfoColl;

  FitGeoParams _CreateGeoParams(const FitGeoParams &origGeoParams);
  FitChannelInfo _CreateChannelInfo(const FitGeoParams &digGeoParams);
  std::vector<float> _GenerateSipmGains(const FitGeoParams &digGeoParams);

  SidesArray<std::vector<std::vector<float>>> m_sipmGains;

  // Parameters
  bool _publishChannelInfo = true;
  bool _simulateNoise = true;
  bool _simulateIntraChannelXTalk = true;
  unsigned int _nTopMatsPerLayer = 11;
  unsigned int _nSideShortMatsPerLayer = 10;
  unsigned int _nSideLongMatsPerLayer = 7;

  double _sipmGainVariation = 0.0125; // relative fluctuation of the gain
  double _adcNoise = 0.1;             // Sigma of the noise in the electronics ADC conversion

  // utility tools
  std::unique_ptr<FitDigitizationTools::PhotonTool> m_photonTool;
  std::unique_ptr<FitDigitizationTools::DistributionFiberTool> m_distrFiberTool;
  std::unique_ptr<FitDigitizationTools::SiPMTool> m_sipmTool;

  std::mt19937 m_generator;
  std::normal_distribution<double> m_gauss_dist;
};

} // namespace Herd

#endif /* HERD_PSDGEOMETRICDIGITIZERALGO_H_ */