MeanNpeFiller.cc

Go to the documentation of this file.

/*
  For the simulation of threshold values. This module fills all the things FdElectronicsSimulator needs to run into an event.
  \author Joachim Debatin
*/

#include "MeanNpeFiller.h"

#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <utl/TabulatedFunction.h>
#include <utl/MathConstants.h>
#include <utl/String.h>

#include <fwk/CentralConfig.h>

#include <evt/Event.h>

#include <evt/DefaultShowerGeometryProducer.h>

#include <fevt/FEvent.h>
#include <fevt/Eye.h>
#include <fevt/Telescope.h>
#include <fevt/TelescopeSimData.h>
#include <fevt/PixelSimData.h>
#include <fevt/Pixel.h>

#include <det/Detector.h>

#include <fdet/FDetector.h>
#include <fdet/Eye.h>
#include <fdet/Telescope.h>
#include <fdet/Pixel.h>
#include <fdet/Camera.h>

using namespace MeanNpeFillerKG;
using namespace std;
using namespace utl;
using namespace fwk;
using namespace fevt;
using namespace det;
using namespace fdet;


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

  fVerbosityLevel = 0;
  if (topB.GetChild("verbosityLevel"))
    topB.GetChild("verbosityLevel").GetData(fVerbosityLevel);

  if (topB.GetChild("Npe")) {
    topB.GetChild("Npe").GetData(fNpe);
  } else {
    ostringstream info;
    info <<"Npelevels must be set! Abort Simulation.";
    ERROR(info);
    return eFailure;
  }

  if (topB.GetChild("NRuns"))
    topB.GetChild("NRuns").GetData(fMaxRunNumber);
  fRunNumber = 1;

  if (topB.GetChild("Eye"))
    topB.GetChild("Eye").GetData(fEye);

  if (topB.GetChild("Tel"))
    topB.GetChild("Tel").GetData(fTel);

  // the model config signature taken from TelescopeSimulatorKG
  ostringstream configSS;
  configSS << "TelescopeSimulatorKG";
  if (false)//fDoNoShadow || fDoNoShadowSupport || fHasNoMercedes)
    configSS << "(camera-shadow / support-shadow / mercedes) = ("
             << (true ? "yes" : "no") << " / "
             << (true ? "yes" : "no") << " / "
             << (true ? "yes" : "no") << ")";
  fGeneralConfigSignature = configSS.str();

  // info output
  ostringstream info;
  if (fVerbosityLevel) {
    info << "verbosity level: " << fVerbosityLevel << "\n"
            "Used Npe values: " << Join(" ", fNpe);
    INFO(info);
  }

  Detector::GetInstance().Update(UTCDateTime(2012,1,1,0,0,0).GetTimeStamp());

  return eSuccess;
}


VModule::ResultFlag
MeanNpeFiller::Run(evt::Event& event)
{
  if (fRunNumber > fMaxRunNumber){
    ostringstream msg;
    msg << "Reached maximum run number. Exit!";
    INFO(msg);
    return eBreakLoop;
  }

  const double meanPe = fNpe[(fRunNumber-1) % fNpe.size()];

  if (fVerbosityLevel >= 1) {
    ostringstream msg;
    msg << "This is run number: " << fRunNumber << " meanPe =  " << meanPe;
    INFO(msg);
  }

  ++fRunNumber;

  if (!event.HasSimShower())
    event.MakeSimShower(evt::DefaultShowerGeometryProducer());  // ... pretty irrelevant

  if(!event.HasFEvent())
    event.MakeFEvent();

  const Detector& detector = Detector::GetInstance();
  const FDetector& detFD = detector.GetFDetector();

  FEvent& fEvent = event.GetFEvent();

  if (fVerbosityLevel >= 1) {
    ostringstream msg;
    msg << "Eye: "<< fEye <<" Tel: "<< fTel;
    INFO(msg);
  }

  fEvent.MakeEye(fEye);
  fEvent.GetEye(fEye).MakeTelescope(fTel);
  fevt::Telescope& telEvent = fEvent.GetEye(fEye).GetTelescope(fTel);

  // create configSignature see Init()
  telEvent.MakeSimData();
  fevt::TelescopeSimData& telSim = telEvent.GetSimData();
  const fdet::Telescope& detTel = detFD.GetTelescope(telEvent);
  telSim.SetConfigSignatureStr(detTel.GetConfigSignatureStr(fGeneralConfigSignature));

  //Needed to correct for scaling in ElectronicsSimulator
  const double normWavelength = detFD.GetReferenceLambda();

  //create pixels and pixelSimData
  for (unsigned int pixelId = 1; pixelId <= 440; ++pixelId) {
    telEvent.MakePixel(pixelId);
    telEvent.GetPixel(pixelId).MakeSimData();
    fevt::Pixel& pix = telEvent.GetPixel(pixelId);
    fevt::PixelSimData& pixSimData = pix.GetSimData();
    const fdet::Pixel& detPixel = detFD.GetPixel(pix);
    const TabulatedFunction& qEff = detPixel.GetQEfficiency();
    const double qEffAtNormWavelength = qEff.Y(normWavelength);

    // Correct for scaling in ElectronicsSimulator
    const double backphotons = meanPe / qEffAtNormWavelength;
    pixSimData.SetMeanBgPhotonFlux(backphotons);
  }

  return eSuccess;
}