KLMLikelihoodParametersImporter Class Reference

This class imports KLMLikelihoodParameters into the database. More...

#include <KLMLikelihoodParametersImporter.h>

Public Member Functions
	KLMLikelihoodParametersImporter ()=default
	Constructor.
	~KLMLikelihoodParametersImporter ()=default
	Destructor.
void	writeLikelihoodParameters ()
	Write KLMLikelihoodParameters into the database. More...
void	readLikelihoodParameters ()
	Read KLMLikelihoodParameters from the database.

Detailed Description

This class imports KLMLikelihoodParameters into the database.

Definition at line 20 of file KLMLikelihoodParametersImporter.h.

Member Function Documentation

writeLikelihoodParameters()

void writeLikelihoodParameters

(

)

Write KLMLikelihoodParameters into the database.

cppcheck-suppress unreachableCode

Definition at line 30 of file KLMLikelihoodParametersImporter.cc.

{
  B2WARNING("The method KLMLikelihoodParametersImporter::writeMuidParameters() is temporary unavailable, sorry! :(");
  return;
  DBImportObjPtr<KLMLikelihoodParameters> likelihoodParameters;
  likelihoodParameters.construct();
  std::vector<std::string> const hypotheses = {"Positron", "Electron", "Deuteron", "Antideuteron", "Proton", "Antiproton", "PionPlus", "PionMinus", "KaonPlus", "KaonMinus", "MuonPlus", "MuonMinus" };
  for (unsigned int hypothesis = 0; hypothesis < hypotheses.size(); hypothesis++) {
    GearDir content("/Detector/Muid/MuidParameters//Experiment[@exp=\"0\"]/");
    content.append(hypotheses[hypothesis]);
    for (int outcome = 1; outcome <= MuidElementNumbers::getMaximalOutcome(); ++outcome) {
      GearDir outcomeContent(content);
      outcomeContent.append((boost::format("/LayerProfile/Outcome[@outcome=\"%1%\"]/") % (outcome)).str());
      for (int lastLayer = 0; lastLayer <= MuidElementNumbers::getMaximalBarrelLayer(); ++lastLayer) {
        if (!(MuidElementNumbers::checkExtrapolationOutcome(outcome, lastLayer)))
          break;
        std::vector<double> layerPDF = outcomeContent.getArray((boost::format("LastLayer[@layer=\"%1%\"]") % (lastLayer)).str());
        likelihoodParameters->setLongitudinalPDF(hypothesis, outcome, lastLayer, layerPDF);
      }
    }
    for (int detector = 0; detector <= MuidElementNumbers::getMaximalDetector(); ++detector) {
      GearDir detectorContent(content);
      if (detector == 0)
        detectorContent.append("/TransversePDF/BarrelAndEndcap");
      if (detector == 1)
        detectorContent.append("/TransversePDF/BarrelOnly");
      if (detector == 2)
        detectorContent.append("/TransversePDF/EndcapOnly");
      for (int halfNdof = 1; halfNdof <= MuidElementNumbers::getMaximalHalfNdof(); ++halfNdof) {
        double reducedChiSquaredThreshold = detectorContent.getDouble((boost::format("DegreesOfFreedom[@ndof=\"%1%\"]/Tail/Threshold") %
                                            (2 * halfNdof)).str());
        double reducedChiSquaredScaleY = detectorContent.getDouble((boost::format("DegreesOfFreedom[@ndof=\"%1%\"]/Tail/ScaleY") %
                                                                    (2 * halfNdof)).str());
        double reducedChiSquaredScaleX = detectorContent.getDouble((boost::format("DegreesOfFreedom[@ndof=\"%1%\"]/Tail/ScaleX") %
                                                                    (2 * halfNdof)).str());
        std::vector<double> reducedChiSquaredPDF = detectorContent.getArray((boost::format("DegreesOfFreedom[@ndof=\"%1%\"]/Histogram") %
                                                   (2 * halfNdof)).str());
        likelihoodParameters->setTransversePDF(hypothesis, detector, halfNdof * 2, reducedChiSquaredPDF);
        likelihoodParameters->setTransverseThreshold(hypothesis, detector, halfNdof * 2, reducedChiSquaredThreshold);
        likelihoodParameters->setTransverseScaleY(hypothesis, detector, halfNdof * 2, reducedChiSquaredScaleY);
        likelihoodParameters->setTransverseScaleX(hypothesis, detector, halfNdof * 2, reducedChiSquaredScaleX);
      }
    }
  }
  IntervalOfValidity Iov(0, 0, -1, -1);
  likelihoodParameters.import(Iov);
}

---

The documentation for this class was generated from the following files:

• klm/calibration/include/KLMLikelihoodParametersImporter.h
• klm/calibration/src/KLMLikelihoodParametersImporter.cc