Inheritance diagram for trkf::SpacePointAna:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	SpacePointAna (fhicl::ParameterSet const &pset)

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	analyze (const art::Event &evt) override

void	bookHistograms (bool mc)

Private Attributes
const SpacePointAlg	fSptalgTime

const SpacePointAlg	fSptalgSep

const SpacePointAlg	fSptalgDefault

std::string	fHitModuleLabel

bool	fUseClusterHits

std::string	fClusterModuleLabel

bool	fUseMC

double	fMinX

double	fMaxX

double	fMinY

double	fMaxY

double	fMinZ

double	fMaxZ

bool	fBooked {false}

TH1F *	fHDTUE {nullptr}

TH1F *	fHDTVE {nullptr}

TH1F *	fHDTWE {nullptr}

TH1F *	fHDTUPull {nullptr}

TH1F *	fHDTVPull {nullptr}

TH1F *	fHDTWPull {nullptr}

TH1F *	fHDTUV {nullptr}

TH1F *	fHDTVW {nullptr}

TH1F *	fHDTWU {nullptr}

TH2F *	fHDTUVU {nullptr}

TH2F *	fHDTUVV {nullptr}

TH2F *	fHDTVWV {nullptr}

TH2F *	fHDTVWW {nullptr}

TH2F *	fHDTWUW {nullptr}

TH2F *	fHDTWUU {nullptr}

TH1F *	fHS {nullptr}

TH1F *	fHchisq {nullptr}

TH1F *	fHx {nullptr}

TH1F *	fHy {nullptr}

TH1F *	fHz {nullptr}

TH1F *	fHAmpU {nullptr}

TH1F *	fHAmpV {nullptr}

TH1F *	fHAmpW {nullptr}

TH1F *	fHAreaU {nullptr}

TH1F *	fHAreaV {nullptr}

TH1F *	fHAreaW {nullptr}

TH1F *	fHSumU {nullptr}

TH1F *	fHSumV {nullptr}

TH1F *	fHSumW {nullptr}

TH1F *	fHMCdx {nullptr}

TH1F *	fHMCdy {nullptr}

TH1F *	fHMCdz {nullptr}

TH1F *	fHMCxpull {nullptr}

TH1F *	fHMCypull {nullptr}

TH1F *	fHMCzpull {nullptr}

int	fNumEvent {0}

Detailed Description

Definition at line 51 of file SpacePointAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 122 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"))
     , fMinX(pset.get<double>("MinX", -1.e10))
     , fMaxX(pset.get<double>("MaxX", 1.e10))
     , fMinY(pset.get<double>("MinY", -1.e10))
     , fMaxY(pset.get<double>("MaxY", 1.e10))
     , fMinZ(pset.get<double>("MinZ", -1.e10))
     , fMaxZ(pset.get<double>("MaxZ", 1.e10))
   {
     mf::LogInfo("SpacePointAna") << "SpacePointAna configured with the following parameters:\n"
                                  << "  HitModuleLabel = " << fHitModuleLabel << "\n"
                                  << "  UseClusterHits = " << fUseClusterHits << "\n"
                                  << "  ClusterModuleLabel = " << fClusterModuleLabel << "\n"
                                  << "  UseMC = " << fUseMC;
   }

Member Function Documentation

void trkf::SpacePointAna::analyze ( const art::Event & evt )

overrideprivate

Definition at line 224 of file SpacePointAna_module.cc.

References bookHistograms(), trkf::SpacePointAlg::correctedTime(), geo::CryostatID::Cryostat, fClusterModuleLabel, fHAmpU, fHAmpV, fHAmpW, fHAreaU, fHAreaV, fHAreaW, fHchisq, fHDTUE, fHDTUPull, fHDTUV, fHDTUVU, fHDTUVV, fHDTVE, fHDTVPull, fHDTVW, fHDTVWV, fHDTVWW, fHDTWE, fHDTWPull, fHDTWU, fHDTWUU, fHDTWUW, fHitModuleLabel, fHMCdx, fHMCdy, fHMCdz, fHMCxpull, fHMCypull, fHMCzpull, fHS, fHSumU, fHSumV, fHSumW, fHx, fHy, fHz, fMaxX, fMaxY, fMaxZ, fMinX, fMinY, fMinZ, fNumEvent, fSptalgDefault, fSptalgSep, fSptalgTime, fUseClusterHits, fUseMC, trkf::SpacePointAlg::getAssociatedHits(), art::ProductRetriever::getByLabel(), hits(), cheat::BackTrackerService::HitToXYZ(), recob::Hit::Integral(), art::Event::isRealData(), art::Handle< T >::isValid(), geo::kU, geo::kV, geo::kZ, trkf::SpacePointAlg::makeMCTruthSpacePoints(), trkf::SpacePointAlg::makeSpacePoints(), trkf::SpacePointAlg::merge(), MF_LOG_DEBUG, recob::Hit::PeakAmplitude(), recob::Hit::PeakTime(), geo::PlaneID::Plane, art::PtrVector< T >::push_back(), art::PtrVector< T >::reserve(), recob::Hit::ROISummedADC(), trkf::SpacePointAlg::separation(), recob::Hit::SigmaPeakTime(), art::PtrVector< T >::size(), cheat::BackTrackerService::SpacePointHitsToWeightedXYZ(), t1, t2, lar::to_element, geo::TPCID::TPC, lar::dump::vector(), recob::Hit::View(), geo::WireID::Wire, recob::Hit::WireID(), geo::WireID::WireID(), and x.

   {
     ++fNumEvent;
 
     // Make sure histograms are booked.
 
     bool mc = !evt.isRealData();
     bookHistograms(mc);
 
     // Get Services.
 
     art::ServiceHandle<geo::Geometry const> geom;
 
     // 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.
 
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
     auto const detProp =
       art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
 
     if (mc && fUseMC) {
 
       art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
 
       // Loop over hits and fill hit-electron time difference histogram.
 
       for (auto const& hitPtr : hits) {
         const recob::Hit& hit = *hitPtr;
 
         double tpeak = hit.PeakTime();
         double terr = hit.SigmaPeakTime();
 
         bool tav_ok = true;
         double tav = 0.;
         try {
           std::vector<double> hitxyz = bt_serv->HitToXYZ(clockData, hit);
           tav = detProp.ConvertXToTicks(
             hitxyz[0], hit.WireID().Plane, hit.WireID().TPC, hit.WireID().Cryostat);
         }
         catch (cet::exception& x) {
           tav_ok = false;
         }
         if (tav_ok) {
           if (hit.View() == geo::kU) {
             fHDTUE->Fill(tpeak - tav);
             fHDTUPull->Fill((tpeak - tav) / terr);
           }
           else if (hit.View() == geo::kV) {
             fHDTVE->Fill(tpeak - tav);
             fHDTVPull->Fill((tpeak - tav) / terr);
           }
           else if (hit.View() == geo::kZ) {
             fHDTWE->Fill(tpeak - tav);
             fHDTWPull->Fill((tpeak - tav) / terr);
           }
           else
             throw cet::exception("SpacePointAna") << "Bad view = " << hit.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(clockData, detProp, hits, spts1);
       else
         fSptalgTime.makeSpacePoints(clockData, detProp, hits, spts1);
 
       // Report number of space points.
 
       MF_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(clockData, detProp, hits, spts2);
       else
         fSptalgSep.makeSpacePoints(clockData, detProp, hits, spts2);
 
       // Report number of space points.
 
       MF_LOG_DEBUG("SpacePointAna")
         << "Found " << spts2.size() << " space points using special seperation cut.";
     }
 
     // Make space points using default cuts.
 
     if (mc && fUseMC)
       fSptalgDefault.makeMCTruthSpacePoints(clockData, detProp, hits, spts3);
     else
       fSptalgDefault.makeSpacePoints(clockData, detProp, hits, spts3);
 
     // Report number of space points.
 
     MF_LOG_DEBUG("SpacePointAna") << "Found " << spts3.size()
                                   << " space points using default cuts.";
 
     if (!fSptalgTime.merge()) {
 
       // Loop over space points and fill time histograms.
 
       for (auto const& spt : spts1) {
         if (spt.XYZ()[0] < fMinX || spt.XYZ()[0] > fMaxX || spt.XYZ()[1] < fMinY ||
             spt.XYZ()[1] > fMaxY || spt.XYZ()[2] < fMinZ || spt.XYZ()[2] > fMaxZ)
           continue;
 
         // 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 (auto const& hit1 : spthits | transform(to_element)) {
           geo::WireID hit1WireID = hit1.WireID();
           unsigned int tpc1 = hit1WireID.TPC;
           unsigned int plane1 = hit1WireID.Plane;
           unsigned int wire1 = hit1WireID.Wire;
 
           geo::View_t view1 = hit1.View();
           double t1 = fSptalgTime.correctedTime(detProp, hit1);
 
           for (auto const& hit2 : spthits | transform(to_element)) {
             geo::WireID hit2WireID = hit2.WireID();
             unsigned int tpc2 = hit2WireID.TPC;
             unsigned int plane2 = hit2WireID.Plane;
             unsigned int wire2 = hit2WireID.Wire;
 
             // Require same tpc, different view.
 
             if (tpc1 == tpc2 && plane1 != plane2) {
 
               geo::View_t view2 = hit2.View();
               double t2 = fSptalgTime.correctedTime(detProp, hit2);
 
               if (view1 == geo::kU) {
                 if (view2 == geo::kV) {
                   fHDTUV->Fill(t1 - t2);
                   fHDTUVU->Fill(double(wire1), t1 - t2);
                   fHDTUVV->Fill(double(wire2), t1 - t2);
                 }
                 if (view2 == geo::kZ) {
                   fHDTWU->Fill(t2 - t1);
                   fHDTWUW->Fill(double(wire2), t2 - t1);
                   fHDTWUU->Fill(double(wire1), t2 - t1);
                 }
               }
               if (view1 == geo::kV) {
                 if (view2 == geo::kZ) {
                   fHDTVW->Fill(t1 - t2);
                   fHDTVWV->Fill(double(wire1), t1 - t2);
                   fHDTVWW->Fill(double(wire2), t1 - t2);
                 }
                 if (view2 == geo::kU) {
                   fHDTUV->Fill(t2 - t1);
                   fHDTUVU->Fill(double(wire2), t2 - t1);
                   fHDTUVV->Fill(double(wire1), t2 - t1);
                 }
               }
               if (view1 == geo::kZ) {
                 if (view2 == geo::kU) {
                   fHDTWU->Fill(t1 - t2);
                   fHDTWUW->Fill(double(wire1), t1 - t2);
                   fHDTWUU->Fill(double(wire2), t1 - t2);
                 }
                 if (view2 == geo::kV) {
                   fHDTVW->Fill(t2 - t1);
                   fHDTVWV->Fill(double(wire2), t2 - t1);
                   fHDTVWW->Fill(double(wire1), t2 - t1);
                 }
               }
             }
           }
         }
       }
 
       // Loop over space points and fill seperation histograms.
 
       for (auto const& spt : spts2) {
         if (spt.XYZ()[0] < fMinX || spt.XYZ()[0] > fMaxX || spt.XYZ()[1] < fMinY ||
             spt.XYZ()[1] > fMaxY || spt.XYZ()[2] < fMinZ || spt.XYZ()[2] > fMaxZ)
           continue;
 
         // 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 (auto const& spt : spts3) {
       if (spt.XYZ()[0] < fMinX || spt.XYZ()[0] > fMaxX || spt.XYZ()[1] < fMinY ||
           spt.XYZ()[1] > fMaxY || spt.XYZ()[2] < fMinZ || spt.XYZ()[2] > fMaxZ)
         continue;
 
       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 (auto const& hitPtr : spthits) {
         const recob::Hit& hit = *hitPtr;
 
         geo::View_t view = hit.View();
 
         if (view == geo::kU) {
           fHAmpU->Fill(hit.PeakAmplitude());
           fHAreaU->Fill(hit.Integral());
           fHSumU->Fill(hit.ROISummedADC());
         }
         if (view == geo::kV) {
           fHAmpV->Fill(hit.PeakAmplitude());
           fHAreaV->Fill(hit.Integral());
           fHSumV->Fill(hit.ROISummedADC());
         }
         if (view == geo::kZ) {
           fHAmpW->Fill(hit.PeakAmplitude());
           fHAreaW->Fill(hit.Integral());
           fHSumW->Fill(hit.ROISummedADC());
         }
       }
 
       if (mc && fUseMC) {
 
         art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
 
         try {
           std::vector<double> mcxyz = bt_serv->SpacePointHitsToWeightedXYZ(clockData, 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&) {
         }
       }
     }
   }

void trkf::SpacePointAna::bookHistograms ( bool mc )

private

Definition at line 145 of file SpacePointAna_module.cc.

References dir, fBooked, fHAmpU, fHAmpV, fHAmpW, fHAreaU, fHAreaV, fHAreaW, fHchisq, fHDTUE, fHDTUPull, fHDTUV, fHDTUVU, fHDTUVV, fHDTVE, fHDTVPull, fHDTVW, fHDTVWV, fHDTVWW, fHDTWE, fHDTWPull, fHDTWU, fHDTWUU, fHDTWUW, fHMCdx, fHMCdy, fHMCdz, fHMCxpull, fHMCypull, fHMCzpull, fHS, fHSumU, fHSumV, fHSumW, fHx, fHy, fHz, fSptalgTime, fUseMC, Get, geo::kU, geo::kV, geo::kZ, trkf::SpacePointAlg::merge(), and geo::GeometryCore::TPC().

Referenced by analyze().

   {
     if (!fBooked) {
       fBooked = true;
 
       art::ServiceHandle<geo::Geometry const> geom;
       auto const& tpc = geom->TPC({0, 0});
       art::ServiceHandle<art::TFileService> tfs;
       art::TFileDirectory dir = tfs->mkdir("sptana", "SpacePointAna histograms");
 
       unsigned int nwiresU = 0, nwiresV = 0, nwiresW = 0;
 
       // Figure out the number of wires in U, V, and W planes.
 
       // Loop over cryostats, tpcs, and planes.
 
       auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout>()->Get();
       for (auto const& pgeom : wireReadoutGeom.Iterate<geo::PlaneGeo>()) {
         unsigned int const nwires = pgeom.Nwires();
         switch (pgeom.View()) {
         case geo::kU: nwiresU = nwires; break;
         case geo::kV: nwiresV = nwires; break;
         case geo::kZ: nwiresW = nwires; break;
         default: {
         }
         }
       }
 
       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.);
         fHDTUVU = dir.make<TH2F>(
           "DTUVU", "U-V time difference vs. U", 100, 0., double(nwiresU), 100, -20., 20.);
         fHDTUVV = dir.make<TH2F>(
           "DTUVV", "U-V time difference vs. V", 100, 0., double(nwiresV), 100, -20., 20.);
         fHDTVWV = dir.make<TH2F>(
           "DTVWV", "V-W time difference vs. V", 100, 0., double(nwiresV), 100, -20., 20.);
         fHDTVWW = dir.make<TH2F>(
           "DTVWW", "V-W time difference vs. W", 100, 0., double(nwiresW), 100, -20., 20.);
         fHDTWUW = dir.make<TH2F>(
           "DTWUW", "W-U time difference vs. W", 100, 0., double(nwiresW), 100, -20., 20.);
         fHDTWUU = dir.make<TH2F>(
           "DTWUU", "W-U time difference vs. U", 100, 0., double(nwiresU), 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. * tpc.HalfWidth());
       fHy = dir.make<TH1F>("ypos", "Y Position", 100, -tpc.HalfHeight(), tpc.HalfHeight());
       fHz = dir.make<TH1F>("zpos", "Z Position", 100, 0., tpc.Length());
       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.);
       }
     }
   }

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

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

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

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

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

bool trkf::SpacePointAna::fBooked {false}

private

Definition at line 78 of file SpacePointAna_module.cc.

Referenced by bookHistograms().

std::string trkf::SpacePointAna::fClusterModuleLabel

private

Definition at line 67 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

TH1F* trkf::SpacePointAna::fHAmpU {nullptr}

private

Definition at line 99 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAmpV {nullptr}

private

Definition at line 100 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAmpW {nullptr}

private

Definition at line 101 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaU {nullptr}

private

Definition at line 102 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaV {nullptr}

private

Definition at line 103 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaW {nullptr}

private

Definition at line 104 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHchisq {nullptr}

private

Definition at line 95 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUE {nullptr}

private

Definition at line 79 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUPull {nullptr}

private

Definition at line 82 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUV {nullptr}

private

Definition at line 85 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTUVU {nullptr}

private

Definition at line 88 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTUVV {nullptr}

private

Definition at line 89 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVE {nullptr}

private

Definition at line 80 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVPull {nullptr}

private

Definition at line 83 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVW {nullptr}

private

Definition at line 86 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTVWV {nullptr}

private

Definition at line 90 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTVWW {nullptr}

private

Definition at line 91 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWE {nullptr}

private

Definition at line 81 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWPull {nullptr}

private

Definition at line 84 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWU {nullptr}

private

Definition at line 87 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTWUU {nullptr}

private

Definition at line 93 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTWUW {nullptr}

private

Definition at line 92 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

std::string trkf::SpacePointAna::fHitModuleLabel

private

Definition at line 65 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

TH1F* trkf::SpacePointAna::fHMCdx {nullptr}

private

Definition at line 108 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCdy {nullptr}

private

Definition at line 109 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCdz {nullptr}

private

Definition at line 110 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCxpull {nullptr}

private

Definition at line 111 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCypull {nullptr}

private

Definition at line 112 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCzpull {nullptr}

private

Definition at line 113 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHS {nullptr}

private

Definition at line 94 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumU {nullptr}

private

Definition at line 105 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumV {nullptr}

private

Definition at line 106 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumW {nullptr}

private

Definition at line 107 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHx {nullptr}

private

Definition at line 96 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHy {nullptr}

private

Definition at line 97 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHz {nullptr}

private

Definition at line 98 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

double trkf::SpacePointAna::fMaxX

private

Definition at line 70 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMaxY

private

Definition at line 72 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMaxZ

private

Definition at line 74 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMinX

private

Definition at line 69 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMinY

private

Definition at line 71 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMinZ

private

Definition at line 73 of file SpacePointAna_module.cc.

Referenced by analyze().

int trkf::SpacePointAna::fNumEvent {0}

private

Definition at line 117 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgDefault

private

Definition at line 64 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgSep

private

Definition at line 63 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgTime

private

Definition at line 62 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

bool trkf::SpacePointAna::fUseClusterHits

private

Definition at line 66 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

bool trkf::SpacePointAna::fUseMC

private

Definition at line 68 of file SpacePointAna_module.cc.

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

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

larreco/v10_01_07/source/larreco/TrackFinder/SpacePointAna_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