#include <TimeDet.h>

Inheritance diagram for TimeDet:

Collaboration diagram for TimeDet:

[legend]

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

	TimeDet ()

virtual	~TimeDet ()

virtual void	Initialize ()

virtual Bool_t	ProcessHits (FairVolume *v=0)

virtual void	Register ()

virtual TClonesArray *	GetCollection (Int_t iColl) const

virtual void	Reset ()

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 ()

TimeDetPoint *	AddHit (Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgCode, TVector3 Lpos, TVector3 Lmom)

virtual void	EndOfEvent ()

virtual void	FinishPrimary ()

virtual void	FinishRun ()

virtual void	BeginPrimary ()

virtual void	PostTrack ()

virtual void	PreTrack ()

virtual void	BeginEvent ()

Private Member Functions
	TimeDet (const TimeDet &)

TimeDet &	operator= (const TimeDet &)

Int_t	InitMedium (const char *name)

Private Attributes
Int_t	fTrackID

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

Double_t	fzPos
	energy loss

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.

TGeoVolume *	fDetector
	Overlap along y.

TClonesArray *	fTimeDetPointCollection

Detailed Description

Definition at line 14 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 73 of file TimeDet.cxx.

  : FairDetector(name, active, kTimeDet),
    fTrackID(-1),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fELoss(-1),
    //
    fzPos(0),
    fxSize(500),
    fySize(1000),
    fxBar(168),
    fyBar(6),
    fzBar(1),
    fdzBarCol(2.4),
    fdzBarRow(1.2),
    fNCol(3),
    fNRow(182),
    fxCenter(0),
    fyCenter(0),
    //
    fDetector(0),
    //
    fTimeDetPointCollection(new TClonesArray("TimeDetPoint"))
{
  fNBars = fNCol * fNRow;
  if(fNCol>1) fxOv = (fxBar*fNCol - fxSize) / (double)(fNCol-1); else fxOv = 0;
  if(fNRow>1) fyOv = (fyBar*fNRow - fySize) / (double)(fNRow-1); else fyOv = 0;
 
}

◆ TimeDet() [2/3]

TimeDet::TimeDet ( )

default constructor

Definition at line 39 of file TimeDet.cxx.

  : FairDetector("TimeDet", kTRUE, kTimeDet),
    fTrackID(-1),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fELoss(-1),
    //
    fzPos(0),
    fxSize(500),
    fySize(1000),
    fxBar(168),
    fyBar(6),
    fzBar(1),
    fdzBarCol(2.4),
    fdzBarRow(1.2),
    fNCol(3),
    fNRow(182),
    fxCenter(0),
    fyCenter(0),
    //
    fDetector(0),
    //
    fTimeDetPointCollection(new TClonesArray("TimeDetPoint"))
{
  fNBars = fNCol * fNRow;
  if(fNCol>1) fxOv = (fxBar*fNCol - fxSize) / (double)(fNCol-1); else fxOv = 0;
  if(fNRow>1) fyOv = (fyBar*fNRow - fySize) / (double)(fNRow-1); else fyOv = 0;
}

◆ ~TimeDet()

TimeDet::~TimeDet ( )

virtual

destructor

Definition at line 112 of file TimeDet.cxx.

{
  if (fTimeDetPointCollection) {
    fTimeDetPointCollection->Delete();
    delete fTimeDetPointCollection;
  }
}

◆ TimeDet() [3/3]

TimeDet::TimeDet ( const TimeDet & )

private

Member Function Documentation

◆ AddHit()

TimeDetPoint * TimeDet::AddHit	(	Int_t	trackID,
		Int_t	detID,
		TVector3	pos,
		TVector3	mom,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss,
		Int_t	pdgCode,
		TVector3	Lpos,
		TVector3	Lmom
	)

This method is an example of how to add your own point of type TimeDetPoint to the clones array

Definition at line 283 of file TimeDet.cxx.

{
  TClonesArray& clref = *fTimeDetPointCollection;
  Int_t size = clref.GetEntriesFast();
  // cout << "veto hit called "<< pos.z()<<endl;
  return new(clref[size]) TimeDetPoint(trackID, detID, pos, mom,
                 time, length, eLoss, pdgCode,Lpos,Lmom);
}

◆ BeginEvent()

virtual void TimeDet::BeginEvent ( )

inlinevirtual

Definition at line 81 of file TimeDet.h.

81{;}

◆ BeginPrimary()

virtual void TimeDet::BeginPrimary ( )

inlinevirtual

Definition at line 78 of file TimeDet.h.

78{;}

◆ ConstructGeometry()

void TimeDet::ConstructGeometry ( )

Create the detector geometry

Definition at line 244 of file TimeDet.cxx.

{
  TGeoVolume *top = gGeoManager->GetTopVolume();
  
  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;
}

◆ EndOfEvent()

void TimeDet::EndOfEvent ( )

virtual

Definition at line 207 of file TimeDet.cxx.

{
  fTimeDetPointCollection->Clear();
}

◆ FinishPrimary()

virtual void TimeDet::FinishPrimary ( )

inlinevirtual

Definition at line 76 of file TimeDet.h.

76{;}

◆ FinishRun()

virtual void TimeDet::FinishRun ( )

inlinevirtual

Definition at line 77 of file TimeDet.h.

77{;}

◆ GetBarRowCol()

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

Definition at line 296 of file TimeDet.cxx.

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

◆ GetCollection()

TClonesArray * TimeDet::GetCollection ( Int_t iColl ) const

virtual

Gets the produced collections

Definition at line 229 of file TimeDet.cxx.

{
  if (iColl == 0) { return fTimeDetPointCollection; }
  else { return NULL; }
}

◆ GetXCol()

double TimeDet::GetXCol ( int ic ) const

Definition at line 304 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 312 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 320 of file TimeDet.cxx.

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

◆ Initialize()

void TimeDet::Initialize ( )

virtual

Initialization of the detector is done here

Definition at line 106 of file TimeDet.cxx.

{
  FairDetector::Initialize();
}

◆ InitMedium()

Int_t TimeDet::InitMedium ( const char * name )

private

Definition at line 122 of file TimeDet.cxx.

{
  
   static FairGeoLoader *geoLoad=FairGeoLoader::Instance();
   static FairGeoInterface *geoFace=geoLoad->getGeoInterface();
   static FairGeoMedia *media=geoFace->getMedia();
   static FairGeoBuilder *geoBuild=geoLoad->getGeoBuilder();
 
   FairGeoMedium *ShipMedium=media->getMedium(name);
 
   if (!ShipMedium)
   {
     Fatal("InitMedium","Material %s not defined in media file.", name);
     return -1111;
   }
   TGeoMedium* medium=gGeoManager->GetMedium(name);
   if (medium!=NULL)
     return ShipMedium->getMediumIndex();
 
   return geoBuild->createMedium(ShipMedium);
  
  return 0;
}

◆ operator=()

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

private

◆ PostTrack()

virtual void TimeDet::PostTrack ( )

inlinevirtual

Definition at line 79 of file TimeDet.h.

79{;}

◆ PreTrack()

virtual void TimeDet::PreTrack ( )

inlinevirtual

Definition at line 80 of file TimeDet.h.

80{;}

◆ ProcessHits()

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

virtual

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

This method is called from the MC stepping

Definition at line 148 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();
 
    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(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 = (ShipStack*) gMC->GetStack();
    stack->AddPoint(kTimeDet);
  }
  
  return kTRUE;
}

◆ Register()

void TimeDet::Register ( )

virtual

Registers the produced collections in FAIRRootManager.

This will create a branch in the output tree called TimeDetPoint, setting the last parameter to kFALSE means: this collection will not be written to the file, it will exist only during the simulation.

Definition at line 214 of file TimeDet.cxx.

{
 
  FairRootManager::Instance()->Register("TimeDetPoint", "TimeDet",
                                        fTimeDetPointCollection, kTRUE);
}

◆ Reset()

void TimeDet::Reset ( )

virtual

has to be called after each event to reset the containers

Definition at line 237 of file TimeDet.cxx.

{
  fTimeDetPointCollection->Clear();
}

◆ SetBarZ()

void TimeDet::SetBarZ ( Double_t dz )

inline

Definition at line 55 of file TimeDet.h.

55{fzBar = dz;}

◆ SetBarZspacing()

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

inline

Definition at line 50 of file TimeDet.h.

    {
       fdzBarRow = row;
       fdzBarCol = column;
    }

◆ SetSizeX()

void TimeDet::SetSizeX ( Double_t x )

inline

Definition at line 56 of file TimeDet.h.

56{fxSize = x;}

◆ SetSizeY()

void TimeDet::SetSizeY ( Double_t y )

inline

Definition at line 57 of file TimeDet.h.

57{fySize = y;}

◆ SetZposition()

void TimeDet::SetZposition ( Double_t z )

inline

Sets detector position along z

Definition at line 49 of file TimeDet.h.

49{fzPos = z;}

Member Data Documentation

◆ fDetector

TGeoVolume* TimeDet::fDetector

private

Overlap along y.

Definition at line 118 of file TimeDet.h.

◆ fdzBarCol

Double_t TimeDet::fdzBarCol

private

depth of the bar

Definition at line 105 of file TimeDet.h.

◆ fdzBarRow

Double_t TimeDet::fdzBarRow

private

z-distance between columns

Definition at line 106 of file TimeDet.h.

◆ fELoss

Double_t TimeDet::fELoss

private

length

Definition at line 93 of file TimeDet.h.

◆ fLength

Double_t TimeDet::fLength

private

time

Definition at line 92 of file TimeDet.h.

◆ fMom

TLorentzVector TimeDet::fMom

private

position at entrance

Definition at line 90 of file TimeDet.h.

◆ fNBars

Int_t TimeDet::fNBars

private

y-position of the detector center

Definition at line 114 of file TimeDet.h.

◆ fNCol

Int_t TimeDet::fNCol

private

z-distance between rows

Definition at line 108 of file TimeDet.h.

◆ fNRow

Int_t TimeDet::fNRow

private

Number of columns.

Definition at line 109 of file TimeDet.h.

◆ fPos

TLorentzVector TimeDet::fPos

private

volume id

Definition at line 89 of file TimeDet.h.

◆ fTime

Double_t TimeDet::fTime

private

momentum at entrance

Definition at line 91 of file TimeDet.h.

◆ fTimeDetPointCollection

TClonesArray* TimeDet::fTimeDetPointCollection

private

container for data points

Definition at line 121 of file TimeDet.h.

◆ fTrackID

Int_t TimeDet::fTrackID

private

Track information to be stored until the track leaves the active volume.

Definition at line 87 of file TimeDet.h.

◆ fVolumeID

Int_t TimeDet::fVolumeID

private

track index

Definition at line 88 of file TimeDet.h.

◆ fxBar

Double_t TimeDet::fxBar

private

height of the detector

Definition at line 101 of file TimeDet.h.

◆ fxCenter

Double_t TimeDet::fxCenter

private

Number of rows.

Definition at line 111 of file TimeDet.h.

◆ fxOv

Double_t TimeDet::fxOv

private

Number of bars.

Definition at line 115 of file TimeDet.h.

◆ fxSize

Double_t TimeDet::fxSize

private

z-position of veto station

Definition at line 98 of file TimeDet.h.

◆ fyBar

Double_t TimeDet::fyBar

private

length of the bar

Definition at line 102 of file TimeDet.h.

◆ fyCenter

Double_t TimeDet::fyCenter

private

x-position of the detector center

Definition at line 112 of file TimeDet.h.

◆ fyOv

Double_t TimeDet::fyOv

private

Overlap along x.

Definition at line 116 of file TimeDet.h.

◆ fySize

Double_t TimeDet::fySize

private

width of the detector

Definition at line 99 of file TimeDet.h.

◆ fzBar

Double_t TimeDet::fzBar

private

width of the bar

Definition at line 103 of file TimeDet.h.

◆ fzPos

Double_t TimeDet::fzPos

private

energy loss

Detector parameters.

Definition at line 96 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

Detailed Description

Constructor & Destructor Documentation

◆ TimeDet() [1/3]

◆ TimeDet() [2/3]

◆ ~TimeDet()

◆ TimeDet() [3/3]

Member Function Documentation

◆ AddHit()

◆ BeginEvent()

◆ BeginPrimary()

◆ ConstructGeometry()

◆ EndOfEvent()

◆ FinishPrimary()

◆ FinishRun()

◆ GetBarRowCol()

◆ GetCollection()

◆ GetXCol()

◆ GetYRow()

◆ GetZBar()

◆ Initialize()

◆ InitMedium()

◆ operator=()

◆ PostTrack()

◆ PreTrack()

◆ ProcessHits()

◆ Register()

◆ Reset()

◆ SetBarZ()

◆ SetBarZspacing()

◆ SetSizeX()

◆ SetSizeY()

◆ SetZposition()

Member Data Documentation

◆ fDetector

◆ fdzBarCol

◆ fdzBarRow

◆ fELoss

◆ fLength

◆ fMom

◆ fNBars

◆ fNCol

◆ fNRow

◆ fPos

◆ fTime

◆ fTimeDetPointCollection

◆ fTrackID

◆ fVolumeID

◆ fxBar

◆ fxCenter

◆ fxOv

◆ fxSize

◆ fyBar

◆ fyCenter

◆ fyOv

◆ fySize

◆ fzBar

◆ fzPos