1 #ifndef __JFIT_JREGRESSORHELPER__
2 #define __JFIT_JREGRESSORHELPER__
20 namespace JPP {
using namespace JFIT; }
36 template<
class JModel_t>
38 template<
class JModel_t>
40 template<
class JModel_t>
47 template<
class JModel_t,
template<
class>
class JMinimiser_t>
70 static const int NUMBER_OF_PDFS = 6;
98 template<
class JModel_t,
template<
class>
class JMinimiser_t>
132 const double epsilon = 1.0e-10)
137 const JPDF_t::JSupervisor supervisor(
new JPDF_t::JDefaultResult(
JMATH::zero));
139 for (
int i = 0;
i != NUMBER_OF_PDFS; ++
i) {
141 const string file_name =
getFilename(fileDescriptor, pdf_t[
i]);
143 _pdf[
i].load(file_name.c_str());
145 _pdf[
i].setExceptionHandler(supervisor);
150 for (
int i = 1;
i < NUMBER_OF_PDFS;
i += 2) {
152 _pdf[
i ].add(_pdf[
i-1]);
156 _pdf[
i-1].swap(buffer);
196 for (
int i = 0;
i != NUMBER_OF_PDFS; ++
i) {
197 _pdf[
i].transform(transformer);
198 _npe[
i].transform(transformer);
212 for (
int i = 0;
i != NUMBER_OF_PDFS; ++
i) {
213 if (!_pdf[
i].empty()) {
214 _pdf[
i].compress(range);
JRegressorStorage()
Default constructor.
void transform(const transformer_type &transformer)
Transform PDFs and NPEs.
Template data structure for PDF tables.
const JNPEs_t & getNPE() const
Get NPEs.
direct light from EM showers
void blur(const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10, const double quantile=0.99)
Blur PDF.
JPDF_t::transformer_type transformer_type
Data structure for fit of straight line in positive z-direction.
static const JZero zero
Function object to assign zero value.
Template data structure for storage for PDF tables.
Various implementations of functional maps.
Numbering scheme for PDF types.
JTOOLS::JSplineFunction1S_t JFunction1D_t
scattered light from muon
JPHYSICS::JNPETable< double, double, JNPEMaplist_t > JNPE_t
time integrated PDF
const JPDFs_t & getPDF() const
Get PDFs.
scattered light from delta-rays
JRegressorStorage(const std::string &fileDescriptor, const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10)
Constructor.
scattered light from EM showers
std::array< JNPE_t, NUMBER_OF_PDFS > JNPEs_t
NPEs.
transformablemultifunction_type::transformer_type transformer_type
Fit method based on the Levenberg-Marquardt method.
direct light from delta-rays
z range($ZMAX-$ZMIN)< $MINIMAL_DZ." fi fi typeset -Z 4 STRING typeset -Z 2 FLOOR JPlot1D -f $
JPHYSICS::JPDFTable< JFunction1D_t, JPDFMaplist_t > JPDF_t
time dependent PDF
Auxiliary class to define a range between two values.
Simple fit method based on Powell's algorithm, see reference: Numerical Recipes in C++...
std::string getFilename(const std::string &file_name)
Get file name part, i.e. part after last JEEP::PATHNAME_SEPARATOR if any.
JTOOLS::JMAPLIST< JTOOLS::JPolint1FunctionalMap, JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JPolint0FunctionalGridMap >::maplist JPDFMaplist_t
std::array< JPDF_t, NUMBER_OF_PDFS > JPDFs_t
PDFs.
JTOOLS::JMAPLIST< JTOOLS::JPolint1FunctionalMapH, JTOOLS::JPolint1FunctionalGridMap, JTOOLS::JPolint1FunctionalGridMap >::maplist JNPEMaplist_t
void setRmax(const double Rmax)
Set maximal road width of PDF.