SiPM.cc

Go to the documentation of this file.

#include <mdet/SiPM.h>

#include <mdet/MHierarchyInfo.h>
//
#include <utl/RandomEngine.h>
#include <utl/MathConstants.h>
//
#include <fwk/RandomEngineRegistry.h>
//
#include <cmath>
#include <CLHEP/Random/RandGauss.h>
//
//#include <CLHEP/RandomObjects/RandMultiGauss.h>
using CLHEP::RandGauss;

namespace mdet {

  const char* const
  SiPM::kComponentName = MHierarchyInfo::kComponentsNames[10];

  const char* const
  SiPM::kComponentId   = MHierarchyInfo::kComponentsIds[10];


  SiPM::SiPM
    (int pId, const det::VManager::IndexMap& parentMap, const SiPMArray& parent) :
      MDetectorComponent<SiPM>::Type(pId, parentMap),
      fSiPMArray(parent)
  {
    Register(fPulseAmplitude1Mean);
    Register(fPulseAmplitude1StdDev);
    Register(fPulseAmplitude2Mean);
    Register(fPulseAmplitude2StdDev);
    Register(fPulseAmplitude3Mean);
    Register(fPulseAmplitude3StdDev);
    Register(fPulseTime1Mean);
    Register(fPulseTime1StdDev);
    Register(fPulseTime2Mean);
    Register(fPulseTime2StdDev);
    Register(fPulseTime3Mean);
    Register(fPulseTime3StdDev);
    Register(fPulseTime4Mean);
    Register(fPulseTime4StdDev);
    Register(fPulseRelevantWidth);
  }


  const IdsContainer&
  SiPM::GetChannelSiPMScintillatorFiberIds()
    const
  {
    // see the other getters.
    return GetData<IdsContainer,utl::ThrowOnZeroDereference,utl::ShadowPtr>(fChannelSiPMScintillatorFiberIds, "channelScintillatorFiberIds");
  }


  double
  SiPM::GetPulseAmplitude1Mean()
    const
  {
    return GetData(fPulseAmplitude1Mean, "pulseAmplitude1Mean");
  }

  double
  SiPM::GetPulseAmplitude1StdDev()
    const
  {
    return GetData(fPulseAmplitude1StdDev, "pulseAmplitude1StdDev");
  }

  double
  SiPM::GetPulseAmplitude2Mean()
    const
  {
    return GetData(fPulseAmplitude2Mean, "pulseAmplitude2Mean");
  }

  double
  SiPM::GetPulseAmplitude2StdDev()
    const
  {
    return GetData(fPulseAmplitude2StdDev, "pulseAmplitude2StdDev");
  }

  double
  SiPM::GetPulseAmplitude3Mean()
    const
  {
    return GetData(fPulseAmplitude3Mean, "pulseAmplitude3Mean");
  }

  double
  SiPM::GetPulseAmplitude3StdDev()
    const
  {
    return GetData(fPulseAmplitude3StdDev, "pulseAmplitude3StdDev");
  }

  double
  SiPM::GetPulseTime1Mean()
    const
  {
    return GetData(fPulseTime1Mean, "pulseTime1Mean");
  }

  double
  SiPM::GetPulseTime1StdDev()
    const
  {
    return GetData(fPulseTime1StdDev, "pulseTime1StdDev");
  }

  double
  SiPM::GetPulseTime2Mean()
    const
  {
    return GetData(fPulseTime2Mean, "pulseTime2Mean");
  }

  double
  SiPM::GetPulseTime2StdDev()
    const
  {
    return GetData(fPulseTime2StdDev, "pulseTime2StdDev");
  }

  double
  SiPM::GetPulseTime3Mean()
    const
  {
    return GetData(fPulseTime3Mean, "pulseTime3Mean");
  }

  double
  SiPM::GetPulseTime3StdDev()
    const
  {
    return GetData(fPulseTime3StdDev, "pulseTime3StdDev");
  }

  double
  SiPM::GetPulseTime4Mean()
    const
  {
    return GetData(fPulseTime4Mean, "pulseTime4Mean");
  }

  double
  SiPM::GetPulseTime4StdDev()
    const
  {
    return GetData(fPulseTime4StdDev, "pulseTime4StdDev");
  }

  double
  SiPM::GetPulseRelevantWidth()
    const
  {
    return GetData(fPulseRelevantWidth, "pulseRelevantWidth");
  }


  /*If correlation between parameters wants to be considered, it is necessary to load the parameters in PC space,
   * sort the values and the anti transform to the parameters space
   * (only the eigen vector of the covariance matrix are needed). Not rocket science, but a waste of resources,
   * as it has been shown, considering the correlations or not, give the same description of data.
   * */
  SiPM::PE
  SiPM::MakePEAt(double t)
    const
  {
    utl::RandomEngine& rnd = fwk::RandomEngineRegistry::GetInstance().Get(fwk::RandomEngineRegistry::eDetector);

    double *parameterMeans = new double[7];
    double *parameterSigmas = new double[7];

    parameterMeans[0]  = GetPulseAmplitude1Mean(); //A1
    parameterSigmas[0] = GetPulseAmplitude1StdDev();
    parameterMeans[1]  = GetPulseTime1Mean();// t1
    parameterSigmas[1] = GetPulseTime1StdDev();
    parameterMeans[2]  = GetPulseAmplitude2Mean(); //A2
    parameterSigmas[2] = GetPulseAmplitude2StdDev();
    parameterMeans[3]  = GetPulseTime2Mean();//t2
    parameterSigmas[3] = GetPulseTime2StdDev();
    parameterMeans[4]  = GetPulseAmplitude3Mean(); //A3
    parameterSigmas[4] = GetPulseAmplitude3StdDev();
    parameterMeans[5]  = GetPulseTime3Mean();//t3
    parameterSigmas[5] = GetPulseTime3StdDev();
    parameterMeans[6]  = GetPulseTime4Mean();//t4
    parameterSigmas[6] = GetPulseTime4StdDev();

    double *parameters = new double[7];
    for (int i = 0; i < 7; i++){
        if(parameterSigmas[i]!=0.0){
                parameters[i] = CLHEP::RandGauss::shoot(& rnd.GetEngine(), parameterMeans[i], parameterSigmas[i]);
        }else{
                parameters[i] = parameterMeans[i];
        }
    }

 //   P2X(principalComponents, parameters, 5); //get correlated parameters


    return PE(*this, parameters[0], parameters[2], parameters[4], parameters[1], parameters[3], parameters[5], parameters[6], t);
  }

  void
  SiPM::Update(bool invalidateData, bool invalidateComponents)
  {
    MDetectorComponent<SiPM>::Type::Update(invalidateData, invalidateComponents);
  }

}