#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 36 of file FlashHypothesisAnaAlg.h.

Member Typedef Documentation

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

Definition at line 40 of file FlashHypothesisAnaAlg.h.

Constructor & Destructor Documentation

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

inline

Definition at line 42 of file FlashHypothesisAnaAlg.h.

References FillOpDetPositions(), RunComparison(), and SetOutputObjects().

                                                    :
     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 30 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(), and FlashHypothesisAnaAlg().

 {
 
   fOpDetPositions_Y.resize(geom.NOpDets());
   fOpDetPositions_Z.resize(geom.NOpDets());
 
   double xyz[3];
   for(size_t i_opdet=0; i_opdet<geom.NOpDets(); i_opdet++){
     geom.Cryostat(0).OpDet(i_opdet).GetCenter(xyz);
     fOpDetPositions_Y[i_opdet] = (float)xyz[1];
     fOpDetPositions_Z[i_opdet] = (float)xyz[2];
   }
 
 }

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 45 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(), and FlashHypothesisAnaAlg().

 {
   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 13 of file FlashHypothesisAnaAlg.cxx.

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

Referenced by opdet::FlashHypothesisAna::beginJob(), and FlashHypothesisAnaAlg().

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