SiPM.h

Go to the documentation of this file.

#ifndef _mdet_SiPM_h
#define _mdet_SiPM_h

#include <mdet/MDetectorComponent.h>
//
#include <utl/ShadowPtr.h>
#include <utl/Validated.h>
#include <utl/Math.h>
#include <utl/MathConstants.h>
#include <utl/ConsecutiveEnumFactory.h>
//
#include <complex>
#include <cmath>


namespace det {
  struct ParentCreator;
  struct ComponentUpdater;
}

namespace mdet {
  typedef std::vector<int> IdsContainer;


  class SiPMArray;

  class SiPM :
   public MDetectorComponent<SiPM>::Type {
  public:

    static const char* const kComponentName;

    static const char* const kComponentId;

    const SiPMArray& GetSiPMArray() const {
      return fSiPMArray;
    }
    class PE {
    public:
      double TimeSpaceEvaluate(double t) const {
        /*
         * ns
         */

        double x = 0;
        if(t>=fStartTime){
                //x = fPulseShape.Eval(t);
                x = -1*fAmplitude1*(1-exp(-(t-fStartTime)/fTime1))*(fAmplitude2*exp(-(t-fStartTime)/fTime2)+
                                        fAmplitude3*exp(-(t-fStartTime)/fTime3)+exp(-(t-fStartTime)/fTime4));
                                        
                if(x>0) x = 0; //Positive values mean we are out of the signal range.
        }
        return kOverallSign * x;
      }
      double UpperLimit() const {
        return GetStartTime() + fSiPM.GetPulseRelevantWidth();
      }
      double LowerLimit() const {
        return GetStartTime();
      }

      double GetStartTime() const {
          return fStartTime;
      }

      double GetAmplitude1() const{
          return fAmplitude1;
      }

      double GetAmplitude2() const{
          return fAmplitude2;
      }

      double GetAmplitude3() const{
          return fAmplitude3;
      }

      double GetRiseTime() const{
          return fTime1;
      }

      double GetFallTime1() const{
          return fTime2;
      }

      double GetFallTime2() const{
          return fTime3;
      }

      double GetFallTime3() const{
          return fTime4;
      }

    private:
      static const int kOverallSign = 1;
      PE(const SiPM& sipm, double a1, double a2, double a3, double t1, double t2, double t3, double t4, double t0) :
        fSiPM(sipm), fAmplitude1(a1), fAmplitude2(a2), fAmplitude3(a3), fTime1(t1), fTime2(t2), fTime3(t3), fTime4(t4), fStartTime(t0) { }

      const SiPM& fSiPM;
      double fAmplitude1;
      double fAmplitude2;
      double fAmplitude3;
      double fTime1;
      double fTime2;
      double fTime3;
      double fTime4;
      double fStartTime;
      //TFormula fPulseShape;
      //double fStartTime;
      // Allow construction.
      friend class SiPM;
    };
    PE MakePEAt(double t) const;
  private:
    void Update(bool invalidateData, bool invalidateComponents);
    double GetPulseAmplitude1Mean() const;
    double GetPulseAmplitude1StdDev() const;
    double GetPulseAmplitude2Mean() const;
    double GetPulseAmplitude2StdDev() const;

    double GetPulseAmplitude3Mean() const;
    double GetPulseAmplitude3StdDev() const;
    double GetPulseTime1Mean() const;
    double GetPulseTime1StdDev() const;
    double GetPulseTime2Mean() const;
    double GetPulseTime2StdDev() const;
    double GetPulseTime3Mean() const;
    double GetPulseTime3StdDev() const;

    double GetPulseTime4Mean() const;
    double GetPulseTime4StdDev() const;

    double GetPulseRelevantWidth() const;

   // void P2X(double *p, double *x, int nTest);
    SiPM(int pId, const det::VManager::IndexMap& parentMap, const SiPMArray& parent);
    ~SiPM() { }
    const IdsContainer& GetChannelSiPMScintillatorFiberIds() const;
    friend struct det::ParentCreator;
    template<class T>
    friend void boost::checked_delete(T*) BOOST_NOEXCEPT;
    friend class SiPMArray;
    friend class Module;
    friend struct det::ComponentUpdater;

    typedef const char* PulseParametrizationForConfig;

    mutable utl::Validated<PulseParametrizationForConfig> fPulseParametrization;
    mutable utl::ShadowPtr<IdsContainer> fChannelSiPMScintillatorFiberIds;

    mutable utl::Validated<double> fPulseAmplitude1Mean;

    mutable utl::Validated<double> fPulseAmplitude1StdDev;

    mutable utl::Validated<double> fPulseAmplitude2Mean;

    mutable utl::Validated<double> fPulseAmplitude2StdDev;

    mutable utl::Validated<double> fPulseAmplitude3Mean;

    mutable utl::Validated<double> fPulseAmplitude3StdDev;

    mutable utl::Validated<double> fPulseTime1Mean;

    mutable utl::Validated<double> fPulseTime1StdDev;

    mutable utl::Validated<double> fPulseTime2Mean;

    mutable utl::Validated<double> fPulseTime2StdDev;

    mutable utl::Validated<double> fPulseTime3Mean;

    mutable utl::Validated<double> fPulseTime3StdDev;

    mutable utl::Validated<double> fPulseTime4Mean;

    mutable utl::Validated<double> fPulseTime4StdDev;

    mutable utl::Validated<double> fPulseRelevantWidth;

    const SiPMArray& fSiPMArray;
  };

}

#endif // _mdet_SiPM_h