ginternalDigitization.h

Go to the documentation of this file.

#pragma once
 
#include "gdynamicdigitization.h"  // Base class for dynamic digitization.
#include <map>
#include <vector>
#include <string>
 
class GFluxDigitization : public GDynamicDigitization {
public:
    // inherit the base (const std::shared_ptr<GOptions>&) ctor
    using GDynamicDigitization::GDynamicDigitization;
    bool defineReadoutSpecsImpl() override;
 
    std::unique_ptr<GDigitizedData> digitizeHitImpl(GHit* ghit, size_t hitn) override;
};
 
class GParticleCounterDigitization : public GDynamicDigitization {
public:
    // inherit the base (const std::shared_ptr<GOptions>&) ctor
    using GDynamicDigitization::GDynamicDigitization;
    bool defineReadoutSpecsImpl() override;
 
    std::unique_ptr<GDigitizedData> digitizeHitImpl(GHit* ghit, size_t hitn) override;
};
 
class GDosimeterDigitization : public GDynamicDigitization {
public:
    // inherit the base (const std::shared_ptr<GOptions>&) ctor
    using GDynamicDigitization::GDynamicDigitization;
    bool defineReadoutSpecsImpl() override;
 
    std::unique_ptr<GDigitizedData> digitizeHitImpl(GHit* ghit, size_t hitn) override;
 
    bool loadConstantsImpl(int runno, std::string const &variation) override;
 
private:
    // Calibration maps (keys are particle IDs).
    // Currently:
    // 11 (electrons), 211 (pions), 2112 (neutrons), 2212 (protons)
    std::map<int, std::vector<double>> nielfactorMap;
    std::map<int, std::vector<double>> E_nielfactorMap;
    std::map<int, double> pMassMeV; 
 
    double getNielFactorForParticleAtEnergy(int pid, double energy);
};