49 const double epsilon = 1.0e-10)
61 const int N =
sizeof(pdf_t) /
sizeof(pdf_t[0]);
65 const JPDF_t::JSupervisor supervisor(
new JPDF_t::JDefaultResult(
zero));
67 for (
int i = 0; i !=
N; ++i) {
69 const string file_name =
getFilename(fileDescriptor, pdf_t[i]);
71 cout <<
"loading input from file " << file_name <<
"... " << flush;
73 pdf[i].load(file_name.c_str());
75 pdf[i].setExceptionHandler(supervisor);
82 cout <<
"adding PDFs... " << flush;
92 cout <<
"bluring PDFs... " << flush;
100 }
else if (TTS < 0.0) {
102 THROW(JValueOutOfRange,
"Illegal value of TTS [ns]: " << TTS);
124 const double t1)
const
129 pdfB(R, theta, phi, t1) * E +
183 const double epsilon = 1.0e-10)
191 const int N =
sizeof(pdf_t) /
sizeof(pdf_t[0]);
193 const JPDF_t::JSupervisor supervisor(
new JPDF_t::JDefaultResult(
zero));
195 for (
int i = 0; i !=
N; ++i) {
197 const string file_name =
getFilename(fileDescriptor, pdf_t[i]);
199 cout <<
"loading input from file " << file_name <<
"... " << flush;
203 pdf.
load(file_name.c_str());
212 cout <<
"OK" << endl;
217 cout <<
"bluring PDFs... " << flush;
221 cout <<
"OK" << endl;
223 }
else if (TTS < 0.0) {
225 THROW(JValueOutOfRange,
"Illegal value of TTS [ns]: " << TTS);
249 const double t1)
const
void load(const char *file_name)
Load from input file.
do echo Generating $dir eval D
JPP::JPDFTable< JFunction1D_t, JPDFMaplist_t > JPDF_t
JPP::JMAPLIST< JPP::JPolint1FunctionalMap, JPP::JPolint1FunctionalMap, JPP::JPolint0FunctionalGridMap, JPP::JPolint0FunctionalGridMap >::maplist JPDFMaplist_t
The elements in a collection are sorted according to their abscissa values and a given distance opera...
This include file containes various data structures that can be used as specific return types for the...
scattered light from EM shower
JPP::JSplineFunction1D< JPP::JSplineElement2S< double, double >, JPP::JCollection, JPP::JResultPDF< double > > JFunction1D_t
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
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 pdfA
PDF for shower.
then for HISTOGRAM in h0 h1
static const JZero zero
Function object to assign zero value.
Definition of zero value for any class.
Auxiliary data structure for shower PDF.
JPP::JMAPLIST< JPP::JPolint1FunctionalMap, JPP::JPolint0FunctionalGridMap, JPP::JPolint0FunctionalGridMap >::maplist JPDFMaplist_t
JMuonPDF_t(const std::string &fileDescriptor, const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10)
Constructor.
Numbering scheme for PDF types.
result_type calculate(const double E, const double R, const double theta, const double phi, const double t1) const
Get PDF.
scattered light from muon
JFunction1D_t::result_type result_type
JFunction1D_t::result_type result_type
JPDF_t pdfB
PDF for average energy losses.
result_type calculate(const double E, const double D, const double cd, const double theta, const double phi, const double t1) const
Get PDF.
scattered light from delta-rays
direct light from EM shower
scattered light from EM showers
General purpose class for a collection of sorted elements.
then usage $script[distance] fi case set_variable R
JPDF_t pdfC
PDF for delta-rays.
direct light from delta-rays
double getDeltaRaysFromMuon(const double E)
Equivalent EM-shower energy due to delta-rays per unit muon track length.
std::string getFilename(const std::string &file_name)
Get file name part, i.e. part after last JEEP::PATHNAME_SEPARATOR if any.
JShowerPDF_t(const std::string &fileDescriptor, const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10)
Constructor.
JPDF_t pdfA
PDF for minimum ionisong particle.
JPP::JSplineFunction1D< JPP::JSplineElement2S< double, double >, JPP::JCollection, JPP::JResultPDF< double > > JFunction1D_t
then usage $script[input file[working directory[option]]] nWhere option can be N
then usage $script[input file[working directory[option]]] nWhere option can be E
JPP::JPDFTable< JFunction1D_t, JPDFMaplist_t > JPDF_t