Inheritance diagram for ClusteringValidation::ClusteringValidation:

Public Member Functions
	ClusteringValidation (fhicl::ParameterSet const &p)

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 Member Functions
void	analyze (art::Event const &evt)

void	beginJob ()

void	endJob ()

Private Attributes
std::vector< std::string >	fClusterModuleLabels

std::string	fHitsModuleLabel

int	fMinHitsInPlane

TCanvas *	fCanvas

std::map< std::string, std::unique_ptr< ClusterAnalyser > >	clusterAnalysis

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 644 of file ClusteringValidation_module.cc.

Constructor & Destructor Documentation

ClusteringValidation::ClusteringValidation::ClusteringValidation ( fhicl::ParameterSet const & p )

explicit

Definition at line 668 of file ClusteringValidation_module.cc.

   : EDAnalyzer(pset)
 {
   fMinHitsInPlane = pset.get<int>("MinHitsInPlane");
   fClusterModuleLabels = pset.get<std::vector<std::string>>("ClusterModuleLabels");
   fHitsModuleLabel = pset.get<std::string>("HitsModuleLabel");
 
   fCanvas = new TCanvas("fCanvas", "", 800, 600);
   gStyle->SetOptStat(0);
 }

Member Function Documentation

void ClusteringValidation::ClusteringValidation::analyze ( art::Event const & evt )

privatevirtual

Implements art::EDAnalyzer.

Definition at line 680 of file ClusteringValidation_module.cc.

 {
   // Get the hits from the event
   art::Handle<std::vector<recob::Hit>> hitHandle;
   std::vector<art::Ptr<recob::Hit>> hits;
   if (evt.getByLabel(fHitsModuleLabel, hitHandle)) art::fill_ptr_vector(hits, hitHandle);
 
   // Get clustering information from event
   // and give to the ClusterAnalyser to analyse
   for (auto clustering : fClusterModuleLabels) {
 
     // Get the clusters from the event
     art::Handle<std::vector<recob::Cluster>> clusterHandle;
     std::vector<art::Ptr<recob::Cluster>> clusters;
     if (evt.getByLabel(clustering, clusterHandle)) art::fill_ptr_vector(clusters, clusterHandle);
 
     // Find the associations (the hits) for the clusters
     art::FindManyP<recob::Hit> fmh(clusterHandle, evt, clustering);
 
     // Analyse this particular clustering
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
     clusterAnalysis.at(clustering)->Analyse(clockData, hits, clusters, fmh, fMinHitsInPlane);
   }
 }

void ClusteringValidation::ClusteringValidation::beginJob ( )

privatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 706 of file ClusteringValidation_module.cc.

 {
 
   // Construct the cluster analysers here
   // -- one for each of the clustering types to compare
   for (auto clustering : fClusterModuleLabels)
     clusterAnalysis[clustering] = (std::unique_ptr<ClusterAnalyser>)new ClusterAnalyser(clustering);
 }

void ClusteringValidation::ClusteringValidation::endJob ( )

privatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 716 of file ClusteringValidation_module.cc.

 {
 
   // Make a map of all the histograms for each clustering
   std::map<std::string, TObjArray> allHistograms;
   for (auto clustering : fClusterModuleLabels)
     allHistograms[clustering] = clusterAnalysis.at(clustering)->GetHistograms();
 
   // Write histograms to file
   TFile* file = TFile::Open("validationHistograms.root", "RECREATE");
   for (auto clustering : fClusterModuleLabels) {
     file->cd();
     file->mkdir(clustering.c_str());
     file->cd(clustering.c_str());
     for (int histIt = 0; histIt < allHistograms.begin()->second.GetEntriesFast(); ++histIt)
       allHistograms.at(clustering).At(histIt)->Write();
   }
 
   // Write images of overlaid histograms
   for (int histIt = 0; histIt < allHistograms.begin()->second.GetEntriesFast(); ++histIt) {
     fCanvas->cd();
     fCanvas->Clear();
     const char* name = allHistograms.begin()->second.At(histIt)->GetName();
     TLegend* l = new TLegend(0.6, 0.8, 0.8, 0.9, name, "brNDC");
     int clusterings = 1;
     for (std::map<std::string, TObjArray>::iterator clusteringIt = allHistograms.begin();
          clusteringIt != allHistograms.end();
          ++clusteringIt, ++clusterings) {
       TH1* h = (TH1*)allHistograms.at(clusteringIt->first).At(histIt);
       h->SetLineColor(clusterings);
       h->SetMarkerColor(clusterings);
       if (clusterings == 1)
         h->Draw();
       else
         h->Draw("same");
       l->AddEntry(h, clusteringIt->first.c_str(), "lp");
     }
     l->Draw("same");
     //fCanvas->SaveAs(name+TString(".png"));
     file->cd();
     fCanvas->Write(name);
   }
 
   file->Close();
   delete file;
 
   if (clusterAnalysis.find("blurredclusteringdc") != clusterAnalysis.end())
     clusterAnalysis.at("blurredclusteringdc")->WriteFile();
 }