MicrobooneOpDetResponse_service.cc

Go to the documentation of this file.

//
//  \file MicrobooneOpDetResponse_service.cc
//

#include "art/Framework/Services/Registry/ServiceDefinitionMacros.h"

#include "larana/OpticalDetector/OpDetResponseInterface.h"

#include "lardataobj/Simulation/SimPhotons.h"

namespace opdet {
  class MicrobooneOpDetResponse : public OpDetResponseInterface {
  public:
    MicrobooneOpDetResponse(fhicl::ParameterSet const& pset);

  private:
    void doReconfigure(fhicl::ParameterSet const& p) override;
    bool doDetected(int OpChannel, const sim::OnePhoton& Phot, int& newOpChannel) const override;
    bool doDetectedLite(int OpChannel, int& newOpChannel) const override;

    float fQE; // Quantum efficiency of tube

    float fWavelengthCutLow;  // Sensitive wavelength range
    float fWavelengthCutHigh; //
  };                          // class MicrobooneOpDetResponse
}

DECLARE_ART_SERVICE_INTERFACE_IMPL(opdet::MicrobooneOpDetResponse,
                                   opdet::OpDetResponseInterface,
                                   LEGACY)

namespace opdet {
  //--------------------------------------------------------------------
  MicrobooneOpDetResponse::MicrobooneOpDetResponse(fhicl::ParameterSet const& pset)
  {
    this->doReconfigure(pset);
  }

  //--------------------------------------------------------------------
  void MicrobooneOpDetResponse::doReconfigure(fhicl::ParameterSet const& pset)
  {
    fQE = pset.get<double>("QuantumEfficiency");
    //double tempfQE=            pset.get<double>("QuantumEfficiency");
    fWavelengthCutLow = pset.get<double>("WavelengthCutLow");
    fWavelengthCutHigh = pset.get<double>("WavelengthCutHigh");

  }

  //--------------------------------------------------------------------
  bool MicrobooneOpDetResponse::doDetected(int OpChannel,
                                           const sim::OnePhoton& Phot,
                                           int& newOpChannel) const
  {

    newOpChannel = OpChannel;

    double wavel = wavelength(Phot.Energy);
    // Check wavelength acceptance
    if (wavel < fWavelengthCutLow) return false;
    if (wavel > fWavelengthCutHigh) return false;

    return true;
  }

  //--------------------------------------------------------------------
  bool MicrobooneOpDetResponse::doDetectedLite(int OpChannel, int& newOpChannel) const
  {
    newOpChannel = OpChannel;

    return true;
  }

} // namespace

DEFINE_ART_SERVICE_INTERFACE_IMPL(opdet::MicrobooneOpDetResponse, opdet::OpDetResponseInterface)