Inheritance diagram for PhotonPropogationICARUS:

Public Member Functions
	PhotonPropogationICARUS (fhicl::ParameterSet const &pset)

virtual void	configure (fhicl::ParameterSet const &pset)

virtual void	produce (art::Event &e) override

virtual void	beginJob () override
	Begin job method. More...

virtual void	endJob () override
	End job method. More...

Private Attributes
art::InputTag	fSimPhotonModuleLabel
	The full collection of SimPhotons. More...

unsigned int	fNumEvent = 0
	Number of events seen. More...

geo::GeometryCore const *	fGeometry = nullptr
	Pointer to Geometry service. More...

CLHEP::HepRandomEngine &	fPhotonEngine

Static Private Attributes
static constexpr unsigned int	MaxPhotons = 10000000U
	We don't keep more than this number of photons per `sim::SimPhoton`. More...

Detailed Description

Definition at line 54 of file PhotonPropogationICARUS_module.cc.

Constructor & Destructor Documentation

PhotonPropogationICARUS::PhotonPropogationICARUS ( fhicl::ParameterSet const & pset )

explicit

Constructor.

Arguments:

pset - Fcl parameters.

Definition at line 93 of file PhotonPropogationICARUS_module.cc.

   : EDProducer{pset}, fGeometry(lar::providerFrom<geo::Geometry>())
   , fPhotonEngine(art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this, "HepJamesRandom", "icarusphoton", pset, "SeedPhoton"))
 //  , fDetectorProperties(lar::providerFrom<detinfo::DetectorPropertiesService>())
 {
     configure(pset);
     
     produces<std::vector<sim::SimPhotons>>();
 
     // Report.
     mf::LogDebug("PhotonPropogationICARUS") << "PhotonPropogationICARUS configured";
 } // PhotonPropagationICARUS::PhotonPropagationICARUS()

Member Function Documentation

void PhotonPropogationICARUS::beginJob ( )

overridevirtual

Begin job method.

Definition at line 120 of file PhotonPropogationICARUS_module.cc.

 {
     // Access ART's TFileService, which will handle creating and writing
     // histograms and n-tuples for us.
 //    art::ServiceHandle<art::TFileService> tfs;
     
 //    art::TFileDirectory dir = tfs->mkdir(Form("PhotonPropogation"));
 
     return;
 }

void PhotonPropogationICARUS::configure ( fhicl::ParameterSet const & pset )

virtual

Reconfigure method.

Arguments:

pset - Fcl parameter set.

Definition at line 113 of file PhotonPropogationICARUS_module.cc.

 {
     fSimPhotonModuleLabel = pset.get<art::InputTag>("SimPhotonsModuleLabel");
 }

void PhotonPropogationICARUS::endJob ( )

overridevirtual

End job method.

Definition at line 243 of file PhotonPropogationICARUS_module.cc.

 {
     mf::LogInfo("PhotonPropogationICARUS") << "Looked at " << fNumEvent << " events" << std::endl;
 }

void PhotonPropogationICARUS::produce ( art::Event & event )

overridevirtual

Produce method.

Arguments:

evt - Art event.

This is the primary method.

Definition at line 140 of file PhotonPropogationICARUS_module.cc.

 {
     ++fNumEvent;
     
     /*
      * In LArSoft, detected photons are organised optical channel by channel,
      * each channel with its own sim::SimPhotons.
      * Each channel contains individual information for each of the photons
      * it has detected.
      * Here we go through all detected photons from all the channels,
      * and for each we assign a delay to the start time, parametrized by the
      * distance between its generation point and the PMT.
      *
      */
     // Agreed convention is to ALWAYS output to the event store so get a pointer to our collection
     auto simPhotonVec = std::make_unique<std::vector<sim::SimPhotons>>();
     
     // Read in the digit List object(s).
     auto const& srcSimPhotons
       = *(event.getValidHandle< std::vector<sim::SimPhotons> >(fSimPhotonModuleLabel));
        
     if (srcSimPhotons.empty()) {
         mf::LogWarning("PhotonPropagationICARUS") << "No photons! Nice event you have here.";
         event.put(std::move(simPhotonVec));
         return;
     }
 
     // get hold of all needed services
 //    auto const& pvs = *(art::ServiceHandle<phot::PhotonVisibilityService>());
     CLHEP::RandLandau landauGen(fPhotonEngine);
 
     // Loop through the input photons (this might need to be more complicated...)
     for(const auto& simPhoton : srcSimPhotons)
     {
         auto const channel = simPhoton.OpChannel();
         if (simPhoton.empty()) continue;
         
         auto const& PMTcenter = fGeometry->OpDetGeoFromOpChannel(channel).GetCenter();
        
         // TODO restore the "MF_LOG_TRACE" line for normal operations
         // (MF_LOG_TRACE will print only in debug mode, with debug qualifiers and proper messagefacility settings)
       //  MF_LOG_TRACE("LightPropagationICARUS")
 //        mf::LogVerbatim("LightPropagationICARUS")
 //          << "Processing photon channel #" << channel << ", detector center: " << PMTcenter;
         // 
         // fix the arrival time
         //
         sim::SimPhotons localPhoton = simPhoton; // modify a copy of the original photon
         unsigned int photonNo = 0;
         for (auto& onePhoton: localPhoton) 
         {
              //
              // sanity check: if there are too many photons we are in trouble
              // (like in: not enough computing resources)
              //
              if (photonNo++ >= MaxPhotons) {
                mf::LogError("LightPropagationICARUS")
                  << "Too many photons to handle! only " << (photonNo - 1) << " saved.";
                break;
              }
                
              geo::Point_t const position = geo::vect::toPoint(onePhoton.InitialPosition);
              double const dis = (position - PMTcenter).R();
 
              double mean  = 0.18*dis; // TODO
              double sigma = 0.75*dis; // TODO
              
              //double time_plus = -1;
 
              //while (time_plus<dis/21.74){time_plus=mean + landauGen.fire() * sigma; }// TODO
 
             const double minPropTime = dis / 21.74; // d / (c/n) in [ns]
             //const double maxPropTime = 2000 / 21.74; // dimension / (c/n) in [ns]
             double time_plus;
             do {
             time_plus = mean + landauGen.fire() * sigma; // TODO
                 //time_plus = landauGen.fire(mean,sigma); // TODO
             } while (time_plus < minPropTime);
              
              // TODO restore the "MF_LOG_TRACE" line for normal operations
              MF_LOG_TRACE("LightPropagationICARUS");
 //             mf::LogVerbatim("LightPropagationICARUS")
 //               << "Photon #" << photonNo
 //               << " (at " << position << ", " << dis << " cm far from PMT) given offset "
 //               << time_plus;
              //double time_plus = dis/21.74; 
              //TRandom r;
                 
              onePhoton.Time += time_plus;
                 //onePhoton.Time = time_plus;
              
         } // for photons in simPhoton
         // move the photon with the new arrival time into the new output collection
         simPhotonVec->push_back(std::move(localPhoton));
     } // for all simPhoton channels
     
     // Add tracks and associations to event.
     event.put(std::move(simPhotonVec));
     
 } // LightPropagationICARUS::produce()