trkf::SpacePointAna Class Reference

Inheritance diagram for trkf::SpacePointAna:

Public Types
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer

Public Member Functions
	SpacePointAna (fhicl::ParameterSet const &pset)
virtual	~SpacePointAna ()
void	bookHistograms (bool mc)
void	analyze (const art::Event &evt)
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 Attributes
const SpacePointAlg	fSptalgTime
const SpacePointAlg	fSptalgSep
const SpacePointAlg	fSptalgDefault
std::string	fHitModuleLabel
bool	fUseClusterHits
std::string	fClusterModuleLabel
bool	fUseMC
bool	fBooked
TH1F *	fHDTUE
TH1F *	fHDTVE
TH1F *	fHDTWE
TH1F *	fHDTUPull
TH1F *	fHDTVPull
TH1F *	fHDTWPull
TH1F *	fHDTUV
TH1F *	fHDTVW
TH1F *	fHDTWU
TH1F *	fHS
TH1F *	fHchisq
TH1F *	fHx
TH1F *	fHy
TH1F *	fHz
TH1F *	fHAmpU
TH1F *	fHAmpV
TH1F *	fHAmpW
TH1F *	fHAreaU
TH1F *	fHAreaV
TH1F *	fHAreaW
TH1F *	fHSumU
TH1F *	fHSumV
TH1F *	fHSumW
TH1F *	fHMCdx
TH1F *	fHMCdy
TH1F *	fHMCdz
TH1F *	fHMCxpull
TH1F *	fHMCypull
TH1F *	fHMCzpull
int	fNumEvent

Detailed Description

Definition at line 50 of file SpacePointAna_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

trkf::SpacePointAna::SpacePointAna ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 119 of file SpacePointAna_module.cc.

References fClusterModuleLabel, fHitModuleLabel, fUseClusterHits, and fUseMC.

  : EDAnalyzer(pset)
   , fSptalgTime(pset.get<fhicl::ParameterSet>("SpacePointAlgTime"))
   , fSptalgSep(pset.get<fhicl::ParameterSet>("SpacePointAlgSep"))
   , fSptalgDefault(pset.get<fhicl::ParameterSet>("SpacePointAlgDefault"))
   , fHitModuleLabel(pset.get<std::string>("HitModuleLabel"))
   , fUseClusterHits(pset.get<bool>("UseClusterHits"))
   , fClusterModuleLabel(pset.get<std::string>("ClusterModuleLabel"))
   , fUseMC(pset.get<bool>("UseMC"))
   , fBooked(false)
   , fHDTUE(0)
   , fHDTVE(0)
   , fHDTWE(0)
   , fHDTUPull(0)
   , fHDTVPull(0)
   , fHDTWPull(0)
   , fHDTUV(0)
   , fHDTVW(0)
   , fHDTWU(0)
   , fHS(0)
   , fHchisq(0)
   , fHx(0)
   , fHy(0)
   , fHz(0)
   , fHAmpU(0)
   , fHAmpV(0)
   , fHAmpW(0)
   , fHAreaU(0)
   , fHAreaV(0)
   , fHAreaW(0)
   , fHSumU(0)
   , fHSumV(0)
   , fHSumW(0)
   , fHMCdx(0)
   , fHMCdy(0)
   , fHMCdz(0)
   , fHMCxpull(0)
   , fHMCypull(0)
   , fHMCzpull(0)
   , fNumEvent(0)
  {

    // Report.

    mf::LogInfo("SpacePointAna") 
      << "SpacePointAna configured with the following parameters:\n"
      << "  HitModuleLabel = " << fHitModuleLabel << "\n"
      << "  UseClusterHits = " << fUseClusterHits << "\n"
      << "  ClusterModuleLabel = " << fClusterModuleLabel << "\n"
      << "  UseMC = " << fUseMC;
  }

trkf::SpacePointAna::~SpacePointAna ( )

virtual

Definition at line 176 of file SpacePointAna_module.cc.

  {}

Member Function Documentation

void trkf::SpacePointAna::analyze

(

const art::Event &

evt

)

Definition at line 236 of file SpacePointAna_module.cc.

References art::PtrVector< T >::begin(), bookHistograms(), recob::SpacePoint::Chisq(), detinfo::DetectorProperties::ConvertXToTicks(), trkf::SpacePointAlg::correctedTime(), art::PtrVector< T >::end(), recob::SpacePoint::ErrXYZ(), fClusterModuleLabel, fHAmpU, fHAmpV, fHAmpW, fHAreaU, fHAreaV, fHAreaW, fHchisq, fHDTUE, fHDTUPull, fHDTUV, fHDTVE, fHDTVPull, fHDTVW, fHDTWE, fHDTWPull, fHDTWU, fHitModuleLabel, fHMCdx, fHMCdy, fHMCdz, fHMCxpull, fHMCypull, fHMCzpull, fHS, fHSumU, fHSumV, fHSumW, fHx, fHy, fHz, fNumEvent, fSptalgDefault, fSptalgSep, fSptalgTime, fUseClusterHits, fUseMC, trkf::SpacePointAlg::getAssociatedHits(), art::DataViewImpl::getByLabel(), hits(), cheat::BackTrackerService::HitToXYZ(), recob::Hit::Integral(), art::Event::isRealData(), art::Handle< T >::isValid(), geo::kU, geo::kV, geo::kZ, LOG_DEBUG, trkf::SpacePointAlg::makeMCTruthSpacePoints(), trkf::SpacePointAlg::makeSpacePoints(), trkf::SpacePointAlg::merge(), recob::Hit::PeakAmplitude(), recob::Hit::PeakTime(), geo::PlaneID::Plane, trkf::SpacePointAlg::separation(), recob::Hit::SigmaPeakTime(), cheat::BackTrackerService::SpacePointHitsToWeightedXYZ(), recob::Hit::SummedADC(), t1, t2, geo::TPCID::TPC, lar::dump::vector(), recob::Hit::View(), recob::Hit::WireID(), x, and recob::SpacePoint::XYZ().

  {
    ++fNumEvent;

    // Make sure histograms are booked.

    bool mc = !evt.isRealData();
    bookHistograms(mc);

    // Get Services.

    art::ServiceHandle<geo::Geometry> geom;
    const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();

    // Get Hits.

    art::PtrVector<recob::Hit> hits;

    if(fUseClusterHits) {

      // Get clusters.

      art::Handle< std::vector<recob::Cluster> > clusterh;
      evt.getByLabel(fClusterModuleLabel, clusterh);

      // Get hits from all clusters.
      art::FindManyP<recob::Hit> fm(clusterh, evt, fClusterModuleLabel);

      if(clusterh.isValid()) {
        int nclus = clusterh->size();

        for(int i = 0; i < nclus; ++i) {
          art::Ptr<recob::Cluster> pclus(clusterh, i);
          std::vector< art::Ptr<recob::Hit> > clushits = fm.at(i);
          int nhits = clushits.size();
          hits.reserve(hits.size() + nhits);

          for(std::vector< art::Ptr<recob::Hit> >::const_iterator ihit = clushits.begin(); ihit != clushits.end(); ++ihit) {
            hits.push_back(*ihit);
          }
        }
      }
    }
    else {

      // Get unclustered hits.

      art::Handle< std::vector<recob::Hit> > hith;
      evt.getByLabel(fHitModuleLabel, hith);
      if(hith.isValid()) {
        int nhits = hith->size();
        hits.reserve(nhits);

        for(int i = 0; i < nhits; ++i)
          hits.push_back(art::Ptr<recob::Hit>(hith, i));
      }
    }

    // Fill histograms that don't depend on space points.

    if(mc && fUseMC) {

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

      // Loop over hits and fill hit-electron time difference histogram.

      for(art::PtrVector<recob::Hit>::const_iterator ihit = hits.begin();
          ihit != hits.end(); ++ihit) {
        const recob::Hit& hit = **ihit;

        //unsigned int channel = hit.Channel();
        //geo::View_t geo_view = geom->View(channel);
        //geo::View_t hit_view = hit.View();
        //assert(geo_view == hit_view);
        double tpeak = hit.PeakTime();
        double terr = hit.SigmaPeakTime();

        //assert(channel == hit.Channel());

        // Loop over electrons associated with this hit/channel and fill
        // hit-electron time difference histograms.

        // loop over the map of TDC to sim::IDE to get the TDC for each energy dep
        // Find the average time in ticks for this hit.

        //double sumw = 0.;
        //double sumt = 0.;

        //const std::map<unsigned short, std::vector<sim::IDE> > &idemap = simchan.TDCIDEMap();
        //std::map<unsigned short, std::vector<sim::IDE> >::const_iterator mitr = idemap.begin();
        //for(mitr = idemap.begin(); mitr != idemap.end(); mitr++) {
        //  double tdc = double((*mitr).first);
        //  if(tdc >= tstart && tdc <= tend) {
        //    const std::vector<sim::IDE>& idevec = (*mitr).second;
        //    for(std::vector<sim::IDE>::const_iterator iide=idevec.begin();
        //      iide != idevec.end(); ++iide) {
        //      const sim::IDE& ide = *iide;
        //      double w = ide.numElectrons;
        //      sumw += w;
        //      sumt += w*tdc;
        //    }
        //  }
        //}
        //double tav = 0.;
        //if(sumw != 0.)
        //  tav = sumt / sumw;

        bool tav_ok = true;
        double tav = 0.;
        try {
          std::vector<double> hitxyz = bt_serv->HitToXYZ(*ihit);
          tav = detprop->ConvertXToTicks(hitxyz[0], (*ihit)->WireID().Plane, (*ihit)->WireID().TPC, (*ihit)->WireID().Cryostat);
        }
        catch(cet::exception& x) {
          tav_ok = false;
        }
        if(tav_ok) {
          if((*ihit)->View() == geo::kU) {
            fHDTUE->Fill(tpeak - tav);
            fHDTUPull->Fill((tpeak - tav) / terr);
          }
          else if((*ihit)->View() == geo::kV) {
            fHDTVE->Fill(tpeak - tav);
            fHDTVPull->Fill((tpeak - tav) / terr);
          }
          else if((*ihit)->View() == geo::kZ) {
            fHDTWE->Fill(tpeak - tav);
            fHDTWPull->Fill((tpeak - tav) / terr);
          }
          else
            throw cet::exception("SpacePointAna") << "Bad view = " << (*ihit)->View() << "\n";
        }
      }
    }
    
    std::vector<recob::SpacePoint> spts1;  // For time histograms.
    std::vector<recob::SpacePoint> spts2;  // For separation histogram.
    std::vector<recob::SpacePoint> spts3;  // Default cuts.

    // If nonzero time cut is specified, make space points using that
    // time cut (for time histograms).

    if(!fSptalgTime.merge()) {
      if(mc && fUseMC)
        fSptalgTime.makeMCTruthSpacePoints(hits, spts1);
      else
        fSptalgTime.makeSpacePoints(hits, spts1);

      // Report number of space points.

      LOG_DEBUG("SpacePointAna") << "Found " << spts1.size() 
                                    << " space points using special time cut.";
    }

    // If nonzero separation cut is specified, make space points using that 
    // separation cut (for separation histogram).

    if(!fSptalgSep.merge()) {
      if(mc && fUseMC)
        fSptalgSep.makeMCTruthSpacePoints(hits, spts2);
      else
        fSptalgSep.makeSpacePoints(hits, spts2);

      // Report number of space points.

      LOG_DEBUG("SpacePointAna") << "Found " << spts2.size() 
                                    << " space points using special seperation cut.";
    }

    // Make space points using default cuts.

    if(mc && fUseMC)
      fSptalgDefault.makeMCTruthSpacePoints(hits, spts3);
    else
      fSptalgDefault.makeSpacePoints(hits, spts3);

    // Report number of space points.

    LOG_DEBUG("SpacePointAna") << "Found " << spts3.size() 
                                  << " space points using default cuts.";

    if(!fSptalgTime.merge()) {

      // Loop over space points and fill time histograms.

      for(std::vector<recob::SpacePoint>::const_iterator i = spts1.begin(); 
          i != spts1.end(); ++i) {
        const recob::SpacePoint& spt = *i;

        // Get hits associated with this SpacePoint.

        const art::PtrVector<recob::Hit>& spthits = fSptalgTime.getAssociatedHits(spt);

        // Make a double loop over hits and fill hit time difference histograms.

        for(art::PtrVector<recob::Hit>::const_iterator ihit = spthits.begin();
            ihit != spthits.end(); ++ihit) {
          const recob::Hit& hit1 = **ihit;

          geo::WireID hit1WireID = hit1.WireID();
          unsigned int tpc1, plane1;
          tpc1 = hit1WireID.TPC;
          plane1 = hit1WireID.Plane;
                
          geo::View_t view1 = hit1.View();
          double t1 = fSptalgTime.correctedTime(hit1);

          for(art::PtrVector<recob::Hit>::const_iterator jhit = spthits.begin();
              jhit != spthits.end(); ++jhit) {
            const recob::Hit& hit2 = **jhit;

            geo::WireID hit2WireID = hit2.WireID();
            unsigned int tpc2, plane2;
            tpc2 = hit2WireID.TPC;
            plane2 = hit2WireID.Plane;

            // Require same tpc, different view.

            if(tpc1 == tpc2 && plane1 != plane2) {

              geo::View_t view2 = hit2.View();
              double t2 = fSptalgTime.correctedTime(hit2);

              if(view1 == geo::kU) {
                if(view2 == geo::kV)
                  fHDTUV->Fill(t1-t2);
                if(view2 == geo::kZ)
                  fHDTWU->Fill(t2-t1);
              }
              if(view1 == geo::kV) {
                if(view2 == geo::kZ)
                  fHDTVW->Fill(t1-t2);
                if(view2 == geo::kU)
                  fHDTUV->Fill(t2-t1);
              }
              if(view1 == geo::kZ) {
                if(view2 == geo::kU)
                  fHDTWU->Fill(t1-t2);
                if(view2 == geo::kV)
                  fHDTVW->Fill(t2-t1);
              }
            }
          }
        }
      }

      // Loop over space points and fill seperation histograms.

      for(std::vector<recob::SpacePoint>::const_iterator i = spts2.begin(); 
          i != spts2.end(); ++i) {
        const recob::SpacePoint& spt = *i;

        // Get hits associated with this SpacePoint.

        const art::PtrVector<recob::Hit>& spthits = fSptalgSep.getAssociatedHits(spt);

        // Fill separation histogram.

        double sep = fSptalgSep.separation(spthits);
        fHS->Fill(sep);
      }
    }

    // Loop over default space points and fill histograms.

    for(std::vector<recob::SpacePoint>::const_iterator i = spts3.begin();
        i != spts3.end(); ++i) {
      const recob::SpacePoint& spt = *i;

      fHchisq->Fill(spt.Chisq());
      fHx->Fill(spt.XYZ()[0]);
      fHy->Fill(spt.XYZ()[1]);
      fHz->Fill(spt.XYZ()[2]);

      // Get hits associated with this SpacePoint.

      std::vector<art::Ptr<recob::Hit> > spthits;
      const art::PtrVector<recob::Hit>& av_spthits = fSptalgDefault.getAssociatedHits(spt);
      for( auto const& ptr : av_spthits ){ spthits.push_back(ptr);}

      // Fill single hit histograms.

      for(art::PtrVector<recob::Hit>::const_iterator ihit = spthits.begin();
          ihit != spthits.end(); ++ihit) {
        const recob::Hit& hit = **ihit;

        geo::View_t view = hit.View();

        if(view == geo::kU) {
          fHAmpU->Fill(hit.PeakAmplitude());
          fHAreaU->Fill(hit.Integral());
          fHSumU->Fill(hit.SummedADC());
        }
        if(view == geo::kV) {
          fHAmpV->Fill(hit.PeakAmplitude());
          fHAreaV->Fill(hit.Integral());
          fHSumV->Fill(hit.SummedADC());
        }
        if(view == geo::kZ) {
          fHAmpW->Fill(hit.PeakAmplitude());
          fHAreaW->Fill(hit.Integral());
          fHSumW->Fill(hit.SummedADC());
        }
      }

      if(mc && fUseMC) {

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

        try {
          std::vector<double> mcxyz = bt_serv->SpacePointHitsToWeightedXYZ(spthits);
          fHMCdx->Fill(spt.XYZ()[0] - mcxyz[0]);
          fHMCdy->Fill(spt.XYZ()[1] - mcxyz[1]);
          fHMCdz->Fill(spt.XYZ()[2] - mcxyz[2]);
          if(spt.ErrXYZ()[0] > 0.)
            fHMCxpull->Fill((spt.XYZ()[0] - mcxyz[0]) / std::sqrt(spt.ErrXYZ()[0]));
          if(spt.ErrXYZ()[2] > 0.)
            fHMCypull->Fill((spt.XYZ()[1] - mcxyz[1]) / std::sqrt(spt.ErrXYZ()[2]));
          if(spt.ErrXYZ()[5] > 0.)
            fHMCzpull->Fill((spt.XYZ()[2] - mcxyz[2]) / std::sqrt(spt.ErrXYZ()[5]));
        }
        catch(cet::exception& x) {}
      }
    }
  }

void trkf::SpacePointAna::bookHistograms

(

bool

mc

)

Definition at line 182 of file SpacePointAna_module.cc.

References geo::GeometryCore::DetHalfHeight(), geo::GeometryCore::DetHalfWidth(), geo::GeometryCore::DetLength(), dir, fBooked, fHAmpU, fHAmpV, fHAmpW, fHAreaU, fHAreaV, fHAreaW, fHchisq, fHDTUE, fHDTUPull, fHDTUV, fHDTVE, fHDTVPull, fHDTVW, fHDTWE, fHDTWPull, fHDTWU, fHMCdx, fHMCdy, fHMCdz, fHMCxpull, fHMCypull, fHMCzpull, fHS, fHSumU, fHSumV, fHSumW, fHx, fHy, fHz, fSptalgTime, fUseMC, art::TFileDirectory::make(), trkf::SpacePointAlg::merge(), and art::TFileDirectory::mkdir().

Referenced by analyze().

  {
    if(!fBooked) {
      fBooked = true;

      art::ServiceHandle<geo::Geometry> geom;
      art::ServiceHandle<art::TFileService> tfs;
      art::TFileDirectory dir = tfs->mkdir("sptana", "SpacePointAna histograms");

      if(mc && fUseMC) {
        fHDTUE = dir.make<TH1F>("MCDTUE", "U-Drift Electrons Time Difference", 100, -5., 5.);
        fHDTVE = dir.make<TH1F>("MCDTVE", "V-Drift Electrons Time Difference", 100, -5., 5.);
        fHDTWE = dir.make<TH1F>("MCDTWE", "W-Drift Electrons Time Difference", 100, -5., 5.);
        fHDTUPull = dir.make<TH1F>("MCDTUPull", "U-Drift Electrons Time Pull", 100, -50., 50.);
        fHDTVPull = dir.make<TH1F>("MCDTVPull", "V-Drift Electrons Time Pull", 100, -50., 50.);
        fHDTWPull = dir.make<TH1F>("MCDTWPull", "W-Drift Electrons Time Pull", 100, -50., 50.);
      }
      if(!fSptalgTime.merge()) {
        fHDTUV = dir.make<TH1F>("DTUV", "U-V time difference", 100, -20., 20.);
        fHDTVW = dir.make<TH1F>("DTVW", "V-W time difference", 100, -20., 20.);
        fHDTWU = dir.make<TH1F>("DTWU", "W-U time difference", 100, -20., 20.);
        fHS = dir.make<TH1F>("DS", "Spatial Separation", 100, -2., 2.);
      }
      fHchisq = dir.make<TH1F>("chisq", "Chisquare", 100, 0., 20.);

      fHx = dir.make<TH1F>("xpos", "X Position",
                           100, 0., 2.*geom->DetHalfWidth());
      fHy = dir.make<TH1F>("ypos", "Y Position",
                           100, -geom->DetHalfHeight(), geom->DetHalfHeight());
      fHz = dir.make<TH1F>("zpos", "Z Position",
                           100, 0., geom->DetLength());
      fHAmpU = dir.make<TH1F>("ampU", "U Hit Amplitude", 50, 0., 50.);
      fHAmpV = dir.make<TH1F>("ampV", "V Hit Amplitude", 50, 0., 50.);
      fHAmpW = dir.make<TH1F>("ampW", "W Hit Amplitude", 50, 0., 50.);
      fHAreaU = dir.make<TH1F>("areaU", "U Hit Area", 100, 0., 500.);
      fHAreaV = dir.make<TH1F>("areaV", "V Hit Area", 100, 0., 500.);
      fHAreaW = dir.make<TH1F>("areaW", "W Hit Area", 100, 0., 500.);
      fHSumU = dir.make<TH1F>("sumU", "U Hit Sum ADC", 100, 0., 500.);
      fHSumV = dir.make<TH1F>("sumV", "V Hit Sum ADC", 100, 0., 500.);
      fHSumW = dir.make<TH1F>("sumW", "W Hit Sum ADC", 100, 0., 500.);
      if(mc && fUseMC) {
        fHMCdx = dir.make<TH1F>("MCdx", "X MC Residual", 100, -1., 1.);
        fHMCdy = dir.make<TH1F>("MCdy", "Y MC Residual", 100, -1., 1.);
        fHMCdz = dir.make<TH1F>("MCdz", "Z MC Residual", 100, -1., 1.);
        fHMCxpull = dir.make<TH1F>("MCxpull", "X MC Pull", 100, -50., 50.);
        fHMCypull = dir.make<TH1F>("MCypull", "Y MC Pull", 100, -50., 50.);
        fHMCzpull = dir.make<TH1F>("MCzpull", "Z MC Pull", 100, -50., 50.);
      }
    }
  }

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

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

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

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.

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

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 art::EventObserverBase::registerProducts ( MasterProductRegistry &  , ProductDescriptions &  , ModuleDescription const &    )

inlineinherited

Definition at line 33 of file EventObserverBase.h.

  {}

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

bool trkf::SpacePointAna::fBooked

private

Definition at line 81 of file SpacePointAna_module.cc.

Referenced by bookHistograms().

std::string trkf::SpacePointAna::fClusterModuleLabel

private

Definition at line 76 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

TH1F* trkf::SpacePointAna::fHAmpU

private

Definition at line 96 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAmpV

private

Definition at line 97 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAmpW

private

Definition at line 98 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaU

private

Definition at line 99 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaV

private

Definition at line 100 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaW

private

Definition at line 101 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHchisq

private

Definition at line 92 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUE

private

Definition at line 82 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUPull

private

Definition at line 85 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUV

private

Definition at line 88 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVE

private

Definition at line 83 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVPull

private

Definition at line 86 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVW

private

Definition at line 89 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWE

private

Definition at line 84 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWPull

private

Definition at line 87 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWU

private

Definition at line 90 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

std::string trkf::SpacePointAna::fHitModuleLabel

private

Definition at line 74 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

TH1F* trkf::SpacePointAna::fHMCdx

private

Definition at line 105 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCdy

private

Definition at line 106 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCdz

private

Definition at line 107 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCxpull

private

Definition at line 108 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCypull

private

Definition at line 109 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCzpull

private

Definition at line 110 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHS

private

Definition at line 91 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumU

private

Definition at line 102 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumV

private

Definition at line 103 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumW

private

Definition at line 104 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHx

private

Definition at line 93 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHy

private

Definition at line 94 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHz

private

Definition at line 95 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

int trkf::SpacePointAna::fNumEvent

private

Definition at line 114 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgDefault

private

Definition at line 73 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgSep

private

Definition at line 72 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgTime

private

Definition at line 71 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

bool trkf::SpacePointAna::fUseClusterHits

private

Definition at line 75 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

bool trkf::SpacePointAna::fUseMC

private

Definition at line 77 of file SpacePointAna_module.cc.

Referenced by analyze(), bookHistograms(), and SpacePointAna().

---

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

• larreco/v06_64_02/source/larreco/TrackFinder/SpacePointAna_module.cc