Contains a tracked event, i.e. all information from hits to tracks. More...

#include <QwEvent.h>

Inheritance diagram for QwEvent:

Collaboration diagram for QwEvent:

Public Member Functions
	QwEvent ()
	Default constructor. More...

virtual	~QwEvent ()
	Virtual destructor. More...

const QwEventHeader *	GetEventHeader () const

void	SetEventHeader (const QwEventHeader &eventheader)

void	SetEventHeader (const QwEventHeader *eventheader)

void	LoadBeamProperty (const TString &map)
	Load the beam properties from a map file. More...

double	EnergyLossHydrogen (const double vertex_z)
	Calculate the energy loss in the hydrogen target. More...

void	Clear (Option_t *option="")

void	Reset (Option_t *option="")

void	CalculateKinematics (const QwTrack *track)
	Calculate the kinematic variables for a given track. More...

void	Print (Option_t *option="") const
	Print the event. More...

const TVector3 &	GetVertexPosition () const

const TVector3 &	GetVertexMomentum () const

Double_t	GetScatteringAngle () const

Double_t	GetScatteringVertexZ () const

	ClassDef (QwEvent, 3)

Hit list maintenance for output to ROOT files
QwHit *	CreateNewHit ()
	Create a new hit. More...

void	AddHit (const QwHit *hit)
	Add an existing hit as a copy. More...

void	ClearHits (Option_t *option="")
	Clear the list of hits. More...

void	ResetHits (Option_t *option="")
	Reset the list of hits. More...

void	AddHitContainer (const QwHitContainer *hitlist)
	Add the hits in a hit container as a copy. More...

QwHitContainer *	GetHitContainer ()
	Get the list of hits as a hit container. More...

Int_t	GetNumberOfHits () const
	Get the number of hits. More...

const std::vector< QwHit * > &	GetListOfHits () const
	Get the list of hits. More...

const QwHit *	GetHit (const int hit) const
	Get the specified hit. More...

void	PrintHits (Option_t *option="") const
	Print the list of hits. More...

Tree line list maintenance for output to ROOT files
QwTreeLine *	CreateNewTreeLine ()
	Create a new tree line. More...

void	AddTreeLine (const QwTreeLine *treeline)
	Add an existing tree line as a copy. More...

void	AddTreeLineList (const QwTreeLine *treelinelist)
	Add a list of existing tree lines as a copy. More...

void	AddTreeLineList (const std::vector< QwTreeLine * > &treelinelist)
	Add a list of existing tree lines as a copy. More...

void	ClearTreeLines (Option_t *option="")
	Clear the list of tree lines. More...

void	ResetTreeLines (Option_t *option="")
	Reset the list of tree lines. More...

Int_t	GetNumberOfTreeLines () const
	Get the number of tree lines. More...

const std::vector< QwTreeLine * > &	GetListOfTreeLines () const
	Get the list of tree lines. More...

const std::vector< QwTreeLine * >	GetListOfTreeLines (EQwRegionID region, EQwDirectionID direction) const
	Get the list of tree lines in region and direction. More...

const QwTreeLine *	GetTreeLine (const int tl) const
	Get the specified tree line. More...

void	PrintTreeLines (Option_t *option="") const
	Print the list of tree lines. More...

Partial track list maintenance for output to ROOT files
QwPartialTrack *	CreateNewPartialTrack ()
	Create a new partial track. More...

void	AddPartialTrack (const QwPartialTrack *partialtrack)
	Add an existing partial track as a copy. More...

void	AddPartialTrackList (const std::vector< QwPartialTrack * > &partialtracklist)
	Add a list of existing partial tracks as a copy. More...

void	ClearPartialTracks (Option_t *option="")
	Clear the list of partial tracks. More...

void	ResetPartialTracks (Option_t *option="")
	Reset the list of partial tracks. More...

Int_t	GetNumberOfPartialTracks () const
	Get the number of partial tracks. More...

const std::vector < QwPartialTrack * >	GetListOfPartialTracks (EQwRegionID region, EQwDetectorPackage package) const
	Get the list of partial tracks in region and direction. More...

const std::vector < QwPartialTrack * > &	GetListOfPartialTracks () const
	Get the list of partial tracks. More...

const QwPartialTrack *	GetPartialTrack (const int pt) const
	Get the specified partial track. More...

void	PrintPartialTracks (Option_t *option="") const
	Print the list of partial tracks. More...

Track list maintenance for output to ROOT files
QwTrack *	CreateNewTrack ()
	Create a new track. More...

void	AddTrack (const QwTrack *track)
	Add an existing track as a copy. More...

void	AddTrackList (const std::vector< QwTrack * > &tracklist)
	Add a list of existing partial tracks as a copy. More...

void	ClearTracks (Option_t *option="")
	Clear the list of tracks. More...

void	ResetTracks (Option_t *option="")
	Reset the list of tracks. More...

Int_t	GetNumberOfTracks () const
	Get the number of tracks. More...

const std::vector< QwTrack * > &	GetListOfTracks () const
	Get the list of tracks. More...

const QwTrack *	GetTrack (const int t) const
	Get the specified track. More...

void	PrintTracks (Option_t *option="") const
	Print the list of tracks. More...

Public Member Functions inherited from QwObjectCounter< QwEvent >
	QwObjectCounter ()
	Default constructor. More...

	QwObjectCounter (const QwObjectCounter &)
	Copy constructor. More...

virtual	~QwObjectCounter ()
	Destructor. More...

Data Fields
QwEventHeader *	fEventHeader

Int_t	fNQwHits
	Number of QwHits in the array. More...

std::vector< QwHit * >	fQwHits
	Array of QwHits. More...

Int_t	fNQwTreeLines
	Number of QwTreeLines in the array. More...

std::vector< QwTreeLine * >	fQwTreeLines
	Array of QwTreeLines. More...

Int_t	fNQwPartialTracks
	Number of QwPartialTracks in the array. More...

std::vector< QwPartialTrack * >	fQwPartialTracks
	Array of QwPartialTracks. More...

Int_t	fNQwTracks
	Number of QwTracks in the array. More...

std::vector< QwTrack * >	fQwTracks
	Array of QwTracks. More...

Double_t	fPrimaryQ2
	< Momentum transfer Q^2 assuming elastic scattering with hydrogen energy loss More...

QwTreeLine *	fTreeLine [kNumPackages][kNumRegions][kNumDirections]

Main detector light yield
std::vector< float >	fMD_LeftNbOfPEs

std::vector< float >	fMD_RightNbOfPEs

std::vector< float >	fMD_TotalNbOfPEs

Generic kinematic variables
double	fHydrogenEnergyLoss
	Pre-scattering target energy loss assuming LH2 target. More...

double	fScatteringAngle
	Scattering angle. More...

double	fScatteringVertexZ
	Scattering vertex z position. More...

double	fScatteringVertexR
	Scattering vertex radial distance. More...

double	fCrossSection

double	fPreScatteringEnergy

double	fOriginVertexEnergy

TVector3	fVertexPosition
	Vertex position. More...

TVector3	fVertexMomentum
	Vertex momentum. More...

Kinematics under various assumptions
QwKinematics	fKin
	Inclusive scattering. More...

QwKinematics	fKinWithLoss
	Scattering with hydrogen energy loss. More...

QwKinematics	fKinElastic
	Scattering assuming elastic reaction. More...

QwKinematics	fKinElasticWithLoss
	Scattering assuming elastic reaction and hydrogen energy loss. More...

Static Public Attributes
static double	fBeamEnergy = 1.165 * Qw::GeV
	Electron beam energy. More...

Additional Inherited Members
Static Public Member Functions inherited from QwObjectCounter< QwEvent >
static size_t	GetObjectsCreated ()
	Get number of objects ever created. More...

static size_t	GetObjectsAlive ()
	Get number of objects still alive. More...

Detailed Description

Contains a tracked event, i.e. all information from hits to tracks.

A QwEvent contains all event information, from hits over partial track to complete tracks. It serves as the final product of the tracking code.

Definition at line 156 of file QwEvent.h.

Constructor & Destructor Documentation

QwEvent::QwEvent ( )

Default constructor.

Definition at line 20 of file QwEvent.cc.

References Clear(), fEventHeader, fNQwHits, fNQwPartialTracks, fNQwTracks, fNQwTreeLines, fTreeLine, kNumDirections, kNumPackages, and kNumRegions.

 : fScatteringAngle(0.0),fScatteringVertexZ(0.0),fScatteringVertexR(0.0),
   fCrossSection(-1.0),fPreScatteringEnergy(0.0),fOriginVertexEnergy(0.0),
   fPrimaryQ2(0.0)
 {
   // Reset the event header
   fEventHeader = 0;
 
   fNQwHits = 0;
   fNQwTreeLines = 0;
   fNQwPartialTracks = 0;
   fNQwTracks = 0;
 
   // Loop over all pointer objects
   for (int i = 0; i < kNumPackages; ++i) {
     for (int j = 0; j < kNumRegions; ++j) {
       // Null the treeline pointers
       for (int k = 0; k < kNumDirections; ++k) {
         fTreeLine[i][j][k] = 0;
       } // end of loop over directions
     } // end of loop over regions
   } // end of loop over packages
 
   // Clear the event
   Clear();
 }

Here is the call graph for this function:

QwEvent::~QwEvent ( )

virtual

Virtual destructor.

Definition at line 48 of file QwEvent.cc.

References Clear(), fEventHeader, fTreeLine, kNumDirections, kNumPackages, kNumRegions, and QwTreeLine::next.

 {
   // Clear the event
   Clear();
   
   // Loop over all allocated objects
   for (int i = 0; i < kNumPackages; ++i) {
 
     for (int j = 0; j < kNumRegions; ++j) {
       
         // Delete all those treelines
          for (int k = 0; k < kNumDirections; ++k) {
            QwTreeLine* tl = fTreeLine[i][j][k];
            while (tl) {
              QwTreeLine* tl_next = tl->next;
              delete tl;
              tl = tl_next;
            }
         } // end of loop over directions
     } // end of loop over regions
   } // end of loop over packages
 
  
   // Delete the event header
   if (fEventHeader) delete fEventHeader;
 }

Here is the call graph for this function:

Member Function Documentation

void QwEvent::AddHit ( const QwHit * hit )

Add an existing hit as a copy.

Definition at line 281 of file QwEvent.cc.

References QwLog::endl(), fNQwHits, fQwHits, and QwError.

Referenced by AddHitContainer(), and CreateNewHit().

 {
   if (hit) fQwHits.push_back(new QwHit(hit));
   else QwError << "trying to add null hit" << QwLog::endl;
   fNQwHits++;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::AddHitContainer ( const QwHitContainer * hitlist )

Add the hits in a hit container as a copy.

Definition at line 322 of file QwEvent.cc.

References AddHit().

Referenced by QwTreeEventBuffer::GetSpecificEvent().

 {
   for (QwHitContainer::const_iterator hit = hitlist->begin();
        hit != hitlist->end(); hit++) {
     const QwHit* p = &(*hit);
     AddHit(p);
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::AddPartialTrack ( const QwPartialTrack * partialtrack )

Add an existing partial track as a copy.

Definition at line 442 of file QwEvent.cc.

References fNQwPartialTracks, and fQwPartialTracks.

Referenced by AddPartialTrackList(), CreateNewPartialTrack(), and QwTreeEventBuffer::GetSpecificEvent().

 {
   fQwPartialTracks.push_back(new QwPartialTrack(partialtrack));
   fNQwPartialTracks++;
 }

Here is the caller graph for this function:

void QwEvent::AddPartialTrackList ( const std::vector< QwPartialTrack * > & partialtracklist )

Add a list of existing partial tracks as a copy.

Definition at line 449 of file QwEvent.cc.

References AddPartialTrack().

 {
   for (size_t i = 0; i < partialtracklist.size(); i++)
     if (partialtracklist[i]->IsValid())
       AddPartialTrack(partialtracklist[i]);
 }

Here is the call graph for this function:

void QwEvent::AddTrack ( const QwTrack * track )

Add an existing track as a copy.

Definition at line 515 of file QwEvent.cc.

References fNQwTracks, and fQwTracks.

Referenced by AddTrackList(), and CreateNewTrack().

 {
   fQwTracks.push_back(new QwTrack(track));
   ++fNQwTracks;
 }

Here is the caller graph for this function:

void QwEvent::AddTrackList ( const std::vector< QwTrack * > & tracklist )

Add a list of existing partial tracks as a copy.

Definition at line 522 of file QwEvent.cc.

References AddTrack().

 {
   for (size_t i = 0; i < tracklist.size(); i++)
     //if (tracklist[i]->IsValid()) // TODO
       AddTrack(tracklist[i]);
 }

Here is the call graph for this function:

void QwEvent::AddTreeLine ( const QwTreeLine * treeline )

Add an existing tree line as a copy.

Definition at line 354 of file QwEvent.cc.

References fNQwTreeLines, and fQwTreeLines.

Referenced by AddTreeLineList(), CreateNewTreeLine(), and QwTreeEventBuffer::GetSpecificEvent().

 {
   fQwTreeLines.push_back(new QwTreeLine(treeline));
   fNQwTreeLines++;
 }

Here is the caller graph for this function:

void QwEvent::AddTreeLineList ( const QwTreeLine * treelinelist )

Add a list of existing tree lines as a copy.

Definition at line 361 of file QwEvent.cc.

References AddTreeLine(), and QwTreeLine::next.

 {
   for (const QwTreeLine *treeline = treelinelist;
          treeline; treeline = treeline->next){
     if (treeline->IsValid()){
        AddTreeLine(treeline);
     }
   } 
 }

Here is the call graph for this function:

void QwEvent::AddTreeLineList ( const std::vector< QwTreeLine * > & treelinelist )

Add a list of existing tree lines as a copy.

Definition at line 372 of file QwEvent.cc.

References AddTreeLine().

 {
   for (size_t i = 0; i < treelinelist.size(); i++)
     if (treelinelist[i]->IsValid())
       AddTreeLine(treelinelist[i]);
 }

Here is the call graph for this function:

void QwEvent::CalculateKinematics ( const QwTrack * track )

Calculate the kinematic variables for a given track.

Calculate the kinematic variables for a given track

Parameters

track Reconstructed track

Definition at line 161 of file QwEvent.cc.

References Qw::cm, Qw::deg, QwLog::endl(), EnergyLossHydrogen(), fBeamEnergy, QwTrack::fFront, fHydrogenEnergyLoss, fKin, fKinElastic, fKinElasticWithLoss, fKinWithLoss, QwTrack::fMomentum, QwKinematics::fNu, QwKinematics::fP0, QwKinematics::fPp, fPrimaryQ2, QwKinematics::fQ2, fScatteringAngle, fScatteringVertexR, fScatteringVertexZ, fVertexMomentum, fVertexPosition, QwKinematics::fW2, QwKinematics::fX, QwKinematics::fY, QwPartialTrack::GetMomentumDirection(), QwPartialTrack::GetMomentumDirectionTheta(), QwPartialTrack::GetPosition(), QwPartialTrack::GetVertexZ(), Qw::GeV, Qw::GeV2, Qw::MeV, Qw::Mp, and QwWarning.

 {
   if (! track) {
     QwWarning << "QwEvent::CalculateKinematics: "
         << "called with null track pointer." << QwLog::endl;
     return;
   }
 
   if (! (track->fFront)) {
     QwWarning << "QwEvent::CalculateKinematics: "
         << "full track with null front track pointer." << QwLog::endl;
     return;
   }
 
   // Beam energy
   Double_t energy = fBeamEnergy;
 
   // Longitudinal vertex in target from front partial track
   Double_t vertex_z = track->fFront->GetVertexZ();
   fVertexPosition = track->fFront->GetPosition(vertex_z);
   fVertexMomentum = track->fFront->GetMomentumDirection() * track->fMomentum;
   fScatteringVertexZ = fVertexPosition.Z() / Qw::cm;
   fScatteringVertexR = fVertexPosition.Perp() / Qw::cm;
 
   // Scattering angle from front partial track
   Double_t theta = track->fFront->GetMomentumDirectionTheta();
   Double_t cos_theta = cos(theta);
   fScatteringAngle = theta / Qw::deg;
 
   // Prescattering energy loss
   Double_t pre_loss = EnergyLossHydrogen(vertex_z);
   fHydrogenEnergyLoss = pre_loss;
 
   // Generic scattering without energy loss
   Double_t P0 = energy;
   Double_t PP = track->fMomentum;
   Double_t Q2 = 2.0 * P0 * PP * (1 - cos_theta);
   fKin.fP0 = P0 / Qw::GeV;
   fKin.fPp = PP / Qw::GeV;
   fKin.fQ2 = Q2 / Qw::GeV2;
   fKin.fNu = (P0 - PP) / Qw::GeV;
   fKin.fW2 = (Qw::Mp * Qw::Mp + 2.0 * Qw::Mp * (P0 - PP) - Q2) / Qw::GeV2;
   fKin.fX = Q2 / (2.0 * Qw::Mp * (P0 - PP));
   fKin.fY = (P0 - PP) / P0;
 
   // Generic scattering with energy loss
   P0 = energy - pre_loss;
   PP = track->fMomentum;
   Q2 = 2.0 * P0 * PP * (1 - cos_theta);
   fKinWithLoss.fP0 = P0 / Qw::GeV;
   fKinWithLoss.fPp = PP / Qw::GeV;
   fKinWithLoss.fQ2 = Q2 / Qw::GeV2;
   fKinWithLoss.fNu = (P0 - PP) / Qw::GeV;
   fKinWithLoss.fW2 = (Qw::Mp * Qw::Mp + 2.0 * Qw::Mp * (P0 - PP) - Q2) / Qw::GeV2;
   fKinWithLoss.fX = Q2 / (2.0 * Qw::Mp * (P0 - PP));
   fKinWithLoss.fY = (P0 - PP) / P0;
 
   // Elastic scattering without energy loss
   P0 = energy * Qw::MeV;
   PP = Qw::Mp * P0 / (Qw::Mp + P0 * (1 - cos_theta));
   Q2 = 2.0 * P0 * PP * (1 - cos_theta);
   fKinElastic.fP0 = P0 / Qw::GeV;
   fKinElastic.fPp = PP / Qw::GeV;
   fKinElastic.fQ2 = Q2 / Qw::GeV2;
   fKinElastic.fNu = (P0 - PP) / Qw::GeV;
   fKinElastic.fW2 = (Qw::Mp * Qw::Mp + 2.0 * Qw::Mp * (P0 - PP) - Q2) / Qw::GeV2;
   fKinElastic.fX = Q2 / (2.0 * Qw::Mp * (P0 - PP));
   fKinElastic.fY = (P0 - PP) / P0;
 
   // Elastic scattering with energy loss
   P0 = (energy - pre_loss) * Qw::MeV;
   PP = Qw::Mp * P0 / (Qw::Mp + P0 * (1 - cos_theta));
   Q2 = 2.0 * P0 * PP * (1 - cos_theta);
   fKinElasticWithLoss.fP0 = P0 / Qw::GeV;
   fKinElasticWithLoss.fPp = PP / Qw::GeV;
   fKinElasticWithLoss.fQ2 = Q2 / Qw::GeV2;
   fKinElasticWithLoss.fNu = (P0 - PP) / Qw::GeV;
   fKinElasticWithLoss.fW2 = (Qw::Mp * Qw::Mp + 2.0 * Qw::Mp * (P0 - PP) - Q2) / Qw::GeV2;
   fKinElasticWithLoss.fX = Q2 / (2.0 * Qw::Mp * (P0 - PP));
   fKinElasticWithLoss.fY = (P0 - PP) / P0;
   
   fPrimaryQ2 = fKinElasticWithLoss.fQ2;
 }

Here is the call graph for this function:

QwEvent::ClassDef	(	QwEvent	,
		3
	)

void QwEvent::Clear ( Option_t * option = "" )

Definition at line 77 of file QwEvent.cc.

References ClearHits(), ClearPartialTracks(), ClearTracks(), and ClearTreeLines().

Referenced by QwEvent(), and ~QwEvent().

 {
   ClearHits(option);
   ClearTreeLines(option);
   ClearPartialTracks(option);
   ClearTracks(option);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::ClearHits ( Option_t * option = "" )

Clear the list of hits.

Definition at line 289 of file QwEvent.cc.

References fNQwHits, and fQwHits.

Referenced by Clear(), and ResetHits().

 {
   for (size_t i = 0; i < fQwHits.size(); i++)
     delete fQwHits.at(i);
   fQwHits.clear();
   fNQwHits = 0;
 }

Here is the caller graph for this function:

void QwEvent::ClearPartialTracks ( Option_t * option = "" )

Clear the list of partial tracks.

Definition at line 457 of file QwEvent.cc.

References fNQwPartialTracks, and fQwPartialTracks.

Referenced by Clear(), and ResetPartialTracks().

 {
   for (size_t i = 0; i < fQwPartialTracks.size(); i++)
     delete fQwPartialTracks.at(i);
   fQwPartialTracks.clear();
   fNQwPartialTracks = 0;
 }

Here is the caller graph for this function:

void QwEvent::ClearTracks ( Option_t * option = "" )

Clear the list of tracks.

Definition at line 530 of file QwEvent.cc.

References fNQwTracks, and fQwTracks.

Referenced by Clear(), and ResetTracks().

 {
   for (size_t i = 0; i < fQwTracks.size(); i++)
     delete fQwTracks.at(i);
   fQwTracks.clear();
   fNQwTracks = 0;
 }

Here is the caller graph for this function:

void QwEvent::ClearTreeLines ( Option_t * option = "" )

Clear the list of tree lines.

Definition at line 380 of file QwEvent.cc.

References fNQwTreeLines, and fQwTreeLines.

Referenced by Clear(), and ResetTreeLines().

 {
   for (size_t i = 0; i < fQwTreeLines.size(); i++) {
     QwTreeLine* tl = fQwTreeLines.at(i);
     delete tl;
   }
   fQwTreeLines.clear();
   fNQwTreeLines = 0;
 }

Here is the caller graph for this function:

QwHit * QwEvent::CreateNewHit ( )

Create a new hit.

Definition at line 273 of file QwEvent.cc.

References AddHit().

 {
   QwHit* hit = new QwHit();
   AddHit(hit);
   return hit;
 }

Here is the call graph for this function:

QwPartialTrack * QwEvent::CreateNewPartialTrack ( )

Create a new partial track.

Definition at line 434 of file QwEvent.cc.

References AddPartialTrack().

 {
   QwPartialTrack* partialtrack = new QwPartialTrack();
   AddPartialTrack(partialtrack);
   return partialtrack;
 }

Here is the call graph for this function:

QwTrack * QwEvent::CreateNewTrack ( )

Create a new track.

Definition at line 507 of file QwEvent.cc.

References AddTrack().

 {
   QwTrack* track = new QwTrack();
   AddTrack(track);
   return track;
 }

Here is the call graph for this function:

QwTreeLine * QwEvent::CreateNewTreeLine ( )

Create a new tree line.

Definition at line 346 of file QwEvent.cc.

References AddTreeLine().

 {
   QwTreeLine* treeline = new QwTreeLine();
   AddTreeLine(treeline);
   return treeline;
 }

Here is the call graph for this function:

double QwEvent::EnergyLossHydrogen ( const double vertex_z )

Calculate the energy loss in the hydrogen target.

Calculate the energy loss in the hydrogen target

Parameters

vertex_z Longitudinal position of vertex (absolute coordinates)

Returns: Energy loss

Definition at line 127 of file QwEvent.cc.

References Qw::MeV.

Referenced by CalculateKinematics().

 {
   // Target parameters
   double target_z_length = 34.35;       // Target Length (cm)
   double target_z_position= -652.67;    // Target center position (cm) in Z
   double target_z_front = target_z_position - 0.5 * target_z_length;
   double target_z_back  = target_z_position + 0.5 * target_z_length;
 
   // Energy loss in hydrogen target
   double pre_loss = 0.0;
   if (vertex_z < target_z_front) {
     // Before target
     pre_loss = 0.05 * Qw::MeV;
 
   } else if (vertex_z >= target_z_front && vertex_z <= target_z_back) {
     // Inside target
     double depth = vertex_z - target_z_front;
     //pre_loss = 20.08 * Qw::MeV * depth/target_z_length; // linear approximation
     pre_loss = 0.05 + depth*0.6618 - 0.003462*depth*depth; // quadratic fit approximation
 
   } else {
     // After target
     pre_loss = 18.7 * Qw::MeV;  // averaged total Eloss, full target (excluding downstream Al window)
   }
 
   return pre_loss;
 }

Here is the caller graph for this function:

const QwEventHeader* QwEvent::GetEventHeader ( ) const

inline

Definition at line 192 of file QwEvent.h.

References fEventHeader.

192 { return fEventHeader; };

QwEvent::fEventHeader

QwEventHeader * fEventHeader

Definition: QwEvent.h:161

const QwHit * QwEvent::GetHit ( const int hit ) const

Get the specified hit.

Definition at line 304 of file QwEvent.cc.

References fQwHits, and GetNumberOfHits().

 {
   if (hit < 0 || hit > GetNumberOfHits()) return 0;
   return fQwHits.at(hit);
 }

Here is the call graph for this function:

QwHitContainer * QwEvent::GetHitContainer ( )

Get the list of hits as a hit container.

Definition at line 332 of file QwEvent.cc.

References QwLog::endl(), fQwHits, and QwError.

Referenced by QwTreeEventBuffer::GetHitContainer().

 {
   QwHitContainer* hitlist = new QwHitContainer();
   for (std::vector<QwHit*>::const_iterator hit = fQwHits.begin();
       hit != fQwHits.end(); hit++)
     if (*hit)
       hitlist->push_back(*hit);
     else
       QwError << "null hit in hit list" << QwLog::endl;
   return hitlist;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

const std::vector<QwHit*>& QwEvent::GetListOfHits ( ) const

inline

Get the list of hits.

Definition at line 232 of file QwEvent.h.

References fQwHits.

232 { return fQwHits; };

QwEvent::fQwHits

std::vector< QwHit * > fQwHits

Array of QwHits.

Definition: QwEvent.h:169

const std::vector< QwPartialTrack * > QwEvent::GetListOfPartialTracks	(	EQwRegionID	region,
		EQwDetectorPackage	package
	)		const

Get the list of partial tracks in region and direction.

Definition at line 488 of file QwEvent.cc.

References fQwPartialTracks.

 {
   std::vector<QwPartialTrack*> partialtracks;
   for (std::vector<QwPartialTrack*>::const_iterator
       pt  = fQwPartialTracks.begin();
       pt != fQwPartialTracks.end(); pt++) {
     if ((*pt)->GetRegion() == region) {
       if ((*pt)->GetPackage() == package) {
         partialtracks.push_back(*pt);
       }
     }
   }
   return partialtracks;
 }

const std::vector<QwPartialTrack*>& QwEvent::GetListOfPartialTracks ( ) const

inline

Get the list of partial tracks.

Definition at line 282 of file QwEvent.h.

References fQwPartialTracks.

282 { return fQwPartialTracks; };

QwEvent::fQwPartialTracks

std::vector< QwPartialTrack * > fQwPartialTracks

Array of QwPartialTracks.

Definition: QwEvent.h:177

const std::vector<QwTrack*>& QwEvent::GetListOfTracks ( ) const

inline

Get the list of tracks.

Definition at line 304 of file QwEvent.h.

References fQwTracks.

304 { return fQwTracks; };

QwEvent::fQwTracks

std::vector< QwTrack * > fQwTracks

Array of QwTracks.

Definition: QwEvent.h:181

const std::vector<QwTreeLine*>& QwEvent::GetListOfTreeLines ( ) const

inline

Get the list of tree lines.

Definition at line 256 of file QwEvent.h.

References fQwTreeLines.

256 { return fQwTreeLines; };

QwEvent::fQwTreeLines

std::vector< QwTreeLine * > fQwTreeLines

Array of QwTreeLines.

Definition: QwEvent.h:173

const std::vector< QwTreeLine * > QwEvent::GetListOfTreeLines	(	EQwRegionID	region,
		EQwDirectionID	direction
	)		const

Get the list of tree lines in region and direction.

Definition at line 417 of file QwEvent.cc.

References fQwTreeLines.

 {
   std::vector<QwTreeLine*> treelines;
   for (std::vector<QwTreeLine*>::const_iterator
       tl  = fQwTreeLines.begin();
       tl != fQwTreeLines.end(); tl++) {
     if ((*tl)->GetRegion() == region) {
       if ((*tl)->GetDirection() == direction) {
         treelines.push_back(*tl);
       }
     }
   }
   return treelines;
 }

Int_t QwEvent::GetNumberOfHits ( ) const

inline

Get the number of hits.

Definition at line 230 of file QwEvent.h.

References fNQwHits.

Referenced by GetHit().

230 { return fNQwHits; };

QwEvent::fNQwHits

Int_t fNQwHits

Number of QwHits in the array.

Definition: QwEvent.h:168

Here is the caller graph for this function:

Int_t QwEvent::GetNumberOfPartialTracks ( ) const

inline

Get the number of partial tracks.

Definition at line 278 of file QwEvent.h.

References fNQwPartialTracks.

Referenced by GetPartialTrack().

278 { return fNQwPartialTracks; };

QwEvent::fNQwPartialTracks

Int_t fNQwPartialTracks

Number of QwPartialTracks in the array.

Definition: QwEvent.h:176

Here is the caller graph for this function:

Int_t QwEvent::GetNumberOfTracks ( ) const

inline

Get the number of tracks.

Definition at line 302 of file QwEvent.h.

References fNQwTracks.

Referenced by GetTrack().

302 { return fNQwTracks; };

QwEvent::fNQwTracks

Int_t fNQwTracks

Number of QwTracks in the array.

Definition: QwEvent.h:180

Here is the caller graph for this function:

Int_t QwEvent::GetNumberOfTreeLines ( ) const

inline

Get the number of tree lines.

Definition at line 254 of file QwEvent.h.

References fNQwTreeLines.

Referenced by GetTreeLine(), and QwTrackingWorker::ProcessEvent().

254 { return fNQwTreeLines; };

QwEvent::fNQwTreeLines

Int_t fNQwTreeLines

Number of QwTreeLines in the array.

Definition: QwEvent.h:172

Here is the caller graph for this function:

const QwPartialTrack * QwEvent::GetPartialTrack ( const int pt ) const

Get the specified partial track.

Definition at line 472 of file QwEvent.cc.

References fQwPartialTracks, and GetNumberOfPartialTracks().

 {
   if (pt < 0 || pt > GetNumberOfPartialTracks()) return 0;
   return fQwPartialTracks.at(pt);
 }

Here is the call graph for this function:

Double_t QwEvent::GetScatteringAngle ( ) const

inline

Definition at line 319 of file QwEvent.h.

References fScatteringAngle.

319 { return fScatteringAngle;};

QwEvent::fScatteringAngle

double fScatteringAngle

Scattering angle.

Definition: QwEvent.h:336

Double_t QwEvent::GetScatteringVertexZ ( ) const

inline

Definition at line 320 of file QwEvent.h.

References fScatteringVertexZ.

320 { return fScatteringVertexZ;};

QwEvent::fScatteringVertexZ

double fScatteringVertexZ

Scattering vertex z position.

Definition: QwEvent.h:337

const QwTrack * QwEvent::GetTrack ( const int t ) const

Get the specified track.

Definition at line 545 of file QwEvent.cc.

References fQwTracks, and GetNumberOfTracks().

Referenced by DrawTreeLineVDC(), and QwTrackingWorker::ProcessEvent().

 {
   if (t < 0 || t > GetNumberOfTracks()) return 0;
   return fQwTracks.at(t);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

const QwTreeLine * QwEvent::GetTreeLine ( const int tl ) const

Get the specified tree line.

Definition at line 397 of file QwEvent.cc.

References fQwTreeLines, and GetNumberOfTreeLines().

 {
   if (tl < 0 || tl > GetNumberOfTreeLines()) return 0;
   return fQwTreeLines.at(tl);
 }

Here is the call graph for this function:

const TVector3& QwEvent::GetVertexMomentum ( ) const

inline

Definition at line 318 of file QwEvent.h.

References fVertexMomentum.

318 { return fVertexMomentum;};

QwEvent::fVertexMomentum

TVector3 fVertexMomentum

Vertex momentum.

Definition: QwEvent.h:343

const TVector3& QwEvent::GetVertexPosition ( ) const

inline

Definition at line 317 of file QwEvent.h.

References fVertexPosition.

317 { return fVertexPosition;};

QwEvent::fVertexPosition

TVector3 fVertexPosition

Vertex position.

Definition: QwEvent.h:342

void QwEvent::LoadBeamProperty ( const TString & map )

Load the beam properties from a map file.

Parameters

map	Name of map file

Definition at line 99 of file QwEvent.cc.

References QwLog::endl(), fBeamEnergy, Qw::GeV, Qw::MeV, QwMessage, and QwParameterFile::TrimComment().

 {
   QwParameterFile mapstr(map.Data());
   while (mapstr.ReadNextLine())
   {
     mapstr.TrimComment();       // Remove everything after a comment character.
     mapstr.TrimWhitespace();    // Get rid of leading and trailing spaces.
     if (mapstr.LineIsEmpty())  continue;
 
     TString varname, varvalue;
     if (mapstr.HasVariablePair("=",varname,varvalue)) {
       //  This is a declaration line.  Decode it.
       varname.ToLower();
       if (varname == "energy") {
         fBeamEnergy = atof(varvalue.Data()) * Qw::MeV;
         QwMessage << "Beam energy set to " << fBeamEnergy / Qw::GeV << " GeV" << QwLog::endl;
       }
     }
   }
 }

Here is the call graph for this function:

void QwEvent::Print ( Option_t * option = "" ) const

Print the event.

Definition at line 246 of file QwEvent.cc.

References QwLog::endl(), fEventHeader, fKin, QwKinematics::fP0, QwKinematics::fPp, QwKinematics::fQ2, PrintHits(), PrintPartialTracks(), PrintTracks(), PrintTreeLines(), and QwOut.

 {
   // Event header
   QwOut << *fEventHeader << QwLog::endl;
   // Event kinematics
   QwOut << "P0 = " << fKin.fP0 << " GeV/c" << QwLog::endl;
   QwOut << "PP = " << fKin.fPp << " GeV/c" << QwLog::endl;
   QwOut << "Q^2 = " << fKin.fQ2 << " (GeV/c)^2" << QwLog::endl;
   //std::cout << "weight = " << fCrossSectionWeight << std::endl;
   //std::cout << "energy = " << fTotalEnergy/Qw::MeV << " MeV" << std::endl;
   //std::cout << "momentum = " << fMomentum / Qw::MeV << " MeV" << std::endl;
   //std::cout << "vertex position = " << fVertexPosition.Z()/Qw::cm << " cm" << std::endl;
   //std::cout << "vertex momentum = " << fVertexMomentum.Z()/Qw::MeV << " MeV" << std::endl;
 
   // Event content
   QwOut << "Hits in this event:" << QwLog::endl;
   PrintHits(option);
   QwOut << "Tree lines in this event:" << QwLog::endl;
   PrintTreeLines(option);
   QwOut << "Partial tracks in this event:" << QwLog::endl;
   PrintPartialTracks(option);
   QwOut << "Tracks in this event:" << QwLog::endl;
   PrintTracks(option);
 }

Here is the call graph for this function:

void QwEvent::PrintHits ( Option_t * option = "" ) const

Print the list of hits.

Definition at line 311 of file QwEvent.cc.

References fQwHits.

Referenced by Print().

 {
   for (std::vector<QwHit*>::const_iterator hit = fQwHits.begin();
       hit != fQwHits.end(); ++hit) {
     if ((*hit)->GetDriftDistance() != -5)
       std::cout << **hit << std::endl;
   }
 }

Here is the caller graph for this function:

void QwEvent::PrintPartialTracks ( Option_t * option = "" ) const

Print the list of partial tracks.

Definition at line 479 of file QwEvent.cc.

References fQwPartialTracks.

Referenced by Print().

 {
   for (std::vector<QwPartialTrack*>::const_iterator parttrack =
       fQwPartialTracks.begin(); parttrack != fQwPartialTracks.end();
       ++parttrack)
     std::cout << **parttrack << std::endl;
 }

Here is the caller graph for this function:

void QwEvent::PrintTracks ( Option_t * option = "" ) const

Print the list of tracks.

Definition at line 552 of file QwEvent.cc.

References fQwTracks.

Referenced by Print().

 {
   for (std::vector<QwTrack*>::const_iterator track = fQwTracks.begin();
       track != fQwTracks.end(); track++)
     std::cout << **track << std::endl;
 }

Here is the caller graph for this function:

void QwEvent::PrintTreeLines ( Option_t * option = "" ) const

Print the list of tree lines.

Definition at line 404 of file QwEvent.cc.

References fQwTreeLines, and QwTreeLine::next.

Referenced by Print().

 {
   for (std::vector<QwTreeLine*>::const_iterator treeline = fQwTreeLines.begin();
       treeline != fQwTreeLines.end(); treeline++) {
     std::cout << **treeline << std::endl;
     QwTreeLine* tl = (*treeline)->next;
     while (tl) {
       std::cout << *tl << std::endl;
       tl = tl->next;
     }
   }
 }

Here is the caller graph for this function:

void QwEvent::Reset ( Option_t * option = "" )

Definition at line 86 of file QwEvent.cc.

References ResetHits(), ResetPartialTracks(), ResetTracks(), and ResetTreeLines().

Referenced by main().

 {
   ResetHits(option);
   ResetTreeLines(option);
   ResetPartialTracks(option);
   ResetTracks(option);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::ResetHits ( Option_t * option = "" )

Reset the list of hits.

Definition at line 298 of file QwEvent.cc.

References ClearHits().

Referenced by Reset().

 {
   ClearHits();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::ResetPartialTracks ( Option_t * option = "" )

Reset the list of partial tracks.

Definition at line 466 of file QwEvent.cc.

References ClearPartialTracks().

Referenced by Reset().

 {
   ClearPartialTracks();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::ResetTracks ( Option_t * option = "" )

Reset the list of tracks.

Definition at line 539 of file QwEvent.cc.

References ClearTracks().

Referenced by Reset().

 {
   ClearTracks();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::ResetTreeLines ( Option_t * option = "" )

Reset the list of tree lines.

Definition at line 391 of file QwEvent.cc.

References ClearTreeLines().

Referenced by Reset().

 {
   ClearTreeLines();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwEvent::SetEventHeader ( const QwEventHeader & eventheader )

inline

Definition at line 193 of file QwEvent.h.

References fEventHeader.

Referenced by QwTreeEventBuffer::GetSpecificEvent().

                                                           {
       if (fEventHeader) delete fEventHeader;
       fEventHeader = new QwEventHeader(eventheader);
     };

Here is the caller graph for this function:

void QwEvent::SetEventHeader ( const QwEventHeader * eventheader )

inline

Definition at line 197 of file QwEvent.h.

References fEventHeader.

                                                           {
       if (fEventHeader) delete fEventHeader;
       fEventHeader = new QwEventHeader(*eventheader);
     };