Deprecated/UpgradeASCIITests/SdSimulationCalibrationFillerASCII/SdSimulationCalibrationFiller.cc

Go to the documentation of this file.

#include <det/Detector.h>

#include <sdet/SDetector.h>

#include <evt/Event.h>

#include <fwk/CentralConfig.h>

#include <sevt/SEvent.h>
#include <sevt/Station.h>
#include <sevt/StationSimData.h>
#include <sevt/StationCalibData.h>
#include <sevt/PMTCalibData.h>
#include <sevt/PMTSimData.h>
#include <sevt/PMT.h>

#include <sdet/PMT.h>
#include <sdet/PMTConstants.h>

#include <utl/Reader.h>

#include "SdSimulationCalibrationFiller.h"

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

using namespace SdSimulationCalibrationFillerASCII;


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

   topB.GetChild("calibrationVersion").GetData(fcalibVersion);

   Branch calibBranch =
     CentralConfig::GetInstance()->GetTopBranch("SdCalibrator");

   calibBranch.GetChild("peakVEMConversionFactor").
     GetData(fPeakVEMConversionFactor);
   calibBranch.GetChild("chargeVEMConversionFactor").
     GetData(fChargeVEMConversionFactor);
   calibBranch.GetChild("onlineChargeVEMFactor").
     GetData(fOnlineChargeVEMFactor);

   cout << "got magic numbers " << fPeakVEMConversionFactor << ' ' <<
     fChargeVEMConversionFactor << ' ' << fOnlineChargeVEMFactor << endl;

  return eSuccess;
}


VModule::ResultFlag
SdSimulationCalibrationFiller::Run(Event& event)
{
  if (!event.HasSEvent())
    event.MakeSEvent();

  SEvent& sevent = event.GetSEvent();

  for (SEvent::CandidateStationIterator seIt = sevent.CandidateStationsBegin();
       seIt != sevent.CandidateStationsEnd(); ++seIt) {

    string signature;
    if (seIt->HasSimData()) {
      signature = seIt->GetSimData().GetSimulatorSignature();

      // FIXME ( for some reason the SdSimCalibrationManager does not load the values in the SdSimCalibrationConstants.xml )
      if (signature == "G4TankSimulatorFullASCII")
        signature="G4TankSimulatorFullOG";
      else if (signature == "G4TankSimulatorASCII")
        signature="G4TankSimulatorOG";
    } else
      ERROR("No SimData found, though it should exist at this point");

    if (!seIt->HasCalibData())
      seIt->MakeCalibData();

    sevt::StationCalibData& statCalib = seIt->GetCalibData();

    statCalib.SetVersion(fCalibVersion);
    statCalib.SetStartSecond(1);
    statCalib.SetEndSecond(3);
    statCalib.SetNT1(100);
    statCalib.SetNT2(20);
    statCalib.SetNTot(1);

    for (sevt::Station::PMTIterator pmtIt = seIt->PMTsBegin(sdet::PMTConstants::eAnyType);
         pmtIt != seIt->PMTsEnd(sdet::PMTConstants::eAnyType); ++pmtIt) {

      if (!pmtIt->HasCalibData())
        pmtIt->MakeCalibData();

      try {
        PMTCalibData& pmtCal = pmtIt->GetCalibData();

        // Currently, there is no way to make a new PMT in SDetector (and it is forbidden to get a
        // PMT with an index higher than those required by the WCD PMTs). This will be resolved,
        // but until now, the sdet PMT contents of WCD PMT1 are called here.
        const sdet::PMT& detPMT =
          det::Detector::GetInstance().GetSDetector().GetStation(*seIt).GetPMT((pmtIt->GetType() == sdet::PMTConstants::eWaterCherenkovLarge) ? pmtIt->GetId() : 1);

        // apply calibration magic numbers from Laura (in reverse)
        pmtCal.SetVEMPeak(detPMT.GetVEMPeak(signature) * fPeakVEMConversionFactor);
        pmtCal.SetVEMCharge(detPMT.GetVEMCharge(signature) * fChargeVEMConversionFactor/fOnlineChargeVEMFactor);

        pmtCal.SetBaseline(detPMT.GetBaseline(sdet::PMTConstants::eHighGain),
                           detPMT.GetBaselineRMS(sdet::PMTConstants::eHighGain),
                           sevt::PMTConstants::eHighGain);
        pmtCal.SetBaseline(detPMT.GetBaseline(sdet::PMTConstants::eLowGain),
                           detPMT.GetBaselineRMS(sdet::PMTConstants::eLowGain),
                           sevt::PMTConstants::eLowGain);

        pmtCal.SetIsTubeOk(true);
        pmtCal.SetRate(7e3);
        pmtCal.SetDynodeAnodeRatio(detPMT.GetDA(), detPMT.GetDARMS());

        pmtCal.SetEvolution(0);
        pmtCal.SetNumberTDA(0);
        pmtCal.SetDynodeAnodeDelay(0, 0);
        pmtCal.SetDynodeAnodeDelayChi2(0);
      } catch (utl::NonExistentComponentException ex) {
        ostringstream msg;
        msg << "Skipping station " << seIt->GetId() << ", pmt " << pmtIt->GetId();
        INFO(msg);
        continue;
      }
    }
  }

  return eSuccess;
}


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