Inheritance diagram for shwf::ShowerReco:

Public Types
using	ModuleType = EDProducer

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Member Functions
	ShowerReco (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)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

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
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Private Member Functions
void	beginJob ()

void	beginRun (art::Run &run)

void	produce (art::Event &evt)

void	GetVertexAndAnglesFromCluster (art::Ptr< recob::Cluster > clust, unsigned int plane)

void	LongTransEnergy (geo::GeometryCore const *geom, geo::WireReadoutGeom const &wireReadoutGeom, detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, unsigned int set, std::vector< art::Ptr< recob::Hit >> hitlist)

void	ClearandResizeVectors (unsigned int nPlanes)

Private Attributes
int	fRun

int	fEvent

int	fSubRun

float	slope [3]

float	angle [3]

std::string	fClusterModuleLabel

float	ftimetick

double	fMean_wire_pitch

fhicl::ParameterSet	fCaloPSet

std::vector< double >	fRMS_2cm

std::vector< int >	fNpoints_2cm

std::vector< double >	fCorr_MeV_2cm

std::vector< double >	fCorr_Charge_2cm

std::vector< int >	fNpoints_corr_ADC_2cm

std::vector< int >	fNpoints_corr_MeV_2cm

std::vector< double >	fTotChargeADC

std::vector< double >	fTotChargeMeV

std::vector< double >	fTotChargeMeV_MIPs

std::vector< double >	fChargeADC_2cm

std::vector< double >	fChargeMeV_2cm

std::vector< double >	fChargeMeV_2cm_refined

std::vector< double >	fChargeMeV_2cm_axsum

std::vector< std::vector< double > >	fDistribChargeADC

std::vector< std::vector< double > >	fDistribChargeMeV

std::vector< std::vector< double > >	fDistribHalfChargeMeV

std::vector< std::vector< double > >	fDistribChargeposition

std::vector< std::vector< double > >	fSingleEvtAngle

std::vector< std::vector< double > >	fSingleEvtAngleVal

std::vector< unsigned int >	fWire_vertex

std::vector< double >	fTime_vertex

std::vector< double >	fWire_vertexError

std::vector< double >	fTime_vertexError

std::vector< unsigned int >	fWire_last

std::vector< double >	fTime_last

std::vector< double >	xyz_vertex_fit

std::vector< std::vector< double > >	fNPitch

float	Kin_En

std::vector< float >	vdEdx

std::vector< float >	vresRange

std::vector< float >	vdQdx

std::vector< float >	deadwire

float	Trk_Length

float	fTrkPitchC

float	fdEdxlength

float	fcalodEdxlength

bool	fUseArea

double	xphi

double	xtheta

unsigned int	fNPlanes

unsigned int	fNAngles

TTree *	ftree_shwf

double	fWirePitch

double	fTimeTick

double	fDriftVelocity

double	fWireTimetoCmCm

std::vector< int >	fNhitsperplane

std::vector< double >	fTotADCperplane

Detailed Description

Definition at line 61 of file ShowerReco_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

shwf::ShowerReco::ShowerReco ( fhicl::ParameterSet const & pset )

explicit

Definition at line 166 of file ShowerReco_module.cc.

References fcalodEdxlength, fCaloPSet, fClusterModuleLabel, fdEdxlength, fUseArea, and fhicl::ParameterSet::get().

                                                       : EDProducer{pset}
   {
     fClusterModuleLabel = pset.get<std::string>("ClusterModuleLabel");
     fCaloPSet = pset.get<fhicl::ParameterSet>("CaloAlg");
 
     fdEdxlength =
       pset.get<double>("dEdxlength"); // distance that gets used to determine e/gamma separation
     fcalodEdxlength =
       pset.get<double>("calodEdxlength"); // cutoff distance for hits saved to the calo object.
     fUseArea = pset.get<bool>("UseArea");
 
     produces<std::vector<recob::Shower>>();
     produces<art::Assns<recob::Shower, recob::Cluster>>();
     produces<art::Assns<recob::Shower, recob::Hit>>();
     produces<std::vector<anab::Calorimetry>>();
     produces<art::Assns<recob::Shower, anab::Calorimetry>>();
   }

Member Function Documentation

void shwf::ShowerReco::beginJob ( )

privatevirtual

Todo:: the call to geo->Nplanes() assumes this is a single cryostat and single TPC detector; need to generalize to multiple cryostats and TPCs

Get TFileService and define output Histograms

All-knowing tree with reconstruction information

Reimplemented from art::EDProducer.

Definition at line 185 of file ShowerReco_module.cc.

References fChargeADC_2cm, fChargeMeV_2cm, fChargeMeV_2cm_axsum, fChargeMeV_2cm_refined, fDistribChargeADC, fDistribChargeMeV, fDistribChargeposition, fDistribHalfChargeMeV, fEvent, fMean_wire_pitch, fNAngles, fNhitsperplane, fNPitch, fNPlanes, fNpoints_2cm, fRMS_2cm, fRun, fSubRun, ftimetick, fTotADCperplane, fTotChargeADC, fTotChargeMeV, fTotChargeMeV_MIPs, ftree_shwf, Get, detinfo::sampling_rate(), xphi, xtheta, and xyz_vertex_fit.

   {
     auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
     fNPlanes = wireReadoutGeom.Nplanes();
     fMean_wire_pitch = wireReadoutGeom.Plane({0, 0, 0}).WirePitch(); // wire pitch in cm
 
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
     ftimetick = sampling_rate(clockData) / 1000.;
 
     art::ServiceHandle<art::TFileService const> tfs;
 
     ftree_shwf = tfs->make<TTree>("ShowerReco",
                                   "Results"); 
     ftree_shwf->Branch("run", &fRun, "run/I");
     ftree_shwf->Branch("subrun", &fSubRun, "subrun/I");
     ftree_shwf->Branch("event", &fEvent, "event/I");
     ftree_shwf->Branch("nplanes", &fNPlanes, "nplanes/I");
     ftree_shwf->Branch("nangles", &fNAngles, "nangles/I");
     ftree_shwf->Branch("xtheta", &xtheta, "xtheta/D");
     ftree_shwf->Branch("xphi", &xphi, "xphi/D");
     ftree_shwf->Branch("ftotChargeADC", "std::vector<double>", &fTotChargeADC);
     ftree_shwf->Branch("ftotChargeMeV", "std::vector<double>", &fTotChargeMeV);
     ftree_shwf->Branch("fTotChargeMeV_MIPs", "std::vector<double>", &fTotChargeMeV_MIPs);
     ftree_shwf->Branch("NPitch", "std::vector< std::vector<double> >", &fNPitch);
 
     // this should be temporary - until the omega is sorted out.
     ftree_shwf->Branch("RMS_2cm", "std::vector<double>", &fRMS_2cm);
     ftree_shwf->Branch("Npoints_2cm", "std::vector<int>", &fNpoints_2cm);
     ftree_shwf->Branch("ChargeADC_2cm", "std::vector<double>", &fChargeADC_2cm);
     ftree_shwf->Branch("ChargeMeV_2cm", "std::vector<double>", &fChargeMeV_2cm);
     ftree_shwf->Branch("ChargeMeV_2cm_refined", "std::vector<double>", &fChargeMeV_2cm_refined);
     ftree_shwf->Branch("ChargeMeV_2cm_axsum", "std::vector<double>", &fChargeMeV_2cm_axsum);
     ftree_shwf->Branch("fNhitsperplane", "std::vector<int>", &fNhitsperplane);
     ftree_shwf->Branch("fTotADCperplane", "std::vector<double>", &fTotADCperplane);
     ftree_shwf->Branch(
       "ChargedistributionADC", "std::vector<std::vector<double>>", &fDistribChargeADC);
     ftree_shwf->Branch(
       "ChargedistributionMeV", "std::vector<std::vector<double>>", &fDistribChargeMeV);
     ftree_shwf->Branch(
       "DistribHalfChargeMeV", "std::vector<std::vector<double>>", &fDistribHalfChargeMeV);
     ftree_shwf->Branch(
       "ChargedistributionPosition", "std::vector<std::vector<double>>", &fDistribChargeposition);
     ftree_shwf->Branch("xyz_vertex_fit", "std::vector<double>", &xyz_vertex_fit);
   }

void shwf::ShowerReco::beginRun ( art::Run & run )

privatevirtual

Reimplemented from art::EDProducer.

Definition at line 233 of file ShowerReco_module.cc.

References fChargeADC_2cm, fChargeMeV_2cm, fChargeMeV_2cm_axsum, fChargeMeV_2cm_refined, fCorr_Charge_2cm, fCorr_MeV_2cm, fDistribChargeADC, fDistribChargeMeV, fDistribChargeposition, fDistribHalfChargeMeV, fDriftVelocity, fNAngles, fNhitsperplane, fNPitch, fNPlanes, fNpoints_2cm, fNpoints_corr_ADC_2cm, fNpoints_corr_MeV_2cm, fRMS_2cm, fTime_last, fTime_vertex, fTime_vertexError, fTimeTick, fTotADCperplane, fTotChargeADC, fTotChargeMeV, fTotChargeMeV_MIPs, fWire_last, fWire_vertex, fWire_vertexError, fWirePitch, fWireTimetoCmCm, Get, detinfo::sampling_rate(), vdEdx, vdQdx, vresRange, and xyz_vertex_fit.

   {
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
     auto const detProp =
       art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataForJob();
 
     fWirePitch = art::ServiceHandle<geo::WireReadout>()
                    ->Get()
                    .Plane({0, 0, 0})
                    .WirePitch(); // wire pitch in cm
     fTimeTick = sampling_rate(clockData) / 1000.;
     fDriftVelocity = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
     fWireTimetoCmCm = (fTimeTick * fDriftVelocity) / fWirePitch;
   }

void shwf::ShowerReco::ClearandResizeVectors ( unsigned int nPlanes )

private

Referenced by produce().

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::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::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 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

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

void shwf::ShowerReco::GetVertexAndAnglesFromCluster	(	art::Ptr< recob::Cluster >	clust,
		unsigned int	plane
	)

private

Actual routine that reconstruct the shower

Definition at line 806 of file ShowerReco_module.cc.

References angle, DEFINE_ART_MODULE, recob::Cluster::EndTick(), recob::Cluster::EndWire(), fTime_last, fTime_vertex, fTime_vertexError, fWire_last, fWire_vertex, fWire_vertexError, recob::Cluster::SigmaStartTick(), recob::Cluster::SigmaStartWire(), slope, recob::Cluster::StartAngle(), recob::Cluster::StartTick(), and recob::Cluster::StartWire().

Referenced by produce().

   {
     // convert to cm/cm units needed in the calculation
     angle[plane] = clust->StartAngle();
     slope[plane] = std::tan(clust->StartAngle());
     fWire_vertex[plane] = clust->StartWire();
     fTime_vertex[plane] = clust->StartTick();
 
     fWire_vertexError[plane] = clust->SigmaStartWire(); // wire coordinate of vertex for each plane
     fTime_vertexError[plane] = clust->SigmaStartTick(); // time coordinate of vertex for each plane
 
     fWire_last[plane] = clust->EndWire(); // wire coordinate of last point for each plane
     fTime_last[plane] = clust->EndTick();
   }

void shwf::ShowerReco::LongTransEnergy	(	geo::GeometryCore const *	geom,
		geo::WireReadoutGeom const &	wireReadoutGeom,
		detinfo::DetectorClocksData const &	clockData,
		detinfo::DetectorPropertiesData const &	detProp,
		unsigned int	set,
		std::vector< art::Ptr< recob::Hit >>	hitlist
	)

private

third loop to get only points inside of 1RMS of value.

Definition at line 607 of file ShowerReco_module.cc.

References tca::dEdx(), calo::CalorimetryAlg::dEdx_AMP(), calo::CalorimetryAlg::dEdx_AREA(), fcalodEdxlength, fCaloPSet, fChargeMeV_2cm, fChargeMeV_2cm_refined, fCorr_MeV_2cm, fdEdxlength, fDistribChargeMeV, fDistribChargeposition, fNPitch, fNPlanes, fNpoints_2cm, fNpoints_corr_MeV_2cm, fRMS_2cm, fTime_vertex, fTotChargeADC, fTotChargeMeV, fTrkPitchC, fUseArea, fWire_vertex, fWireTimetoCmCm, geo::kCollection, Kin_En, geo::kMysteryType, geo::WireReadoutGeom::SignalType(), slope, lar::to_element, Trk_Length, vdEdx, vdQdx, vresRange, xphi, and xtheta.

Referenced by produce().

   {
     // alogorithm for energy vs dx of the shower (roto-translation) COLLECTION
     // VIEW
 
     calo::CalorimetryAlg calalg(fCaloPSet);
 
     double totCnrg = 0,
            totCnrg_corr = 0; // tot enegry of the shower in collection
 
     double time;
     unsigned int wire = 0, plane = fNPlanes - 1;
 
     double mevav2cm = 0.;
     double npoints_calo = 0;
 
     int direction = -1;
 
     //override direction if phi (XZ angle) is less than 90 degrees
     if (fabs(xphi) < 90) direction = 1;
 
     //variables to check whether a hit is close to the shower axis.
     double ortdist, linedist;
     double wire_on_line, time_on_line;
 
     //get effective pitch using 3D angles
     util::GeometryUtilities const gser{*geom, wireReadoutGeom, clockData, detProp};
     double newpitch = gser.PitchInView(plane, xphi, xtheta);
 
     using lar::to_element;
     using ranges::views::transform;
     for (auto const& hit : hitlist | transform(to_element)) {
       time = hit.PeakTime();
       wire = hit.WireID().Wire;
       plane = hit.WireID().Plane;
 
       double dEdx_new;
 
       if (fUseArea) { dEdx_new = calalg.dEdx_AREA(clockData, detProp, hit, newpitch); }
       else // this will hopefully go away, once all of the calibration factors
            // are calculated.
       {
         dEdx_new = calalg.dEdx_AMP(clockData, detProp, hit, newpitch);
       }
 
       //calculate total energy.
       totCnrg_corr += dEdx_new;
 
       // calculate the wire,time coordinates of the hit projection on to the 2D shower axis
       gser.GetPointOnLine(slope[plane] / fWireTimetoCmCm,
                           fWire_vertex[plane],
                           fTime_vertex[plane],
                           wire,
                           time,
                           wire_on_line,
                           time_on_line);
       linedist =
         gser.Get2DDistance(wire_on_line, time_on_line, fWire_vertex[plane], fTime_vertex[plane]);
       ortdist = gser.Get2DDistance(wire_on_line, time_on_line, wire, time);
 
       //calculate the distance from the vertex using the effective pitch metric
       double wdist = (((double)wire - (double)fWire_vertex[plane]) * newpitch) *
                      direction; //wdist is always positive
 
       if ((wdist < fcalodEdxlength) && (wdist > 0.2)) {
 
         vdEdx.push_back(dEdx_new);
         vresRange.push_back(fabs(wdist));
         vdQdx.push_back(hit.PeakAmplitude() / newpitch);
         Trk_Length = wdist;
         fTrkPitchC = fNPitch[set][plane];
         Kin_En += dEdx_new * newpitch;
         npoints_calo++;
 
         if (wdist < fdEdxlength &&
             ((direction == 1 && wire > fWire_vertex[plane]) // take no hits before vertex
                                                             // (depending on direction)
              || (direction == -1 && wire < fWire_vertex[plane])) &&
             ortdist < 4.5 && linedist < fdEdxlength) {
           fChargeMeV_2cm[set] += dEdx_new;
           fNpoints_2cm[set]++;
         }
 
         // fill out for 4cm preshower
 
         fDistribChargeMeV[set].push_back(dEdx_new); // vector with the first De/Dx points
         fDistribChargeposition[set].push_back(
           wdist); // vector with the first De/Dx points' positions
 
       } // end inside range if statement
 
     } // end first loop on hits.
 
     auto const signalType =
       hitlist.empty() ? geo::kMysteryType : wireReadoutGeom.SignalType(hitlist.front()->WireID());
 
     if (signalType == geo::kCollection) {
       fTotChargeADC[set] = totCnrg * newpitch;
       fTotChargeMeV[set] = totCnrg_corr * newpitch;
     }
 
     // calculate average dE/dx
     if (fNpoints_2cm[set] > 0) { mevav2cm = fChargeMeV_2cm[set] / fNpoints_2cm[set]; }
 
     // second loop to calculate RMS
     for (auto const& hit : hitlist | transform(to_element)) {
       time = hit.PeakTime();
       wire = hit.WireID().Wire;
       plane = hit.WireID().Plane;
       double dEdx = 0;
 
       if (fUseArea) { dEdx = calalg.dEdx_AREA(clockData, detProp, hit, newpitch); }
       else // this will hopefully go away, once all of the calibration factors
            // are calculated.
       {
         dEdx = calalg.dEdx_AMP(clockData, detProp, hit, newpitch);
       }
 
       gser.GetPointOnLine(slope[plane] / fWireTimetoCmCm,
                           fWire_vertex[plane],
                           fTime_vertex[plane],
                           wire,
                           time,
                           wire_on_line,
                           time_on_line);
       linedist =
         gser.Get2DDistance(wire_on_line, time_on_line, fWire_vertex[plane], fTime_vertex[plane]);
       ortdist = gser.Get2DDistance(wire_on_line, time_on_line, wire, time);
 
       double wdist = (((double)wire - (double)fWire_vertex[plane]) * newpitch) * direction;
 
       if ((wdist < fcalodEdxlength) && (wdist > 0.2)) {
         if (wdist < fdEdxlength &&
             ((direction == 1 && wire > fWire_vertex[plane]) ||
              (direction == -1 && wire < fWire_vertex[plane])) &&
             ortdist < 4.5 && linedist < fdEdxlength) {
           fRMS_2cm[set] += (dEdx - mevav2cm) * (dEdx - mevav2cm);
         }
 
       } // end if on correct hits.
     }   // end RMS_calculating loop.
 
     if (fNpoints_2cm[set] > 0) { fRMS_2cm[set] = TMath::Sqrt(fRMS_2cm[set] / fNpoints_2cm[set]); }
 
 
     for (auto const& hit : hitlist | transform(to_element)) {
       time = hit.PeakTime();
       wire = hit.WireID().Wire;
       plane = hit.WireID().Plane;
 
       double dEdx = 0;
       if (fUseArea) { dEdx = calalg.dEdx_AREA(clockData, detProp, hit, newpitch); }
       else // this will hopefully go away, once all of the calibration factors
            // are calculated.
       {
         dEdx = calalg.dEdx_AMP(clockData, detProp, hit, newpitch);
       }
 
       gser.GetPointOnLine(slope[plane] / fWireTimetoCmCm,
                           fWire_vertex[plane],
                           fTime_vertex[plane],
                           wire,
                           time,
                           wire_on_line,
                           time_on_line);
       linedist =
         gser.Get2DDistance(wire_on_line, time_on_line, fWire_vertex[plane], fTime_vertex[plane]);
       ortdist = gser.Get2DDistance(wire_on_line, time_on_line, wire, time);
 
       double wdist = (((double)wire - (double)fWire_vertex[plane]) * newpitch) * direction;
 
       if ((wdist < fcalodEdxlength) && (wdist > 0.2 &&
                                         ((direction == 1 && wire > fWire_vertex[plane]) ||
                                          (direction == -1 && wire < fWire_vertex[plane])) &&
                                         ortdist < 4.5 && linedist < fdEdxlength)) {
         if (wdist < fdEdxlength) {
           if (((dEdx > (mevav2cm - fRMS_2cm[set])) && (dEdx < (mevav2cm + fRMS_2cm[set]))) ||
               (newpitch > 0.3 * fdEdxlength)) {
             fCorr_MeV_2cm[set] += dEdx;
             fNpoints_corr_MeV_2cm[set]++;
           }
 
         } // end if on good hits
       }
     } // end of third loop on hits
 
     if (fNpoints_corr_MeV_2cm[set] > 0) {
       fCorr_MeV_2cm[set] /= fNpoints_corr_MeV_2cm[set];
       fChargeMeV_2cm_refined[set] = fCorr_MeV_2cm[set];
     }
   }

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

void shwf::ShowerReco::produce ( art::Event & evt )

privatevirtual

Get Clusters

Todo:: really need to determine the values of the arguments of the recob::Shower ctor

Implements art::EDProducer.

Definition at line 352 of file ShowerReco_module.cc.

References angle, ClearandResizeVectors(), util::CreateAssn(), deadwire, e, art::EventID::event(), fChargeADC_2cm, fChargeMeV_2cm, fClusterModuleLabel, fCorr_Charge_2cm, fCorr_MeV_2cm, fDriftVelocity, fEvent, fNAngles, fNhitsperplane, fNPlanes, fNpoints_2cm, fNpoints_corr_ADC_2cm, fNpoints_corr_MeV_2cm, fRMS_2cm, fRun, fSubRun, fTime_vertex, fTime_vertexError, fTimeTick, fTotADCperplane, fTotChargeADC, fTotChargeMeV, fTotChargeMeV_MIPs, ftree_shwf, fTrkPitchC, fWire_last, fWire_vertex, fWire_vertexError, Get, art::ProductRetriever::getByLabel(), art::ProductRetriever::getProductID(), GetVertexAndAnglesFromCluster(), hits(), art::Event::id(), util::kBogusD, art::Ptr< T >::key(), Kin_En, LongTransEnergy(), geo::PlaneID::Plane, art::ProductRetriever::productGetter(), art::PtrVector< T >::push_back(), art::Event::put(), art::EventID::run(), detinfo::sampling_rate(), recob::Shower::set_direction(), recob::Shower::set_direction_err(), art::PtrVector< T >::size(), art::EventID::subRun(), detinfo::trigger_offset(), Trk_Length, vdEdx, vdQdx, vresRange, X, xphi, xtheta, and xyz_vertex_fit.

   {
     auto const* geom = lar::providerFrom<geo::Geometry>();
     auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
     auto const detProp =
       art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
 
     util::GeometryUtilities const gser{*geom, wireReadoutGeom, clockData, detProp};
     constexpr geo::TPCID tpcid{0, 0};
     fNPlanes = wireReadoutGeom.Nplanes(tpcid);
     auto Shower3DVector = std::make_unique<std::vector<recob::Shower>>();
     auto cassn = std::make_unique<art::Assns<recob::Shower, recob::Cluster>>();
     auto hassn = std::make_unique<art::Assns<recob::Shower, recob::Hit>>();
     auto calorimetrycol = std::make_unique<std::vector<anab::Calorimetry>>();
     auto calassn = std::make_unique<art::Assns<anab::Calorimetry, recob::Shower>>();
 
     art::Handle<std::vector<recob::Cluster>> clusterListHandle;
     evt.getByLabel(fClusterModuleLabel, clusterListHandle);
 
     art::Handle<std::vector<art::PtrVector<recob::Cluster>>> clusterAssociationHandle;
     evt.getByLabel(fClusterModuleLabel, clusterAssociationHandle);
 
     art::FindManyP<recob::Hit> fmh(clusterListHandle, evt, fClusterModuleLabel);
 
     fRun = evt.id().run();
     fSubRun = evt.id().subRun();
     fEvent = evt.id().event();
 
     // find all the hits associated to all the clusters (once and for all);
     // the index of the query matches the index of the cluster in the collection
     // (conveniently carried around in its art pointer)
     art::FindManyP<recob::Hit> ClusterHits(clusterListHandle, evt, fClusterModuleLabel);
 
     std::vector<art::PtrVector<recob::Cluster>>::const_iterator clusterSet =
       clusterAssociationHandle->begin();
 
     // loop over vector of vectors (each size of NPlanes) and reconstruct showers from each of those
     for (size_t iClustSet = 0; iClustSet < clusterAssociationHandle->size(); iClustSet++) {
 
       const art::PtrVector<recob::Cluster> CurrentClusters = (*(clusterSet++));
 
       // do some error checking - i.e. are the clusters themselves present.
       if (clusterListHandle->size() < 2 || CurrentClusters.size() < 2) {
         ftree_shwf->Fill();
         return;
       }
 
       ClearandResizeVectors(fNPlanes);
 
       std::vector<std::vector<art::Ptr<recob::Hit>>> hitlist_all;
       hitlist_all.resize(fNPlanes);
 
       for (size_t iClust = 0; iClust < CurrentClusters.size(); iClust++) {
         art::Ptr<recob::Cluster> const& pclust = CurrentClusters[iClust];
 
         // get all the hits for this cluster;
         // pclust is a art::Ptr to the original cluster collection stored in the event;
         // its key corresponds to its index in the collection
         // (and therefore to the index in the query)
         std::vector<art::Ptr<recob::Hit>> const& hitlist = ClusterHits.at(pclust.key());
 
         unsigned int p(0); //c=channel, p=plane, w=wire
 
         if (hitlist.size() == 0) continue;
 
         p = (*hitlist.begin())->WireID().Plane;
         // get vertex position and slope information to start with - ii is the
         // posistion of the correct cluster:
         GetVertexAndAnglesFromCluster(pclust, p);
 
         double ADCcharge = 0;
         //loop over cluster hits
         for (art::Ptr<recob::Hit> const& hit : hitlist) {
           p = hit->WireID().Plane;
           hitlist_all[p].push_back(hit);
           ADCcharge += hit->PeakAmplitude();
         }
         fNhitsperplane[p] = hitlist_all[p].size();
         fTotADCperplane[p] = ADCcharge;
       } // End loop on clusters.
       // Now I have the Hitlists and the relevent clusters parameters saved.
 
       // find best set:
       unsigned int bp1 = 0, bp2 = 0;
       double minerror1 = 99999999, minerror2 = 9999999;
       for (unsigned int ii = 0; ii < fNPlanes; ++ii) {
         double locerror =
           fWire_vertexError[ii] * fWire_vertexError[ii] +
           fTime_vertexError[ii] * fTime_vertexError[ii]; // time coordinate of vertex for each plane
 
         if (minerror1 >= locerror) // the >= sign is to favorize collection
         {
           minerror1 = locerror;
           bp1 = ii;
         }
       }
       for (unsigned int ij = 0; ij < fNPlanes; ++ij) {
         double locerror =
           fWire_vertexError[ij] * fWire_vertexError[ij] +
           fTime_vertexError[ij] * fTime_vertexError[ij]; // time coordinate of vertex for each plane
 
         if (minerror2 >= locerror && ij != bp1) {
           minerror2 = locerror;
           bp2 = ij;
         }
       }
 
       gser.Get3DaxisN(bp1, bp2, angle[bp1], angle[bp2], xphi, xtheta);
 
       std::vector<geo::Point_t> position;
       position.reserve(fNPlanes);
       // get starting positions for all planes -- FIXME: only position[0] is used.
       for (auto const& plane : wireReadoutGeom.Iterate<geo::PlaneGeo>(tpcid)) {
         position.push_back(plane.GetBoxCenter());
       }
 
       // Assuming there is no problem ( and we found the best pair that comes
       // close in time ) we try to get the Y and Z coordinates for the start of
       // the shower.
       double fTimeTick = sampling_rate(clockData) / 1000.;
       double fDriftVelocity = detProp.DriftVelocity(detProp.Efield(), detProp.Temperature());
       try {
         int chan1 = wireReadoutGeom.PlaneWireToChannel({0, 0, bp1, fWire_vertex[bp1]});
         int chan2 = wireReadoutGeom.PlaneWireToChannel({0, 0, bp2, fWire_vertex[bp2]});
 
         auto const intersection = wireReadoutGeom.ChannelsIntersect(chan1, chan2).value();
         xyz_vertex_fit[1] = intersection.y;
         xyz_vertex_fit[2] = intersection.z;
         xyz_vertex_fit[0] =
           (fTime_vertex[bp1] - trigger_offset(clockData)) * fDriftVelocity * fTimeTick +
           position[0].X();
       }
       catch (cet::exception const& e) {
         mf::LogWarning("ShowerReco") << "caught exception \n" << e;
         xyz_vertex_fit[1] = 0;
         xyz_vertex_fit[2] = 0;
         xyz_vertex_fit[0] = 0;
       }
 
       // if collection is not best plane, project starting point from that
       if (bp1 != fNPlanes - 1 && bp2 != fNPlanes - 1) {
         geo::PlaneID const lastPlaneID{0, 0, fNPlanes - 1};
         auto const& plane = wireReadoutGeom.Plane(lastPlaneID);
         auto pos = plane.GetBoxCenter();
         pos.SetY(xyz_vertex_fit[1]);
         pos.SetZ(xyz_vertex_fit[2]);
         auto const wirevertex = plane.NearestWireID(pos).Wire;
 
         double drifttick =
           (xyz_vertex_fit[0] / detProp.DriftVelocity(detProp.Efield(), detProp.Temperature())) *
           (1. / fTimeTick);
         fWire_vertex[fNPlanes - 1] = wirevertex; // wire coordinate of vertex for each plane
         fTime_vertex[fNPlanes - 1] =
           drifttick -
           (pos.X() / detProp.DriftVelocity(detProp.Efield(), detProp.Temperature())) *
             (1. / fTimeTick) +
           trigger_offset(clockData);
       }
 
       if (fabs(xphi) < 5.) {
         xtheta = gser.Get3DSpecialCaseTheta(
           bp1, bp2, fWire_last[bp1] - fWire_vertex[bp1], fWire_last[bp2] - fWire_vertex[bp2]);
       }
 
       // zero the arrays just to make sure
       for (unsigned int i = 0; i < fNAngles; ++i) {
         fTotChargeADC[i] = 0;
         fTotChargeMeV[i] = 0;
         fTotChargeMeV_MIPs[i] = 0;
         fNpoints_corr_ADC_2cm[i] = 0;
         fNpoints_corr_MeV_2cm[i] = 0;
 
         fRMS_2cm[i] = 0;
         fNpoints_2cm[i] = 0;
 
         fCorr_MeV_2cm[i] = 0;
         fCorr_Charge_2cm[i] = 0;
 
         fChargeADC_2cm[i] = 0; //Initial charge in ADC/cm for each plane angle calculation 4cm
         fChargeMeV_2cm[i] = 0; //initial charge in MeV/cm for each angle calculation first 4cm
       }
 
       // do loop and choose Collection. With ful calorimetry can do all.
       if (!(fabs(xphi) > 89 && fabs(xphi) < 91)) // do not calculate pitch for extreme angles
         LongTransEnergy(geom,
                         wireReadoutGeom,
                         clockData,
                         detProp,
                         0,
                         hitlist_all[fNPlanes - 1]); // temporary only plane 2.
 
 
       // make an art::PtrVector of the clusters
       art::PtrVector<recob::Cluster> prodvec;
       for (unsigned int i = 0; i < clusterListHandle->size(); ++i) {
         art::Ptr<recob::Cluster> prod(clusterListHandle, i);
         prodvec.push_back(prod);
       }
 
       //create a singleSpacePoint at vertex.
       std::vector<recob::SpacePoint> spcpts;
 
       //get direction cosines and set them for the shower
       // TBD determine which angle to use for the actual shower
       double fPhi = xphi;
       double fTheta = xtheta;
 
       TVector3 dcosVtx(std::cos(fPhi * TMath::Pi() / 180) * std::sin(fTheta * TMath::Pi() / 180),
                        std::cos(fTheta * TMath::Pi() / 180),
                        std::sin(fPhi * TMath::Pi() / 180) * std::sin(fTheta * TMath::Pi() / 180));
       // fill with bogus values for now
       TVector3 dcosVtxErr(util::kBogusD, util::kBogusD, util::kBogusD);
       recob::Shower singShower;
       singShower.set_direction(dcosVtx);
       singShower.set_direction_err(dcosVtxErr);
 
       Shower3DVector->push_back(singShower);
       // associate the shower with its clusters
       util::CreateAssn(evt, *Shower3DVector, prodvec, *cassn);
 
       // get the hits associated with each cluster and associate those with the shower
       for (size_t p = 0; p < prodvec.size(); ++p) {
         std::vector<art::Ptr<recob::Hit>> hits = fmh.at(p);
         util::CreateAssn(evt, *Shower3DVector, hits, *hassn);
       }
 
       geo::PlaneID planeID(0, 0, fNPlanes - 1);
       calorimetrycol->emplace_back(
         Kin_En, vdEdx, vdQdx, vresRange, deadwire, Trk_Length, fTrkPitchC, planeID);
 
       art::PtrVector<recob::Shower> ssvec;
 
       art::ProductID aid = evt.getProductID<std::vector<recob::Shower>>();
       art::Ptr<recob::Shower> aptr(aid, 0, evt.productGetter(aid));
       ssvec.push_back(aptr);
 
       util::CreateAssn(evt, *calorimetrycol, ssvec, *calassn);
       ftree_shwf->Fill();
     } // end loop on Vectors of "Associated clusters"
 
     evt.put(std::move(Shower3DVector));
     evt.put(std::move(cassn));
     evt.put(std::move(hassn));
     evt.put(std::move(calorimetrycol));
     evt.put(std::move(calassn));
   }

void art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

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

Member Data Documentation

float shwf::ShowerReco::angle[3]

private

Definition at line 85 of file ShowerReco_module.cc.

Referenced by GetVertexAndAnglesFromCluster(), and produce().

std::vector<float> shwf::ShowerReco::deadwire

private

Definition at line 142 of file ShowerReco_module.cc.

Referenced by produce().

float shwf::ShowerReco::fcalodEdxlength

private

Definition at line 147 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and ShowerReco().

fhicl::ParameterSet shwf::ShowerReco::fCaloPSet

private

Definition at line 92 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and ShowerReco().

std::vector<double> shwf::ShowerReco::fChargeADC_2cm

private

Definition at line 106 of file ShowerReco_module.cc.

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

std::vector<double> shwf::ShowerReco::fChargeMeV_2cm

private

Definition at line 107 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fChargeMeV_2cm_axsum

private

Definition at line 110 of file ShowerReco_module.cc.

Referenced by beginJob(), and beginRun().

std::vector<double> shwf::ShowerReco::fChargeMeV_2cm_refined

private

Definition at line 109 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

std::string shwf::ShowerReco::fClusterModuleLabel

private

Definition at line 87 of file ShowerReco_module.cc.

Referenced by produce(), and ShowerReco().

std::vector<double> shwf::ShowerReco::fCorr_Charge_2cm

private

Definition at line 97 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

std::vector<double> shwf::ShowerReco::fCorr_MeV_2cm

private

Definition at line 96 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

float shwf::ShowerReco::fdEdxlength

private

Definition at line 145 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and ShowerReco().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribChargeADC

private

Definition at line 112 of file ShowerReco_module.cc.

Referenced by beginJob(), and beginRun().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribChargeMeV

private

Definition at line 114 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribChargeposition

private

Definition at line 117 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribHalfChargeMeV

private

Definition at line 115 of file ShowerReco_module.cc.

Referenced by beginJob(), and beginRun().

double shwf::ShowerReco::fDriftVelocity

private

Definition at line 158 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

int shwf::ShowerReco::fEvent

private

Definition at line 82 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

double shwf::ShowerReco::fMean_wire_pitch

private

Definition at line 91 of file ShowerReco_module.cc.

Referenced by beginJob().

unsigned int shwf::ShowerReco::fNAngles

private

Definition at line 152 of file ShowerReco_module.cc.

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

std::vector<int> shwf::ShowerReco::fNhitsperplane

private

Definition at line 161 of file ShowerReco_module.cc.

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

std::vector<std::vector<double> > shwf::ShowerReco::fNPitch

private

Definition at line 135 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

unsigned int shwf::ShowerReco::fNPlanes

private

Definition at line 151 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<int> shwf::ShowerReco::fNpoints_2cm

private

Definition at line 95 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<int> shwf::ShowerReco::fNpoints_corr_ADC_2cm

private

Definition at line 99 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

std::vector<int> shwf::ShowerReco::fNpoints_corr_MeV_2cm

private

Definition at line 100 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fRMS_2cm

private

Definition at line 94 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

int shwf::ShowerReco::fRun

private

Definition at line 82 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

std::vector<std::vector<double> > shwf::ShowerReco::fSingleEvtAngle

private

Definition at line 119 of file ShowerReco_module.cc.

std::vector<std::vector<double> > shwf::ShowerReco::fSingleEvtAngleVal

private

Definition at line 120 of file ShowerReco_module.cc.

int shwf::ShowerReco::fSubRun

private

Definition at line 82 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

std::vector<double> shwf::ShowerReco::fTime_last

private

Definition at line 129 of file ShowerReco_module.cc.

Referenced by beginRun(), and GetVertexAndAnglesFromCluster().

std::vector<double> shwf::ShowerReco::fTime_vertex

private

Definition at line 123 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fTime_vertexError

private

Definition at line 126 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), and produce().

float shwf::ShowerReco::ftimetick

private

Definition at line 89 of file ShowerReco_module.cc.

Referenced by beginJob().

double shwf::ShowerReco::fTimeTick

private

Definition at line 157 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

std::vector<double> shwf::ShowerReco::fTotADCperplane

private

Definition at line 162 of file ShowerReco_module.cc.

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

std::vector<double> shwf::ShowerReco::fTotChargeADC

private

Definition at line 102 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fTotChargeMeV

private

Definition at line 103 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fTotChargeMeV_MIPs

private

Definition at line 104 of file ShowerReco_module.cc.

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

TTree* shwf::ShowerReco::ftree_shwf

private

Definition at line 153 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

float shwf::ShowerReco::fTrkPitchC

private

Definition at line 144 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and produce().

bool shwf::ShowerReco::fUseArea

private

Definition at line 148 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and ShowerReco().

std::vector<unsigned int> shwf::ShowerReco::fWire_last

private

Definition at line 128 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), and produce().

std::vector<unsigned int> shwf::ShowerReco::fWire_vertex

private

Definition at line 122 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fWire_vertexError

private

Definition at line 125 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), and produce().

double shwf::ShowerReco::fWirePitch

private

Definition at line 156 of file ShowerReco_module.cc.

Referenced by beginRun().

double shwf::ShowerReco::fWireTimetoCmCm

private

Definition at line 159 of file ShowerReco_module.cc.

Referenced by beginRun(), and LongTransEnergy().

float shwf::ShowerReco::Kin_En

private

Definition at line 138 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and produce().

float shwf::ShowerReco::slope[3]

private

Definition at line 84 of file ShowerReco_module.cc.

Referenced by GetVertexAndAnglesFromCluster(), and LongTransEnergy().

float shwf::ShowerReco::Trk_Length

private

Definition at line 143 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and produce().

std::vector<float> shwf::ShowerReco::vdEdx

private

Definition at line 139 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

std::vector<float> shwf::ShowerReco::vdQdx

private

Definition at line 141 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

std::vector<float> shwf::ShowerReco::vresRange

private

Definition at line 140 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

double shwf::ShowerReco::xphi

private

Definition at line 150 of file ShowerReco_module.cc.

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

double shwf::ShowerReco::xtheta

private

Definition at line 150 of file ShowerReco_module.cc.

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

std::vector<double> shwf::ShowerReco::xyz_vertex_fit

private

Definition at line 132 of file ShowerReco_module.cc.

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

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

larreco/v10_01_07/source/larreco/ShowerFinder/ShowerReco_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