MicrobooneOpDetResponse_service.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
//
//  \file MicrobooneOpDetResponse_service.cc
//
////////////////////////////////////////////////////////////////////////

#include "larana/OpticalDetector/MicrobooneOpDetResponse.h"
#include "lardataobj/Simulation/SimPhotons.h"

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

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");

        /**
         * Don't apply QE here.  It is applied in the uboone
         * electronics simulation.
         **
        // Correct out the prescaling applied during simulation
        auto const* LarProp = lar::providerFrom<detinfo::LArPropertiesService>();
        fQE = tempfQE / LarProp->ScintPreScale();

        if (fQE > 1.0001 ) {
            mf::LogWarning("MicrobooneOpDetResponse_service") << "Quantum efficiency set in OpDetResponse_service, " << tempfQE
                                                              << " is too large.  It is larger than the prescaling applied during simulation, "
                                                              << LarProp->ScintPreScale()
                                                              << ".  Final QE must be equalt to or smaller than the QE applied at simulation time.";
            assert(false);
        }
        **/
    }


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

        newOpChannel = OpChannel;

        /**
         * Don't apply QE here.  It is applied in the uboone
         * electronics simulation.
         **
        // Check QE
        if ( CLHEP::RandFlat::shoot(1.0) > fQE ) return false;
        **/

        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;

        /**
         * Don't apply QE here.  It is applied in the uboone
         * electronics simulation.
         **
        // Check QE
        if ( CLHEP::RandFlat::shoot(1.0) > fQE ) return false;
        **/

        return true;
    }



} // namespace

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