#include "FlashHypothesisAnaAlg.h"

Public Types
using	Providers_t = FlashHypothesisCreator::Providers_t

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

void	SetOutputObjects (TTree , TH1F , TH1F , TH1F , TH1F , TH1F , TH1F , TH1F , TH1F , TH1F , geo::Geometry const &)

void	FillOpDetPositions (geo::Geometry const &)

void	RunComparison (const unsigned int run, const unsigned int event, std::vector< sim::MCTrack > const &, std::vector< sim::SimPhotons > const &, Providers_t providers, opdet::OpDigiProperties const &opdigip, phot::PhotonVisibilityService const &pvs)

Private Attributes
unsigned int	fCounterIndex

float	fdEdx

float	fXOffset

TTree *	fTree

FlashHypothesisCreator	fFHCreator

SimPhotonCounterAlg	fSPCAlg

FlashHypothesisComparison	fFHCompare

sim::MCTrackCollectionAnaAlg	fMCTAlg

std::vector< float >	fOpDetPositions_Y

std::vector< float >	fOpDetPositions_Z

Detailed Description

Definition at line 43 of file FlashHypothesisAnaAlg.h.

Member Typedef Documentation

using opdet::FlashHypothesisAnaAlg::Providers_t = FlashHypothesisCreator::Providers_t

Definition at line 46 of file FlashHypothesisAnaAlg.h.

Constructor & Destructor Documentation

opdet::FlashHypothesisAnaAlg::FlashHypothesisAnaAlg ( fhicl::ParameterSet const & p )

inline

Definition at line 48 of file FlashHypothesisAnaAlg.h.

       : fCounterIndex(p.get<unsigned int>("SimPhotonCounterIndex", 0))
       , fdEdx(p.get<float>("dEdx", 2.1))
       , fXOffset(p.get<float>("HypothesisXOffset", 0.0))
       , fSPCAlg(p.get<fhicl::ParameterSet>("SimPhotonCounterAlgParams"))
     {}

Member Function Documentation

void opdet::FlashHypothesisAnaAlg::FillOpDetPositions ( geo::Geometry const & geom )

Definition at line 45 of file FlashHypothesisAnaAlg.cxx.

References geo::GeometryCore::Cryostat(), fOpDetPositions_Y, fOpDetPositions_Z, geo::OpDetGeo::GetCenter(), geo::GeometryCore::NOpDets(), and geo::CryostatGeo::OpDet().

Referenced by opdet::FlashHypothesisAna::beginJob().

 {
 
   fOpDetPositions_Y.resize(geom.NOpDets());
   fOpDetPositions_Z.resize(geom.NOpDets());
 
   for (size_t i_opdet = 0; i_opdet < geom.NOpDets(); i_opdet++) {
     auto const xyz = geom.Cryostat().OpDet(i_opdet).GetCenter();
     fOpDetPositions_Y[i_opdet] = (float)xyz.Y();
     fOpDetPositions_Z[i_opdet] = (float)xyz.Z();
   }
 }

void opdet::FlashHypothesisAnaAlg::RunComparison	(	const unsigned int	run,
		const unsigned int	event,
		std::vector< sim::MCTrack > const &	mctrackVec,
		std::vector< sim::SimPhotons > const &	simPhotonsVec,
		Providers_t	providers,
		opdet::OpDigiProperties const &	opdigip,
		phot::PhotonVisibilityService const &	pvs
	)

Definition at line 58 of file FlashHypothesisAnaAlg.cxx.

References opdet::SimPhotonCounterAlg::AddSimPhotonsVector(), fCounterIndex, fdEdx, fFHCompare, fFHCreator, sim::MCTrackCollectionAnaAlg::FillTree(), fMCTAlg, fOpDetPositions_Y, fOpDetPositions_Z, fSPCAlg, fTree, fXOffset, lar::ProviderPack< Providers >::get(), opdet::FlashHypothesisCreator::GetFlashHypothesisCollection(), opdet::SimPhotonCounterAlg::GetSimPhotonCounter(), opdet::SimPhotonCounterAlg::InitializeCounters(), and opdet::FlashHypothesisComparison::RunComparison().

Referenced by opdet::FlashHypothesisAna::analyze().

 {
   auto const* geom = providers.get<geo::GeometryCore>();
 
   FlashHypothesisCollection fhc(geom->NOpDets());
   for (auto const& mctrack : mctrackVec) {
     if (mctrack.size() == 0) continue;
     std::vector<float> dEdxVector(mctrack.size() - 1, fdEdx);
     fhc = fhc + fFHCreator.GetFlashHypothesisCollection(
                   mctrack, dEdxVector, providers, pvs, opdigip, fXOffset);
   }
 
   fSPCAlg.InitializeCounters(*geom, opdigip);
   fSPCAlg.AddSimPhotonsVector(simPhotonsVec);
 
   fFHCompare.RunComparison(run,
                            event,
                            fhc,
                            fSPCAlg.GetSimPhotonCounter(fCounterIndex),
                            fOpDetPositions_Y,
                            fOpDetPositions_Z);
 
   fMCTAlg.FillTree(run, event, mctrackVec);
 
   fTree->Fill();
 }

void opdet::FlashHypothesisAnaAlg::SetOutputObjects	(	TTree *	tree,
		TH1F *	h_h_p,
		TH1F *	h_s_p,
		TH1F *	h_c_p,
		TH1F *	h_h_l,
		TH1F *	h_s_l,
		TH1F *	h_c_l,
		TH1F *	h_h_t,
		TH1F *	h_s_t,
		TH1F *	h_c_t,
		geo::Geometry const &	geom
	)

Title: FlashHypothesisAnaAlg Class Author: Wes Ketchum (wketc.nosp@m.hum@.nosp@m.lanl..nosp@m.gov)

Description: Alg that compares the flash hypotheses with truth photons and stores the results in a TTree.

Definition at line 26 of file FlashHypothesisAnaAlg.cxx.

References fFHCompare, fMCTAlg, fTree, geo::GeometryCore::NOpDets(), opdet::FlashHypothesisComparison::SetOutputObjects(), and sim::MCTrackCollectionAnaAlg::SetOutputTree().

Referenced by opdet::FlashHypothesisAna::beginJob().

 {
   fTree = tree;
   fFHCompare.SetOutputObjects(
     tree, h_h_p, h_s_p, h_c_p, h_h_l, h_s_l, h_c_l, h_h_t, h_s_t, h_c_t, geom.NOpDets(), false);
 
   fMCTAlg.SetOutputTree(tree, false);
 }