Inheritance diagram for ems::MultiEMShowers:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

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

	MultiEMShowers (MultiEMShowers const &)=delete

	MultiEMShowers (MultiEMShowers &&)=delete

MultiEMShowers &	operator= (MultiEMShowers const &)=delete

MultiEMShowers &	operator= (MultiEMShowers &&)=delete

void	beginJob () override

void	endJob () override

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

void	reconfigure (fhicl::ParameterSet const &p)

std::string	workerType () const

bool	modifiesEvent () const

void	registerProducts (MasterProductRegistry &, ProductDescriptions &, ModuleDescription const &)

std::string const &	processName () const

bool	wantAllEvents () const

bool	wantEvent (Event const &e)

fhicl::ParameterSetID	selectorConfig () const

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

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)

base_engine_t &	createEngine (seed_t seed)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)

seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const

detail::CachedProducts &	cachedProducts ()

void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Member Functions
bool	convCluster (art::Event const &evt)

double	getMinDist (std::vector< art::Ptr< recob::Hit > > const &v, TVector3 const &convmc, size_t view, size_t tpc, size_t cryo)

Private Attributes
int	fConvGood

int	fConvWrong

int	fConvBothGood

int	fGammasInside

TTree *	fEvTree

int	fEvNumber

int	fNGroups

double	fPi0mom

double	fGmom1

double	fGmom2

double	fMcth

int	fNgammas

int	fEvFidVol

int	fEvComp

int	fEvGMomCut

int	fEvInput

TVector3	fGdir1

TVector3	fGdir2

TVector3	fPrimary

int	fEvReco

int	fEv2Groups

int	fEv2Good

int	fCountph

int	fCountreco

TTree *	fShTree

TTree *	fRecoTree

double	fStartX

double	fStartY

double	fStartZ

double	fDedxZ

double	fDedxV

double	fDedxU

double	fMCrecovtx

double	fMCrecoTh

double	fMCrecovtxgood

double	fMCrecoThgood

double	fRecth

double	fRecthgood

double	fDistConvrecomc1

double	fDistConvrecomc2

double	fGdirmcreco1

double	fGdirmcreco2

double	fGdirmcreco1good

double	fGdirmcreco2good

art::InputTag	fHitsModuleLabel

art::InputTag	fCluModuleLabel

art::InputTag	fTrk3DModuleLabel

art::InputTag	fVtxModuleLabel

art::InputTag	fShsModuleLabel

Detailed Description

Definition at line 238 of file MultiEMShowers_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 39 of file EDAnalyzer.h.

using art::EDAnalyzer::WorkerType = WorkerT<EDAnalyzer>

inherited

Definition at line 38 of file EDAnalyzer.h.

Constructor & Destructor Documentation

ems::MultiEMShowers::MultiEMShowers ( fhicl::ParameterSet const & p )

explicit

Definition at line 309 of file MultiEMShowers_module.cc.

References fConvBothGood, fConvGood, fConvWrong, fEv2Good, fEv2Groups, fEvComp, fEvFidVol, fEvGMomCut, fEvInput, fEvReco, and reconfigure().

   :
   EDAnalyzer(p)  
 {
         fConvGood = 0;
         fConvWrong = 0;
         fConvBothGood = 0;
         fEvFidVol = 0;
         fEvComp = 0;
         fEvGMomCut = 0;
         fEvReco = 0;
         fEvInput = 0;
         fEv2Groups = 0;
         fEv2Good = 0;
         reconfigure(p);
 }

ems::MultiEMShowers::MultiEMShowers ( MultiEMShowers const & )

delete

ems::MultiEMShowers::MultiEMShowers ( MultiEMShowers && )

delete

Member Function Documentation

void ems::MultiEMShowers::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 395 of file MultiEMShowers_module.cc.

 {
         fEvNumber = e.id().event();
         fNGroups = 0; 
         fStartX = 0.0; fStartY = 0.0; fStartZ = 0.0;
         fPi0mom = 0.0; fNgammas = 0;
         fDistConvrecomc1 = 0.0; fDistConvrecomc2 = 0.0;
         fMCrecovtx = -400.0;    
         fMCrecovtxgood = -400.0;
         fRecth = -400.0;
         fRecthgood = -400.0;
         fMCrecoTh = -400.0;
         fMCrecoThgood = -400.0;
         fGammasInside = 0;
         fCountph = 0;
         fCountreco = 0;
         fGdirmcreco1 = 0.0;
         fGdirmcreco2 = 0.0;
         fGdirmcreco1good = 0.0;
         fGdirmcreco2good = 0.0;
         fDedxZ = 0.0; fDedxV = 0.0; fDedxU = 0.0;
 
         ems::MCinfo mc(e);
         fPrimary = mc.GetPrimary();
         fPi0mom = mc.GetMompi0();
         fGmom1 = mc.GetMomGamma1();
         fGmom2 = mc.GetMomGamma2();
         fGdir1 = mc.GetDirgamma1();
         fGdir2 = mc.GetDirgamma2();
         fNgammas = mc.GetNgammas();
         TVector3 pospi0 = mc.GetPospi0();
 
         double cosinemc = mc.GetCosine();
         fMcth = 180.0F * (std::acos(cosinemc)) / TMath::Pi();
         TVector3 convp[2]; 
         convp[0] = mc.GetPosgamma1();
         convp[1] = mc.GetPosgamma2();
         const double maxdist = 2.0; //cm
 
         // check whether two photons are inside fid vol
         if (mc.IsInside1() && mc.IsInside2()) 
         {
                 fGammasInside = 1;
                 fEvFidVol++;
         }
 
         if ((fGmom1 > 0.1) && (fGmom2 > 0.1)) fEvGMomCut++;
 
         if (mc.IsCompton()) fEvComp++;
 
         art::Handle< std::vector<recob::Shower> > shsListHandle;
         art::Handle< std::vector<recob::Track> > trkListHandle;
         art::Handle< std::vector<recob::Vertex> > vtxListHandle;
         art::Handle< std::vector<recob::Cluster> > cluListHandle;
         art::Handle< std::vector<recob::Hit> > hitListHandle;
         if (e.getByLabel(fShsModuleLabel, shsListHandle) &&
                 e.getByLabel(fTrk3DModuleLabel, trkListHandle) &&
             e.getByLabel(fVtxModuleLabel, vtxListHandle) &&
             e.getByLabel(fCluModuleLabel, cluListHandle) &&
             e.getByLabel(fHitsModuleLabel, hitListHandle))
         {
                 art::FindManyP< recob::Cluster > cluFromShs(shsListHandle, e, fShsModuleLabel);
                 art::FindManyP< recob::Cluster > cluFromTrk(trkListHandle, e, fTrk3DModuleLabel);
                 art::FindManyP< recob::Vertex > vtxFromTrk(trkListHandle, e, fVtxModuleLabel);
                 art::FindManyP< recob::Hit > hitFromClu(cluListHandle, e, fCluModuleLabel);
 
                 fNGroups = shsListHandle->size();       
 
                 fCountph = 0;
                 if (fNgammas == 2) // pi0
                 {
                         int idph = -1; 
                         for (size_t s = 0; s < shsListHandle->size(); ++s)
                         {
                                 const recob::Shower& sh = (*shsListHandle)[s];
                                 double mindist = maxdist; bool found = false; 
                         
                                 for (int i = 0; i < fNgammas; ++i)
                                 {
                                         double dist = sqrt(pma::Dist2(sh.ShowerStart(), convp[i]));
                                         if ((dist < mindist) && (idph != i)) 
                                         { mindist =  dist; idph = i; found = true; }
                                 }
                                 if (found) { fConvGood++; fCountph++; }
                                 else { fConvWrong++; }
                         }
                         if (fCountph == 2) fConvBothGood++;
                 
                         // plot a few variables if there are 2 showers
                         if (fCountph == 2)
                                 for (size_t s = 0; s < shsListHandle->size(); ++s)      
                                 {
                                         const recob::Shower& sh = (*shsListHandle)[s];
                                         TVector3 pos = sh.ShowerStart(); 
                                         fStartX = pos.X(); fStartY = pos.Y(); fStartZ = pos.Z();
                                         std::vector<double> const& vecdedx = sh.dEdx();
                                         
                                         if (vecdedx.size() == 3)
                                         {
                                                 fDedxZ = vecdedx[0]; fDedxV = vecdedx[1]; fDedxU = vecdedx[2];
                                         }
                                         
                                         fShTree->Fill();        
                                 }
                 }
                 else  // other than pi0
                 {
                         for (size_t s = 0; s < shsListHandle->size(); ++s)
                         {
                                 const recob::Shower& sh = (*shsListHandle)[s];
                                 double mindist = maxdist; 
                         
                                 double dist = sqrt(pma::Dist2(sh.ShowerStart(), fPrimary));
                                 if (dist < mindist) 
                                 { 
                                         TVector3 pos = sh.ShowerStart(); 
                                         fStartX = pos.X(); fStartY = pos.Y(); fStartZ = pos.Z();
                                         std::vector<double> vecdedx = sh.dEdx();
                                         if (vecdedx.size() == 3)
                                         {
                                                 fDedxZ = vecdedx[0]; fDedxV = vecdedx[1]; fDedxU = vecdedx[2];
                                         }
                                 }
                                         
                                 fShTree->Fill();        
                         }
                 }
                 // compute the crossing point
 
                 //cut from mc and clusters
         
                 if (mc.IsInside1() && mc.IsInside2() && (fGmom1 > 0.1) && (fGmom2 > 0.1) && (!mc.IsCompton()) && convCluster(e))
                 {
                         fCountreco = 1;
                         if (fNGroups == 2) fEv2Groups++;        
                         if ((fNGroups == 2) && (fCountph == 2)) fEv2Good++;
                         // cut from reco
                         //if (countph == 2)
                         if (fNGroups == 2)
                         {
                                 std::vector< std::pair<TVector3, TVector3> > lines;
                                 const recob::Shower& sh1 = (*shsListHandle)[0];
                                 const recob::Shower& sh2 = (*shsListHandle)[1];
 
                                 std::pair<TVector3, TVector3> frontback1(sh1.ShowerStart(), sh1.ShowerStart() + sh1.Direction());
                                 std::pair<TVector3, TVector3> frontback2(sh2.ShowerStart(), sh2.ShowerStart() + sh2.Direction());
                                 lines.push_back(frontback1); lines.push_back(frontback2);
 
                                 TVector3 result;
                                 pma::SolveLeastSquares3D(lines, result); // mse.
 
                                 double dist1_0 = pma::Dist2(result, sh1.ShowerStart());
                                 double dist2_0 = pma::Dist2(result, sh1.ShowerStart() + sh1.Direction());
                                 double dist1_1 = pma::Dist2(result, sh2.ShowerStart());
                                 double dist2_1 = pma::Dist2(result, sh2.ShowerStart() + sh2.Direction());
                                 if ((dist1_0 > dist2_0) || (dist1_1 > dist2_1)) {}
                                 else
                                 {
                                         fMCrecovtx = std::sqrt(pma::Dist2(pospi0, result));
 
                                         if (fCountph == 2) fMCrecovtxgood = fMCrecovtx;
 
                                         double cosine_reco = sh1.Direction() * sh2.Direction();
                                         fRecth = 180.0F * (std::acos(cosine_reco)) / TMath::Pi();
 
                                         fGdirmcreco1 = fGdir1 * sh1.Direction();
                                         fGdirmcreco2 = fGdir2 * sh2.Direction();
                                         if (fCountph == 2)
                                         {
                                                 fGdirmcreco1good = fGdirmcreco1;
                                                 fGdirmcreco2good = fGdirmcreco2;
                                         }
 
                                         if (fCountph == 2) fRecthgood = fRecth;
 
                                         fMCrecoTh = fRecth - fMcth;
 
                                         if (fCountph == 2) fMCrecoThgood = fMCrecoTh;                                   
 
                                         fEvReco++;
                                         fRecoTree->Fill();
                                 }
                         }
                         fEvInput++;
                         //fRecoTree->Fill();
                 }
         }
 
         fEvTree->Fill();                
 }

void ems::MultiEMShowers::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 337 of file MultiEMShowers_module.cc.

References fConvBothGood, fConvGood, fConvWrong, fCountph, fCountreco, fDedxU, fDedxV, fDedxZ, fEvNumber, fEvTree, fGammasInside, fGdirmcreco1, fGdirmcreco1good, fGdirmcreco2, fGdirmcreco2good, fGmom1, fGmom2, fMCrecoTh, fMCrecoThgood, fMCrecovtx, fMCrecovtxgood, fMcth, fNgammas, fNGroups, fPi0mom, fRecoTree, fRecth, fRecthgood, fShTree, fStartX, fStartY, fStartZ, and art::TFileDirectory::make().

 {
         art::ServiceHandle<art::TFileService> tfs;
 
         fEvTree = tfs->make<TTree>("MultiShowers", "showers3d");
         fEvTree->Branch("fEvNumber", &fEvNumber, "fEvNumber/I");
         fEvTree->Branch("fNGroups", &fNGroups, "fNGroups/I");
         fEvTree->Branch("fPi0mom", &fPi0mom, "fPi0mom/D");
         fEvTree->Branch("fNgammas", &fNgammas, "fNgammas/I");
         fEvTree->Branch("fGmom1", &fGmom1, "fGmom1/D");
         fEvTree->Branch("fGmom2", &fGmom2, "fGmom2/D");
         fEvTree->Branch("fMcth", &fMcth, "fMcth/D");
         fEvTree->Branch("fRecth", &fRecth, "fRecth/D");
         fEvTree->Branch("fMCrecovtx", &fMCrecovtx, "fMCrecovtx/D");
         fEvTree->Branch("fMCrecoTh", &fMCrecoTh, "fMCrecoTh/D");
         fEvTree->Branch("fConvGood", &fConvGood, "fConvGood/I");
         fEvTree->Branch("fConvWrong", &fConvWrong, "fConvWrong/I");
         fEvTree->Branch("fConvBothGood", &fConvBothGood, "fConvBothGood/I");
         fEvTree->Branch("fGammasInside", &fGammasInside, "fGammasInside/I");
         fEvTree->Branch("fCountph", &fCountph, "fCountph/I");
         fEvTree->Branch("fCountreco", &fCountreco, "fCountreco/I");
 
         fRecoTree = tfs->make<TTree>("Cascades", "conv points");
         fRecoTree->Branch("fRecth", &fRecth, "fRecth/D");
         fRecoTree->Branch("fMCrecovtx", &fMCrecovtx, "fMCrecovtx/D");
         fRecoTree->Branch("fMCrecoTh", &fMCrecoTh, "fMCrecoTh/D");
         fRecoTree->Branch("fRecthgood", &fRecthgood, "fRecthgood/D");
         fRecoTree->Branch("fMCrecovtxgood", &fMCrecovtxgood, "fMCrecovtxgood/D");
         fRecoTree->Branch("fMCrecoThgood", &fMCrecoThgood, "fMCrecoThgood/D");
         fRecoTree->Branch("fGdirmcreco1", &fGdirmcreco1, "fGdirmcreco1/D");
         fRecoTree->Branch("fGdirmcreco2", &fGdirmcreco2, "fGdirmcreco2/D");
         fRecoTree->Branch("fGdirmcreco1good", &fGdirmcreco1good, "fGdirmcreco1good/D");
         fRecoTree->Branch("fGdirmcreco2good", &fGdirmcreco2good, "fGdirmcreco2good/D");
 
         fShTree = tfs->make<TTree>("Shower", "conv point");
         
         fShTree->Branch("fStartX", &fStartX, "fStartX/D");
         fShTree->Branch("fStartY", &fStartY, "fStartY/D");
         fShTree->Branch("fStartZ", &fStartZ, "fStartZ/D");
         fShTree->Branch("fDedxZ", &fDedxZ, "fDedxZ/D");
         fShTree->Branch("fDedxV", &fDedxV, "fDedxV/D");
         fShTree->Branch("fDedxU", &fDedxU, "fDedxU/D");
 }

detail::CachedProducts& art::EventObserverBase::cachedProducts ( )

inlineprotectedinherited

Definition at line 79 of file EventObserverBase.h.

References art::EventObserverBase::selectors_.

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

   {
     return selectors_;
   }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const & it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const & it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool ems::MultiEMShowers::convCluster ( art::Event const & evt )

private

Definition at line 587 of file MultiEMShowers_module.cc.

References geo::GeometryCore::Cryostat(), fCluModuleLabel, fHitsModuleLabel, geo::GeometryCore::FindCryostatAtPosition(), geo::GeometryCore::FindTPCAtPosition(), art::DataViewImpl::getByLabel(), getMinDist(), ems::MCinfo::GetPosgamma1(), ems::MCinfo::GetPosgamma2(), geo::GeometryCore::HasTPC(), hits(), geo::CryostatGeo::MaxPlanes(), and geo::TPCID::TPC.

Referenced by analyze().

 {
         ems::MCinfo mc(evt);
         TVector3 convp[2]; 
         convp[0] = mc.GetPosgamma1();
         convp[1] = mc.GetPosgamma2();
 
         double vtx[3] = {convp[0].X(), convp[0].Y(), convp[0].Z()};
 
         art::ServiceHandle<geo::Geometry> geom;
         geo::TPCID idtpc = geom->FindTPCAtPosition(vtx);
         size_t cryoid = geom->FindCryostatAtPosition(vtx);
 
         art::Handle< std::vector<recob::Hit> > hitListHandle;
         art::Handle< std::vector<recob::Cluster> > cluListHandle;
 
         //map: conversion point, vec of id clusters in each view
         std::map < size_t, std::vector< size_t > > used;
 
 
         art::FindManyP< recob::Hit > fbc(cluListHandle, evt, fCluModuleLabel);
 
         double maxdist = 1.0; // 1 cm
         if (geom->HasTPC(idtpc))
         {
                 const geo::CryostatGeo& cryostat = geom->Cryostat(cryoid);
                 if (evt.getByLabel(fHitsModuleLabel, hitListHandle) &&
                         evt.getByLabel(fCluModuleLabel, cluListHandle))
                 {
                         size_t conv = 0;
                         while (conv < 2)
                         {
                                 for (size_t view = 0; view < cryostat.MaxPlanes(); view++)
                                 {
 
                                         double mindist = maxdist; int clid = 0;
                                         for (size_t c = 0; c < cluListHandle->size(); ++c)
                                         {
                                         
                                                 bool exist = false;
                                                 for (auto const & ids : used)
                                                         for (auto i : ids.second)
                                                                 if (i == c) exist = true;
                                                 if (exist) continue;
 
                                                 std::vector< art::Ptr<recob::Hit> > hits = fbc.at(c);
                                                 if (hits.size() < 20) continue;
                                                 if (hits[0]->WireID().Plane != view) continue;
                                 
                                                 double dist = getMinDist(hits, convp[conv], view, idtpc.TPC, cryoid);
                                                 if (dist < mindist)
                                                 {
                                                         mindist = dist;
                                                         clid = c;
                                                 }                                       
                                         }
                                         if (mindist < maxdist) used[conv].push_back(clid);                      
                                 }
                                 conv++;
                         }
                 }
         }
         bool result = false;
         
         if (used.size() > 1)     
                 for (auto const & ids : used)
                 {
                         if (ids.second.size() > 1) result = true;
                         else {result = false; break;}
                 }
         
         return result;
 }

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make
	)

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make,
		label_t const &	engine_label
	)

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

CurrentProcessingContext const * art::EDAnalyzer::currentContext ( ) const

protectedinherited

Definition at line 114 of file EDAnalyzer.cc.

References art::EDAnalyzer::current_context_.

Referenced by art::EDAnalyzer::workerType().

   {
     return current_context_.get();
   }

void ems::MultiEMShowers::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 381 of file MultiEMShowers_module.cc.

References fEv2Good, fEv2Groups, fEvComp, fEvFidVol, fEvGMomCut, fEvInput, and fEvReco.

 {
   mf::LogInfo log("MultiEMShower");
   log << "******************** fEvFidVol =  " << fEvFidVol << "\n";
   log << "******************** fEvGMomCut = " << fEvGMomCut << "\n";
   log << "******************** fEvComp =    " << fEvComp << "\n";
   log << "******************** fEvReco =    " << fEvReco << "\n";
   log << "******************** fEvInput =   " << fEvInput << "\n";
   log << "******************** fEv2Groups = " << fEv2Groups << "\n";
   log << "******************** fEv2Good =   " << fEv2Good << "\n";
   if (fEvInput)
     log << "******************** reco %  =    " << double(fEvReco)/double(fEvInput) << "\n"; 
 }

EngineCreator::seed_t EngineCreator::get_seed_value	(	fhicl::ParameterSet const &	pset,
		char const	key[] = `"seed"`,
		seed_t const	implicit_seed = `-1`
	)

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

double ems::MultiEMShowers::getMinDist	(	std::vector< art::Ptr< recob::Hit > > const &	v,
		TVector3 const &	convmc,
		size_t	view,
		size_t	tpc,
		size_t	cryo
	)

private

Definition at line 661 of file MultiEMShowers_module.cc.

References DEFINE_ART_MODULE, pma::Dist2(), pma::GetProjectionToPlane(), proj, and pma::WireDriftToCm().

Referenced by convCluster().

 {
         double mindist = 9999;
         // MC vertex projected to view
         TVector2 proj = pma::GetProjectionToPlane(convmc, view, tpc, cryo);
         
         // loop over hits to find the closest to MC 2d vtx
         for (size_t h = 0; h < v.size(); ++h)
         {
                 if ((v[h]->WireID().Plane == view) &&
                         (v[h]->WireID().TPC == tpc))
                 {
                         TVector2 hpoint = pma::WireDriftToCm(v[h]->WireID().Wire, v[h]->PeakTime(), view, tpc, cryo);
                         
                         double dist = pma::Dist2(proj, hpoint);
                         if (dist < mindist) 
                         {
                                 mindist = dist;
                         }
                 }
         }
         
         return mindist;
 }

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

inlineinherited

Definition at line 61 of file EventObserverBase.h.

References art::detail::CachedProducts::getOneTriggerResults(), and art::EventObserverBase::selectors_.

Referenced by art::OutputModule::doWriteEvent().

   {
     return selectors_.getOneTriggerResults(e);
   }

template<typename T , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const & it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

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

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const & it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool art::EventObserverBase::modifiesEvent ( ) const

inlineinherited

Definition at line 25 of file EventObserverBase.h.

   {
     return false;
   }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

MultiEMShowers& ems::MultiEMShowers::operator= ( MultiEMShowers const & )

delete

MultiEMShowers& ems::MultiEMShowers::operator= ( MultiEMShowers && )

delete

void art::Consumer::prepareForJob ( fhicl::ParameterSet const & pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

std::string const& art::EventObserverBase::processName ( ) const

inlineinherited

Definition at line 41 of file EventObserverBase.h.

References art::EventObserverBase::process_name_.

Referenced by art::FileDumperOutput::printPrincipal(), and art::RootOutput::RootOutput().

   {
     return process_name_;
   }

void ems::MultiEMShowers::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 326 of file MultiEMShowers_module.cc.

References fCluModuleLabel, fHitsModuleLabel, fShsModuleLabel, fTrk3DModuleLabel, fVtxModuleLabel, and fhicl::ParameterSet::get().

Referenced by MultiEMShowers().

 {
   fHitsModuleLabel = p.get< art::InputTag >("HitsModuleLabel");
   fCluModuleLabel = p.get< art::InputTag >("ClustersModuleLabel");
   fTrk3DModuleLabel = p.get< art::InputTag >("Trk3DModuleLabel");
   fVtxModuleLabel = p.get< art::InputTag >("VtxModuleLabel");
   fShsModuleLabel = p.get< art::InputTag >("ShsModuleLabel");
   
   return;
 }

void art::EventObserverBase::registerProducts	(	MasterProductRegistry &	,
		ProductDescriptions &	,
		ModuleDescription const &
	)

inlineinherited

Definition at line 33 of file EventObserverBase.h.

36 {}

fhicl::ParameterSetID art::EventObserverBase::selectorConfig ( ) const

inlineinherited

Definition at line 56 of file EventObserverBase.h.

References art::EventObserverBase::selector_config_id_.

Referenced by art::RootOutputFile::writeOne().

   {
     return selector_config_id_;
   }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

void art::Consumer::validateConsumedProduct	(	BranchType const	bt,
		ProductInfo const &	pi
	)

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }

bool art::EventObserverBase::wantAllEvents ( ) const

inlineinherited

Definition at line 46 of file EventObserverBase.h.

References art::EventObserverBase::wantAllEvents_.

Referenced by art::EDAnalyzer::doEvent(), art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), art::RootOutput::RootOutput(), and art::OutputWorker::wantAllEvents().

   {
     return wantAllEvents_;
   }

bool art::EventObserverBase::wantEvent ( Event const & e )

inlineinherited

Definition at line 51 of file EventObserverBase.h.

References art::EventObserverBase::selectors_, and art::detail::CachedProducts::wantEvent().

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

   {
     return selectors_.wantEvent(e);
   }

std::string art::EDAnalyzer::workerType ( ) const

inlineinherited

Definition at line 109 of file EDAnalyzer.h.

References art::EDAnalyzer::currentContext().

     {
       return "WorkerT<EDAnalyzer>";
     }

Member Data Documentation

art::InputTag ems::MultiEMShowers::fCluModuleLabel

private

Definition at line 302 of file MultiEMShowers_module.cc.

Referenced by analyze(), convCluster(), and reconfigure().

int ems::MultiEMShowers::fConvBothGood

private

Definition at line 261 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fConvGood

private

Definition at line 259 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fConvWrong

private

Definition at line 260 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fCountph

private

Definition at line 287 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fCountreco

private

Definition at line 288 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDedxU

private

Definition at line 293 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDedxV

private

Definition at line 293 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDedxZ

private

Definition at line 293 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDistConvrecomc1

private

Definition at line 297 of file MultiEMShowers_module.cc.

Referenced by analyze().

double ems::MultiEMShowers::fDistConvrecomc2

private

Definition at line 297 of file MultiEMShowers_module.cc.

Referenced by analyze().

int ems::MultiEMShowers::fEv2Good

private

Definition at line 286 of file MultiEMShowers_module.cc.

Referenced by analyze(), endJob(), and MultiEMShowers().

int ems::MultiEMShowers::fEv2Groups

private

Definition at line 285 of file MultiEMShowers_module.cc.

Referenced by analyze(), endJob(), and MultiEMShowers().

int ems::MultiEMShowers::fEvComp

private

Definition at line 276 of file MultiEMShowers_module.cc.

Referenced by analyze(), endJob(), and MultiEMShowers().

int ems::MultiEMShowers::fEvFidVol

private

Definition at line 275 of file MultiEMShowers_module.cc.

Referenced by analyze(), endJob(), and MultiEMShowers().

int ems::MultiEMShowers::fEvGMomCut

private

Definition at line 277 of file MultiEMShowers_module.cc.

Referenced by analyze(), endJob(), and MultiEMShowers().

int ems::MultiEMShowers::fEvInput

private

Definition at line 278 of file MultiEMShowers_module.cc.

Referenced by analyze(), endJob(), and MultiEMShowers().

int ems::MultiEMShowers::fEvNumber

private

Definition at line 266 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fEvReco

private

Definition at line 284 of file MultiEMShowers_module.cc.

Referenced by analyze(), endJob(), and MultiEMShowers().

TTree* ems::MultiEMShowers::fEvTree

private

Definition at line 265 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fGammasInside

private

Definition at line 262 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

TVector3 ems::MultiEMShowers::fGdir1

private

Definition at line 279 of file MultiEMShowers_module.cc.

Referenced by analyze().

TVector3 ems::MultiEMShowers::fGdir2

private

Definition at line 280 of file MultiEMShowers_module.cc.

Referenced by analyze().

double ems::MultiEMShowers::fGdirmcreco1

private

Definition at line 298 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGdirmcreco1good

private

Definition at line 299 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGdirmcreco2

private

Definition at line 298 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGdirmcreco2good

private

Definition at line 299 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGmom1

private

Definition at line 271 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGmom2

private

Definition at line 272 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

art::InputTag ems::MultiEMShowers::fHitsModuleLabel

private

Definition at line 301 of file MultiEMShowers_module.cc.

Referenced by analyze(), convCluster(), and reconfigure().

double ems::MultiEMShowers::fMCrecoTh

private

Definition at line 294 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMCrecoThgood

private

Definition at line 295 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMCrecovtx

private

Definition at line 294 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMCrecovtxgood

private

Definition at line 295 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMcth

private

Definition at line 273 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fNgammas

private

Definition at line 274 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fNGroups

private

Definition at line 267 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fPi0mom

private

Definition at line 270 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

TVector3 ems::MultiEMShowers::fPrimary

private

Definition at line 281 of file MultiEMShowers_module.cc.

Referenced by analyze().

TTree* ems::MultiEMShowers::fRecoTree

private

Definition at line 291 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fRecth

private

Definition at line 296 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fRecthgood

private

Definition at line 296 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

art::InputTag ems::MultiEMShowers::fShsModuleLabel

private

Definition at line 305 of file MultiEMShowers_module.cc.

Referenced by analyze(), and reconfigure().

TTree* ems::MultiEMShowers::fShTree

private

Definition at line 290 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fStartX

private

Definition at line 292 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fStartY

private

Definition at line 292 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fStartZ

private

Definition at line 292 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

art::InputTag ems::MultiEMShowers::fTrk3DModuleLabel

private

Definition at line 303 of file MultiEMShowers_module.cc.

Referenced by analyze(), and reconfigure().

art::InputTag ems::MultiEMShowers::fVtxModuleLabel

private

Definition at line 304 of file MultiEMShowers_module.cc.

Referenced by analyze(), and reconfigure().

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

larreco/v06_64_02/source/larreco/DirOfGamma/MultiEMShowers_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation