46 int main(
int argc,
char **argv)
66 JParser<> zap(
"Auxiliary program to draw PDF of Cherenkov light from EM-shower including shower profile.");
73 zap[
'E'] =
make_field(E,
"shower energy [GeV]");
75 zap[
'c'] =
make_field(cd,
"cosine emission angle");
76 zap[
'D'] =
make_field(dir,
"(theta, phi) of PMT [rad]");
83 catch(
const exception &error) {
84 FATAL(error.what() << endl);
104 for (
double dt; cin >> dt; ) {
108 cout << setw(2) << *F <<
' '
110 <<
FIXED(5,1) << D <<
' '
111 <<
FIXED(5,2) << cd <<
' '
114 <<
FIXED(5,1) << dt <<
' '
148 for (
int i = 1; i <= h0.GetNbinsX(); ++i) {
150 const double dt = h0.GetBinCenter(i) - t0;
158 h0.SetBinContent(i, value);
Properties of Antares PMT and deep-sea water.
double getAngularAcceptance(const double x)
Angular acceptence of PMT.
int main(int argc, char **argv)
double getAbsorptionLength(const double lambda)
double getScatteringLength(const double lambda)
double absorptionLengthFactor
Scaling of absorption and scattering length.
double scatteringLengthFactor
General purpose messaging.
#define DEBUG(A)
Message macros.
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
I/O formatting auxiliaries.
Properties of KM3NeT PMT and deep-sea water.
Data structure for angles in three dimensions.
double getTheta() const
Get theta angle.
double getPhi() const
Get phi angle.
Utility class to parse command line options.
Probability Density Functions of the time response of a PMT with an implementation of the JAbstractPM...
double getLightFromEMshower(const int type, const double D_m, const double cd, const double theta, const double phi, const double t_ns) const
Probability density function for light from EM-shower.
double getAmbientPressure()
Get ambient pressure.
double getScatteringProbability(const double x)
Function to describe light scattering in water.
double getPhotocathodeArea()
Get photo-cathode area of PMT.
const JPolynome f1(1.0, 2.0, 3.0)
Function.
JAbstractHistogram< double > JHistogram_t
Type definition for scan along axis.
double getQE(const double R, const double mu)
Get QE for given ratio of hit probabilities and expectation value of the number of photo-electrons.
double getIndexOfRefraction()
Get average index of refraction of water corresponding to group velocity.
double getMinimalWavelength()
Get minimal wavelength for PDF evaluations.
@ DIRECT_LIGHT_FROM_EMSHOWER
direct light from EM shower
double getMaximalWavelength()
Get maximal wavelength for PDF evaluations.
static const double C
Physics constants.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Auxiliary data structure for floating point format specification.
Auxiliary data structure for floating point format specification.