#include <SiliconTarget.h>

Inheritance diagram for SiliconTarget:

[legend]

Collaboration diagram for SiliconTarget:

[legend]

Public Member Functions
	SiliconTarget (const char name, Bool_t Active, const char Title="")

	SiliconTarget ()

void	SetSiliconTargetParameters (Double_t targetWidth, Double_t targetHeight, Double_t sensorWidth, Double_t sensorLength, Int_t nLayers, Double_t zPosition, Double_t targetThickness, Double_t targetSpacing, Double_t moduleOffset)

TGeoVolume *	CreateSiliconPlanes (const char name, Double_t sensorWidth, Double_t sensorLength, Double_t planeSpacing, TGeoMedium material, Int_t layerId)

void	ConstructGeometry () override

Bool_t	ProcessHits (FairVolume *vol=0) override

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

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

	Detector (const char *Name, Bool_t Active)

	~Detector () override=default

SiliconTargetPoint *	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
	SiliconTarget (const SiliconTarget &)=delete

SiliconTarget &	operator= (const SiliconTarget &)=delete

Private Attributes
Double_t	fTargetWidth

Double_t	fTargetHeight

Double_t	fSensorWidth

Double_t	fSensorLength

Int_t	fLayers

Double_t	fZPosition

Double_t	fTargetThickness

Double_t	fTargetSpacing

Double_t	fModuleOffset

Additional Inherited Members
Protected Attributes inherited from SHiP::Detector< SiliconTargetPoint >
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< SiliconTargetPoint > *	fDetPoints
	energy loss

TGeoVolume *	fDetector

Detailed Description

Definition at line 14 of file SiliconTarget.h.

Constructor & Destructor Documentation

◆ SiliconTarget() [1/3]

SiliconTarget::SiliconTarget	(	const char *	name,
		Bool_t	Active,
		const char *	Title = `""`
	)

Definition at line 31 of file SiliconTarget.cxx.

33 : Detector(name, Active, kSiliconTarget) {}

kSiliconTarget

@ kSiliconTarget

Definition: ShipDetectorList.h:22

SHiP::Detector< SiliconTargetPoint >::Detector

Detector()=default

◆ SiliconTarget() [2/3]

SiliconTarget::SiliconTarget ( )

Definition at line 28 of file SiliconTarget.cxx.

29 : Detector("SiliconTarget", kTRUE, kSiliconTarget) {}

◆ SiliconTarget() [3/3]

SiliconTarget::SiliconTarget ( const SiliconTarget & )

privatedelete

Member Function Documentation

◆ ConstructGeometry()

void SiliconTarget::ConstructGeometry ( )

overridevirtual

Create the detector geometry

Implements SHiP::Detector< SiliconTargetPoint >.

Definition at line 103 of file SiliconTarget.cxx.

                                      {
  ShipGeo::InitMedium("tungstensifon");
  TGeoMedium* tungsten = gGeoManager->GetMedium("tungstensifon");
  ShipGeo::InitMedium("air");
  TGeoMedium* air = gGeoManager->GetMedium("air");
  ShipGeo::InitMedium("silicon");
  TGeoMedium* Silicon = gGeoManager->GetMedium("silicon");
 
  Double_t totalLength = fLayers * fTargetSpacing;
 
  // --- Create an envelope volume for the detector (green, semi-transparent)
  // ---
  auto envBox = new TGeoBBox("SiliconTarget_env", fTargetWidth / 2.,
                             fTargetHeight / 2., totalLength / 2.);
  auto envVol = new TGeoVolume("SiliconTarget", envBox, air);
  envVol->SetLineColor(kGreen);
  envVol->SetTransparency(50);
 
  auto target = new TGeoBBox("Target", fTargetWidth / 2., fTargetHeight / 2.,
                             fTargetThickness / 2.);
  auto targetVol = new TGeoVolume("TargetVol", target, tungsten);
  targetVol->SetLineColor(kGray);
  targetVol->SetTransparency(40);
 
  for (Int_t i = 0; i < fLayers; i++) {
    // Compute the center position (z) for the current W layer
    Double_t zPos = -totalLength / 2 + i * fTargetSpacing;
 
    // Place the tungsten layer
    envVol->AddNode(targetVol, i,
                    new TGeoTranslation(0, 0, zPos + fTargetThickness / 2.));
 
    TGeoVolume* siliconPlanes = CreateSiliconPlanes(
        "TrackerPlane", fSensorWidth, fSensorLength,
        fTargetSpacing - fTargetThickness - 2. * fModuleOffset, Silicon, i);
    envVol->AddNode(
        siliconPlanes, i,
        new TGeoTranslation(0, 0, zPos + fTargetThickness + fModuleOffset));
  }
 
  // Finally, add the envelope to the top volume with the global z offset
  // fZCenter
  gGeoManager->GetTopVolume()->AddNode(envVol, 1,
                                       new TGeoTranslation(0, 0, fZPosition));
}

◆ CreateSiliconPlanes()

TGeoVolume * SiliconTarget::CreateSiliconPlanes	(	const char *	name,
		Double_t	sensorWidth,
		Double_t	sensorLength,
		Double_t	planeSpacing,
		TGeoMedium *	material,
		Int_t	layerId
	)

Definition at line 50 of file SiliconTarget.cxx.

                                                              {
  Double_t strip_pitch = 75.5 * um;
  Int_t nRows = 2;
  Int_t nCols = 4;
  Int_t nPlanes = 2;
 
  TGeoBBox* SensorShape =
      new TGeoBBox("SensorShape", width / 2, length / 2, 0.3 * mm / 2);
  TGeoVolume* SensorVolume =
      new TGeoVolume("SensorVolume", SensorShape, silicon);
  SensorVolume->SetLineColor(kRed);
  SensorVolume->SetTransparency(40);
 
  auto* Strips =
      SensorVolume->Divide("SLICEX", 1, 1298, -width / 2, strip_pitch);
  AddSensitiveVolume(Strips);
 
  TGeoVolumeAssembly* trackingStation =
      new TGeoVolumeAssembly("TrackingStation");
  // Each tracking station consists of X and Y planes
  for (Int_t plane = 0; plane < nPlanes; plane++) {
    TGeoVolumeAssembly* trackerPlane = new TGeoVolumeAssembly("TrackerPlane");
    // Each plane consists of 8 modules
    for (Int_t column = 0; column < nCols; column++) {
      for (Int_t row = 0; row < nRows; row++) {
        Int_t sensor_id =
            (layerId << 5) + (plane << 4) + (column << 2) + (row << 1);
        // Add 1mm gap between sensors for realistic placement
        trackerPlane->AddNode(
            SensorVolume, sensor_id,
            new TGeoTranslation(
                (-1.5 * width - 2. * mm) + (column * width + column * 1 * mm),
                (length / 2. + 0.5 * mm) - row * (length + 1. * mm), 0));
      }
    }
    if (plane == 0) {
      trackingStation->AddNode(trackerPlane, plane);
    } else if (plane == 1) {
      // Rotate the second plane by 90 degrees and translate such that the
      // planes are evenly spaced.
      trackingStation->AddNode(
          trackerPlane, plane,
          new TGeoCombiTrans(TGeoTranslation(0, 0, spacing),
                             TGeoRotation("y_rot", 0, 0, 90)));
    }
  }
 
  return trackingStation;
}

◆ operator=()

SiliconTarget & SiliconTarget::operator= ( const SiliconTarget & )

privatedelete

◆ ProcessHits()

Bool_t SiliconTarget::ProcessHits ( FairVolume * vol = 0 )

override

This method is called from the MC stepping

Definition at line 149 of file SiliconTarget.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 SiliconTargetPoint at exit of active volume
  if (gMC->IsTrackExiting() || gMC->IsTrackStop() ||
      gMC->IsTrackDisappeared()) {
    if (fELoss == 0.) {
      return kFALSE;
    }
 
    TParticle* p = gMC->GetStack()->GetCurrentTrack();
    fTrackID = gMC->GetStack()->GetCurrentTrackNumber();
    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.;
 
    Int_t strip_id = 0;
    Int_t sensor_id = 0;
    gMC->CurrentVolID(strip_id);
    gMC->CurrentVolOffID(1, sensor_id);
    fVolumeID = (sensor_id << 12) + strip_id;
 
    AddHit(fTrackID, fVolumeID, TVector3(xmean, ymean, zmean),
           TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength, fELoss,
           pdgCode);
 
    ShipStack* stack = dynamic_cast<ShipStack*>(gMC->GetStack());
    stack->AddPoint(kSiliconTarget);
  }
  return kTRUE;
}

◆ SetSiliconTargetParameters()

void SiliconTarget::SetSiliconTargetParameters	(	Double_t	targetWidth,
		Double_t	targetHeight,
		Double_t	sensorWidth,
		Double_t	sensorLength,
		Int_t	nLayers,
		Double_t	zPosition,
		Double_t	targetThickness,
		Double_t	targetSpacing,
		Double_t	moduleOffset
	)

Definition at line 35 of file SiliconTarget.cxx.

                                                                             {
  fTargetWidth = targetWidth;
  fTargetHeight = targetHeight;
  fSensorWidth = sensorWidth;
  fSensorLength = sensorLength;
  fLayers = nLayers;
  fZPosition = zPosition;
  fTargetThickness = targetThickness;
  fTargetSpacing = targetSpacing;
  fModuleOffset = moduleOffset;
}

Member Data Documentation

◆ fLayers

Int_t SiliconTarget::fLayers

private

Definition at line 39 of file SiliconTarget.h.

◆ fModuleOffset

Double_t SiliconTarget::fModuleOffset

private

Definition at line 43 of file SiliconTarget.h.

◆ fSensorLength

Double_t SiliconTarget::fSensorLength

private

Definition at line 38 of file SiliconTarget.h.

◆ fSensorWidth

Double_t SiliconTarget::fSensorWidth

private

Definition at line 37 of file SiliconTarget.h.

◆ fTargetHeight

Double_t SiliconTarget::fTargetHeight

private

Definition at line 36 of file SiliconTarget.h.

◆ fTargetSpacing

Double_t SiliconTarget::fTargetSpacing

private

Definition at line 42 of file SiliconTarget.h.

◆ fTargetThickness

Double_t SiliconTarget::fTargetThickness

private

Definition at line 41 of file SiliconTarget.h.

◆ fTargetWidth

Double_t SiliconTarget::fTargetWidth

private

Definition at line 35 of file SiliconTarget.h.

◆ fZPosition

Double_t SiliconTarget::fZPosition

private

Definition at line 40 of file SiliconTarget.h.

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

SND/SiliconTarget/SiliconTarget.h
SND/SiliconTarget/SiliconTarget.cxx

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiliconTarget() [1/3]

◆ SiliconTarget() [2/3]

◆ SiliconTarget() [3/3]

Member Function Documentation

◆ ConstructGeometry()

◆ CreateSiliconPlanes()

◆ operator=()

◆ ProcessHits()

◆ SetSiliconTargetParameters()

Member Data Documentation

◆ fLayers

◆ fModuleOffset

◆ fSensorLength

◆ fSensorWidth

◆ fTargetHeight

◆ fTargetSpacing

◆ fTargetThickness

◆ fTargetWidth

◆ fZPosition