LaserGenerator.cc

Go to the documentation of this file.

#include "LaserGenerator.h"

#include <fwk/CentralConfig.h>
#include <fwk/RandomEngineRegistry.h>

#include <evt/Event.h>
#include <evt/ShowerSimData.h>
#include <evt/LaserData.h>
#include <evt/DefaultShowerGeometryProducer.h>

#include <utl/AugerCoordinateSystem.h>
#include <utl/AugerUnits.h>
#include <utl/CoordinateSystem.h>
#include <utl/ErrorLogger.h>
#include <utl/RandomEngine.h>
#include <utl/Reader.h>
#include <utl/Point.h>
#include <utl/Particle.h>
#include <utl/UTMPoint.h>
#include <utl/Vector.h>
#include <utl/PhysicalConstants.h>

#include <CLHEP/Random/RandGauss.h>

#include <cstddef>
#include <iostream>
#include <string>
#include <sstream>

using namespace LaserGeneratorNA;
using CLHEP::RandGauss;

using namespace fwk;
using namespace evt;
using namespace utl;
using namespace std;


LaserGenerator::LaserGenerator() : 
  fTheta(0.),
  fPhi(0.),
  fEnergy(0.),
  fEnergyRMS(0.),
  fDoEnergySmearing(false),
  fWavelength(0.),
  fRandomEngine(0)
{
}

LaserGenerator::~LaserGenerator()
{
}

VModule::ResultFlag
LaserGenerator::Init()
{
  Branch topB = CentralConfig::GetInstance()->GetTopBranch("LaserGeneratorNA");

  topB.GetChild("theta").GetData(fTheta);
  topB.GetChild("phi").GetData(fPhi);

  Branch laserB = topB.GetChild("laserEnergy");
  laserB.GetChild("energyMean").GetData(fEnergy);
  laserB.GetChild("energyRMS").GetData(fEnergyRMS);
  laserB.GetChild("doSmearing").GetData(fDoEnergySmearing);

  topB.GetChild("wavelength").GetData(fWavelength);

  fRandomEngine =
    &RandomEngineRegistry::GetInstance().Get(RandomEngineRegistry::eDetector);

  ostringstream info;
  info << "Version: " << GetVersionInfo(VModule::eRevisionNumber) << "\n"
       << " Parameters:\n"
       << "            theta: " << fTheta / degree << " degree\n"
       << "              phi: " << fPhi / degree << " degree\n"
       << "           energy: " << fEnergy / (1e-3 * joule) << " mJ\n"
       << "       energy RMS: " << fEnergyRMS / perCent << "%\n"
       << "  energy smearing: " << (fDoEnergySmearing ? "yes\n" : "no\n")
       << "       wavelength: " << fWavelength / nanometer << " nm\n";
  INFO(info);
 
  return eSuccess;
}

VModule::ResultFlag
LaserGenerator::Run(evt::Event& theEvent)
{
  INFO("LaserGenerator");
  
  if (theEvent.HasSimShower()) {
    ERROR("Shower sim data already exists. Cannot produce Laser simulation.");
    return eFailure;
  }
  theEvent.MakeSimShower(evt::DefaultShowerGeometryProducer());

  evt::ShowerSimData& theShower = theEvent.GetSimShower();

  // NO primary particle --> LASER event
  theShower.SetPrimaryParticle (Particle::eUndefined);

  theShower.SetGroundParticleCoordinateSystemAzimuth(fPhi);
  theShower.SetGroundParticleCoordinateSystemZenith(fTheta);

  double energy = fEnergy;
  if (fDoEnergySmearing) {
    energy = RandGauss::shoot(&fRandomEngine->GetEngine(),
                              fEnergy,
                              fEnergyRMS * fEnergy);
    ostringstream info;
    info << "Laser energy: " << energy / (1e-3 * joule) << " mJ";
    INFO(info);
  }
  theShower.SetEnergy(energy);

  if (!theShower.HasLaserData()) 
    theShower.MakeLaserData();

  evt::LaserData& theLaser = theShower.GetLaserData();
  
  theLaser.SetLaserWavelength(fWavelength);
  theLaser.SetPolarizationStatus(false);

  return eSuccess;
}

VModule::ResultFlag
LaserGenerator::Finish()
{
  return eSuccess;
}