Inheritance diagram for ems::MultiEMShowers:

Public Types
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	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	beginJob () override

void	endJob () override

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

bool	convCluster (art::Event const &evt)

double	getMinDist (detinfo::DetectorPropertiesData const &detProp, 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 227 of file MultiEMShowers_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 308 of file MultiEMShowers_module.cc.

References fCluModuleLabel, fConvBothGood, fConvGood, fConvWrong, fEv2Good, fEv2Groups, fEvComp, fEvFidVol, fEvGMomCut, fEvInput, fEvReco, fHitsModuleLabel, fShsModuleLabel, fTrk3DModuleLabel, fVtxModuleLabel, and fhicl::ParameterSet::get().

                                                             : EDAnalyzer(p)
 {
   fConvGood = 0;
   fConvWrong = 0;
   fConvBothGood = 0;
   fEvFidVol = 0;
   fEvComp = 0;
   fEvGMomCut = 0;
   fEvReco = 0;
   fEvInput = 0;
   fEv2Groups = 0;
   fEv2Good = 0;
   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");
 }

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

delete

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

delete

Member Function Documentation

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

overrideprivatevirtual

Implements art::EDAnalyzer.

Definition at line 385 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;
   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 ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 327 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, and fStartZ.

 {
   art::ServiceHandle<art::TFileService const> 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");
 }

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

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

private

Definition at line 591 of file MultiEMShowers_module.cc.

References geo::CryostatID::Cryostat, geo::GeometryCore::Cryostat(), larg4::dist(), fCluModuleLabel, fHitsModuleLabel, geo::GeometryCore::FindTPCAtPosition(), art::ProductRetriever::getByLabel(), getMinDist(), ems::MCinfo::GetPosgamma1(), ems::MCinfo::GetPosgamma2(), geo::GeometryCore::HasTPC(), hits(), geo::CryostatGeo::MaxPlanes(), geo::GeometryCore::PositionToCryostatID(), and geo::TPCID::TPC.

Referenced by analyze().

 {
   ems::MCinfo mc;
   TVector3 const convp[2]{mc.GetPosgamma1(), mc.GetPosgamma2()};
 
   geo::Point_t const vtx{convp[0].X(), convp[0].Y(), convp[0].Z()};
 
   art::ServiceHandle<geo::Geometry const> geom;
   geo::TPCID idtpc = geom->FindTPCAtPosition(vtx);
   size_t cryoid = geom->PositionToCryostatID(vtx).Cryostat;
 
   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);
 
   auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt);
   double maxdist = 1.0; // 1 cm
   if (geom->HasTPC(idtpc)) {
     const geo::CryostatGeo& cryostat = geom->Cryostat(geo::CryostatID(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(detProp, 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;
 }

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

void ems::MultiEMShowers::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 371 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";
 }

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

double ems::MultiEMShowers::getMinDist	(	detinfo::DetectorPropertiesData const &	detProp,
		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, larg4::dist(), 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(detProp, v[h]->WireID().Wire, v[h]->PeakTime(), view, tpc, cryo);
 
       double dist = pma::Dist2(proj, hpoint);
       if (dist < mindist) { mindist = dist; }
     }
   }
 
   return mindist;
 }

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

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";
   }

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

delete

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

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 art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

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

art::InputTag ems::MultiEMShowers::fCluModuleLabel

private

Definition at line 302 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fConvBothGood

private

Definition at line 251 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fConvGood

private

Definition at line 249 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fConvWrong

private

Definition at line 250 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fCountph

private

Definition at line 277 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fCountreco

private

Definition at line 278 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDedxU

private

Definition at line 287 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDedxV

private

Definition at line 286 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDedxZ

private

Definition at line 285 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fDistConvrecomc1

private

Definition at line 294 of file MultiEMShowers_module.cc.

Referenced by analyze().

double ems::MultiEMShowers::fDistConvrecomc2

private

Definition at line 295 of file MultiEMShowers_module.cc.

Referenced by analyze().

int ems::MultiEMShowers::fEv2Good

private

Definition at line 276 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fEv2Groups

private

Definition at line 275 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fEvComp

private

Definition at line 266 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fEvFidVol

private

Definition at line 265 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fEvGMomCut

private

Definition at line 267 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fEvInput

private

Definition at line 268 of file MultiEMShowers_module.cc.

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

int ems::MultiEMShowers::fEvNumber

private

Definition at line 256 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fEvReco

private

Definition at line 274 of file MultiEMShowers_module.cc.

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

TTree* ems::MultiEMShowers::fEvTree

private

Definition at line 255 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fGammasInside

private

Definition at line 252 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

TVector3 ems::MultiEMShowers::fGdir1

private

Definition at line 269 of file MultiEMShowers_module.cc.

Referenced by analyze().

TVector3 ems::MultiEMShowers::fGdir2

private

Definition at line 270 of file MultiEMShowers_module.cc.

Referenced by analyze().

double ems::MultiEMShowers::fGdirmcreco1

private

Definition at line 296 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGdirmcreco1good

private

Definition at line 298 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGdirmcreco2

private

Definition at line 297 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 261 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fGmom2

private

Definition at line 262 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 MultiEMShowers().

double ems::MultiEMShowers::fMCrecoTh

private

Definition at line 289 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMCrecoThgood

private

Definition at line 291 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMCrecovtx

private

Definition at line 288 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMCrecovtxgood

private

Definition at line 290 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fMcth

private

Definition at line 263 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fNgammas

private

Definition at line 264 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

int ems::MultiEMShowers::fNGroups

private

Definition at line 257 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fPi0mom

private

Definition at line 260 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

TVector3 ems::MultiEMShowers::fPrimary

private

Definition at line 271 of file MultiEMShowers_module.cc.

Referenced by analyze().

TTree* ems::MultiEMShowers::fRecoTree

private

Definition at line 281 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fRecth

private

Definition at line 292 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fRecthgood

private

Definition at line 293 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 MultiEMShowers().

TTree* ems::MultiEMShowers::fShTree

private

Definition at line 280 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fStartX

private

Definition at line 282 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fStartY

private

Definition at line 283 of file MultiEMShowers_module.cc.

Referenced by analyze(), and beginJob().

double ems::MultiEMShowers::fStartZ

private

Definition at line 284 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 MultiEMShowers().

art::InputTag ems::MultiEMShowers::fVtxModuleLabel

private

Definition at line 304 of file MultiEMShowers_module.cc.

Referenced by analyze(), and MultiEMShowers().

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

larreco/v09_25_00/source/larreco/DirOfGamma/MultiEMShowers_module.cc

Public Types

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