DefaultShowerGeometryProducer.cc

Go to the documentation of this file.

#include <evt/DefaultShowerGeometryProducer.h>

#include <det/Detector.h>
#include <fwk/MagneticFieldModel.h>

#include <utl/Point.h>
#include <utl/CoordinateSystemPtr.h>

#include <fwk/LocalCoordinateSystem.h>

#include <sstream>

using namespace evt;


utl::CoordinateSystemPtr
DefaultShowerGeometryProducer::MakeGroundParticleCoordinateSystem(const utl::Point& anyPointOnShowerAxis)
{
  utl::CoordinateSystemPtr cs = fwk::LocalCoordinateSystem::Create(anyPointOnShowerAxis);
  double north = 0;
  if (fCorrectForRealisticMagneticFieldDeclination) {
    // magnetic field declination
    north = det::Detector::GetInstance().GetGeoMagneticField().GetDeclination(anyPointOnShowerAxis);
    std::ostringstream info;
    info << "Use magnetic field declination: " << north/utl::degree << " deg";
    INFO(info);
  }
  return utl::CoordinateSystem::RotationZ(fRotZ + north, cs);
}