CastorSimParameters.cc

Go to the documentation of this file.

#include "CLHEP/Random/RandGaussQ.h"
#include "CondFormats/CastorObjects/interface/CastorGain.h"
#include "CondFormats/CastorObjects/interface/CastorGainWidth.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "DataFormats/HcalDetId/interface/HcalGenericDetId.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "SimCalorimetry/CastorSim/interface/CastorSimParameters.h"
#include <cassert>

CastorSimParameters::CastorSimParameters(double simHitToPhotoelectrons,
                                         double photoelectronsToAnalog,
                                         double samplingFactor,
                                         double timePhase,
                                         bool syncPhase)
    : CaloSimParameters(
          simHitToPhotoelectrons, photoelectronsToAnalog, samplingFactor, timePhase, 6, 4, false, syncPhase),
      theDbService(nullptr),
      theSamplingFactor(samplingFactor),
      nominalfCperPE(1) {}

CastorSimParameters::CastorSimParameters(const edm::ParameterSet &p)
    : CaloSimParameters(p),
      theDbService(nullptr),
      theSamplingFactor(p.getParameter<double>("samplingFactor")),
      nominalfCperPE(p.getParameter<double>("photoelectronsToAnalog")) {}

double CastorSimParameters::getNominalfCperPE() const {
  // return the nominal PMT gain value of CASTOR from the config file.
  return nominalfCperPE;
}

double CastorSimParameters::photoelectronsToAnalog(const DetId &detId) const {
  // calculate factor (PMT gain) using sampling factor value & available
  // electron gain
  return theSamplingFactor / fCtoGeV(detId);
}

double CastorSimParameters::fCtoGeV(const DetId &detId) const {
  assert(theDbService != nullptr);
  HcalGenericDetId hcalGenDetId(detId);
  const CastorGain *gains = theDbService->getGain(hcalGenDetId);
  const CastorGainWidth *gwidths = theDbService->getGainWidth(hcalGenDetId);
  double result = 0.0;
  if (!gains || !gwidths) {
    edm::LogError("CastorAmplifier") << "Could not fetch HCAL conditions for channel " << hcalGenDetId;
  } else {
    // only one gain will be recorded per channel, so just use capID 0 for now
    result = gains->getValue(0);
  }
  return result;
}