ShowerQuality Class Reference

Inheritance diagram for ShowerQuality:

Classes
struct	TreeParams_t
	For convenience: struct to define a set of parameters per shower to be stored in TTree. More...

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	ShowerQuality (fhicl::ParameterSet const &p)
	ShowerQuality (ShowerQuality const &)=delete
	ShowerQuality (ShowerQuality &&)=delete
ShowerQuality &	operator= (ShowerQuality const &)=delete
ShowerQuality &	operator= (ShowerQuality &&)=delete
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)
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)
std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
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 ()

Private Member Functions
void	analyze (art::Event const &e) override
void	beginJob () override
void	SetShowerProducer (const std::string name)
void	SetMCShowerProducer (const std::string name)
void	SetSimChannelProducer (const std::string name)
void	SetMaxEnergyCut (const double energy)
	Set maximum energy for MCShowers to be considered. More...
void	SetMinEnergyCut (const double energy)
	Set minimum energy for MCShowers to be considered. More...
void	InitializeAnaTree ()
	Function to prepare TTree. More...

Private Attributes
::btutil::MCMatchAlg	fBTAlg
	Shower back tracking algorithm. More...
double	_mc_energy_min
	Minimum MC shower energy cut. More...
double	_mc_energy_max
	Maximum MC shower energy cut. More...
std::string	fShowerProducer
	Shower Producer's Name. More...
std::string	fMCShowerProducer
	MCShower Producer's Name. More...
std::string	fSimChannelProducer
	SimChannel Producer's Name. More...
TH1D *	hMatchCorrectness
	Matching correctness. More...
TH1D *	hVtxDX
	X difference (reco-MC) in cm. More...
TH1D *	hVtxDY
	Y difference (reco-MC) in cm. More...
TH1D *	hVtxDZ
	Z difference (reco-MC) in cm. More...
TH1D *	hVtxDR
	3D vtx distance between reco to MC in cm More...
TH1D *	hDCosX
	Direction unit vector X component difference. More...
TH1D *	hDCosY
	Direction unit vector Y component difference. More...
TH1D *	hDCosZ
	Direction unit vector Z component difference. More...
TH1D *	h3DAngleDiff
	Opening angle between reco & MC 3D direction. More...
TH2D *	hEnergyCorr
	Energy correlation reco (x) vs. MC (y) More...
TH1D *	hEnergyAssym
	Energy assym. parameter: (reco E - MC E) / (reco E + MC E) * 2. More...
TH1D *	hEnergyDiff
	Energy difference: reco E - MC E. More...
TH1D *	hMatchedClusterEff
	Matched 3D shower's cluster efficiency (combined across planes) More...
TH1D *	hMatchedClusterPur
	Matched 3D shower's cluster purity (combined across planes) More...
std::map< int, TH1D * >	mDEDX
	dEdx per particle per PDG code More...
TH1D *	hBestPlane
	Best plane id. More...
struct ShowerQuality::TreeParams_t	fTreeParams
TTree *	fTree
	Analysis TTree. More...

Detailed Description

Definition at line 34 of file ShowerQuality_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 146 of file ShowerQuality_module.cc.

References fTree, fhicl::ParameterSet::get(), h3DAngleDiff, hBestPlane, hDCosX, hDCosY, hDCosZ, hEnergyAssym, hEnergyCorr, hEnergyDiff, hMatchCorrectness, hMatchedClusterEff, hMatchedClusterPur, hVtxDR, hVtxDX, hVtxDY, hVtxDZ, mDEDX, SetMaxEnergyCut(), SetMCShowerProducer(), SetMinEnergyCut(), SetShowerProducer(), and SetSimChannelProducer().

                                                       : EDAnalyzer(p) // ,
// More initializers here.
{

  //fShowerProducer   = "";
  //fMCShowerProducer = "";
  //fSimChannelProducer = "";
  SetShowerProducer(p.get<std::string>("ShowerProducer"));
  SetMCShowerProducer(p.get<std::string>("MCShowerProducer"));
  SetSimChannelProducer(p.get<std::string>("SimChannelProducer"));
  SetMinEnergyCut(p.get<double>("MCShowerEnergyMin"));
  SetMaxEnergyCut(p.get<double>("MCShowerEnergyMax"));

  hMatchCorrectness = nullptr;

  hVtxDX = nullptr;
  hVtxDY = nullptr;
  hVtxDZ = nullptr;
  hVtxDR = nullptr;

  hDCosX = nullptr;
  hDCosY = nullptr;
  hDCosZ = nullptr;
  h3DAngleDiff = nullptr;

  hEnergyCorr = nullptr;
  hEnergyAssym = nullptr;
  hEnergyDiff = nullptr;

  hMatchedClusterPur = nullptr;
  hMatchedClusterEff = nullptr;

  mDEDX.clear();
  hBestPlane = nullptr;

  fTree = nullptr;
}

ShowerQuality::ShowerQuality ( ShowerQuality const &  )

delete

ShowerQuality::ShowerQuality ( ShowerQuality &&  )

delete

Member Function Documentation

void ShowerQuality::analyze ( art::Event const &  e )

overrideprivatevirtual

Implements art::EDAnalyzer.

Definition at line 314 of file ShowerQuality_module.cc.

References _mc_energy_max, _mc_energy_min, ShowerQuality::TreeParams_t::best_plane_id, btutil::MCMatchAlg::BTAlg(), btutil::MCMatchAlg::BuildMap(), ShowerQuality::TreeParams_t::cluster_eff, ShowerQuality::TreeParams_t::cluster_pur, btutil::MCMatchAlg::ClusterEP(), energy, fBTAlg, fMCShowerProducer, fShowerProducer, fSimChannelProducer, fTree, fTreeParams, art::ProductRetriever::getValidHandle(), h3DAngleDiff, hBestPlane, hDCosX, hDCosY, hDCosZ, hEnergyAssym, hEnergyCorr, hEnergyDiff, hMatchCorrectness, hMatchedClusterEff, hMatchedClusterPur, hVtxDR, hVtxDX, hVtxDY, hVtxDZ, ShowerQuality::TreeParams_t::match_correctness, ShowerQuality::TreeParams_t::mc_containment, ShowerQuality::TreeParams_t::mc_dcosx, ShowerQuality::TreeParams_t::mc_dcosy, ShowerQuality::TreeParams_t::mc_dcosz, ShowerQuality::TreeParams_t::mc_energy, ShowerQuality::TreeParams_t::mc_pdgid, ShowerQuality::TreeParams_t::mc_reco_anglediff, ShowerQuality::TreeParams_t::mc_reco_dist, ShowerQuality::TreeParams_t::mc_x, ShowerQuality::TreeParams_t::mc_y, ShowerQuality::TreeParams_t::mc_z, btutil::MCBTAlg::MCQ(), mDEDX, art::Handle< T >::provenance(), ShowerQuality::TreeParams_t::reco_dcosx, ShowerQuality::TreeParams_t::reco_dcosy, ShowerQuality::TreeParams_t::reco_dcosz, ShowerQuality::TreeParams_t::reco_dedx, ShowerQuality::TreeParams_t::reco_energy, ShowerQuality::TreeParams_t::reco_x, ShowerQuality::TreeParams_t::reco_y, ShowerQuality::TreeParams_t::reco_z, btutil::MCMatchAlg::ShowerCorrectness(), and lar::dump::vector().

{
  // Retrieve mcshower data product
  auto mcsHandle = e.getValidHandle<std::vector<sim::MCShower>>(fMCShowerProducer);
  auto resHandle = e.getValidHandle<std::vector<recob::Shower>>(fShowerProducer);
  auto schHandle = e.getValidHandle<std::vector<sim::SimChannel>>(fSimChannelProducer);
  const std::vector<sim::MCShower>& ev_mcs(*mcsHandle);
  const std::vector<recob::Shower>& ev_shower(*resHandle);
  const std::vector<sim::SimChannel>& ev_simch(*schHandle);

  if (ev_shower.empty()) return;

  // Get the whole clusters + associated clusters
  art::Handle<std::vector<recob::Cluster>> clsHandle;
  art::FindManyP<recob::Cluster> cluster_m(resHandle, e, fShowerProducer);
  const std::vector<recob::Cluster>& ev_cluster = cluster_m.at(0).front().parentAs<std::vector>();

  // Make clusters in terms of hit vector to feed into BT algorithm
  art::FindManyP<recob::Hit> hit_m(clsHandle, e, clsHandle.provenance()->moduleLabel());
  std::vector<std::vector<art::Ptr<recob::Hit>>> ev_cluster_hit;
  ev_cluster_hit.reserve(clsHandle->size());
  std::map<art::Ptr<recob::Cluster>, size_t> cluster_ptr_map;
  for (size_t i = 0; i < ev_cluster.size(); ++i) {
    const art::Ptr<recob::Cluster> cluster_ptr(clsHandle, i);
    cluster_ptr_map[cluster_ptr] = ev_cluster_hit.size();
    ev_cluster_hit.push_back(hit_m.at(i));
  }

  // Create ass_cluster_v index vector
  std::vector<std::vector<unsigned int>> ass_cluster_v;
  ass_cluster_v.reserve(ev_shower.size());
  for (size_t shower_index = 0; shower_index < ev_shower.size(); ++shower_index) {
    ass_cluster_v.push_back(std::vector<unsigned int>());
    for (auto const& p : cluster_m.at(shower_index))
      ass_cluster_v.back().push_back(cluster_ptr_map[p]);
  }

  // Create G4 track ID vector for which we are interested in
  std::vector<std::vector<unsigned int>> g4_trackid_v;
  std::vector<unsigned int> mc_index_v;
  g4_trackid_v.reserve(ev_mcs.size());
  for (size_t mc_index = 0; mc_index < ev_mcs.size(); ++mc_index) {
    auto const& mcs = ev_mcs[mc_index];
    double energy = mcs.DetProfile().E();
    std::vector<unsigned int> id_v;
    id_v.reserve(mcs.DaughterTrackID().size());
    if (_mc_energy_min < energy && energy < _mc_energy_max) {
      for (auto const& id : mcs.DaughterTrackID()) {
        if (id == mcs.TrackID()) continue;
        id_v.push_back(id);
      }
      id_v.push_back(mcs.TrackID());
      g4_trackid_v.push_back(id_v);
      mc_index_v.push_back(mc_index);
    }
  }

  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
  if (!fBTAlg.BuildMap(clockData, g4_trackid_v, ev_simch, ev_cluster_hit)) {
    std::cerr << "\033[93m[ERROR]\033[00m <<ShowerQuality::analyze>> Failed to "
                 "build back-tracking map for MC..."
              << std::endl;
    return;
  }

  // Find the best-representative reco-ed Shower given an MCShower
  std::vector<std::vector<double>> shower_mcq_vv(ev_shower.size(),
                                                 std::vector<double>(mc_index_v.size(), 0));

  for (size_t shower_index = 0; shower_index < ass_cluster_v.size(); ++shower_index) {

    auto const& ass_cluster = ass_cluster_v[shower_index];

    std::vector<::btutil::WireRange_t> w_v;

    for (auto const& cluster_index : ass_cluster) {

      auto const& ass_hit = ev_cluster_hit[cluster_index];

      w_v.reserve(ass_hit.size() + w_v.size());

      for (auto const& hit_ptr : ass_hit) {
        w_v.emplace_back(hit_ptr->Channel(), hit_ptr->StartTick(), hit_ptr->EndTick());
      }
    }

    auto mcq_v = fBTAlg.BTAlg().MCQ(clockData, w_v);

    auto& shower_mcq_v = shower_mcq_vv[shower_index];

    for (size_t mcs_index = 0; mcs_index < (mcq_v.size() - 1); ++mcs_index) {

      shower_mcq_v[mcs_index] = mcq_v[mcs_index];
    }
  }

  // Loop over MCShower and inspect corresponding shower quality
  for (size_t mcs_index = 0; mcs_index < mc_index_v.size(); ++mcs_index) {

    auto const& mc_shower = ev_mcs[mc_index_v[mcs_index]];

    // Search for the best representative shower
    size_t best_shower_index = shower_mcq_vv.size();
    double max_mcq = 0;
    for (size_t shower_index = 0; shower_index < shower_mcq_vv.size(); ++shower_index) {

      if (shower_mcq_vv[shower_index][mcs_index] > max_mcq) best_shower_index = shower_index;
    }

    if (best_shower_index == shower_mcq_vv.size()) {
      std::string msg;
      std::cerr << "\033[93m[ERROR]\033[00m <<ShowerQuality::analyze>> "
                << "Failed to find a corresponding shower for MCShower " << mc_index_v[mcs_index]
                << std::endl;
      continue;
    }

    auto const& reco_shower = ev_shower[best_shower_index];

    auto res = fBTAlg.ShowerCorrectness(ass_cluster_v[best_shower_index]);

    fTreeParams.match_correctness = res.second;

    if (fTreeParams.match_correctness < 0) {
      std::cerr << "\033[93m[ERROR]\033[00m <<ShowerQuality::analyze>> "
                << "Failed to find a corresponding MCShower for shower " << best_shower_index
                << std::endl;
      continue;
    }

    // MC Info
    fTreeParams.mc_x = mc_shower.DetProfile().X();
    fTreeParams.mc_y = mc_shower.DetProfile().Y();
    fTreeParams.mc_z = mc_shower.DetProfile().Z();

    fTreeParams.mc_energy = mc_shower.DetProfile().E();
    fTreeParams.mc_pdgid = mc_shower.PdgCode();
    fTreeParams.mc_containment = mc_shower.DetProfile().E() / mc_shower.Start().E();

    //fTreeParams.mc_dcosx = mc_shower.DetProfile().Px() / fTreeParams.mc_energy;
    //fTreeParams.mc_dcosy = mc_shower.DetProfile().Py() / fTreeParams.mc_energy;
    //fTreeParams.mc_dcosz = mc_shower.DetProfile().Pz() / fTreeParams.mc_energy;
    fTreeParams.mc_dcosx = mc_shower.Start().Px() / mc_shower.Start().E();
    fTreeParams.mc_dcosy = mc_shower.Start().Py() / mc_shower.Start().E();
    fTreeParams.mc_dcosz = mc_shower.Start().Pz() / mc_shower.Start().E();

    // Reco vtx
    fTreeParams.reco_x = reco_shower.ShowerStart()[0];
    fTreeParams.reco_y = reco_shower.ShowerStart()[1];
    fTreeParams.reco_z = reco_shower.ShowerStart()[2];

    // Reco angle
    fTreeParams.reco_dcosx = reco_shower.Direction()[0];
    fTreeParams.reco_dcosy = reco_shower.Direction()[1];
    fTreeParams.reco_dcosz = reco_shower.Direction()[2];

    // Reco - MC angle diff
    fTreeParams.mc_reco_anglediff = acos(fTreeParams.reco_dcosx * fTreeParams.mc_dcosx +
                                         fTreeParams.reco_dcosy * fTreeParams.mc_dcosy +
                                         fTreeParams.reco_dcosz * fTreeParams.mc_dcosz) /
                                    3.14159265359 * 180.;
    // Reco - MC vtx distance
    fTreeParams.mc_reco_dist = sqrt(pow(fTreeParams.reco_x - fTreeParams.mc_x, 2) +
                                    pow(fTreeParams.reco_y - fTreeParams.mc_y, 2) +
                                    pow(fTreeParams.reco_z - fTreeParams.mc_z, 2));

    // Reco cluster efficiency & purity
    fTreeParams.cluster_eff = 1.;
    fTreeParams.cluster_pur = 1.;
    for (auto const& cluster_index : ass_cluster_v[best_shower_index]) {
      auto ep = fBTAlg.ClusterEP(cluster_index, mcs_index);
      if (ep.first == 0 && ep.second == 0) continue;
      fTreeParams.cluster_eff *= ep.first;
      fTreeParams.cluster_pur *= ep.second;
    }

    // Reco energy & dedx info
    fTreeParams.best_plane_id = reco_shower.best_plane();

    /*
      int best_plane_index = -1;

      for(size_t i=0; i < ass_cluster_v[best_shower_index].size(); ++i) {

        size_t cluster_index = ass_cluster_v[best_shower_index][i];
        //std::cout<<best_plane_index<<" : "<<ev_cluster->at(cluster_index).View()<<std::endl;
        if( ev_cluster->at(cluster_index).View() == reco_shower.best_plane() ) {
          best_plane_index = i;
          break;
        }
      }

      if(best_plane_index < 0) {
      throw ::showerreco::ShowerRecoException(Form("Failed to identify the best plane for shower %zu",
      best_shower_index)
      );
      }
    */

    fTreeParams.reco_energy = reco_shower.Energy().at(reco_shower.best_plane());
    //fTreeParams.reco_dedx_U   = reco_shower.dEdx().at(0);
    //fTreeParams.reco_dedx_V   = reco_shower.dEdx().at(1);
    //fTreeParams.reco_dedx_Y   = reco_shower.dEdx().at(2);
    fTreeParams.reco_dedx = reco_shower.dEdx().at(reco_shower.best_plane());

    //
    // Fill histograms
    //
    hMatchCorrectness->Fill(fTreeParams.match_correctness);

    hVtxDX->Fill(fTreeParams.reco_x - fTreeParams.mc_x);
    hVtxDY->Fill(fTreeParams.reco_y - fTreeParams.mc_y);
    hVtxDZ->Fill(fTreeParams.reco_z - fTreeParams.mc_z);
    hVtxDR->Fill(fTreeParams.mc_reco_dist);

    // Angular info
    hDCosX->Fill(fTreeParams.reco_dcosx - fTreeParams.mc_dcosx);
    hDCosY->Fill(fTreeParams.reco_dcosy - fTreeParams.mc_dcosy);
    hDCosZ->Fill(fTreeParams.reco_dcosz - fTreeParams.mc_dcosz);
    h3DAngleDiff->Fill(fTreeParams.mc_reco_anglediff);

    hEnergyCorr->Fill(fTreeParams.reco_energy, fTreeParams.mc_energy);

    hEnergyAssym->Fill((fTreeParams.reco_energy - fTreeParams.mc_energy) /
                       (fTreeParams.reco_energy + fTreeParams.mc_energy) * 2.);

    hEnergyDiff->Fill(fTreeParams.mc_energy - fTreeParams.reco_energy);

    if (mDEDX.find(fTreeParams.mc_pdgid) == mDEDX.end())

      mDEDX.insert(std::make_pair(
        fTreeParams.mc_pdgid,
        new TH1D(Form("hdEdx_PDG_%d", fTreeParams.mc_pdgid),
                 Form("Reco dE/dx for PDG = %d; dE/dx [MeV/cm]; Showers", fTreeParams.mc_pdgid),
                 100,
                 0,
                 50)));

    hMatchedClusterEff->Fill(fTreeParams.cluster_eff);
    hMatchedClusterPur->Fill(fTreeParams.cluster_pur);

    mDEDX[fTreeParams.mc_pdgid]->Fill(fTreeParams.reco_dedx);

    hBestPlane->Fill(fTreeParams.best_plane_id);

    // Fill Tree
    fTree->Fill();
  }
}

void ShowerQuality::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 184 of file ShowerQuality_module.cc.

References fMCShowerProducer, fShowerProducer, fSimChannelProducer, fTree, Get, h3DAngleDiff, hBestPlane, hDCosX, hDCosY, hDCosZ, hEnergyAssym, hEnergyCorr, hEnergyDiff, hMatchCorrectness, hMatchedClusterEff, hMatchedClusterPur, hVtxDR, hVtxDX, hVtxDY, hVtxDZ, and InitializeAnaTree().

{

  if (fShowerProducer.empty() || fMCShowerProducer.empty() || fSimChannelProducer.empty()) {
    std::string msg;
    msg += "\033[93m[ERROR]\033[00m <<";
    msg += __FUNCTION__;
    msg += ">> Producer's name not set!";
    std::cout << msg.c_str() << std::endl;
    throw ::showerreco::ShowerRecoException(msg.c_str());
  }

  art::ServiceHandle<art::TFileService const> tfs;

  if (fTree) delete fTree;
  fTree = tfs->make<TTree>("fShowerQualityTree", "");

  //
  // Matching correctness histogram initialization
  //
  if (hMatchCorrectness) delete hMatchCorrectness;
  hMatchCorrectness =
    tfs->make<TH1D>("hMatchCorrectness",
                    "Shower 2D Cluster Matching Correctness; Correctness; Showers",
                    101,
                    -0.005,
                    1.005);

  //
  // 3D Vtx (start point) MC/Reco comparison histogram initialization
  //
  if (hVtxDX) delete hVtxDX;
  if (hVtxDY) delete hVtxDY;
  if (hVtxDZ) delete hVtxDZ;
  if (hVtxDR) delete hVtxDR;

  hVtxDX = tfs->make<TH1D>(
    "hVtxDX", "Reco - MC Start X [cm] Displacement; #DeltaX [cm]; Showers", 200, -100, 100);

  hVtxDY = tfs->make<TH1D>(
    "hVtxDY", "Reco - MC Start Y [cm] Displacement; #DeltaY [cm]; Showers", 200, -100, 100);

  hVtxDZ = tfs->make<TH1D>(
    "hVtxDZ", "Reco - MC Start Z [cm] Displacement; #DeltaZ [cm]; Showers", 200, -100, 100);

  hVtxDR = tfs->make<TH1D>(
    "hVtxDR", "Reco - MC Start 3D Vtx Displacement; #DeltaR [cm]; Showers", 200, -100, 100);

  //
  // 3D Angular MC/Reco comparison histogram initialization
  //
  if (hDCosX) delete hDCosX;
  if (hDCosY) delete hDCosY;
  if (hDCosZ) delete hDCosZ;
  if (h3DAngleDiff) delete h3DAngleDiff;

  hDCosX = tfs->make<TH1D>(
    "hDCosX", "Direction Unit Vector Reco - MC #DeltaX; #DeltaCosX; Showers", 100, -2, 2);

  hDCosY = tfs->make<TH1D>(
    "hDCosY", "Direction Unit Vector Reco - MC #DeltaY; #DeltaCosY; Showers", 100, -2, 2);

  hDCosZ = tfs->make<TH1D>(
    "hDCosZ", "Direction Unit Vector Reco - MC #DeltaZ; #DeltaCosZ; Showers", 100, -2, 2);

  h3DAngleDiff =
    tfs->make<TH1D>("h3DAngleDiff",
                    "3D Opening Angle Between Reco & MC; Opening Angle [degrees]; Showers",
                    181,
                    -0.5,
                    180.5);

  //
  // Energy MC/Reco comparison histogram initialization
  //
  if (hEnergyCorr) delete hEnergyCorr;
  if (hEnergyAssym) delete hEnergyAssym;
  if (hEnergyDiff) delete hEnergyDiff;

  hEnergyCorr =
    tfs->make<TH2D>("hEnergyCorr",
                    "Reco (x) vs. MC (y) Energy Comparison; Reco Energy [MeV]; MC Energy [MeV]",
                    200,
                    0,
                    1000,
                    200,
                    0,
                    1000);

  hEnergyAssym = tfs->make<TH1D>("hEnergyAssym",
                                 "MC - Reco Energy Fractional Difference; Assymetry; Showers",
                                 201,
                                 -1.005,
                                 1.005);

  hEnergyDiff = tfs->make<TH1D>(
    "hEnergyDiff", "MC - Reco Energy Difference; Energy Difference [MeV]; Showers", 200, 0, 1000);

  //
  // Shower cluster purity & efficiency histograms initialization
  //
  if (hMatchedClusterEff) delete hMatchedClusterEff;
  if (hMatchedClusterPur) delete hMatchedClusterPur;

  hMatchedClusterEff =
    tfs->make<TH1D>("hMatchedClusterEff_PlaneCombo",
                    "Matched Shower Cluster's Charge Efficiency; Efficiency; Clusters",
                    101,
                    -0.005,
                    1.005);

  hMatchedClusterPur = tfs->make<TH1D>("hMatchedClusterPur_PlaneCombo",
                                       "Matched Shower Cluster's Charge Purity; Purity; Clusters",
                                       101,
                                       -0.005,
                                       1.005);

  //
  // Best plane ID histogram initialization
  //
  auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout>()->Get();
  hBestPlane = tfs->make<TH1D>("hBestPlane",
                               "Best Plane (for energy & dE/dx estimate); Plane ID; Showers",
                               wireReadoutGeom.Nplanes(),
                               -0.5,
                               wireReadoutGeom.Nplanes() - 0.5);

  InitializeAnaTree();
}

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const &  tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

  {
    return collector_.consumes<T, BT>(tag);
  }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

  {
    return collector_;
  }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

  {
    collector_.consumesMany<T, BT>();
  }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const &  )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const &  tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

  {
    return collector_.consumesView<T, BT>(tag);
  }

void art::detail::Analyzer::doBeginJob ( SharedResources const &  resources )

inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

  {
    setupQueues(resources);
    ProcessingFrame const frame{ScheduleID{}};
    beginJobWithFrame(frame);
  }

bool art::detail::Analyzer::doBeginRun ( RunPrincipal &  rp, ModuleContext const &  mc  )

inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{mc.scheduleID()};
    beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
    return true;
  }

bool art::detail::Analyzer::doBeginSubRun ( SubRunPrincipal &  srp, ModuleContext const &  mc  )

inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{mc.scheduleID()};
    beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
    return true;
  }

void art::detail::Analyzer::doEndJob ( )

inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{ScheduleID{}};
    endJobWithFrame(frame);
  }

bool art::detail::Analyzer::doEndRun ( RunPrincipal &  rp, ModuleContext const &  mc  )

inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{mc.scheduleID()};
    endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
    return true;
  }

bool art::detail::Analyzer::doEndSubRun ( SubRunPrincipal &  srp, ModuleContext const &  mc  )

inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{mc.scheduleID()};
    endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
    return true;
  }

bool art::detail::Analyzer::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  )

inherited

Definition at line 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

  {
    auto const e = std::as_const(ep).makeEvent(mc);
    if (wantEvent(mc.scheduleID(), e)) {
      ++counts_run;
      ProcessingFrame const frame{mc.scheduleID()};
      analyzeWithFrame(e, frame);
      ++counts_passed;
    }
    return true;
  }

void art::detail::Analyzer::doRespondToCloseInputFile ( FileBlock const &  fb )

inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToCloseInputFileWithFrame(fb, frame);
  }

void art::detail::Analyzer::doRespondToCloseOutputFiles ( FileBlock const &  fb )

inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToCloseOutputFilesWithFrame(fb, frame);
  }

void art::detail::Analyzer::doRespondToOpenInputFile ( FileBlock const &  fb )

inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToOpenInputFileWithFrame(fb, frame);
  }

void art::detail::Analyzer::doRespondToOpenOutputFiles ( FileBlock const &  fb )

inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToOpenOutputFilesWithFrame(fb, frame);
  }

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

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

  {
    return collector_.getConsumables();
  }

Handle< TriggerResults > art::Observer::getTriggerResults ( Event const &  e ) const

protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

  {
    if (selectors_) {
      return selectors_->getOneTriggerResults(e);
    }

    // The following applies for cases where no SelectEvents entries
    // exist.
    Handle<TriggerResults> h;
    if (e.get(empty_process_name, h)) {
      return h;
    }
    return Handle<TriggerResults>{};
  }

void ShowerQuality::InitializeAnaTree ( )

private

Function to prepare TTree.

Definition at line 564 of file ShowerQuality_module.cc.

References ShowerQuality::TreeParams_t::best_plane_id, ShowerQuality::TreeParams_t::cluster_eff, ShowerQuality::TreeParams_t::cluster_pur, DEFINE_ART_MODULE, fTree, fTreeParams, ShowerQuality::TreeParams_t::match_correctness, ShowerQuality::TreeParams_t::mc_containment, ShowerQuality::TreeParams_t::mc_dcosx, ShowerQuality::TreeParams_t::mc_dcosy, ShowerQuality::TreeParams_t::mc_dcosz, ShowerQuality::TreeParams_t::mc_energy, ShowerQuality::TreeParams_t::mc_pdgid, ShowerQuality::TreeParams_t::mc_reco_anglediff, ShowerQuality::TreeParams_t::mc_reco_dist, ShowerQuality::TreeParams_t::mc_x, ShowerQuality::TreeParams_t::mc_y, ShowerQuality::TreeParams_t::mc_z, ShowerQuality::TreeParams_t::reco_dcosx, ShowerQuality::TreeParams_t::reco_dcosy, ShowerQuality::TreeParams_t::reco_dcosz, ShowerQuality::TreeParams_t::reco_dedx, ShowerQuality::TreeParams_t::reco_dedx_U, ShowerQuality::TreeParams_t::reco_dedx_V, ShowerQuality::TreeParams_t::reco_dedx_Y, ShowerQuality::TreeParams_t::reco_energy, ShowerQuality::TreeParams_t::reco_x, ShowerQuality::TreeParams_t::reco_y, and ShowerQuality::TreeParams_t::reco_z.

Referenced by beginJob().

{

  fTree->Branch("reco_x", &fTreeParams.reco_x, "reco_x/D");
  fTree->Branch("reco_y", &fTreeParams.reco_y, "reco_y/D");
  fTree->Branch("reco_z", &fTreeParams.reco_z, "reco_z/D");
  fTree->Branch("reco_dcosx", &fTreeParams.reco_dcosx, "reco_dcosx/D");
  fTree->Branch("reco_dcosy", &fTreeParams.reco_dcosy, "reco_dcosy/D");
  fTree->Branch("reco_dcosz", &fTreeParams.reco_dcosz, "reco_dcosz/D");
  fTree->Branch("reco_energy", &fTreeParams.reco_energy, "reco_energy/D");

  fTree->Branch("best_plane_id", &fTreeParams.best_plane_id, "best_plane_id/i");

  fTree->Branch("mc_x", &fTreeParams.mc_x, "mc_x/D");
  fTree->Branch("mc_y", &fTreeParams.mc_y, "mc_y/D");
  fTree->Branch("mc_z", &fTreeParams.mc_z, "mc_z/D");
  fTree->Branch("mc_dcosx", &fTreeParams.mc_dcosx, "mc_dcosx/D");
  fTree->Branch("mc_dcosy", &fTreeParams.mc_dcosy, "mc_dcosy/D");
  fTree->Branch("mc_dcosz", &fTreeParams.mc_dcosz, "mc_dcosz/D");
  fTree->Branch("mc_energy", &fTreeParams.mc_energy, "mc_energy/D");

  fTree->Branch("reco_dedx", &fTreeParams.reco_dedx, "reco_dedx_/D");
  fTree->Branch("reco_dedx_U", &fTreeParams.reco_dedx_U, "reco_dedx_U/D");
  fTree->Branch("reco_dedx_V", &fTreeParams.reco_dedx_V, "reco_dedx_V/D");
  fTree->Branch("reco_dedx_Y", &fTreeParams.reco_dedx_Y, "reco_dedx_Y/D");
  fTree->Branch("mc_pdgid", &fTreeParams.mc_pdgid, "mc_pdgid/i");

  fTree->Branch("mc_reco_anglediff", &fTreeParams.mc_reco_anglediff, "mc_reco_anglediff/D");
  fTree->Branch("mc_reco_dist", &fTreeParams.mc_reco_dist, "mc_reco_dist/D");

  fTree->Branch("mc_containment", &fTreeParams.mc_containment, "mc_containment/D");

  fTree->Branch("match_correctness", &fTreeParams.match_correctness, "match_correctness/D");
  fTree->Branch("cluster_eff", &fTreeParams.cluster_eff, "cluster_eff/D");
  fTree->Branch("cluster_pur", &fTreeParams.cluster_pur, "cluster_pur/D");
}

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const &  wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

  {
    return doMakeWorker(wp);
  }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const &  tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

  {
    return collector_.mayConsume<T, BT>(tag);
  }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

  {
    collector_.mayConsumeMany<T, BT>();
  }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const &  )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const &  tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

  {
    return collector_.mayConsumeView<T, BT>(tag);
  }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

  {
    if (md_.has_value()) {
      return *md_;
    }

    throw Exception{errors::LogicError,
                    "There was an error while calling moduleDescription().\n"}
      << "The moduleDescription() base-class member function cannot be called\n"
         "during module construction.  To determine which module is "
         "responsible\n"
         "for calling it, find the '<module type>:<module "
         "label>@Construction'\n"
         "tag in the message prefix above.  Please contact artists@fnal.gov\n"
         "for guidance.\n";
  }

ShowerQuality& ShowerQuality::operator= ( ShowerQuality const &  )

delete

ShowerQuality& ShowerQuality::operator= ( ShowerQuality &&  )

delete

string const & art::Observer::processName ( ) const

protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

  {
    return process_name_;
  }

void ShowerQuality::SetMaxEnergyCut ( const double  energy )

inlineprivate

Set maximum energy for MCShowers to be considered.

Definition at line 60 of file ShowerQuality_module.cc.

References _mc_energy_max, and energy.

Referenced by ShowerQuality().

{ _mc_energy_max = energy; }

void ShowerQuality::SetMCShowerProducer ( const std::string  name )

inlineprivate

Definition at line 56 of file ShowerQuality_module.cc.

References fMCShowerProducer.

Referenced by ShowerQuality().

{ fMCShowerProducer = name; }

void ShowerQuality::SetMinEnergyCut ( const double  energy )

inlineprivate

Set minimum energy for MCShowers to be considered.

Definition at line 63 of file ShowerQuality_module.cc.

References _mc_energy_min, and energy.

Referenced by ShowerQuality().

{ _mc_energy_min = energy; }

void art::ModuleBase::setModuleDescription ( ModuleDescription const &  md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

  {
    md_ = md;
  }

void ShowerQuality::SetShowerProducer ( const std::string  name )

inlineprivate

Setter function for a shower producer name.

Definition at line 55 of file ShowerQuality_module.cc.

References fShowerProducer.

Referenced by ShowerQuality().

{ fShowerProducer = name; }

void ShowerQuality::SetSimChannelProducer ( const std::string  name )

inlineprivate

Definition at line 57 of file ShowerQuality_module.cc.

References fSimChannelProducer.

Referenced by ShowerQuality().

{ fSimChannelProducer = name; }

void art::ModuleBase::sortConsumables ( std::string const &  current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

  {
    // Now that we know we have seen all the consumes declarations,
    // sort the results for fast lookup later.
    collector_.sortConsumables(current_process_name);
  }

bool art::Observer::wantAllEvents ( ) const

inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

    {
      return wantAllEvents_;
    }

bool art::Observer::wantEvent ( ScheduleID  id, Event const &  e  ) const

protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

Referenced by art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

  {
    if (wantAllEvents_) {
      return true;
    }
    bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
    bool const reject_event =
      rejectors_ ? rejectors_->matchEvent(id, e) : false;
    return select_event and not reject_event;
  }

Member Data Documentation

double ShowerQuality::_mc_energy_max

private

Maximum MC shower energy cut.

Definition at line 72 of file ShowerQuality_module.cc.

Referenced by analyze(), and SetMaxEnergyCut().

double ShowerQuality::_mc_energy_min

private

Minimum MC shower energy cut.

Definition at line 69 of file ShowerQuality_module.cc.

Referenced by analyze(), and SetMinEnergyCut().

::btutil::MCMatchAlg ShowerQuality::fBTAlg

private

Shower back tracking algorithm.

Definition at line 66 of file ShowerQuality_module.cc.

Referenced by analyze().

std::string ShowerQuality::fMCShowerProducer

private

MCShower Producer's Name.

Definition at line 78 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and SetMCShowerProducer().

std::string ShowerQuality::fShowerProducer

private

Shower Producer's Name.

Definition at line 75 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and SetShowerProducer().

std::string ShowerQuality::fSimChannelProducer

private

SimChannel Producer's Name.

Definition at line 81 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and SetSimChannelProducer().

TTree* ShowerQuality::fTree

private

Analysis TTree.

Definition at line 138 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), InitializeAnaTree(), and ShowerQuality().

struct ShowerQuality::TreeParams_t ShowerQuality::fTreeParams

private

Referenced by analyze(), and InitializeAnaTree().

TH1D* ShowerQuality::h3DAngleDiff

private

Opening angle between reco & MC 3D direction.

Definition at line 94 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hBestPlane

private

Best plane id.

Definition at line 108 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hDCosX

private

Direction unit vector X component difference.

Definition at line 91 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hDCosY

private

Direction unit vector Y component difference.

Definition at line 92 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hDCosZ

private

Direction unit vector Z component difference.

Definition at line 93 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hEnergyAssym

private

Energy assym. parameter: (reco E - MC E) / (reco E + MC E) * 2.

Definition at line 98 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH2D* ShowerQuality::hEnergyCorr

private

Energy correlation reco (x) vs. MC (y)

Definition at line 96 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hEnergyDiff

private

Energy difference: reco E - MC E.

Definition at line 99 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hMatchCorrectness

private

Matching correctness.

Definition at line 84 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hMatchedClusterEff

private

Matched 3D shower's cluster efficiency (combined across planes)

Definition at line 101 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hMatchedClusterPur

private

Matched 3D shower's cluster purity (combined across planes)

Definition at line 102 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDR

private

3D vtx distance between reco to MC in cm

Definition at line 89 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDX

private

X difference (reco-MC) in cm.

Definition at line 86 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDY

private

Y difference (reco-MC) in cm.

Definition at line 87 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

TH1D* ShowerQuality::hVtxDZ

private

Z difference (reco-MC) in cm.

Definition at line 88 of file ShowerQuality_module.cc.

Referenced by analyze(), beginJob(), and ShowerQuality().

std::map<int, TH1D*> ShowerQuality::mDEDX

private

dEdx per particle per PDG code

Definition at line 105 of file ShowerQuality_module.cc.

Referenced by analyze(), and ShowerQuality().

---

The documentation for this class was generated from the following file:

• larreco/v10_01_12/source/larreco/ShowerFinder/ShowerReco3D/ShowerQuality_module.cc