sbn Namespace Reference

This module creates Common Analysis Files. More...

Namespaces
	bits
	Simple utilities to deal with bit enumerators.
	crt
	details
	evwgh
	util

Classes
struct	selTrackInfo
struct	selBeamInfo
struct	selTriggerInfo
struct	selSimTriggerInfo
struct	selLightInfo
struct	selHitInfo
class	TimedTrackSelector
	Selects tracks with time information. More...
class	TimeTrackTreeStorage
	Fills a ROOT tree with track-based triggering information. More...
class	DumpArtDAQfragments
	Dumps on console the content of artdaq::Fragment collections. More...
class	DumpPMTconfiguration
	Dumps on console the content of sbn::PMTconfiguration data product. More...
class	DumpTrigger
	Dumps on console the content of trigger data products. More...
class	DumpTriggerConfiguration
	Dumps on console the content of icarus::TriggerConfiguration data product. More...
struct	OpDetWaveformMeta
	Derivative information from raw::OpDetWaveform data. More...
class	OpDetWaveformMetaMaker
	Converter from raw::OpDetWaveform into sbn::OpDetWaveformMeta. More...
class	TriggerEmulationTree
	Fills a ROOT tree with trigger emulation results. More...
class	DuplicateEventTracker
	Keeps track of duplicate events in the job. More...
class	EventRegistry
	Keeps a record of all registered events and their source. More...
class	BNBRetriever
class	MWRData
class	EXTRetriever
class	NuMIRetriever
class	ITCSSelectionTool
class	TrackCaloSkimmer
class	TrackCaloSkimmerSelectAllTrack
class	TrackCaloSkimmerSelectAnode2CathodeTrack
class	TrackCaloSkimmerSelectStoppingTrack
class	CRUMBS
struct	CRTSipmGeo
struct	CRTStripGeo
struct	CRTModuleGeo
struct	CRTTaggerGeo
class	CRTGeoAlg
class	TPCGeoAlg
class	CRTHitConverter
class	MCSFitAllPID
class	PMTFlashTriggerFilter
class	PMTFlashTriggerMaker
class	RangePAllPID
class	ShowerCosmicDistance
class	ShowerSelectionVars
class	TrackScatterClosestApproach
class	TrackStoppingChi2Alg
class	TrackStoppingChi2Fitter
class	TransferPFParticleFlashMatch
class	TransferPFParticleT0
class	TransferTrackT0
class	Dazzle
class	Razzle
class	PCAngleKinkFinder
class	PCAngleKinkTree
struct	PCAngleInfo
class	PCAnglePlaneMaker
class	NuVertexChargeTree
	Analyzer module for use with sbn::Stub and sbn::VertexHit objects. More...
class	TrackAreaHit
class	TrackHitDumper
class	TrackHitFilter
class	MergedTrackIdentifier
class	TrackSplitter
class	IStubMerge
	Abstract interface intended for art tools which take a list of stubs and return a new list with some of them merged. More...
class	PlaneStubMerge
	Art tool for merging stubs on the same plane. More...
class	StubBuilder
struct	StubInfo
	Internal struct: contains information on stub and other associated data products. More...
class	TwoPlaneStubMerge
	Art tool for merging stubs across planes. More...
class	VertexChargeVacuum
class	VertexStubTracker
struct	Vector3D
struct	WireInfo
struct	HitTruth
struct	HitInfo
struct	TrackHitInfo
struct	MetaInfo
struct	TrueHit
struct	TrueParticle
struct	TrackTruth
struct	TrackInfo
struct	PMTconfiguration
	Class containing configuration for PMT readout. More...
struct	V1730channelConfiguration
	Class containing configuration for a V1730 channel. More...
struct	V1730Configuration
	Class containing configuration for a V1730 board. More...
class	BNBSpillInfo
class	EXTCountInfo
class	NuMISpillInfo
class	CRUMBSResult
struct	FlashTriggerPrimitive
class	MergedTrackInfo
class	MVAPID
class	RangeP
class	ScatterClosestApproach
class	ShowerDensityFit
class	ShowerTrackFit
class	SimpleFlashMatch
class	StoppingChi2Fit
class	StubHit
class	Stub
class	VertexHit
struct	ExtraTriggerInfo
	Additional information on trigger. More...

Enumerations
enum	EDet { kNOTDEFINED, kSBND, kICARUS }

Functions
TimeTrackTreeStorage::TriggerInputSpec_t	convert (TimeTrackTreeStorage::Config::TriggerSpecConfig const &config)
TriggerEmulationTree::TriggerInputSpec_t	convert (TriggerEmulationTree::Config::TriggerSpecConfig const &config)
double	GetPitch (const geo::GeometryCore *geo, const spacecharge::SpaceCharge *sce, geo::Point_t loc, geo::Vector_t dir, geo::View_t view, geo::TPCID tpc, bool correct_sce, bool track_is_sce_corrected, float xsign=1.)
	Computes the track-pitch on a plane given an input direction and location. More...
geo::Point_t	GetLocation (const spacecharge::SpaceCharge *sce, geo::Point_t loc_w, geo::TPCID TPC, float xsign=1.)
	Get the location in the presence of space charge. More...
double	GetEfield (const detinfo::DetectorPropertiesData &dprop, const spacecharge::SpaceCharge *sce, geo::Point_t loc, geo::TPCID TPC, bool correct_loc_sce, float xsign=1.)
	Get the E-Field in the presence of space charge. More...
geo::Point_t	GetLocationAtWires (const spacecharge::SpaceCharge *sce, const geo::GeometryCore *geo, geo::Point_t loc, geo::TPCID TPC, float xsign=1.)
	Get the SCE-distorted location (i.e. the location "seen" by the wireplanes) More...
bool	StubContains (const sbn::StubInfo &A, const sbn::StubInfo &B)
	Returns whether stub A contains stub B. More...
float	StubDirectionDot (const sbn::StubInfo &A, const sbn::StubInfo &B, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
	Computes the dot product of two stubs. More...
float	StubTimeOffset (const sbn::StubInfo &A, const sbn::StubInfo &B, const detinfo::DetectorClocksData &dclock, const detinfo::DetectorPropertiesData &dprop)
geo::Point_t	TwoStubEndPosition (const sbn::StubInfo &A, const sbn::StubInfo &B, const geo::GeometryCore *geo, const spacecharge::SpaceCharge *sce, const detinfo::DetectorPropertiesData &dprop)
	Returns an updated end position of a stub after merging across two planes. More...
float	StubChargeOffset (const sbn::StubInfo &A, const sbn::StubInfo &B)
	Difference of the total charge between two stubs. More...
float	StubPeakChargeOffset (const sbn::StubInfo &A, const sbn::StubInfo &B)
	Difference of the endpoint charge between two stubs. More...
float	StubPeakdQdxOffset (const sbn::StubInfo &A, const sbn::StubInfo &B, const geo::GeometryCore *geo, const spacecharge::SpaceCharge *sce, const detinfo::DetectorPropertiesData &dprop)
	Difference of the endpoint dQ/dx between two stubs. More...
std::ostream &	operator<< (std::ostream &out, sbn::PMTconfiguration const &config)
	Prints the configuration into a stream with default verbosity. More...
std::ostream &	operator<< (std::ostream &out, sbn::V1730channelConfiguration const &config)
	Prints the configuration into a stream with default verbosity. More...
std::ostream &	operator<< (std::ostream &out, sbn::V1730Configuration const &config)
	Prints the configuration into a stream with default verbosity. More...
std::ostream &	operator<< (std::ostream &out, ExtraTriggerInfo const &info)
Creation of `sbn::OpDetWaveformMeta` from `raw::OpDetWaveform`
The creation of summary objects sbn::OpDetWaveformMeta from optical detector waveforms is possible with two options: one shot: call to convert a single waveform (makeOpDetWaveformMeta()); bulk: converter object reused for multiple conversions (sbn::OpDetWaveformMetaMaker). For usage examples, see their respective documentation.
OpDetWaveformMeta	makeOpDetWaveformMeta (raw::OpDetWaveform const &waveform, detinfo::DetectorTimings const &detTimings)
	Creates a sbn::OpDetWaveformMeta out of a raw::OpDetWaveform. More...
OpDetWaveformMeta	makeOpDetWaveformMeta (raw::OpDetWaveform const &waveform, util::quantities::intervals::microseconds opDetTickPeriod)
	Creates a sbn::OpDetWaveformMeta out of a raw::OpDetWaveform. More...

Variables
constexpr int	default_int = std::numeric_limits<int>::max()
constexpr float	default_float = std::numeric_limits<float>::signaling_NaN()

Detailed Description

This module creates Common Analysis Files.

Set of values to identify a beam.

Multivariate Gaussian sampling utilities.

Original author: J. Zennamo.

Note: In larsim EventWeight, these functions are static members of the WeightCalc class.

The constants are used to identify the gate during which a trigger is acquired, in the raw::Trigger::TriggerBits() bit mask.

Enumeration Type Documentation

enum sbn::EDet

Enumerator
kNOTDEFINED
kSBND
kICARUS

Definition at line 70 of file TrackCaloSkimmer.h.

{kNOTDEFINED, kSBND, kICARUS}; 

Function Documentation

TriggerEmulationTree::TriggerInputSpec_t sbn::convert

(

TriggerEmulationTree::Config::TriggerSpecConfig const &

config

)

Definition at line 231 of file TriggerEmulationTree_module.cc.

  {
    return {
        config.Name()       // name
      , config.TriggerTag() // inputTag
      };
  } // convert(sbn::TriggerSpecConfig)

TimeTrackTreeStorage::TriggerInputSpec_t sbn::convert

(

TimeTrackTreeStorage::Config::TriggerSpecConfig const &

config

)

Definition at line 388 of file TimeTrackTreeStorage_module.cc.

  {
    return {
        config.Name()       // name
      , config.TriggerTag() // inputTag
      };
  } // convert(sbn::TriggerSpecConfig)

double sbn::GetEfield	(	const detinfo::DetectorPropertiesData &	dprop,
		const spacecharge::SpaceCharge *	sce,
		geo::Point_t	loc,
		geo::TPCID	TPC,
		bool	correct_loc_sce,
		float	xsign = 1.
	)

Get the E-Field in the presence of space charge.

Definition at line 14 of file StubMergeAlgorithms.cxx.

                                                                                                                                                                    {

  double EField = dprop.Efield();
  if (sce && sce->EnableSimEfieldSCE()) {
    if (correct_loc_sce) loc = sbn::GetLocation(sce, loc, TPC, xsign);

    // Gets relative E field Distortions
    geo::Vector_t EFieldOffsets = sce->GetEfieldOffsets(loc);
    // Add 1 in X direction as this is the direction of the drift field
    EFieldOffsets = EFieldOffsets + geo::Xaxis();
    // Convert to Absolute E Field from relative
    EFieldOffsets = EField * EFieldOffsets;
    // We only care about the magnitude for recombination
    EField = EFieldOffsets.r();
  }

  return EField;
}

geo::Point_t sbn::GetLocation	(	const spacecharge::SpaceCharge *	sce,
		geo::Point_t	loc_w,
		geo::TPCID	TPC,
		float	xsign = 1.
	)

Get the location in the presence of space charge.

Definition at line 3 of file StubMergeAlgorithms.cxx.

                                                                                                          {
  if (sce && sce->EnableCalSpatialSCE()) {
    geo::Vector_t offset = sce->GetCalPosOffsets(loc_w, TPC.TPC);
    loc_w.SetX(loc_w.X() + xsign * offset.X());
    loc_w.SetY(loc_w.Y() + offset.Y());
    loc_w.SetZ(loc_w.Z() + offset.Z());
  }

  return loc_w;
}

geo::Point_t sbn::GetLocationAtWires	(	const spacecharge::SpaceCharge *	sce,
		const geo::GeometryCore *	geo,
		geo::Point_t	loc,
		geo::TPCID	TPC,
		float	xsign = 1.
	)

Get the SCE-distorted location (i.e. the location "seen" by the wireplanes)

Definition at line 33 of file StubMergeAlgorithms.cxx.

                                                                                                                                           { 
  if (sce && sce->EnableCalSpatialSCE()) {
    // Returned X is the drift -- multiply by the drift direction to undo this
    int corr = geo->TPC(TPC).DriftDir()[0];

    geo::Vector_t offset = sce->GetPosOffsets(loc);
    loc.SetX(loc.X() + corr * xsign * offset.X());
    loc.SetY(loc.Y() + offset.Y());
    loc.SetZ(loc.Z() + offset.Z());
  }

  return loc;
}

double sbn::GetPitch	(	const geo::GeometryCore *	geo,
		const spacecharge::SpaceCharge *	sce,
		geo::Point_t	loc,
		geo::Vector_t	dir,
		geo::View_t	view,
		geo::TPCID	tpc,
		bool	correct_sce,
		bool	track_is_sce_corrected,
		float	xsign = 1.
	)

Computes the track-pitch on a plane given an input direction and location.

Parameters

loc	location in world coordinate system [cm]
dir	direction in world coordinate system (assumed to be unit vector)
view	the view to project the direction onto
tpc	ID of the TPC the selected view belongs to
correct_sce	whether to apply space charge correction
track_is_sce_corrected	whether location loc is already corrected for spacecharge effects
xsign	factor for the x (drift) coordinate (-1 flips it)

Returns

the pitch along direction dir [cm]

Is able to handle the presence of space charge (through the SpaceCharge service) and whether or not the input location/direction are already corrected for SCE.

Definition at line 47 of file StubMergeAlgorithms.cxx.

                                                                {

  double angleToVert = geo->WireAngleToVertical(view, tpc) - 0.5*::util::pi<>();

  geo::Vector_t dir_w;

  // "dir_w" should be the direction that the wires see. If the track already has the field
  // distortion corrections applied, then we need to de-apply them to get the direction as
  // seen by the wire planes
  if (correct_sce && track_is_sce_corrected) {
    // compute the dir of the track trajectory
    geo::Point_t loc_mdx = loc - dir * (geo->WirePitch(view) / 2.);
    geo::Point_t loc_pdx = loc + dir * (geo->WirePitch(view) / 2.);

    // map to the wires
    loc_mdx = GetLocationAtWires(sce, geo, loc_mdx, tpc, xsign);
    loc_pdx = GetLocationAtWires(sce, geo, loc_pdx, tpc, xsign);

    dir_w = (loc_pdx - loc_mdx).Unit(); 
  }
  // If there is no space charge or the track is not yet corrected, then the dir
  // is what we want
  else {
    dir_w = dir;
  }

  double cosgamma = std::abs(std::sin(angleToVert)*dir_w.Y() + std::cos(angleToVert)*dir_w.Z());
  double pitch;
  if (cosgamma) {
    pitch = geo->WirePitch(view)/cosgamma;
  }
  else {
    pitch = 0.;
  }

  // now take the pitch computed on the wires and correct it back to the particle trajectory
  geo::Point_t loc_w = loc;
  if (correct_sce && track_is_sce_corrected) {
    loc_w = sbn::GetLocationAtWires(sce, geo, loc, tpc, xsign);
  }

  geo::Point_t locw_traj = (correct_sce) ? sbn::GetLocation(sce, loc_w, tpc, xsign) : loc_w;
  geo::Point_t locw_pdx_traj = (correct_sce) ? sbn::GetLocation(sce, loc_w + pitch * dir_w, tpc, xsign) : (loc_w + pitch * dir_w);

  pitch = (locw_traj - locw_pdx_traj).R();

  return pitch;
}

sbn::OpDetWaveformMeta sbn::makeOpDetWaveformMeta	(	raw::OpDetWaveform const &	waveform,
		detinfo::DetectorTimings const &	detTimings
	)

Creates a sbn::OpDetWaveformMeta out of a raw::OpDetWaveform.

Parameters

waveform	the input waveform
detTimings	timing service provider

Returns

a sbn::OpDetWaveformMeta object with the summary information

Returns a new summary object extracted from the waveform.

The timing information is used to determine whether the global trigger and beam times are included.

Example:

detinfo::DetectorTimings const detTimings = detinfo::makeDetectorTimings
  (art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(event));

sbn::OpDetWaveformMeta info = sbn::makeOpDetWaveformMeta(waveform, detTimings);

Definition at line 83 of file OpDetWaveformMetaUtils.cxx.

  {
  return sbn::OpDetWaveformMetaMaker{ detTimings }.make(waveform);
}

sbn::OpDetWaveformMeta sbn::makeOpDetWaveformMeta	(	raw::OpDetWaveform const &	waveform,
		util::quantities::intervals::microseconds	opDetTickPeriod
	)

Creates a sbn::OpDetWaveformMeta out of a raw::OpDetWaveform.

Parameters

waveform	the input waveform
opDetTickPeriod	period of the optical detector digitizer

Returns

a sbn::OpDetWaveformMeta object with the summary information

Returns a new summary object extracted from the waveform.

Information requiring timing is not saved.

Example:

util::quantities::intervals::nanoseconds const opDetPeriod { 2.0 };

sbn::OpDetWaveformMeta info = sbn::makeOpDetWaveformMeta(waveform, opDetPeriod);

Definition at line 91 of file OpDetWaveformMetaUtils.cxx.

  {
  return sbn::OpDetWaveformMetaMaker{ opDetTickPeriod }.make(waveform);
}

std::ostream & sbn::operator<< ( std::ostream &  out, sbn::PMTconfiguration const &  config  )

inline

Prints the configuration into a stream with default verbosity.

Definition at line 150 of file PMTconfiguration.h.

  { config.dump(out); return out; }

std::ostream & sbn::operator<< ( std::ostream &  out, sbn::V1730Configuration const &  config  )

inline

Prints the configuration into a stream with default verbosity.

Definition at line 233 of file V1730Configuration.h.

  { config.dump(out); return out; }

std::ostream & sbn::operator<< ( std::ostream &  out, sbn::V1730channelConfiguration const &  config  )

inline

Prints the configuration into a stream with default verbosity.

Definition at line 191 of file V1730channelConfiguration.h.

  { config.dump(out); return out; }

std::ostream & sbn::operator<<	(	std::ostream &	out,
		ExtraTriggerInfo const &	info
	)

Definition at line 114 of file ExtraTriggerInfo.cxx.

{
  if (!info.isValid()) {
    out << "<invalid>";
    return out;
  }
  
  // quite a load:
  out
    <<   "trigger ID=" << dumpTriggerID(info.triggerID)
      << " (type " << name(info.triggerType) << ")"
      << " from source " << name(info.sourceType)
      << " at " << dumpTimestamp(info.triggerTimestamp)
      << " on beam gate ID=" << dumpTriggerID(info.gateID)
      << " at " << dumpTimestamp(info.beamGateTimestamp)
      << " (diff: "
      << timestampDiff(info.beamGateTimestamp, info.triggerTimestamp) << " ns)"
    << "\n"
    << "enable gate opened at " << dumpTimestamp(info.enableGateTimestamp)
    << " (" << timestampDiff(info.beamGateTimestamp, info.enableGateTimestamp)
    << " ns before the gate)"
    << "\n"
      << "counts from this source: trigger="
        << dumpTriggerCount(info.triggerCount)
      << " beam=" << dumpTriggerCount(info.gateCount)
    << "\n"
      << "previous trigger from this source at "
      << dumpTimestamp(info.previousTriggerTimestamp)
      << ", triggers since: "
      << dumpTriggerCount(info.anyTriggerCountFromPreviousTrigger)
      << ", gates since: "
      << dumpTriggerCount(info.anyGateCountFromPreviousTrigger)
      << " ("
      << dumpTriggerCount(info.gateCountFromPreviousTrigger)
      << " from this same source)"
    << "\n"
      << "most recent trigger was from source "
      << name(info.anyPreviousTriggerSourceType) << " at "
      << dumpTimestamp(info.anyPreviousTriggerTimestamp)
    ;
  if (sbn::ExtraTriggerInfo::isValidTimestamp(info.anyPreviousTriggerTimestamp)
    && sbn::ExtraTriggerInfo::isValidTimestamp(info.triggerTimestamp))
  {
    out
      << " ("
      << dumpTimestamp(info.triggerTimestamp - info.anyPreviousTriggerTimestamp)
      << " earlier)";
  }
  out
      << ", and "
      << dumpTriggerCount(info.anyGateCountFromAnyPreviousTrigger)
      << " gates from any source have opened since"
    ;
  if (info.WRtimeToTriggerTime != ExtraTriggerInfo::UnknownCorrection) {
    out << "\nCorrection applied to the timestamps: "
      << dumpTimestamp(info.WRtimeToTriggerTime);
  }
  if (info.triggerLocationBits != 0) {
    out << "\nLocation(s) of trigger:";
    for (std::string const& bitName: names(info.triggerLocation()))
      out << " " << bitName;
  }
  out << "\nWest cryostat: "
      << info.cryostats[ExtraTriggerInfo::WestCryostat].triggerCount
      << " triggers";
  if (auto const& cryo = info.cryostats[ExtraTriggerInfo::WestCryostat];
      cryo.hasLVDS()
      ) {
    out
      << "\n  west wall:  "
      << dumpLVDSmask(cryo.LVDSstatus[ExtraTriggerInfo::WestPMTwall])
      << "\n  east wall:  "
      << dumpLVDSmask(cryo.LVDSstatus[ExtraTriggerInfo::EastPMTwall])
      ;
  }

  out << "\nEast cryostat: "
      << info.cryostats[ExtraTriggerInfo::EastCryostat].triggerCount
      << " triggers";
  if (auto const& cryo = info.cryostats[ExtraTriggerInfo::EastCryostat];
      cryo.hasLVDS()
      ) {
    out
      << "\n  west wall:  "
      << dumpLVDSmask(cryo.LVDSstatus[ExtraTriggerInfo::WestPMTwall])
      << "\n  east wall:  "
      << dumpLVDSmask(cryo.LVDSstatus[ExtraTriggerInfo::EastPMTwall])
      ;
  }

  
  return out;
} // sbn::operator<< (ExtraTriggerInfo)

float sbn::StubChargeOffset	(	const sbn::StubInfo &	A,
		const sbn::StubInfo &	B
	)

Difference of the total charge between two stubs.

Definition at line 195 of file StubMergeAlgorithms.cxx.

                                                                      {
  assert(A.stub.hits.size() == 1 && B.stub.hits.size() == 1);

  // Only count charge on the main stub
  float ACharge = A.stub.CoreCharge();
  float BCharge = B.stub.CoreCharge();

  return abs(ACharge - BCharge);
}

bool sbn::StubContains	(	const sbn::StubInfo &	A,
		const sbn::StubInfo &	B
	)

Returns whether stub A contains stub B.

Definition at line 102 of file StubMergeAlgorithms.cxx.

                                                                 {
  if (A.vhit_hit->WireID().Plane != B.vhit_hit->WireID().Plane) return false;

  for (const art::Ptr<recob::Hit> &h: A.hits) {
    if (A.stub.OnCore(h->WireID())) {
      if (h.key() == B.vhit_hit.key()) {
        return true;
      }
    }
  }
  return false;
}

float sbn::StubDirectionDot	(	const sbn::StubInfo &	A,
		const sbn::StubInfo &	B,
		const geo::GeometryCore *	geo,
		const detinfo::DetectorPropertiesData &	dprop
	)

Computes the dot product of two stubs.

Definition at line 115 of file StubMergeAlgorithms.cxx.

                                                {

  // A and B should chare vertex locations
  //
  // We want this to be the vertex as seen by the wires
  float vert_w = A.vhit->vtxw;
  float vert_x = A.vhit->vtxx;;

  float Ahit_w = A.vhit_hit->WireID().Wire * geo->WirePitch() - vert_w;
  float Ahit_x = dprop.ConvertTicksToX(A.vhit_hit->PeakTime(), A.vhit_hit->WireID()) - vert_x;

  float Bhit_w = B.vhit_hit->WireID().Wire * geo->WirePitch() - vert_w;
  float Bhit_x = dprop.ConvertTicksToX(B.vhit_hit->PeakTime(), B.vhit_hit->WireID()) - vert_x;

  float Amag = sqrt(Ahit_w*Ahit_w + Ahit_x*Ahit_x);
  float Bmag = sqrt(Bhit_w*Bhit_w + Bhit_x*Bhit_x);

  return(Ahit_w*Bhit_w + Ahit_x*Bhit_x) / (Amag*Bmag);
}

float sbn::StubPeakChargeOffset	(	const sbn::StubInfo &	A,
		const sbn::StubInfo &	B
	)

Difference of the endpoint charge between two stubs.

Definition at line 205 of file StubMergeAlgorithms.cxx.

                                                                          {
  return abs(A.vhit->charge - B.vhit->charge);
}

float sbn::StubPeakdQdxOffset	(	const sbn::StubInfo &	A,
		const sbn::StubInfo &	B,
		const geo::GeometryCore *	geo,
		const spacecharge::SpaceCharge *	sce,
		const detinfo::DetectorPropertiesData &	dprop
	)

Difference of the endpoint dQ/dx between two stubs.

Definition at line 209 of file StubMergeAlgorithms.cxx.

                                                {

  // stubs should share the vertex
  geo::Point_t vtx = A.stub.vtx; 
  geo::Point_t end = sbn::TwoStubEndPosition(A, B, geo, sce, dprop); 

  geo::Vector_t dir((end-vtx).Unit());

  // Input point, dir are always space charge corrected here
  float Apitch = GetPitch(geo, sce, end, dir, A.vhit_hit->View(), A.vhit_hit->WireID(), true, true);
  float Bpitch = GetPitch(geo, sce, end, dir, B.vhit_hit->View(), B.vhit_hit->WireID(), true, true);

  return abs(A.vhit->charge / Apitch - B.vhit->charge / Bpitch);
}

float sbn::StubTimeOffset	(	const sbn::StubInfo &	A,
		const sbn::StubInfo &	B,
		const detinfo::DetectorClocksData &	dclock,
		const detinfo::DetectorPropertiesData &	dprop
	)

Definition at line 137 of file StubMergeAlgorithms.cxx.

                                                {

  // Convert the ticks to an X-point to undo tick-offsets between planes. And then convert to ticks.
  return abs(dprop.ConvertTicksToX(A.vhit_hit->PeakTime(), A.vhit_hit->WireID()) - dprop.ConvertTicksToX(B.vhit_hit->PeakTime(), B.vhit_hit->WireID())) / dprop.DriftVelocity() / dclock.TPCClock().TickPeriod();
}

geo::Point_t sbn::TwoStubEndPosition	(	const sbn::StubInfo &	A,
		const sbn::StubInfo &	B,
		const geo::GeometryCore *	geo,
		const spacecharge::SpaceCharge *	sce,
		const detinfo::DetectorPropertiesData &	dprop
	)

Returns an updated end position of a stub after merging across two planes.

Definition at line 145 of file StubMergeAlgorithms.cxx.

                                                {

  static constexpr bool VERBOSE = false;

  // intersect the wires with the geometry service
  geo::Point_t yz;
  bool intersects = geo->WireIDsIntersect(A.vhit_hit->WireID(), B.vhit_hit->WireID(), yz);
  (void) intersects; // TODO: Do we care if the wires intersect? Means the end point will be outside the AV
  double y = yz.Y();
  double z = yz.Z();

  // just average the x-pos
  double x = (dprop.ConvertTicksToX(A.vhit_hit->PeakTime(), A.vhit_hit->WireID()) + dprop.ConvertTicksToX(B.vhit_hit->PeakTime(), B.vhit_hit->WireID())) / 2.;

  geo::Point_t pos(x, y, z);

  // Correct for space charge
  pos = GetLocation(sce, pos, A.vhit_hit->WireID());

  if (VERBOSE) {
    std::cout << "A View: " << geo->View(A.vhit_hit->WireID()) << " Angle to Vertical: " << geo->WireAngleToVertical(geo->View(A.vhit_hit->WireID()), A.vhit_hit->WireID()) << std::endl;
    std::cout << "B View: " << geo->View(B.vhit_hit->WireID()) << " Angle to Vertical: " << geo->WireAngleToVertical(geo->View(B.vhit_hit->WireID()), B.vhit_hit->WireID()) << std::endl;
    std::cout << "A End: " << A.stub.end.X() << " " << A.stub.end.Y() << " " << A.stub.end.Z() << std::endl;
    std::cout << "B End: " << B.stub.end.X() << " " << B.stub.end.Y() << " " << B.stub.end.Z() << std::endl;
    std::cout << "A Wire Pos Y: " << geo->Wire(A.vhit_hit->WireID()).GetCenter().Y() << " Z: " << geo->Wire(A.vhit_hit->WireID()).GetCenter().Z() << std::endl;
    std::cout << "B Wire Pos Y: " << geo->Wire(B.vhit_hit->WireID()).GetCenter().Y() << " Z: " << geo->Wire(B.vhit_hit->WireID()).GetCenter().Z() << std::endl;
    
    std::cout << "A Wire: " << A.vhit_hit->WireID().Wire << " Plane: " << A.vhit_hit->WireID().Plane << " TPC: " << A.vhit_hit->WireID().TPC << std::endl;
    std::cout << "B Wire: " << B.vhit_hit->WireID().Wire << " Plane: " << B.vhit_hit->WireID().Plane << " TPC: " << B.vhit_hit->WireID().TPC << std::endl;
    
    std::cout << "Output Y: " << y << " Z: " << z << std::endl;
    
    std::cout << "A Wire coord: " << geo->WireCoordinate(pos, A.vhit_hit->WireID()) << std::endl;
    std::cout << "B Wire coord: " << geo->WireCoordinate(pos, B.vhit_hit->WireID()) << std::endl;
    
    std::cout << "A X: " << dprop.ConvertTicksToX(A.vhit_hit->PeakTime(), A.vhit_hit->WireID()) << std::endl;
    std::cout << "B X: " << dprop.ConvertTicksToX(B.vhit_hit->PeakTime(), B.vhit_hit->WireID()) << std::endl;
    
    std::cout << "Output X: " << x << std::endl;
    
    
    std::cout << "Final x: " << pos.x() << " y: " << pos.y() << " z: " << pos.z() << std::endl;
  }

  return pos;
}

Variable Documentation

constexpr float sbn::default_float = std::numeric_limits<float>::signaling_NaN()

Definition at line 10 of file CRUMBSResult.h.

constexpr int sbn::default_int = std::numeric_limits<int>::max()

Definition at line 9 of file CRUMBSResult.h.