Inheritance diagram for ClusteringValidation::ClusteringValidation:

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

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

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 )

private

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 ( )

private

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 ( )

private

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();
 }