ScdDriftDiffusionAlgo.h

Go to the documentation of this file.

 
/*
 * ScdDriftDiffusionAlgo.h
 *
 *  Created on: 23 Sep 2020
 *      Author: Valerio Formato
 */
 
#ifndef HERD_SCDDRIFTDIFFUSIONALGO_H_
#define HERD_SCDDRIFTDIFFUSIONALGO_H_
 
// HerdSoftware headers
#include "herd/analysis/dataobjects/ScdGeoParamsColl.h"
#include "herd/analysis/dataobjects/ScdHitsColl.h"
#include "herd/analysis/dataobjects/ScdParticleHitsColl.h"
 
// EventAnalysis headers
#include "algorithm/Algorithm.h"
 
namespace Herd {
 
// clang-format off
// clang-format on
 
namespace ScdDigitizationTools {
class ScdChargePacket;
}
 
class ScdDriftDiffusionAlgo : public EA::Algorithm {
public:
  ScdDriftDiffusionAlgo(const std::string &name);
 
  bool Initialize() override;
 
  bool Process() override;
 
  std::vector<std::string> FreeObjects(const std::vector<std::string> &objs, EA::Memory::Status memStatus) override;
 
protected:
  // Parameters
  float m_pitch;
  float m_minStepSize;
 
  // diffusion parameters
  std::string m_diffusionType;
 
  // - Gaussian drift parameters
  float m_posEffDiffCoeff;
  float m_negEffDiffCoeff;
 
  // - non-linear Gaussian drift parameters
  float m_biasVoltage;
  float m_intrinsicVoltage;
 
private:
  enum class DiffusionType { Gaussian, GaussianNonLinear };
  DiffusionType m_diffType = DiffusionType::GaussianNonLinear;
 
  ScdHitsColl m_siliconDetectorDigHitsColl;
  EA::observer_ptr<ScdGeoParamsColl> m_origScdGeoParamsColl;
  std::shared_ptr<ScdGeoParamsColl> m_digScdGeoParamsColl;
 
  EA::StorePtr m_evStore, m_globStore;
 
  // Helper functions
 
  void GenerateHits(const ScdParticleHitsColl &particleHitsColl, const ScdHitsColl &hitsColl, ScdHitsColl &digHitsColl);
 
  ScdGeoParams CreateGeoParams(const ScdGeoParams &origGeoParams);
 
  std::vector<Point> PartitionTrajectory(const ParticleHit &hit);
 
  void AccumulateSteps(std::vector<Hit> &layerHits, const ScdGeoParams &geoParams, unsigned int firstStrip,
                       unsigned int lastStrip, const ScdDigitizationTools::ScdChargePacket &packet, float edep);
 
  // Helper function to generate different kinds of charge packets
  std::unique_ptr<ScdDigitizationTools::ScdChargePacket> CreatePacket(Point position, DiffusionType diffType,
                                                                      const ScdGeoParams &geoParams);
};
 
} // namespace Herd
 
#endif /* HERD_SCDDriftDiffusionALGO_H_ */