JDrawPD0.cc File Reference

Auxiliary program to draw PDF of Cherenkov light from bright point. More...

#include <string>
#include <iostream>
#include <iomanip>
#include "TROOT.h"
#include "TFile.h"
#include "TH1D.h"
#include "JPhysics/JPDF.hh"
#include "JPhysics/Antares.hh"
#include "JPhysics/KM3NeT.hh"
#include "JTools/JSpline.hh"
#include "JTools/JQuantiles.hh"
#include "JTools/JAbstractHistogram.hh"
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions
double	getAbsorptionLength (const double lambda)
double	getScatteringLength (const double lambda)
int	main (int argc, char **argv)

Variables
double	absorptionLengthFactor
	Scaling of absorption and scattering length. More...
double	scatteringLengthFactor

Detailed Description

Auxiliary program to draw PDF of Cherenkov light from bright point.

Author

mdejong

Definition in file JDrawPD0.cc.

Function Documentation

double getAbsorptionLength ( const double  lambda )

inline

Definition at line 27 of file JDrawPD0.cc.

{
  return absorptionLengthFactor * NAMESPACE::getAbsorptionLength(lambda);
}

double getScatteringLength ( const double  lambda )

inline

Definition at line 33 of file JDrawPD0.cc.

{
  return scatteringLengthFactor * NAMESPACE::getScatteringLength(lambda);
}

int main	(	int	argc,
		char **	argv
	)

Definition at line 45 of file JDrawPD0.cc.

{
  using namespace std;
  using namespace JPP;

  typedef JAbstractHistogram<double> JHistogram_t;

  string        outputFile;
  int           numberOfPoints;
  double        epsilon;
  double        E;
  double        D;
  double        ct;
  vector<int>   function;
  JHistogram_t  histogram;
  int           debug;

  try { 

    JParser<> zap("Auxiliary program to draw PDF of Cherenkov light from bright point.");

    zap['o'] = make_field(outputFile)                                           = "pd0.root";
    zap['n'] = make_field(numberOfPoints,         "points for integration")     = 25;
    zap['e'] = make_field(epsilon,                "precision for integration")  = 1.0e-10;
    zap['A'] = make_field(absorptionLengthFactor, "scaling factor")             = 1.0;
    zap['S'] = make_field(scatteringLengthFactor, "scaling factor")             = 1.0;
    zap['E'] = make_field(E,                      "shower energy [GeV]");
    zap['R'] = make_field(D,                      "distance [m]");
    zap['c'] = make_field(ct,                     "cosine PMT angle");
    zap['F'] = make_field(function,               "PDF type");
    zap['H'] = make_field(histogram,              "histogram binning")          = JHistogram_t();
    zap['d'] = make_field(debug)                  = 0;

    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }


  const JPDF_C
    pdf(NAMESPACE::getPhotocathodeArea(),
        NAMESPACE::getQE,
        NAMESPACE::getAngularAcceptance,
        getAbsorptionLength,
        getScatteringLength,
        NAMESPACE::getScatteringProbability,
        NAMESPACE::getAmbientPressure(),
        getMinimalWavelength(),
        getMaximalWavelength(),
        numberOfPoints,
        epsilon);


  if (outputFile == "") {

    for (double dt; cin >> dt; ) {

      for (vector<int>::const_iterator F = function.begin(); F != function.end(); ++F) {

        cout << setw(2)         << *F             << ' '
             << SCIENTIFIC(7,1) << E              << ' '
             << FIXED(5,1)      << D              << ' '
             << FIXED(5,2)      << ct             << ' '
             << FIXED(5,1)      << dt             << ' '
             << SCIENTIFIC(9,3) << pdf.getLightFromBrightPoint(*F, D, ct, dt) * E << endl;
      }
    }

    return 0;
  }


  TFile out(outputFile.c_str(), "recreate");
 
  //const double t0   =  D * getIndexOfRefraction() / C;        // time [ns]
  const double t0   =  0.0;                                   // time [ns]
 
  if (!histogram.is_valid()) {
    
    if (function.size() == 1 && function[0] == DIRECT_LIGHT_FROM_BRIGHT_POINT) {

      histogram = JHistogram_t(t0 - 20.0, t0 + 50.0);

      histogram.setBinWidth(0.1);

    } else {

      histogram = JHistogram_t(t0 - 20.0, t0 + 500.0);

      histogram.setBinWidth(0.5);
    }
  }

  TH1D h0("h0", NULL, histogram.getNumberOfBins(), histogram.getLowerLimit(), histogram.getUpperLimit());
 
  JSplineFunction1D_t f1;
 
  for (int i = 1; i <= h0.GetNbinsX(); ++i) {

    const double dt = h0.GetBinCenter(i) - t0;

    double value = 0.0;

    for (vector<int>::const_iterator F = function.begin(); F != function.end(); ++F) {
      value += pdf.getLightFromBrightPoint(*F, D, ct, dt) * E;
    }

    h0.SetBinContent(i, value);

    f1[dt] = value;
  }

  f1.compile();

  JQuantiles quantiles(f1);

  DEBUG("int  " << quantiles.getIntegral() << endl);
  DEBUG("x    " << quantiles.getX()        << endl);
  DEBUG("y    " << quantiles.getY()        << endl);
  DEBUG("FWHM " << quantiles.getFWHM()     << endl);
  

  out.Write();
  out.Close();
}

Variable Documentation

double absorptionLengthFactor

Scaling of absorption and scattering length.

Definition at line 24 of file JDrawPD0.cc.

double scatteringLengthFactor

Definition at line 25 of file JDrawPD0.cc.