QwAnalysis
VQwBridgingMethod Class Referenceabstract

Interface to the various bridging methods. More...

#include <VQwBridgingMethod.h>

+ Inheritance diagram for VQwBridgingMethod:

Public Member Functions

 VQwBridgingMethod ()
 Default constructor. More...
 
virtual ~VQwBridgingMethod ()
 Destructor. More...
 
virtual const QwTrackBridge (const QwPartialTrack *front, const QwPartialTrack *back)=0
 Bridge from the front to back partial track (pure virtual) More...
 

Protected Member Functions

virtual double EstimateInitialMomentum (const TVector3 &direction) const
 Estimate the momentum based only on the direction. More...
 

Detailed Description

Interface to the various bridging methods.

Definition at line 31 of file VQwBridgingMethod.h.

Constructor & Destructor Documentation

VQwBridgingMethod::VQwBridgingMethod ( )
inline

Default constructor.

Definition at line 36 of file VQwBridgingMethod.h.

36 { };
virtual VQwBridgingMethod::~VQwBridgingMethod ( )
inlinevirtual

Destructor.

Definition at line 38 of file VQwBridgingMethod.h.

38 { };

Member Function Documentation

virtual const QwTrack* VQwBridgingMethod::Bridge ( const QwPartialTrack front,
const QwPartialTrack back 
)
pure virtual

Bridge from the front to back partial track (pure virtual)

Implemented in QwRayTracer, and QwMatrixLookup.

double VQwBridgingMethod::EstimateInitialMomentum ( const TVector3 &  direction) const
inlineprotectedvirtual

Estimate the momentum based only on the direction.

Estimate the momentum for a first attempt at bridging.

This takes into account only the direction of the front partial track and assumes an elastic electron-proton reaction without any energy loss. Any effects of the energy loss should be small enough that they can be recovered from the bridging algorithm in little time, and this avoids any bias of LH2 over Al events (E_loss).

Parameters
directionDirection of front partial track
Returns
Initial momentum

Definition at line 64 of file VQwBridgingMethod.h.

References Qw::GeV, and Qw::Mp.

65 {
66  double Mp = Qw::Mp; // proton mass
67  double e0 = 1.165 * Qw::GeV; // beam energy
68  // \todo Get the beam energy from the event (?)
69 
70  // Kinematics for elastic e+p scattering
71  double cth = direction.CosTheta();
72  return e0 / (1.0 + e0 / Mp * (1 - cth));
73 }
static const double GeV
Definition: QwUnits.h:95
static const double Mp
Mass of the proton.
Definition: QwUnits.h:119

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