sbndcode/sbndcode/gallery/galleryAnalysis/galleryAnalysis.cpp File Reference

#include "TrackAnalysis.h"
#include "HitAnalysisAlg.h"
#include "MCAssociations.h"
#include "sbndcode/gallery/helpers/expandInputFiles.h"
#include "lardataobj/RecoBase/Track.h"
#include "lardataalg/DetectorInfo/DetectorPropertiesStandardTestHelpers.h"
#include "lardataalg/DetectorInfo/DetectorPropertiesStandard.h"
#include "lardataalg/DetectorInfo/DetectorClocksStandardTestHelpers.h"
#include "lardataalg/DetectorInfo/DetectorClocksStandard.h"
#include "lardataalg/DetectorInfo/LArPropertiesStandardTestHelpers.h"
#include "lardataalg/DetectorInfo/LArPropertiesStandard.h"
#include "larcorealg/Geometry/StandaloneGeometrySetup.h"
#include "larcorealg/Geometry/GeometryCore.h"
#include "sbndcode/Geometry/ChannelMapSBNDAlg.h"
#include "larcorealg/Geometry/StandaloneBasicSetup.h"
#include "gallery/Event.h"
#include "canvas/Utilities/InputTag.h"
#include "messagefacility/MessageLogger/MessageLogger.h"
#include "fhiclcpp/ParameterSet.h"
#include "TFile.h"
#include <string>
#include <vector>
#include <memory>
#include <iostream>
#include <algorithm>
#include <iterator>

Go to the source code of this file.

Functions
int	galleryAnalysis (std::string const &configFile, std::vector< std::string > const &inputFiles)
	Runs the analysis macro. More...
int	galleryAnalysis (std::string const &configFile, std::string filename)
	Version with a single input file. More...
int	main (int argc, char **argv)

Function Documentation

int galleryAnalysis	(	std::string const &	configFile,
		std::vector< std::string > const &	inputFiles
	)

Runs the analysis macro.

Parameters

configFile	path to the FHiCL configuration to be used for the services
inputFiles	vector of path of file names

Returns

an integer as exit code (0 means success)

Definition at line 75 of file sbndcode/sbndcode/gallery/galleryAnalysis/galleryAnalysis.cpp.

{
    /*
     * the "test" environment configuration
     */
    // read FHiCL configuration from a configuration file:
    fhicl::ParameterSet config = lar::standalone::ParseConfiguration(configFile);
  
    // set up message facility (always picked from "services.message")
    lar::standalone::SetupMessageFacility(config, "galleryAnalysis");
  
    // configuration from the "analysis" table of the FHiCL configuration file:
    auto const& analysisConfig = config.get<fhicl::ParameterSet>("analysis");
  
    // ***************************************************************************
    // ***  SERVICE PROVIDER SETUP BEGIN  ****************************************
    // ***************************************************************************
    //
    // Uncomment the things you need
    // (and make sure the corresponding headers are also uncommented)
    //
  
    // geometry setup (it's special)
    auto geom = lar::standalone::SetupGeometry<geo::ChannelMapSBNDAlg>(config.get<fhicl::ParameterSet>("services.Geometry"));
  
    // LArProperties setup
    auto larp = testing::setupProvider<detinfo::LArPropertiesStandard>(config.get<fhicl::ParameterSet>("services.LArPropertiesService"));
  
    // DetectorClocks setup
    auto detclk = testing::setupProvider<detinfo::DetectorClocksStandard>(config.get<fhicl::ParameterSet>("services.DetectorClocksService"));
  
    // DetectorProperties setup
    auto detp = testing::setupProvider<detinfo::DetectorPropertiesStandard>(config.get<fhicl::ParameterSet>("services.DetectorPropertiesService"),
                                                                            detinfo::DetectorPropertiesStandard::providers_type{geom.get(),static_cast<detinfo::LArProperties const*>(larp.get())}); // TODO type cast is required until issue #18001 is solved
  
    // ***************************************************************************
    // ***  SERVICE PROVIDER SETUP END    ****************************************
    // ***************************************************************************
  
    /*
     * the preparation of input file list
     */
    std::vector<std::string> const allInputFiles = expandInputFiles(inputFiles);
  
    /*
     * other parameters
     */
    art::InputTag trackTag = analysisConfig.get<art::InputTag>("tracks");
    art::InputTag hitsTag  = analysisConfig.get<art::InputTag>("hits");
  
    /*
     * preparation of histogram output file
     */
    std::unique_ptr<TFile> pHistFile;
    if (analysisConfig.has_key("histogramFile"))
    {
        std::string fileName = analysisConfig.get<std::string>("histogramFile");
        std::cout << "Creating output file: '" << fileName << "'" << std::endl;
        pHistFile = std::make_unique<TFile>(fileName.c_str(), "RECREATE");
    }
  
    /*
     * preparation of the algorithm class
     */
    TrackAnalysis trackAnalysis(analysisConfig.get<fhicl::ParameterSet>("trackAnalysis"));
  
    trackAnalysis.setup(*geom, pHistFile.get());
    trackAnalysis.prepare();
    auto const clock_data = detclk->DataForJob();
    auto const det_prop_data = detp->DataFor(clock_data);
    HitAnalysis::HitAnalysisAlg hitAnalysisAlg(analysisConfig.get<fhicl::ParameterSet>("hitAnalysisAlg"));
    
    hitAnalysisAlg.setup(*geom, pHistFile.get());
    
    MCAssociations mcAssociations(analysisConfig.get<fhicl::ParameterSet>("mcAssociations"));
    
    mcAssociations.setup(*geom, det_prop_data, pHistFile.get());
    mcAssociations.prepare();
    
    int numEvents(0);
  
    /*
     * the event loop
     */
    for (gallery::Event event(allInputFiles); !event.atEnd(); event.next())
    {
        // *************************************************************************
        // ***  SINGLE EVENT PROCESSING BEGIN  *************************************
        // *************************************************************************
    
        mf::LogVerbatim("galleryAnalysis") << "This is event " << event.fileEntry() << "-" << event.eventEntry();
    
        trackAnalysis.processTracks(*(event.getValidHandle<std::vector<recob::Track>>(trackTag)));
        
        hitAnalysisAlg.fillHistograms(*(event.getValidHandle<std::vector<recob::Hit>>(hitsTag)));
        
        mcAssociations.doTrackHitMCAssociations(event);

        numEvents++;
    
        // *************************************************************************
        // ***  SINGLE EVENT PROCESSING END    *************************************
        // *************************************************************************
    
    } // for
  
    trackAnalysis.finish();
    mcAssociations.finish();
    
    hitAnalysisAlg.endJob(numEvents);
  
    return 0;
} // galleryAnalysis()

int galleryAnalysis	(	std::string const &	configFile,
		std::string	filename
	)

Version with a single input file.

Definition at line 191 of file sbndcode/sbndcode/gallery/galleryAnalysis/galleryAnalysis.cpp.

  { return galleryAnalysis(configFile, std::vector<std::string>{ filename }); }

int main	(	int	argc,
		char **	argv
	)

Definition at line 195 of file sbndcode/sbndcode/gallery/galleryAnalysis/galleryAnalysis.cpp.

                                {
  
  char **pParam = argv + 1, **pend = argv + argc;
  if (pParam == pend) {
    std::cerr << "Usage: " << argv[0] << "  configFile [inputFile ...]"
      << std::endl;
    return 1;
  }
  std::string const configFile = *(pParam++);
  std::vector<std::string> fileNames;
  std::copy(pParam, pend, std::back_inserter(fileNames));
  
  return galleryAnalysis(configFile, fileNames);
} // main()