JKM3Sim.cc File Reference

#include <string>
#include <iostream>
#include <iomanip>
#include <sstream>
#include <fstream>
#include <vector>
#include <map>
#include "JLang/JException.hh"
#include "JMath/JConstants.hh"
#include "JPhysics/JConstants.hh"
#include "JPhysics/KM3NeT.hh"
#include "JPhysics/JDispersion.hh"
#include "Jeep/JStreamToolkit.hh"
#include "Jeep/JParser.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 270 of file JKM3Sim.cc.

{
  using namespace std;
  using namespace JPP;
 
  string    inputFile;
  map_type  precision (0.0);
  map_type  correction(1.0);
  int       debug;
 
  precision[JKM3Sim::RINDEX_WATER_t]      =   1.0e-4;
  precision[JKM3Sim::ABSORPTION_WATER_t]  =   1.0e-4;
  precision[JKM3Sim::MIE_WATER_t]         =   0.2;
  precision[JKM3Sim::ANG_ACCEPT_t]        =   0.5;
  precision[JKM3Sim::Q_EFF_t]             =   0.01;
  
  try {
 
    JParser<> zap;
 
    zap['f'] = make_field(inputFile);
    zap['p'] = make_field(precision)   = JPARSER::initialised();
    zap['c'] = make_field(correction)  = JPARSER::initialised();
    zap['d'] = make_field(debug)       = 1;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  JKM3Sim km3sim;
 
  ifstream in(inputFile.c_str());
 
  in >> km3sim;
 
  in.close();
 
  DEBUG(km3sim);
  
  using namespace KM3NET;
 
 
  {
    /*
      const double g = 9.80665;                                                   // [m/s^2]
      const double P = detectorDepth * DENSITY_SEA_WATER * 1.0e-3 * g * 1.0e1;    // [Atm]
    */
    const double P = 315.0;                                                     // [Atm]
 
    DEBUG("Pressure " << P << " [Atm]" << endl);
 
    const JDispersion dispersion(P);
 
    for (size_t i = 0; i != km3sim.ppckov.size(); ++i) {
 
      const double x = km3sim.ppckov[i];
 
      ASSERT(fabs(km3sim.rindex_water[i] - dispersion.getIndexOfRefractionPhase(x) * correction.get(JKM3Sim::RINDEX_WATER_t)) <= precision.get(JKM3Sim::RINDEX_WATER_t),
             "Test of index of refraction at " << FIXED(6,1) << x << ' ' << FIXED(7,4) << km3sim.rindex_water[i] << ' ' << FIXED(7,4) << dispersion.getIndexOfRefractionPhase(x) * correction.get(JKM3Sim::RINDEX_WATER_t));
    }
  }
  {
    for (const auto& i : km3sim.absorption_water) {
      ASSERT(fabs(i.second - getAbsorptionLength(i.first) * correction.get(JKM3Sim::ABSORPTION_WATER_t)) <= precision.get(JKM3Sim::ABSORPTION_WATER_t),
             "Test of absorption length at " << FIXED(6,1) << i.first << ' ' << FIXED(7,3) << i.second << ' ' << FIXED(7,3) << getAbsorptionLength(i.first) * correction.get(JKM3Sim::ABSORPTION_WATER_t));
    }
  }
  {
    for (size_t i = 0; i != km3sim.ppckov.size(); ++i) {
 
      const double x = km3sim.ppckov[i];
 
      ASSERT(fabs(km3sim.mie_water[i] - getScatteringLength(x) * correction.get(JKM3Sim::MIE_WATER_t)) <= precision.get(JKM3Sim::MIE_WATER_t),
             "Test of scattering length at " << FIXED(6,1) << x << ' ' << FIXED(7,3) << km3sim.mie_water[i] << ' ' << FIXED(7,3) << getScatteringLength(x) * correction.get(JKM3Sim::MIE_WATER_t));
    }
  }
  {
    for (size_t i = 0; i != km3sim.ppckov.size(); ++i) {
 
      const double x = km3sim.ppckov[i];
 
      ASSERT(fabs(km3sim.q_eff[i] - getQE(x) * correction.get(JKM3Sim::Q_EFF_t)) <= precision.get(JKM3Sim::Q_EFF_t),
             "Test of QE at " << FIXED(6,1) << x << ' ' << FIXED(7,4) << km3sim.q_eff[i] << ' ' << FIXED(7,4) << getQE(x) * correction.get(JKM3Sim::Q_EFF_t));
    }
  }
  {
    const double U = 1.0e4;                                                     // m^2 -> cm^2
 
    for (size_t i = 0; i != km3sim.cos_angles.size(); ++i) {
 
      const double x = km3sim.cos_angles[i];
 
      ASSERT(fabs(km3sim.ang_accept[i]*km3sim.A - getPhotocathodeArea() * getAngularAcceptance(x) * U * correction.get(km3sim.ANG_ACCEPT_t)) <= precision.get(km3sim.ANG_ACCEPT_t),
             "Test of PMT acceptance at " << FIXED(6,2) << x << ' ' << FIXED(9,5) << km3sim.ang_accept[i]*km3sim.A << ' ' << FIXED(9,5) << getPhotocathodeArea() * getAngularAcceptance(x) * U * correction.get(km3sim.ANG_ACCEPT_t));
    }
  }
}