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

TH1F *	fHDTUE

TH1F *	fHDTVE

TH1F *	fHDTWE

TH1F *	fHDTUPull

TH1F *	fHDTVPull

TH1F *	fHDTWPull

TH1F *	fHDTUV

TH1F *	fHDTVW

TH1F *	fHDTWU

TH2F *	fHDTUVU

TH2F *	fHDTUVV

TH2F *	fHDTVWV

TH2F *	fHDTVWW

TH2F *	fHDTWUW

TH2F *	fHDTWUU

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 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 121 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))
     , fBooked(false)
     , fHDTUE(0)
     , fHDTVE(0)
     , fHDTWE(0)
     , fHDTUPull(0)
     , fHDTVPull(0)
     , fHDTWPull(0)
     , fHDTUV(0)
     , fHDTVW(0)
     , fHDTWU(0)
     , fHDTUVU(0)
     , fHDTUVV(0)
     , fHDTVWV(0)
     , fHDTVWW(0)
     , fHDTWUW(0)
     , fHDTWUU(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)
   {
     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 264 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(), trkf::SpacePointAlg::separation(), recob::Hit::SigmaPeakTime(), art::PtrVector< T >::size(), cheat::BackTrackerService::SpacePointHitsToWeightedXYZ(), recob::Hit::SummedADC(), 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.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 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 186 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, 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, geo::GeometryCore::Iterate(), geo::kU, geo::kV, geo::kZ, and trkf::SpacePointAlg::merge().

Referenced by analyze().

   {
     if (!fBooked) {
       fBooked = true;
 
       art::ServiceHandle<geo::Geometry const> geom;
       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.
 
       for (auto const& pgeom : geom->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. * 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.);
       }
     }
   }

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

private

Definition at line 77 of file SpacePointAna_module.cc.

Referenced by bookHistograms().

std::string trkf::SpacePointAna::fClusterModuleLabel

private

Definition at line 66 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

TH1F* trkf::SpacePointAna::fHAmpU

private

Definition at line 98 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAmpV

private

Definition at line 99 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAmpW

private

Definition at line 100 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaU

private

Definition at line 101 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaV

private

Definition at line 102 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHAreaW

private

Definition at line 103 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHchisq

private

Definition at line 94 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUE

private

Definition at line 78 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUPull

private

Definition at line 81 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTUV

private

Definition at line 84 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTUVU

private

Definition at line 87 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTUVV

private

Definition at line 88 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVE

private

Definition at line 79 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVPull

private

Definition at line 82 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTVW

private

Definition at line 85 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTVWV

private

Definition at line 89 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTVWW

private

Definition at line 90 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWE

private

Definition at line 80 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWPull

private

Definition at line 83 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHDTWU

private

Definition at line 86 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTWUU

private

Definition at line 92 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH2F* trkf::SpacePointAna::fHDTWUW

private

Definition at line 91 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

std::string trkf::SpacePointAna::fHitModuleLabel

private

Definition at line 64 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

TH1F* trkf::SpacePointAna::fHMCdx

private

Definition at line 107 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCdy

private

Definition at line 108 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCdz

private

Definition at line 109 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCxpull

private

Definition at line 110 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCypull

private

Definition at line 111 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHMCzpull

private

Definition at line 112 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHS

private

Definition at line 93 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumU

private

Definition at line 104 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumV

private

Definition at line 105 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHSumW

private

Definition at line 106 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHx

private

Definition at line 95 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHy

private

Definition at line 96 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

TH1F* trkf::SpacePointAna::fHz

private

Definition at line 97 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

double trkf::SpacePointAna::fMaxX

private

Definition at line 69 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMaxY

private

Definition at line 71 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMaxZ

private

Definition at line 73 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMinX

private

Definition at line 68 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMinY

private

Definition at line 70 of file SpacePointAna_module.cc.

Referenced by analyze().

double trkf::SpacePointAna::fMinZ

private

Definition at line 72 of file SpacePointAna_module.cc.

Referenced by analyze().

int trkf::SpacePointAna::fNumEvent

private

Definition at line 116 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgDefault

private

Definition at line 63 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgSep

private

Definition at line 62 of file SpacePointAna_module.cc.

Referenced by analyze().

const SpacePointAlg trkf::SpacePointAna::fSptalgTime

private

Definition at line 61 of file SpacePointAna_module.cc.

Referenced by analyze(), and bookHistograms().

bool trkf::SpacePointAna::fUseClusterHits

private

Definition at line 65 of file SpacePointAna_module.cc.

Referenced by analyze(), and SpacePointAna().

bool trkf::SpacePointAna::fUseMC

private

Definition at line 67 of file SpacePointAna_module.cc.

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

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

larreco/v09_25_00/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