ShipMCTrack Class Reference

#include <ShipMCTrack.h>

Inheritance diagram for ShipMCTrack:

Collaboration diagram for ShipMCTrack:

Public Member Functions
	ShipMCTrack ()
	ShipMCTrack (Int_t pdgCode, Int_t motherID, Double_t px, Double_t py, Double_t pz, Double_t E, Double_t x, Double_t y, Double_t z, Double_t t, Int_t nPoints, Int_t eventID, Int_t trackID, Double_t w)
	ShipMCTrack (const ShipMCTrack &track)
	ShipMCTrack (TParticle *particle)
virtual	~ShipMCTrack ()
void	Print (Int_t iTrack=0) const
Int_t	GetPdgCode () const
Int_t	GetMotherId () const
Double_t	GetPx () const
Double_t	GetPy () const
Double_t	GetPz () const
Double_t	GetStartX () const
Double_t	GetStartY () const
Double_t	GetStartZ () const
Double_t	GetStartT () const
void	SetProcID (Int_t i)
Int_t	GetProcID () const
TString	GetProcName () const
Double_t	GetMass () const
Double_t	GetEnergy () const
Double_t	GetPt () const
Double_t	GetP () const
Double_t	GetRapidity () const
void	MultiplyWeight (Double_t w)
void	SetWeight (Double_t w)
Double_t	GetWeight () const
void	GetMomentum (TVector3 &momentum)
void	Get4Momentum (TLorentzVector &momentum)
void	GetStartVertex (TVector3 &vertex)
Int_t	GetNPoints (DetectorId detId) const
void	SetMotherId (Int_t id)
void	SetNPoints (Int_t iDet, Int_t np)
void	SetEventID (const Int_t &eventID)
void	SetTrackID (const Int_t &trackID)

Private Member Functions
	ClassDef (ShipMCTrack, 9)

Private Attributes
Int_t	fPdgCode
Int_t	fMotherId
Double32_t	fPx
Double32_t	fPy
Double32_t	fPz
Double32_t	fM
Double32_t	fStartX
Double32_t	fStartY
Double32_t	fStartZ
Double32_t	fStartT
Double32_t	fW
Int_t	fProcID
Int_t	fNPoints
Int_t	fEventID
Int_t	fTrackID

Detailed Description

Definition at line 28 of file ShipMCTrack.h.

Constructor & Destructor Documentation

ShipMCTrack() [1/4]

ShipMCTrack::ShipMCTrack

(

)

Default constructor

Definition at line 16 of file ShipMCTrack.cxx.

    : TObject(),
      fPdgCode(0),
      fMotherId(-1),
      fPx(0.),
      fPy(0.),
      fPz(0.),
      fM(-1.),
      fStartX(0.),
      fStartY(0.),
      fStartZ(0.),
      fStartT(0.),
      fW(1.),
      fProcID(44),
      fNPoints(0),
      fEventID(0),
      fTrackID(0) {}

ShipMCTrack() [2/4]

ShipMCTrack::ShipMCTrack ( Int_t  pdgCode, Int_t  motherID, Double_t  px, Double_t  py, Double_t  pz, Double_t  E, Double_t  x, Double_t  y, Double_t  z, Double_t  t, Int_t  nPoints, Int_t  eventID, Int_t  trackID, Double_t  w  )

explicit

Standard constructor

Definition at line 36 of file ShipMCTrack.cxx.

    : TObject(),
      fPdgCode(pdgCode),
      fMotherId(motherId),
      fPx(px),
      fPy(py),
      fPz(pz),
      fM(M),
      fStartX(x),
      fStartY(y),
      fStartZ(z),
      fStartT(t),
      fW(w),
      fProcID(44),
      fNPoints(nPoints),
      fEventID(eventID),
      fTrackID(trackID) {}

ShipMCTrack() [3/4]

ShipMCTrack::ShipMCTrack

(

const ShipMCTrack &

track

)

Copy constructor

Definition at line 59 of file ShipMCTrack.cxx.

    : TObject(track),
      fPdgCode(track.fPdgCode),
      fMotherId(track.fMotherId),
      fPx(track.fPx),
      fPy(track.fPy),
      fPz(track.fPz),
      fM(track.fM),
      fStartX(track.fStartX),
      fStartY(track.fStartY),
      fStartZ(track.fStartZ),
      fStartT(track.fStartT),
      fW(track.GetWeight()),
      fProcID(track.GetProcID()),
      fNPoints(track.fNPoints),
      fEventID(track.fEventID),
      fTrackID(track.fTrackID) {}

ShipMCTrack() [4/4]

ShipMCTrack::ShipMCTrack ( TParticle *  particle )

explicit

Constructor from TParticle

Definition at line 79 of file ShipMCTrack.cxx.

    : TObject(),
      fPdgCode(part->GetPdgCode()),
      fMotherId(part->GetMother(0)),
      fPx(part->Px()),
      fPy(part->Py()),
      fPz(part->Pz()),
      fM([](const TParticle* p) {
        Double_t m2 = p->Energy() * p->Energy() - p->P() * p->P();
        if (m2 >= 0.) return TMath::Sqrt(m2);
        // Tiny negative: numerical noise for on-shell massless particle
        Double_t e2 = p->Energy() * p->Energy();
        if (-m2 < 1e-10 * e2) return 0.;
        // Genuinely spacelike (virtual photon): store negative mass (ROOT
        // convention)
        return -TMath::Sqrt(-m2);
      }(part)),
      fStartX(part->Vx()),
      fStartY(part->Vy()),
      fStartZ(part->Vz()),
      fStartT(part->T() * 1e09),
      fW(part->GetWeight()),
      fProcID(part->GetUniqueID()),
      fNPoints(0),
      fEventID(0),
      fTrackID(0) {}

~ShipMCTrack()

ShipMCTrack::~ShipMCTrack ( )

virtual

Destructor

Definition at line 108 of file ShipMCTrack.cxx.

{}

Member Function Documentation

ClassDef()

ShipMCTrack::ClassDef ( ShipMCTrack  , 9    )

private

Get4Momentum()

void ShipMCTrack::Get4Momentum ( TLorentzVector &  momentum )

inline

Definition at line 139 of file ShipMCTrack.h.

                                                              {
  momentum.SetXYZT(fPx, fPy, fPz, GetEnergy());
}

GetEnergy()

Double_t ShipMCTrack::GetEnergy

(

)

const

Definition at line 124 of file ShipMCTrack.cxx.

                                      {
  if (fM < 0) {
    // older data, mass not made persistent
    Double_t mass = GetMass();
    return TMath::Sqrt(mass * mass + fPx * fPx + fPy * fPy + fPz * fPz);
  } else {
    return TMath::Sqrt(fM * fM + fPx * fPx + fPy * fPy + fPz * fPz);
  }
}

GetMass()

Double_t ShipMCTrack::GetMass

(

)

const

Definition at line 134 of file ShipMCTrack.cxx.

                                    {
  if (fM < 0) {
    // older data, mass not made persistent
    if (TDatabasePDG::Instance()) {
      TParticlePDG* particle = TDatabasePDG::Instance()->GetParticle(fPdgCode);
      if (particle) {
        return particle->Mass();
      } else {
        return 0.;
      }
    }
  }
  return fM;
}

GetMomentum()

void ShipMCTrack::GetMomentum ( TVector3 &  momentum )

inline

Definition at line 135 of file ShipMCTrack.h.

                                                       {
  momentum.SetXYZ(fPx, fPy, fPz);
}

GetMotherId()

Int_t ShipMCTrack::GetMotherId ( ) const

inline

Definition at line 54 of file ShipMCTrack.h.

{ return fMotherId; }

GetNPoints()

Int_t ShipMCTrack::GetNPoints

(

DetectorId

detId

)

const

Accessors to the number of MCPoints in the detectors

Definition at line 163 of file ShipMCTrack.cxx.

                                                    {
  /*  // TODO: Where does this come from
    if      ( detId == kREF  ) { return (  fNPoints &   1); }
    else if ( detId == kTutDet  ) { return ( (fNPoints & ( 7 <<  1) ) >>  1); }
    else if ( detId == kFairRutherford ) { return ( (fNPoints & (31 <<  4) ) >>
    4); } else { LOG(error) << "Unknown detector ID " << detId ; return 0;
    }
  */
  return 0;
}

GetP()

Double_t ShipMCTrack::GetP ( ) const

inline

Definition at line 68 of file ShipMCTrack.h.

                        {
    return TMath::Sqrt(fPx * fPx + fPy * fPy + fPz * fPz);
  }

GetPdgCode()

Int_t ShipMCTrack::GetPdgCode ( ) const

inline

Accessors

Definition at line 53 of file ShipMCTrack.h.

{ return fPdgCode; }

GetProcID()

Int_t ShipMCTrack::GetProcID ( ) const

inline

Definition at line 63 of file ShipMCTrack.h.

{ return fProcID; }

GetProcName()

TString ShipMCTrack::GetProcName ( ) const

inline

Definition at line 64 of file ShipMCTrack.h.

{ return TMCProcessName[fProcID]; }

GetPt()

Double_t ShipMCTrack::GetPt ( ) const

inline

Definition at line 67 of file ShipMCTrack.h.

{ return TMath::Sqrt(fPx * fPx + fPy * fPy); }

GetPx()

Double_t ShipMCTrack::GetPx ( ) const

inline

Definition at line 55 of file ShipMCTrack.h.

{ return fPx; }

GetPy()

Double_t ShipMCTrack::GetPy ( ) const

inline

Definition at line 56 of file ShipMCTrack.h.

{ return fPy; }

GetPz()

Double_t ShipMCTrack::GetPz ( ) const

inline

Definition at line 57 of file ShipMCTrack.h.

{ return fPz; }

GetRapidity()

Double_t ShipMCTrack::GetRapidity

(

)

const

Definition at line 155 of file ShipMCTrack.cxx.

                                        {
  Double_t e = GetEnergy();
  Double_t y = 0.5 * TMath::Log((e + fPz) / (e - fPz));
  return y;
}

GetStartT()

Double_t ShipMCTrack::GetStartT ( ) const

inline

Definition at line 61 of file ShipMCTrack.h.

{ return fStartT; }

GetStartVertex()

void ShipMCTrack::GetStartVertex ( TVector3 &  vertex )

inline

Definition at line 143 of file ShipMCTrack.h.

                                                        {
  vertex.SetXYZ(fStartX, fStartY, fStartZ);
}

GetStartX()

Double_t ShipMCTrack::GetStartX ( ) const

inline

Definition at line 58 of file ShipMCTrack.h.

{ return fStartX; }

GetStartY()

Double_t ShipMCTrack::GetStartY ( ) const

inline

Definition at line 59 of file ShipMCTrack.h.

{ return fStartY; }

GetStartZ()

Double_t ShipMCTrack::GetStartZ ( ) const

inline

Definition at line 60 of file ShipMCTrack.h.

{ return fStartZ; }

GetWeight()

Double_t ShipMCTrack::GetWeight

(

)

const

Definition at line 151 of file ShipMCTrack.cxx.

{ return fW; }

MultiplyWeight()

void ShipMCTrack::MultiplyWeight ( Double_t  w )

inline

Definition at line 72 of file ShipMCTrack.h.

{ fW = fW * w; }

Print()

void ShipMCTrack::Print

(

Int_t

iTrack = 0

)

const

Definition at line 112 of file ShipMCTrack.cxx.

                                           {
  LOG(debug) << "Track " << trackId << ", mother : " << fMotherId << ", Type "
             << fPdgCode << ", momentum (" << fPx << ", " << fPy << ", " << fPz
             << ") GeV";
  /* LOG(DEBUG2) << "       Ref " << GetNPoints(kREF)
               << ", TutDet " << GetNPoints(kTutDet)
               << ", Rutherford " << GetNPoints(kFairRutherford) ;
 */
}

SetEventID()

void ShipMCTrack::SetEventID

(

const Int_t &

eventID

)

Definition at line 200 of file ShipMCTrack.cxx.

{ fEventID = eventID; }

SetMotherId()

void ShipMCTrack::SetMotherId ( Int_t  id )

inline

Modifiers

Definition at line 83 of file ShipMCTrack.h.

{ fMotherId = id; }

SetNPoints()

void ShipMCTrack::SetNPoints	(	Int_t	iDet,
		Int_t	np
	)

Definition at line 176 of file ShipMCTrack.cxx.

                                                      {
  /*
    if ( iDet == kREF ) {
      if      ( nPoints < 0 ) { nPoints = 0; }
      else if ( nPoints > 1 ) { nPoints = 1; }
      fNPoints = ( fNPoints & ( ~ 1 ) )  |  nPoints;
    }
 
    else if ( iDet == kTutDet ) {
      if      ( nPoints < 0 ) { nPoints = 0; }
      else if ( nPoints > 7 ) { nPoints = 7; }
      fNPoints = ( fNPoints & ( ~ (  7 <<  1 ) ) )  |  ( nPoints <<  1 );
    }
 
    else if ( iDet == kFairRutherford ) {
      if      ( nPoints <  0 ) { nPoints =  0; }
      else if ( nPoints > 31 ) { nPoints = 31; }
      fNPoints = ( fNPoints & ( ~ ( 31 <<  4 ) ) )  |  ( nPoints <<  4 );
    }
 
    else LOG(error) << "Unknown detector ID "  << iDet ;
  */
}

SetProcID()

void ShipMCTrack::SetProcID ( Int_t  i )

inline

Definition at line 62 of file ShipMCTrack.h.

{ fProcID = i; }

SetTrackID()

void ShipMCTrack::SetTrackID

(

const Int_t &

trackID

)

Definition at line 202 of file ShipMCTrack.cxx.

{ fTrackID = trackID; }

SetWeight()

void ShipMCTrack::SetWeight ( Double_t  w )

inline

Definition at line 73 of file ShipMCTrack.h.

{ fW = w; }

Member Data Documentation

fEventID

Int_t ShipMCTrack::fEventID

private

Index of the event

Definition at line 125 of file ShipMCTrack.h.

fM

Double32_t ShipMCTrack::fM

private

Definition at line 96 of file ShipMCTrack.h.

fMotherId

Int_t ShipMCTrack::fMotherId

private

Index of mother track. -1 for primary particles.

Definition at line 93 of file ShipMCTrack.h.

fNPoints

Int_t ShipMCTrack::fNPoints

private

Bitvector representing the number of MCPoints for this track in each subdetector. The detectors are represented by REF: Bit 0 (1 bit, max. value 1) MVD: Bit 1 - 3 (3 bits, max. value 7) STS: Bit 4 - 8 (5 bits, max. value 31) RICH: Bit 9 (1 bit, max. value 1) MUCH: Bit 10 - 14 (5 bits, max. value 31) TRD: Bit 15 - 19 (5 bits, max. value 31) TOF: Bit 20 - 23 (4 bits, max. value 15) ECAL: Bit 24 (1 bit, max. value 1) ZDC: Bit 25 (1 bit, max. value 1) The respective point numbers can be accessed and modified with the inline functions. Bits 26-31 are spare for potential additional detectors.

Definition at line 122 of file ShipMCTrack.h.

fPdgCode

Int_t ShipMCTrack::fPdgCode

private

PDG particle code

Definition at line 90 of file ShipMCTrack.h.

fProcID

Int_t ShipMCTrack::fProcID

private

Geant4 process ID which created the particle

Definition at line 105 of file ShipMCTrack.h.

fPx

Double32_t ShipMCTrack::fPx

private

Momentum components at start vertex [GeV]

Definition at line 96 of file ShipMCTrack.h.

fPy

Double32_t ShipMCTrack::fPy

private

Definition at line 96 of file ShipMCTrack.h.

fPz

Double32_t ShipMCTrack::fPz

private

Definition at line 96 of file ShipMCTrack.h.

fStartT

Double32_t ShipMCTrack::fStartT

private

Definition at line 99 of file ShipMCTrack.h.

fStartX

Double32_t ShipMCTrack::fStartX

private

Coordinates of start vertex [cm, ns]

Definition at line 99 of file ShipMCTrack.h.

fStartY

Double32_t ShipMCTrack::fStartY

private

Definition at line 99 of file ShipMCTrack.h.

fStartZ

Double32_t ShipMCTrack::fStartZ

private

Definition at line 99 of file ShipMCTrack.h.

fTrackID

Int_t ShipMCTrack::fTrackID

private

Index of track in the event

Definition at line 128 of file ShipMCTrack.h.

fW

Double32_t ShipMCTrack::fW

private

weight

Definition at line 102 of file ShipMCTrack.h.

---

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

• shipdata/ShipMCTrack.h
• shipdata/ShipMCTrack.cxx