#include <TimeDet.h>

Inheritance diagram for TimeDet:

[legend]

Collaboration diagram for TimeDet:

[legend]

Public Member Functions
	TimeDet (const char *Name, Bool_t Active)

	TimeDet ()

Bool_t	ProcessHits (FairVolume *v=0) override

void	SetZposition (Double_t z)

void	SetBarZspacing (Double_t row, Double_t column)

void	SetBarZ (Double_t dz)

void	SetSizeX (Double_t x)

void	SetSizeY (Double_t y)

double	GetXCol (int ic) const

double	GetYRow (int ir) const

void	GetBarRowCol (int ib, int &irow, int &icol) const

double	GetZBar (int ir, int ic) const

void	ConstructGeometry () override

Public Member Functions inherited from SHiP::Detector< TimeDetPoint >
	Detector ()=default

	Detector (const char *Name, Bool_t Active, Int_t detID)

	Detector (const char *Name, Bool_t Active)

	~Detector () override=default

TimeDetPoint *	AddHit (Args &&... args)

void	ConstructGeometry () override=0

void	Initialize () override

void	Reset () override

void	EndOfEvent () override

void	Register () override

TClonesArray *	GetCollection (Int_t iColl) const override

void	UpdatePointTrackIndices (const std::map< Int_t, Int_t > &indexMap) override
	Update track indices in point collection after track filtering.

void	SetSpecialPhysicsCuts () override

void	FinishPrimary () override

void	FinishRun () override

void	BeginPrimary () override

void	PostTrack () override

void	PreTrack () override

void	BeginEvent () override

void	CopyClones (TClonesArray cl1, TClonesArray cl2, Int_t offset) override

Public Member Functions inherited from ISTLPointContainer
virtual void	UpdatePointTrackIndices (const std::map< Int_t, Int_t > &indexMap)=0
	Update track indices in point collection after track filtering.

virtual	~ISTLPointContainer ()=default

Private Member Functions
	TimeDet (const TimeDet &)=delete
	Overlap along y.

TimeDet &	operator= (const TimeDet &)=delete

Private Attributes
Double_t	fzPos

Double_t	fxSize
	z-position of veto station

Double_t	fySize
	width of the detector

Double_t	fxBar
	height of the detector

Double_t	fyBar
	length of the bar

Double_t	fzBar
	width of the bar

Double_t	fdzBarCol
	depth of the bar

Double_t	fdzBarRow
	z-distance between columns

Int_t	fNCol
	z-distance between rows

Int_t	fNRow
	Number of columns.

Double_t	fxCenter
	Number of rows.

Double_t	fyCenter
	x-position of the detector center

Int_t	fNBars
	y-position of the detector center

Double_t	fxOv
	Number of bars.

Double_t	fyOv
	Overlap along x.

Additional Inherited Members
Protected Attributes inherited from SHiP::Detector< TimeDetPoint >
Int_t	fEventID

Int_t	fTrackID
	event index

Int_t	fVolumeID
	track index

TLorentzVector	fPos
	volume id

TLorentzVector	fMom
	position at entrance

Double_t	fTime
	momentum at entrance

Double_t	fLength
	time

Double_t	fELoss
	length

std::vector< TimeDetPoint > *	fDetPoints
	energy loss

TGeoVolume *	fDetector

Detailed Description

Definition at line 19 of file TimeDet.h.

Constructor & Destructor Documentation

◆ TimeDet() [1/3]

TimeDet::TimeDet	(	const char *	Name,
		Bool_t	Active
	)

Name : Detector Name Active: kTRUE for active detectors (ProcessHits() will be called) kFALSE for inactive detectors

Definition at line 67 of file TimeDet.cxx.

    : Detector(name, active, kTimeDet),
      //
      fzPos(0),
      fxSize(450),
      fySize(650),
      fxBar(168),
      fyBar(6),
      fzBar(1),
      fdzBarCol(2.4),
      fdzBarRow(1.2),
      fNCol(3),
      fNRow(148),
      fxCenter(0),
      fyCenter(0) {
  fNBars = fNCol * fNRow;
  if (fNCol > 1)
    fxOv = (fxBar * fNCol - fxSize) / static_cast<double>(fNCol - 1);
  else
    fxOv = 0;
  if (fNRow > 1)
    fyOv = (fyBar * fNRow - fySize) / static_cast<double>(fNRow - 1);
  else
    fyOv = 0;
}

◆ TimeDet() [2/3]

TimeDet::TimeDet ( )

default constructor

Definition at line 42 of file TimeDet.cxx.

    : Detector("TimeDet", kTRUE, kTimeDet),
      fzPos(0),
      fxSize(450),
      fySize(650),
      fxBar(168),
      fyBar(6),
      fzBar(1),
      fdzBarCol(2.4),
      fdzBarRow(1.2),
      fNCol(3),
      fNRow(148),
      fxCenter(0),
      fyCenter(0) {
  fNBars = fNCol * fNRow;
  if (fNCol > 1)
    fxOv = (fxBar * fNCol - fxSize) / static_cast<double>(fNCol - 1);
  else
    fxOv = 0;
  if (fNRow > 1)
    fyOv = (fyBar * fNRow - fySize) / static_cast<double>(fNRow - 1);
  else
    fyOv = 0;
}

◆ TimeDet() [3/3]

TimeDet::TimeDet ( const TimeDet & )

privatedelete

Overlap along y.

Member Function Documentation

◆ ConstructGeometry()

void TimeDet::ConstructGeometry ( )

overridevirtual

Create the detector geometry

Implements SHiP::Detector< TimeDetPoint >.

Definition at line 151 of file TimeDet.cxx.

                                {
  TGeoVolume* top = gGeoManager->GetTopVolume();
 
  ShipGeo::InitMedium("polyvinyltoluene");
  TGeoMedium* Scint = gGeoManager->GetMedium("polyvinyltoluene");
 
 
  fDetector = new TGeoVolumeAssembly("Timing Detector");
 
  TGeoVolume* plate =
      gGeoManager->MakeBox("TimeDet", Scint, fxBar / 2, fyBar / 2, fzBar / 2);
  plate->SetLineColor(kBlue);
  AddSensitiveVolume(plate);
 
  for (int ib = 0; ib < fNBars; ib++) {
    int irow = 0, icol = 0;
    GetBarRowCol(ib, irow, icol);
 
    double xbar = 0, ybar = 0, zbar = 0;
    xbar = GetXCol(icol);
    ybar = GetYRow(irow);
    zbar = GetZBar(irow, icol);
 
    fDetector->AddNode(plate, ib, new TGeoTranslation(xbar, ybar, zbar));
 
    // printf("%3i  %3i %2i   %8.3f %8.3f %8.3f\n",ib, irow,icol,
    // xbar,ybar,zbar);
  }
 
  top->AddNode(fDetector, 1, new TGeoTranslation(0, 0, fzPos));
 
 
  return;
}

◆ GetBarRowCol()

void TimeDet::GetBarRowCol	(	int	ib,
		int &	irow,
		int &	icol
	)		const

Definition at line 188 of file TimeDet.cxx.

                                                             {
  irow = ib / fNCol;
  icol = (ib % fNCol);
  return;
}

◆ GetXCol()

double TimeDet::GetXCol ( int ic ) const

Definition at line 194 of file TimeDet.cxx.

                                    {
  ic += 1;
  double x = fxBar * ic - fxOv * (ic - 1) - fxBar / 2;
  return x - fxSize / 2 + fxCenter;
}

◆ GetYRow()

double TimeDet::GetYRow ( int ir ) const

Definition at line 200 of file TimeDet.cxx.

                                    {
  ir += 1;
  double y = fyBar * ir - fyOv * (ir - 1) - fyBar / 2;
  return y - fySize / 2 + fyCenter;
}

◆ GetZBar()

double TimeDet::GetZBar	(	int	ir,
		int	ic
	)		const

Definition at line 206 of file TimeDet.cxx.

                                            {
  double z = (ir % 2) * fdzBarRow + (ic % 2) * fdzBarCol;
  return z;
}

◆ operator=()

TimeDet & TimeDet::operator= ( const TimeDet & )

privatedelete

◆ ProcessHits()

Bool_t TimeDet::ProcessHits ( FairVolume * v = 0 )

override

this method is called for each step during simulation (see FairMCApplication::Stepping())

This method is called from the MC stepping

Definition at line 93 of file TimeDet.cxx.

                                           {
  // Set parameters at entrance of volume. Reset ELoss.
  if (gMC->IsTrackEntering()) {
    fELoss = 0.;
    fTime = gMC->TrackTime() * 1.0e09;
    fLength = gMC->TrackLength();
    gMC->TrackPosition(fPos);
    gMC->TrackMomentum(fMom);
  }
 
  // Sum energy loss for all steps in the active volume
  fELoss += gMC->Edep();
 
  // Create vetoPoint at exit of active volume
  if (gMC->IsTrackExiting() || gMC->IsTrackStop() ||
      gMC->IsTrackDisappeared()) {
    if (fELoss == 0.) {
      return kFALSE;
    }
 
    fTrackID = gMC->GetStack()->GetCurrentTrackNumber();
    fEventID = gMC->CurrentEvent();
    Int_t uniqueId;
    gMC->CurrentVolID(uniqueId);
    if (uniqueId > 1000000)  // Solid scintillator case
    {
      Int_t vcpy;
      gMC->CurrentVolOffID(1, vcpy);
      if (vcpy == 5) uniqueId += 4;  // Copy of half
    }
 
    TParticle* p = gMC->GetStack()->GetCurrentTrack();
    Int_t pdgCode = p->GetPdgCode();
    TLorentzVector Pos;
    gMC->TrackPosition(Pos);
    TLorentzVector Mom;
    gMC->TrackMomentum(Mom);
    Double_t xmean = (fPos.X() + Pos.X()) / 2.;
    Double_t ymean = (fPos.Y() + Pos.Y()) / 2.;
    Double_t zmean = (fPos.Z() + Pos.Z()) / 2.;
 
    // cout << uniqueId << " :(" << xmean << ", " << ymean << ", " << zmean <<
    // "): " << gMC->CurrentVolName() << endl;
 
    AddHit(fEventID, fTrackID, uniqueId, TVector3(xmean, ymean, zmean),
           TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength, fELoss,
           pdgCode, TVector3(Pos.X(), Pos.Y(), Pos.Z()),
           TVector3(Mom.Px(), Mom.Py(), Mom.Pz()));
 
    // Increment number of veto det points in TParticle
    ShipStack* stack = dynamic_cast<ShipStack*>(gMC->GetStack());
    stack->AddPoint(kTimeDet);
  }
 
  return kTRUE;
}

◆ SetBarZ()

void TimeDet::SetBarZ ( Double_t dz )

inline

Definition at line 41 of file TimeDet.h.

41{ fzBar = dz; }

◆ SetBarZspacing()

void TimeDet::SetBarZspacing	(	Double_t	row,
		Double_t	column
	)

inline

Definition at line 37 of file TimeDet.h.

                                                     {
    fdzBarRow = row;
    fdzBarCol = column;
  }

◆ SetSizeX()

void TimeDet::SetSizeX ( Double_t x )

inline

Definition at line 42 of file TimeDet.h.

42{ fxSize = x; }

◆ SetSizeY()

void TimeDet::SetSizeY ( Double_t y )

inline

Definition at line 43 of file TimeDet.h.

43{ fySize = y; }

◆ SetZposition()

void TimeDet::SetZposition ( Double_t z )

inline

Sets detector position along z

Definition at line 36 of file TimeDet.h.

36{ fzPos = z; }

Member Data Documentation

◆ fdzBarCol

Double_t TimeDet::fdzBarCol

private

depth of the bar

Definition at line 64 of file TimeDet.h.

◆ fdzBarRow

Double_t TimeDet::fdzBarRow

private

z-distance between columns

Definition at line 65 of file TimeDet.h.

◆ fNBars

Int_t TimeDet::fNBars

private

y-position of the detector center

Definition at line 73 of file TimeDet.h.

◆ fNCol

Int_t TimeDet::fNCol

private

z-distance between rows

Definition at line 67 of file TimeDet.h.

◆ fNRow

Int_t TimeDet::fNRow

private

Number of columns.

Definition at line 68 of file TimeDet.h.

◆ fxBar

Double_t TimeDet::fxBar

private

height of the detector

Definition at line 60 of file TimeDet.h.

◆ fxCenter

Double_t TimeDet::fxCenter

private

Number of rows.

Definition at line 70 of file TimeDet.h.

◆ fxOv

Double_t TimeDet::fxOv

private

Number of bars.

Definition at line 74 of file TimeDet.h.

◆ fxSize

Double_t TimeDet::fxSize

private

z-position of veto station

Definition at line 57 of file TimeDet.h.

◆ fyBar

Double_t TimeDet::fyBar

private

length of the bar

Definition at line 61 of file TimeDet.h.

◆ fyCenter

Double_t TimeDet::fyCenter

private

x-position of the detector center

Definition at line 71 of file TimeDet.h.

◆ fyOv

Double_t TimeDet::fyOv

private

Overlap along x.

Definition at line 75 of file TimeDet.h.

◆ fySize

Double_t TimeDet::fySize

private

width of the detector

Definition at line 58 of file TimeDet.h.

◆ fzBar

Double_t TimeDet::fzBar

private

width of the bar

Definition at line 62 of file TimeDet.h.

◆ fzPos

Double_t TimeDet::fzPos

private

Detector parameters.

Definition at line 55 of file TimeDet.h.

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

TimeDet/TimeDet.h
TimeDet/TimeDet.cxx

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TimeDet() [1/3]

◆ TimeDet() [2/3]

◆ TimeDet() [3/3]

Member Function Documentation

◆ ConstructGeometry()

◆ GetBarRowCol()

◆ GetXCol()

◆ GetYRow()

◆ GetZBar()

◆ operator=()

◆ ProcessHits()

◆ SetBarZ()

◆ SetBarZspacing()

◆ SetSizeX()

◆ SetSizeY()

◆ SetZposition()

Member Data Documentation

◆ fdzBarCol

◆ fdzBarRow

◆ fNBars

◆ fNCol

◆ fNRow

◆ fxBar

◆ fxCenter

◆ fxOv

◆ fxSize

◆ fyBar

◆ fyCenter

◆ fyOv

◆ fySize

◆ fzBar

◆ fzPos