test_FixedDetPedestalService.cxx File Reference

#include "../FixedDetPedestalService.h"
#include "dunecore/ArtSupport/ArtServiceHelper.h"
#include "larcore/Geometry/Geometry.h"
#include <string>
#include <iostream>
#include <fstream>
#include <vector>
#include <iomanip>
#include <cassert>

Go to the source code of this file.

Functions
int	test_FixedDetPedestalService ()
int	main (int argc, char *argv[])

Function Documentation

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 89 of file test_FixedDetPedestalService.cxx.

                                 {
  return test_FixedDetPedestalService();
}

int test_FixedDetPedestalService

(

)

Definition at line 31 of file test_FixedDetPedestalService.cxx.

                                   {
  const string myname = "test_FixedDetPedestalService: ";
#ifdef NDEBUG
  cout << myname << "NDEBUG must be off." << endl;
  abort();
#endif
  string line = "-----------------------------";

  cout << myname << line << endl;
  cout << myname << "Create top-level FCL." << endl;
  std::ostringstream oss;
  oss << "#include \"pedestals_dune.fcl\"" << endl;
  oss << "#include \"geometry_dune.fcl\"" << endl;
  oss << "services.DetPedestalService: @local::dune_fixedpeds" << endl;
  oss << "services.Geometry: @local::dune35t_geo" << endl;
  oss << "services.ExptGeoHelperInterface: @local::dune_geometry_helper" << endl;
  ArtServiceHelper::load_services(oss.str());

  cout << myname << line << endl;
  cout << myname << "Fetch geometry service." << endl;
  ServiceHandle<Geometry> hgeo;
  cout << myname << "Detector: " << hgeo->DetectorName() << endl;

  cout << myname << line << endl;
  cout << myname << "Fetch pedestal service." << endl;
  ServiceHandle<DetPedestalService> hdps;
  const DetPedestalProvider* pdpp = hdps->provider();
  assert( pdpp != nullptr );

  cout << myname << line << endl;
  cout << myname << "Checking pedestals for selected channels." << endl;
  vector<int> chans = {0, 1, 200, 300, 400, 500, 600, 700, 800, 900};
  int hchan = 10;
  int hview = 10;
  int hmean = 12;
  cout << setw(hchan) << "Channel"
       << setw(hview) << "View"
       << setw(hmean) << "Mean"
       << setw(hmean) << "RMS"
       << endl;
  for ( int chan : chans ) {
    cout << setw(hchan) << chan
         << setw(hview) << hgeo->View(chan)
         << setw(hmean) << pdpp->PedMean(chan)
         << setw(hmean) << pdpp->PedRms(chan)
         << endl;
    float meanExp = 1800;
    if ( hgeo->View(chan) == geo::kZ ) {
      meanExp = 500.0;
    }
    assert( pdpp->PedMean(chan) == meanExp );
  }

  cout << myname << line << endl;
  cout << myname << "Done." << endl;
  return 0;
}