ems::MergeEMShower3D Class Reference

Inheritance diagram for ems::MergeEMShower3D:

Public Types
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

Public Member Functions
	MergeEMShower3D (fhicl::ParameterSet const &p)
	MergeEMShower3D (MergeEMShower3D const &)=delete
	MergeEMShower3D (MergeEMShower3D &&)=delete
MergeEMShower3D &	operator= (MergeEMShower3D const &)=delete
MergeEMShower3D &	operator= (MergeEMShower3D &&)=delete
void	beginJob () override
void	produce (art::Event &e) override
void	reconfigure (fhicl::ParameterSet const &p)
void	mcinfo (art::Event &evt)
int	getClusterBestId (const std::vector< art::Ptr< recob::Hit > > &v)
int	getGammaId (art::Event &evt, const size_t t)
double	getCos3D (const TVector3 &p0, const recob::Track &trk)
double	getClusterAdcSum (const std::vector< art::Ptr< recob::Hit > > &v)
std::vector< ems::ShowersCollection >	collectshowers (art::Event &evt, std::vector< ems::ShowerInfo > &showers, bool refpoint)
TVector3	getBestPoint (const std::vector< recob::Track > &tracks, const std::vector< ShowerInfo > &showers, const TVector3 &p0, const TVector3 &p1, double step)
template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const
template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const
bool	modifiesEvent () 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
void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
void	prepareForJob (fhicl::ParameterSet const &pset)
void	showMissingConsumes () const

Private Attributes
TTree *	fEvTree
int	fVtxIndex
int	fNParts [2]
int	fNPMA
int	fNConv
int	fNTot
int	fWhat
int	fEvWhat
int	fNMerged
int	fNCleanMerged
double	fDistvtxmcreco
double	fMcMom
int	fEvNumber
TVector3	fPi0vtx
TVector3	fRefPoint
std::string	fHitsModuleLabel
std::string	fCluModuleLabel
std::string	fTrk3DModuleLabel
std::string	fVtxModuleLabel
double	fNarrowConeAngle
double	fWideConeAngle
shower::ShowerEnergyAlg	fShowerEnergyAlg

Detailed Description

Definition at line 324 of file MergeEMShower3D_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

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

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 389 of file MergeEMShower3D_module.cc.

References reconfigure().

  : fShowerEnergyAlg(p.get<fhicl::ParameterSet>("ShowerEnergyAlg"))
{
        reconfigure(p);

        produces< std::vector<recob::Shower> >();
        produces< std::vector<recob::Vertex> >();

        produces< art::Assns<recob::Shower, recob::Vertex> >(); 
        produces< art::Assns<recob::Shower, recob::Cluster> >();
        produces< art::Assns<recob::Shower, recob::Hit> >();
}

ems::MergeEMShower3D::MergeEMShower3D ( MergeEMShower3D const &  )

delete

ems::MergeEMShower3D::MergeEMShower3D ( MergeEMShower3D &&  )

delete

Member Function Documentation

void ems::MergeEMShower3D::beginJob ( )

overridevirtual

Reimplemented from art::EDProducer.

Definition at line 402 of file MergeEMShower3D_module.cc.

References fDistvtxmcreco, fEvNumber, fEvTree, fEvWhat, fMcMom, fNCleanMerged, fNConv, fNMerged, fNParts, fNPMA, fNTot, and art::TFileDirectory::make().

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

        fEvTree = tfs->make<TTree>("ShowerTestEv", "pi0 tests");
        fEvTree->Branch("fEvNumber", &fEvNumber, "fEvNumber/I");
        fEvTree->Branch("fNPartsA", &fNParts[0], "fNPartsA/I");
        fEvTree->Branch("fNPartsB", &fNParts[1], "fNPartsB/I");
        fEvTree->Branch("fNTot", &fNTot, "fNTot/I");
        fEvTree->Branch("fNPMA", &fNPMA, "fNPMA/I");
        fEvTree->Branch("fNConv", &fNConv, "fNConv/I");
        fEvTree->Branch("fEvWhat", &fEvWhat, "fEvWhat/I");
        fEvTree->Branch("fMcMom", &fMcMom, "fMcMom/D");

        fEvTree->Branch("fDistvtxmcreco", &fDistvtxmcreco, "fDistvtxmcreco/D");
        fEvTree->Branch("fNMerged", &fNMerged, "fNMerged/I");
        fEvTree->Branch("fNCleanMerged", &fNCleanMerged, "fNCleanMerged/I");
}

std::vector< ems::ShowersCollection > ems::MergeEMShower3D::collectshowers	(	art::Event &	evt,
		std::vector< ems::ShowerInfo > &	showers,
		bool	refpoint
	)

Definition at line 616 of file MergeEMShower3D_module.cc.

References geo::GeometryCore::CryostatBoundaries(), pma::Dist2(), fNarrowConeAngle, fNConv, fRefPoint, fTrk3DModuleLabel, getBestPoint(), art::DataViewImpl::getByLabel(), pma::GetProjectionToSegment(), and pma::SolveLeastSquares3D().

Referenced by produce().

{
        std::vector< ems::ShowersCollection > gammawithconv;

        const double smallcone = fNarrowConeAngle; // degree.
        bool merge = true;
        if (refpoint)
        {
                art::ServiceHandle< geo::Geometry > geom;
                art::Handle< std::vector<recob::Track> > trkListHandle;         
                if (evt.getByLabel(fTrk3DModuleLabel, trkListHandle)) 
                {
                        double dsize[6];
                        geom->CryostatBoundaries(dsize, 0);
                        TVector3 geoP0(dsize[0], dsize[2], dsize[4]), geoP1(dsize[1], dsize[3], dsize[5]);
                        TVector3 pov = getBestPoint(*trkListHandle, showers, geoP0, geoP1, 5.0);
                        fRefPoint = pov;

                        for (size_t is = 0; is < showers.size(); is++)
                                showers[is].SetP0Dist(pov);
                        std::sort(showers.begin(), showers.end(), ems::bDistLess());

                        for (size_t is = 0; is < showers.size(); is++)
                        {
                                size_t m_idx = 0;
                                double a, a_min = 360.0;
                                for (size_t m = 0; m < gammawithconv.size(); ++m)
                                {
                                        a = gammawithconv[m].Angle(showers[is].GetFront());
                                        if ((a < fNarrowConeAngle) && (a < a_min))
                                        {
                                                m_idx = m; a_min = a; 
                                        }

                                        std::vector< std::pair<TVector3, TVector3> > lines;
                                        lines.push_back(std::pair<TVector3, TVector3>(
                                                        showers[is].GetFront(), showers[is].GetFront() + showers[is].GetDir()));
                                        lines.push_back(std::pair<TVector3, TVector3>(
                                                        gammawithconv[m].Front(), gammawithconv[m].Front() +  gammawithconv[m].Dir()));

                                        TVector3 isect, pm;
                                        pma::SolveLeastSquares3D(lines, isect);
                                        pm = pma::GetProjectionToSegment(isect,
                                                        gammawithconv[m].Front(), gammawithconv[m].Front() +  gammawithconv[m].Dir());

                                        if (pma::Dist2(pov, gammawithconv[m].Front()) < pma::Dist2(pov, pm))
                                        {
                                                a = gammawithconv[m].Angle(isect);
                                                if ((a < fNarrowConeAngle) && (a < a_min))
                                                {
                                                        m_idx = m; a_min = a;
                                                }
                                        }
                                }
                                if (a_min < fNarrowConeAngle) 
                                {
                                        gammawithconv[m_idx].Merge(showers[is]);
                                }
                                else if (showers[is].HasConPoint())
                                {
                                        ems::ShowersCollection sc(showers[is]);
                                        gammawithconv.push_back(sc);
                                }

                        }
                }
        }
        else
        {
                for (size_t i = 0; i < showers.size(); ++i)
                        if (showers[i].HasConPoint())
                        {
                                ems::ShowersCollection sc(showers[i]);
                                gammawithconv.push_back(sc);
                        }

                fNConv = gammawithconv.size();
                while (merge)
                {
                        merge = false;
                        size_t i = 0;
                        while (i < gammawithconv.size())
                        {
                                size_t best = 0; double a_min = smallcone; 
                                for (size_t j = 0; j < gammawithconv.size(); j++)
                                        if (i != j)
                                        {
                                                double a = gammawithconv[i].MinAngle(gammawithconv[j]);
                                                if (a < a_min)
                                                {
                                                        a_min = a; best = j; merge = true;
                                                }
                                        }
                                if (merge)
                                {
                                        gammawithconv[i].Merge(gammawithconv[best]);
                                        gammawithconv.erase(gammawithconv.begin() + best);
                                        break;
                                }
                                i++;
                        }
                }
        }
        return  gammawithconv;
}

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

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::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

  {
    return current_context_.get();
  }

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

TVector3 ems::MergeEMShower3D::getBestPoint	(	const std::vector< recob::Track > &	tracks,
		const std::vector< ShowerInfo > &	showers,
		const TVector3 &	p0,
		const TVector3 &	p1,
		double	step
	)

Definition at line 806 of file MergeEMShower3D_module.cc.

References pma::Dist2(), f, and getCos3D().

Referenced by collectshowers().

{
        TVector3 best, p;

        double f, fmin = 1.0e10;

        double dx = step; 
        double dy = step; 
        double dz = step; 

        double x0 = p0.X();
        while (x0 < p1.X())
        {
                double y0 = p0.Y();
                while (y0 < p1.Y())
                {
                        double z0 = p0.Z();
                        while (z0 < p1.Z())
                        {
                                f = 0.0;
                                //size_t nOK = 0;
                                TVector3 mid(0., 0., 0.);

                                p.SetXYZ(x0, y0, z0);
                                for (size_t t = 0; t < tracks.size(); t++)
                                //      if (showers[t].OK)
                                {
                                        auto const & trk = tracks[t];

                                        double cos = -acos( getCos3D(p, trk) ) / (0.5*TMath::Pi());
                                        //double cos = getCos3D(p, trk);

                                        if (showers[t].HasConPoint()) cos *= 3.0;
                                        cos *= sqrt( showers[t].GetAdcSum() );

                                        mid += 0.5 * (trk.Vertex() + trk.End());

                                        f += cos;
                                //      nOK++;
                                }
                                //if (!nOK) return best;

                        //      f /= nOK;
                        //      mid *= 1.0 / nOK;

                                f = -f + 0.0001 * sqrt(pma::Dist2(p, mid));
                                if (f < fmin)
                                {
                                        fmin = f; best = p;
                                }

                                z0 += dz;
                        }
                        y0 += dy;
                }
                x0 += dx;
        }
        return best;
}

double ems::MergeEMShower3D::getClusterAdcSum

(

const std::vector< art::Ptr< recob::Hit > > &

v

)

Definition at line 890 of file MergeEMShower3D_module.cc.

References DEFINE_ART_MODULE.

Referenced by produce().

{
        double sum = 0.0;
        for (auto ptr : v)
        {
                sum += ptr->SummedADC();
        }
        return sum;
}

int ems::MergeEMShower3D::getClusterBestId

(

const std::vector< art::Ptr< recob::Hit > > &

v

)

Definition at line 723 of file MergeEMShower3D_module.cc.

References cheat::BackTrackerService::HitToTrackIDEs(), and max.

Referenced by getGammaId().

{
        art::ServiceHandle< cheat::BackTrackerService > bt_serv;

        std::map< int, size_t > ids;
        for (auto ptr : v)
        {
                auto hid = bt_serv->HitToTrackIDEs(ptr);
                if (hid.size()) ids[hid.front().trackID]++;
        }

        int best_id = 9999;
        size_t max = 0;
        for (auto p : ids)
        {
                if ((p.second > (int)(0.8 * v.size())) && (p.second > max))
                {
                        max = p.second; best_id = p.first;
                }
        }

        return best_id;
}

double ems::MergeEMShower3D::getCos3D	(	const TVector3 &	p0,
		const recob::Track &	trk
	)

Definition at line 869 of file MergeEMShower3D_module.cc.

References dir, pma::Dist2(), recob::Track::End(), recob::Track::Vertex(), and recob::Track::VertexDirection().

Referenced by getBestPoint().

{
        TVector3 p1, dir;

        if (pma::Dist2(p0, trk.Vertex()) < pma::Dist2(p0, trk.End()))
                p1 = trk.Vertex();
        else
                p1 = trk.End();

        p1 -= p0;
        double m = p1.Mag();
        if (m > 0.0)
        {
                p1 *= 1.0 / m;
                double c = fabs(p1 * trk.VertexDirection());
                if (c > 1.0) c = 1.0;
                return c;
        }
        else return 0.0;
}

int ems::MergeEMShower3D::getGammaId	(	art::Event &	evt,
		const size_t	t
	)

Definition at line 773 of file MergeEMShower3D_module.cc.

References fCluModuleLabel, fHitsModuleLabel, fTrk3DModuleLabel, fVtxModuleLabel, fWhat, art::DataViewImpl::getByLabel(), and getClusterBestId().

Referenced by produce().

{
        fWhat = 0;
        int cid = 0, gid = 0;
        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 (evt.getByLabel(fTrk3DModuleLabel, trkListHandle) &&
            evt.getByLabel(fVtxModuleLabel, vtxListHandle) &&
            evt.getByLabel(fCluModuleLabel, cluListHandle) &&
            evt.getByLabel(fHitsModuleLabel, hitListHandle))
        {
                art::FindManyP< recob::Cluster > cluFromTrk(trkListHandle, evt, fTrk3DModuleLabel);
                art::FindManyP< recob::Vertex > vtxFromTrk(trkListHandle, evt, fVtxModuleLabel);
                art::FindManyP< recob::Hit > hitFromClu(cluListHandle, evt, fCluModuleLabel);

                auto src_clu_list = cluFromTrk.at(t);
                for (size_t c = 0; c < src_clu_list.size(); ++c)  
                {
                        std::vector< art::Ptr<recob::Hit> > v = hitFromClu.at(src_clu_list[c].key());
                        cid = getClusterBestId(v);

                        if ((fWhat == 0) && (cid == 9999)) fWhat = 1; // confused 2D
                        if ((fWhat == 0) && (gid == 0)) gid = cid;    
                        if ((fWhat == 0) && (cid != gid)) fWhat = 2;  // confused 3D 
                }       
        }

        if (fWhat == 0) return gid;
        else return 9999;
}

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const &  instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

  {
    return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                               instanceName);
  }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID ( ModuleDescription const &  moduleDescription, std::string const &  instanceName  ) const

inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

  {
    auto const& pd =
      get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
    return pd.productID();
  }

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

void ems::MergeEMShower3D::mcinfo

(

art::Event &

evt

)

Definition at line 747 of file MergeEMShower3D_module.cc.

References art::fill_ptr_vector(), fMcMom, fPi0vtx, art::DataViewImpl::getByLabel(), simb::MCParticle::P(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by produce().

{
        int pi0_idx = -1;

        art::Handle< std::vector<simb::MCTruth> > mctruthListHandle;
  std::vector< art::Ptr<simb::MCTruth> > mclist;
  if (evt.getByLabel("generator", mctruthListHandle))
        {
                art::fill_ptr_vector(mclist, mctruthListHandle);
                if (mclist.size())
                {
                        size_t np = mclist[0]->NParticles();
                        for (size_t i = 0; i < np; ++i)
                                if (mclist[0]->GetParticle(i).PdgCode() == 111)
                                        {pi0_idx = i; break;}
                }
        }

        if (pi0_idx < 0) return; 

        const simb::MCParticle & pi0 = mclist[0]->GetParticle(pi0_idx);
        TVector3 pi0_vtx(pi0.Vx(), pi0.Vy(), pi0.Vz());
        fPi0vtx = pi0_vtx;
        fMcMom = pi0.P();
}

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

    {
      return true;
    }

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

staticinherited

Definition at line 76 of file Consumer.cc.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

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

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

delete

MergeEMShower3D& ems::MergeEMShower3D::operator= ( MergeEMShower3D &&  )

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

void ems::MergeEMShower3D::produce ( art::Event &  e )

overridevirtual

Implements art::EDProducer.

Definition at line 434 of file MergeEMShower3D_module.cc.

References collectshowers(), util::CreateAssn(), dir, pma::Dist2(), art::EventID::event(), fCluModuleLabel, fDistvtxmcreco, fEvNumber, fEvTree, fHitsModuleLabel, fMcMom, fNCleanMerged, fNConv, fNMerged, fNParts, fNPMA, fNTot, fPi0vtx, fRefPoint, fShowerEnergyAlg, fTrk3DModuleLabel, fVtxIndex, fVtxModuleLabel, fWhat, fWideConeAngle, art::DataViewImpl::getByLabel(), getClusterAdcSum(), getGammaId(), hits(), art::Event::id(), geo::kU, geo::kV, geo::kZ, mcinfo(), art::Event::put(), recob::Shower::set_total_energy(), recob::Shower::set_total_energy_err(), and shower::ShowerEnergyAlg::ShowerEnergy().

{
        std::unique_ptr< std::vector< recob::Shower > > cascades(new std::vector< recob::Shower >);
        std::unique_ptr< std::vector< recob::Vertex > > vertices(new std::vector< recob::Vertex >);

        std::unique_ptr< art::Assns< recob::Shower, recob::Vertex > > shs2vtx(new art::Assns< recob::Shower, recob::Vertex >);
        std::unique_ptr< art::Assns< recob::Shower, recob::Cluster > > shs2cl(new art::Assns< recob::Shower, recob::Cluster >);
        std::unique_ptr< art::Assns< recob::Shower, recob::Hit > > shs2hit(new art::Assns< recob::Shower, recob::Hit >);

        std::vector< ShowersCollection > gammawithconv;

        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 (evt.getByLabel(fTrk3DModuleLabel, trkListHandle) &&
            evt.getByLabel(fVtxModuleLabel, vtxListHandle) &&
            evt.getByLabel(fCluModuleLabel, cluListHandle) &&
            evt.getByLabel(fHitsModuleLabel, hitListHandle))
        {

                fEvNumber = evt.id().event();
                fNMerged = 0; fNCleanMerged = 0; fNParts[0] = 0; fNParts[1] = 0;
                fNPMA = 0; fNConv = 0; fNTot = 0; fMcMom = 0; fWhat = 0;
                fDistvtxmcreco = -10.0;

                art::FindManyP< recob::Cluster > cluFromTrk(trkListHandle, evt, fTrk3DModuleLabel);
                art::FindManyP< recob::Vertex > vtxFromTrk(trkListHandle, evt, fVtxModuleLabel);
                art::FindManyP< recob::Hit > hitFromClu(cluListHandle, evt, fCluModuleLabel);

                // all reconstructed shower fragments
                std::vector< ShowerInfo > showers;
                for (size_t t = 0; t < trkListHandle->size(); ++t)
                {
                        const recob::Track& trk = (*trkListHandle)[t];

                        auto src_clu_list = cluFromTrk.at(t);
                        double adcsum = 0.0;
                        std::vector<double> ensum(3, 0.0);
                        
                        for (size_t c = 0; c < src_clu_list.size(); ++c)  
                        {
                                std::vector< art::Ptr<recob::Hit> > v = hitFromClu.at(src_clu_list[c].key());
                                adcsum += getClusterAdcSum(v);

                                if (src_clu_list[c]->View() == geo::kU) { ensum[0] += adcsum; }
                                else if (src_clu_list[c]->View() == geo::kV) { ensum[1] += adcsum; }
                                else if (src_clu_list[c]->View() == geo::kZ) { ensum[2] += adcsum; }
                        }

                        auto cnv = vtxFromTrk.at(t);
                        if (cnv.size()) 
                        {
                                int gid = getGammaId(evt, t);
                                ShowerInfo si(t, gid, true, ensum, trk);
                                showers.push_back(si);  
                        }
                        else    
                        {       
                                int gid = getGammaId(evt, t);
                                ShowerInfo si(t, gid, false, ensum, trk);
                                showers.push_back(si);
                        }
                }

                fNTot = showers.size();

                gammawithconv = collectshowers(evt, showers, true); // true switches on refpoint
        
                // with pma segments reconstructed
                // procceed with pma segments when two conversions 
                std::vector< ShowerInfo > pmaseg;
                for (size_t i = 0; i < showers.size(); ++i)
                        if (!showers[i].HasConPoint())
                                pmaseg.push_back(showers[i]);                   
                
                fNPMA = pmaseg.size();

                const double bigcone = fWideConeAngle; // degree.

                if (gammawithconv.size())
                        for (size_t i = 0; i < pmaseg.size(); i++)
                        {
                                double a_min = bigcone; size_t best = 0;
                                for (size_t j = 0; j < gammawithconv.size(); j++)
                                {
                                        TVector3 halfpoint = (pmaseg[i].GetFront() + pmaseg[i].GetEnd()) * 0.5;
                                        double a = gammawithconv[j].first.Angleto(halfpoint);
                                        if (a < a_min)
                                        {
                                                a_min = a; best = j;
                                        }
                                }
                                if (a_min < bigcone)
                                        gammawithconv[best].Merge(pmaseg[i]);           
                        }
                
                mcinfo(evt);
                

                for (size_t i = 0; i < gammawithconv.size(); ++i)
                                if (gammawithconv[i].IsClean()) fNCleanMerged++;
                
                fNMerged = gammawithconv.size(); 
                if (gammawithconv.size() == 2)
                {
                        fNParts[0] = gammawithconv[0].Size();
                        fNParts[1] = gammawithconv[1].Size();
                }

                fDistvtxmcreco = std::sqrt(pma::Dist2(fRefPoint, fPi0vtx));     
                fEvTree->Fill();

                fVtxIndex = 0;
                for (size_t i = 0; i < gammawithconv.size(); ++i)
                {
                        int id = i;
                        TVector3 v0(0., 0., 0.);
                        TVector3 dir = gammawithconv[i].Dir();
                        TVector3 front = gammawithconv[i].Front();
                         
                        std::vector<double> dedx = gammawithconv[i].DeDx();
                        std::vector< double > vd;
                        recob::Shower cas(
                                dir, v0, front, v0,
                                vd, vd, dedx, vd, gammawithconv[i].PlaneId(), id); 


                        std::vector< art::Ptr<recob::Cluster> > cls;
                        std::vector< art::Ptr<recob::Hit> > hits;
                        for (size_t p = 0; p < gammawithconv[i].Size(); p++)
                        {
                                int trkKey = gammawithconv[i].GetParts()[p].GetKey();
                                auto src_clu_list = cluFromTrk.at(trkKey);

                                for (size_t c = 0; c < src_clu_list.size(); c++)
                                {
                                        cls.push_back(src_clu_list[c]);

                                        auto v = hitFromClu.at(src_clu_list[c].key());
                                        for (size_t h = 0; h < v.size(); h++) hits.push_back(v[h]);
                                }

                                auto ver_list = vtxFromTrk.at(trkKey);
                        }

                        std::vector<double> totE;
                        std::vector<double> totEerr;
                        for (int i = 0; i<3; ++i){
                          totE.push_back(fShowerEnergyAlg.ShowerEnergy(hits,i));
                          totEerr.push_back(0);
                        }
                        cas.set_total_energy(totE);
                        cas.set_total_energy_err(totEerr);
                        cascades->push_back(cas);

                        double vtx_pos[3] = {front.X(), front.Y(), front.Z()};
                        vertices->push_back(recob::Vertex(vtx_pos, fVtxIndex));
                        fVtxIndex++;

                        if (vertices->size())
                        {
                                size_t vtx_idx = (size_t)(vertices->size() - 1);
                                util::CreateAssn(*this, evt, *cascades, *vertices, *shs2vtx, vtx_idx, vtx_idx + 1);
                        }

                        util::CreateAssn(*this, evt, *cascades, cls, *shs2cl);
                        util::CreateAssn(*this, evt, *cascades, hits, *shs2hit);

                }

        }

        evt.put(std::move(cascades));
        evt.put(std::move(vertices));

        evt.put(std::move(shs2vtx));
        evt.put(std::move(shs2cl));
        evt.put(std::move(shs2hit));
}

void ems::MergeEMShower3D::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 421 of file MergeEMShower3D_module.cc.

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

Referenced by MergeEMShower3D().

{
        fHitsModuleLabel = p.get< std::string >("HitsModuleLabel");
        fCluModuleLabel = p.get< std::string >("ClustersModuleLabel");
        fTrk3DModuleLabel = p.get< std::string >("Trk3DModuleLabel");
        fVtxModuleLabel = p.get< std::string >("VtxModuleLabel");

        fNarrowConeAngle = p.get< double >("NarrowConeAngle");
        fWideConeAngle = p.get< double >("WideConeAngle");

        return;
}

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

Member Data Documentation

std::string ems::MergeEMShower3D::fCluModuleLabel

private

Definition at line 378 of file MergeEMShower3D_module.cc.

Referenced by getGammaId(), produce(), and reconfigure().

double ems::MergeEMShower3D::fDistvtxmcreco

private

Definition at line 370 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

int ems::MergeEMShower3D::fEvNumber

private

Definition at line 372 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

TTree* ems::MergeEMShower3D::fEvTree

private

Definition at line 359 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

int ems::MergeEMShower3D::fEvWhat

private

Definition at line 367 of file MergeEMShower3D_module.cc.

Referenced by beginJob().

std::string ems::MergeEMShower3D::fHitsModuleLabel

private

Definition at line 377 of file MergeEMShower3D_module.cc.

Referenced by getGammaId(), produce(), and reconfigure().

double ems::MergeEMShower3D::fMcMom

private

Definition at line 371 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), mcinfo(), and produce().

double ems::MergeEMShower3D::fNarrowConeAngle

private

Definition at line 382 of file MergeEMShower3D_module.cc.

Referenced by collectshowers(), and reconfigure().

int ems::MergeEMShower3D::fNCleanMerged

private

Definition at line 368 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

int ems::MergeEMShower3D::fNConv

private

Definition at line 365 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), collectshowers(), and produce().

int ems::MergeEMShower3D::fNMerged

private

Definition at line 368 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

int ems::MergeEMShower3D::fNParts[2]

private

Definition at line 364 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

int ems::MergeEMShower3D::fNPMA

private

Definition at line 365 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

int ems::MergeEMShower3D::fNTot

private

Definition at line 365 of file MergeEMShower3D_module.cc.

Referenced by beginJob(), and produce().

TVector3 ems::MergeEMShower3D::fPi0vtx

private

Definition at line 374 of file MergeEMShower3D_module.cc.

Referenced by mcinfo(), and produce().

TVector3 ems::MergeEMShower3D::fRefPoint

private

Definition at line 375 of file MergeEMShower3D_module.cc.

Referenced by collectshowers(), and produce().

shower::ShowerEnergyAlg ems::MergeEMShower3D::fShowerEnergyAlg

private

Definition at line 385 of file MergeEMShower3D_module.cc.

Referenced by produce().

std::string ems::MergeEMShower3D::fTrk3DModuleLabel

private

Definition at line 379 of file MergeEMShower3D_module.cc.

Referenced by collectshowers(), getGammaId(), produce(), and reconfigure().

int ems::MergeEMShower3D::fVtxIndex

private

Definition at line 362 of file MergeEMShower3D_module.cc.

Referenced by produce().

std::string ems::MergeEMShower3D::fVtxModuleLabel

private

Definition at line 380 of file MergeEMShower3D_module.cc.

Referenced by getGammaId(), produce(), and reconfigure().

int ems::MergeEMShower3D::fWhat

private

Definition at line 367 of file MergeEMShower3D_module.cc.

Referenced by getGammaId(), and produce().

double ems::MergeEMShower3D::fWideConeAngle

private

Definition at line 383 of file MergeEMShower3D_module.cc.

Referenced by produce(), and reconfigure().

---

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

• larreco/v06_64_02/source/larreco/DirOfGamma/MergeEMShower3D_module.cc