Inheritance diagram for phot::PhotonLibraryAnalyzer:

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

void	beginJob ()

void	analyze (const art::Event &evt)

Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDAnalyzer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept

	Analyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	Analyzer (Table< Config > const &config)

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

fhicl::ParameterSetID	selectorConfig () const

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
std::string	fAltXAxis

int	fOpDet

bool	fEachSlice

bool	fEachDetector

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 28 of file PhotonLibraryAnalyzer_module.cc.

Constructor & Destructor Documentation

phot::PhotonLibraryAnalyzer::PhotonLibraryAnalyzer ( fhicl::ParameterSet const & pset )

explicit

Definition at line 53 of file PhotonLibraryAnalyzer_module.cc.

     : EDAnalyzer(pset)
     , fAltXAxis{pset.get<std::string>("alt_x_axis")}
     , fOpDet{pset.get<int>("opdet")}
     , fEachSlice{pset.get<bool>("each_slice")}
     , fEachDetector{pset.get<bool>("each_detector")}
   {
     std::cout<<"Photon library analyzer constructor "<<std::endl;
   }

Member Function Documentation

void phot::PhotonLibraryAnalyzer::analyze ( const art::Event & evt )

Definition at line 252 of file PhotonLibraryAnalyzer_module.cc.

   {
 
   }

void phot::PhotonLibraryAnalyzer::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 64 of file PhotonLibraryAnalyzer_module.cc.

   {
     mf::LogInfo("PhotonLibraryAnalyzer")<<"Analyzing photon library - begin"<< std::endl;
 
 
     art::ServiceHandle<art::TFileService const> tfs;
     art::ServiceHandle<PhotonVisibilityService const> pvs;
     art::ServiceHandle<geo::Geometry const> geom;
 
     int NOpDet = pvs->NOpChannels();
 
     sim::PhotonVoxelDef TheVoxelDef = pvs->GetVoxelDef();
     auto const& UpperCorner = TheVoxelDef.GetRegionUpperCorner();
     auto const& LowerCorner = TheVoxelDef.GetRegionLowerCorner();
 
     mf::LogInfo("PhotonLibraryAnalyzer") << "UpperCorner: " << lar::dump::vector3D(UpperCorner) << "\n"
                                          << "LowerCorner: " << lar::dump::vector3D(LowerCorner);
 
     auto const [ XSteps, YSteps, ZSteps ] = TheVoxelDef.GetSteps(); // unsigned int
 
     // for c2: FullVolume is unused, just call tfs->make
     // TH3D *FullVolume = tfs->make<TH3D>("FullVolume","FullVolume",
     tfs->make<TH3D>("FullVolume","FullVolume",
                     XSteps,LowerCorner.X(),UpperCorner.X(),
                     YSteps,LowerCorner.Y(),UpperCorner.Y(),
                     ZSteps,LowerCorner.Z(),UpperCorner.Z());
 
 
     int reportnum=10000;
 
     int newX, newY;
     if (fAltXAxis == "Z") {
       newX = 2; // Z
       newY = 1; // Y
     }
     else {
       newX = 1; // Y
       newY = 2; // Z
     }
 
 
 
     mf::LogInfo("PhotonLibraryAnalyzer")<<"Analyzing photon library - making historams"<< std::endl;
 
     TH2D* XProjection;
     if (fAltXAxis == "Z") XProjection = tfs->make<TH2D>("XProjection","XProjection",ZSteps,0,ZSteps,YSteps,0,YSteps);
     else                  XProjection = tfs->make<TH2D>("XProjection","XProjection",YSteps,0,YSteps,ZSteps,0,ZSteps);
     TH2D* YProjection = tfs->make<TH2D>("YProjection","YProjection",XSteps,0,XSteps,ZSteps,0,ZSteps);
     TH2D* ZProjection = tfs->make<TH2D>("ZProjection","ZProjection",XSteps,0,XSteps,YSteps,0,YSteps);
 
     //    TH1D * PMTsNoVisibility = tfs->make<TH1D>("PMTsNoVisibility","PMTsNoVisibility", NOpDet,0,NOpDet);
 
     TH1D* VisByN = tfs->make<TH1D>("VisByN","VisByN", NOpDet, 0, NOpDet);
 
     TH2D* XInvisibles;
     if (fAltXAxis == "Z") XInvisibles = tfs->make<TH2D>("XInvisibles","XInvisibles",ZSteps,0,ZSteps,YSteps,0,YSteps);
     else                  XInvisibles = tfs->make<TH2D>("XInvisibles","XInvisibles",YSteps,0,YSteps,ZSteps,0,ZSteps);
     TH2D* YInvisibles = tfs->make<TH2D>("YInvisibles","YInvisibles",XSteps,0,XSteps,ZSteps,0,ZSteps);
     TH2D* ZInvisibles = tfs->make<TH2D>("ZInvisibles","ZInvisibles",XSteps,0,XSteps,YSteps,0,YSteps);
 
 
 
 
 
     std::vector<TH2D*> TheXCrossSections;
     std::vector<TH2D*> TheYCrossSections;
     std::vector<TH2D*> TheZCrossSections;
     if (fEachSlice) {
 
 
       for(unsigned int i=0; i!=XSteps; ++i)
       {
         std::stringstream ss("");
         ss.flush();
         ss<<"projX"<<i;
         if (fAltXAxis == "Z")
           TheXCrossSections.push_back(tfs->make<TH2D>(ss.str().c_str(),ss.str().c_str(), ZSteps, 0,ZSteps, YSteps, 0,YSteps));
         else
           TheXCrossSections.push_back(tfs->make<TH2D>(ss.str().c_str(),ss.str().c_str(), YSteps, 0,YSteps, ZSteps, 0,ZSteps));
 
 
       }
 
       for(unsigned int i=0; i!=YSteps; ++i)
       {
         std::stringstream ss("");
         ss.flush();
         ss<<"projY"<<i;
         TheYCrossSections.push_back(tfs->make<TH2D>(ss.str().c_str(),ss.str().c_str(), XSteps, 0,XSteps, ZSteps, 0, ZSteps));
 
       }
 
       for(unsigned int i=0; i!=ZSteps; ++i)
       {
         std::stringstream ss("");
         ss.flush();
         ss<<"projZ"<<i;
         TheZCrossSections.push_back(tfs->make<TH2D>(ss.str().c_str(),ss.str().c_str(), XSteps, 0,XSteps, YSteps, 0,YSteps));
 
       }
     }
 
 
     std::vector<TH2D*> TheXProjections;
     std::vector<TH2D*> TheYProjections;
     std::vector<TH2D*> TheZProjections;
 
     if (fEachDetector) {
 
       mf::LogInfo("PhotonLibraryAnalyzer")<<"Making projections for each of " << NOpDet << " photon detectors" << std::endl;
 
       for(int i=0; i<NOpDet; ++i)
       {
         char ss[99];
 
         sprintf(ss, "ProjXOpDet%d", i);
         if (fAltXAxis == "Z")
           TheXProjections.push_back(tfs->make<TH2D>(ss, ss, ZSteps, 0,ZSteps, YSteps, 0,YSteps));
         else
           TheXProjections.push_back(tfs->make<TH2D>(ss, ss, YSteps, 0,YSteps, ZSteps, 0,ZSteps));
 
         sprintf(ss, "ProjYOpDet%d", i);
         TheYProjections.push_back(tfs->make<TH2D>(ss, ss, XSteps, 0,XSteps, ZSteps, 0, ZSteps));
 
         sprintf(ss, "ProjZOpDet%d", i);
         TheZProjections.push_back(tfs->make<TH2D>(ss, ss, XSteps, 0,XSteps, YSteps, 0,YSteps));
       }
     }
 
 
     mf::LogInfo("PhotonLibraryAnalyzer")<<"Analyzing photon library - running through voxels "<< std::endl;
 
 
     for(unsigned int i=0; i!=TheVoxelDef.GetNVoxels(); ++i)
     {
       if(i%reportnum==0) std::cout<<"Photon library analyzer at voxel " << i<<std::endl;
 
       auto const Coords = TheVoxelDef.GetVoxelCoords(i);
 
       const float* Visibilities = pvs->GetLibraryEntries(i);
       size_t NOpChannels = pvs->NOpChannels();
 
       float TotalVis=0;
       if (fOpDet < 0) {
         for(size_t ichan=0; ichan!=NOpChannels; ++ichan)
         {
           TotalVis+=Visibilities[ichan];
         }
       }
       else {
         TotalVis = Visibilities[fOpDet];
       }
 
       VisByN->Fill(NOpChannels);
 
       if(TotalVis==0)
       {
         XInvisibles->Fill(Coords[newX],Coords[newY]);
         YInvisibles->Fill(Coords[0],Coords[2]);
         ZInvisibles->Fill(Coords[0],Coords[1]);
       }
 
       if (fEachSlice) {
         TheXCrossSections.at(Coords.at(0))->Fill(Coords[newX],Coords[newY],TotalVis);
         TheYCrossSections.at(Coords.at(1))->Fill(Coords[0],Coords[2],TotalVis);
         TheZCrossSections.at(Coords.at(2))->Fill(Coords[0],Coords[1],TotalVis);
       }
 
       if (fEachDetector) {
         for(size_t ichan=0; ichan!=NOpChannels; ++ichan) {
           TheXProjections.at(ichan)->Fill(Coords[newX],Coords[newY],Visibilities[ichan]);
           TheYProjections.at(ichan)->Fill(Coords[0],Coords[2],Visibilities[ichan]);
           TheZProjections.at(ichan)->Fill(Coords[0],Coords[1],Visibilities[ichan]);
         }
       }
 
       // Always make the summed projections
       XProjection->Fill(Coords[newX], Coords[newY], TotalVis);
       YProjection->Fill(Coords[0], Coords[2], TotalVis);
       ZProjection->Fill(Coords[0], Coords[1], TotalVis);
 
     }
 
     mf::LogInfo("PhotonLibraryAnalyzer")<<"Analyzing photon library - end"<< std::endl;
   }