SctSensorSD Class Reference

#include <SctSensorSD.h>

Inheritance diagram for SctSensorSD:

Collaboration diagram for SctSensorSD:

Public Member Functions
	SctSensorSD (const std::string &name, const std::string &hitCollectionName)
	~SctSensorSD ()
G4bool	ProcessHits (G4Step *, G4TouchableHistory *) override
void	Initialize (G4HCofThisEvent *) override final
template<class... Args>
void	AddHit (Args &&... args)
	Templated method to stuff a single hit into the sensitive detector class. More...

Protected Attributes
SG::WriteHandle< SiHitCollection >	m_HitColl

Private Member Functions
	FRIEND_TEST (SctSensorSDtest, Initialize)
	FRIEND_TEST (SctSensorSDtest, ProcessHits)
	FRIEND_TEST (SctSensorSDtest, indexMethod)
	FRIEND_TEST (SctSensorSDtest, AddHit)
void	indexMethod (const G4TouchableHistory *myTouch, double coord1z, int &brlEcap, int &layerDisk, int &etaMod, int &phiMod, int &side)

Detailed Description

Definition at line 24 of file SctSensorSD.h.

Constructor & Destructor Documentation

SctSensorSD()

SctSensorSD::SctSensorSD	(	const std::string &	name,
		const std::string &	hitCollectionName
	)

Definition at line 31 of file SctSensorSD.cxx.

  : G4VSensitiveDetector( name )
  , m_HitColl( hitCollectionName )
{
}

~SctSensorSD()

SctSensorSD::~SctSensorSD ( )

inline

Definition at line 35 of file SctSensorSD.h.

{ /* If all goes well we do not own myHitColl here */ }

Member Function Documentation

AddHit()

template<class... Args>

void SctSensorSD::AddHit ( Args &&...  args )

inline

Templated method to stuff a single hit into the sensitive detector class.

This could get rather tricky, but the idea is to allow fast simulations to use the very same SD classes as the standard simulation.

Definition at line 46 of file SctSensorSD.h.

{ m_HitColl->Emplace( args... ); }

FRIEND_TEST() [1/4]

SctSensorSD::FRIEND_TEST ( SctSensorSDtest  , AddHit    )

private

FRIEND_TEST() [2/4]

SctSensorSD::FRIEND_TEST ( SctSensorSDtest  , indexMethod    )

private

FRIEND_TEST() [3/4]

SctSensorSD::FRIEND_TEST ( SctSensorSDtest  , Initialize    )

private

FRIEND_TEST() [4/4]

SctSensorSD::FRIEND_TEST ( SctSensorSDtest  , ProcessHits    )

private

indexMethod()

void SctSensorSD::indexMethod ( const G4TouchableHistory *  myTouch, double  coord1z, int &  brlEcap, int &  layerDisk, int &  etaMod, int &  phiMod, int &  side  )

private

Definition at line 113 of file SctSensorSD.cxx.

                                                                                                 {
 
 
 //
  // In the case of the TB the positioning integers won't be initialized
  // and the identifying integer will be zero. There is no need to do
  // anything else
 
  const int sensorCopyNo = myTouch->GetVolume()->GetCopyNo();
 
  if (verboseLevel>5){
    for (int i=0;i<myTouch->GetHistoryDepth();i++){
      std::string detname=myTouch->GetVolume(i)->GetLogicalVolume()->GetName();
      int copyno=myTouch->GetVolume(i)->GetCopyNo();
      G4cout << "Volume "<<detname << " Copy Nr. " << copyno << G4endl;
    }
  }
 
  // Barrel geometry from GeoModel. Flag it with a 1000 in the IdTag
  //
  if(sensorCopyNo/1000) {
    // Barrel
    if(sensorCopyNo == 1000) {
      side = myTouch->GetVolume(1)->GetCopyNo();
      etaMod = myTouch->GetVolume(2)->GetCopyNo();
      phiMod = myTouch->GetVolume(3)->GetCopyNo();
      layerDisk = myTouch->GetVolume(5)->GetCopyNo();
    } else if(sensorCopyNo == 1100) {
      // SLHC stave layout
      int etaModTmp = myTouch->GetVolume(1)->GetCopyNo();
      int sign = (etaModTmp>=0)? +1 : -1;
      side = std::abs(etaModTmp)%2;
      etaMod = sign * std::abs(etaModTmp)/2;
      phiMod = myTouch->GetVolume(2)->GetCopyNo();
      layerDisk = myTouch->GetVolume(4)->GetCopyNo();
    } else if(sensorCopyNo/100 == 12) {
      // Segmented sensor (SLHC only)
      int iSegment = sensorCopyNo%100;
      int sensorEnvCopyNo = myTouch->GetVolume(1)->GetCopyNo();
      if (sensorEnvCopyNo == 1000) {
        side = myTouch->GetVolume(2)->GetCopyNo();
        etaMod = myTouch->GetVolume(3)->GetCopyNo() + iSegment;
        phiMod = myTouch->GetVolume(4)->GetCopyNo();
        layerDisk = myTouch->GetVolume(6)->GetCopyNo();
      } else if  (sensorEnvCopyNo == 1100) {
        // Stave layout
        int etaModTmp = myTouch->GetVolume(2)->GetCopyNo();
        int sign = (etaModTmp>=0)? +1 : -1;
        side = std::abs(etaModTmp)%2;
        etaMod = sign * std::abs(etaModTmp)/2 + iSegment;
        phiMod = myTouch->GetVolume(3)->GetCopyNo();
        layerDisk = myTouch->GetVolume(5)->GetCopyNo();
      } else {
        G4ExceptionDescription description;
        description << "ProcessHits: Unrecognized geometry in SCT sensitive detector. Please contact the maintainer of the SCT Detector Description.";
        G4Exception("SctSensorSD", "UnrecognizedSCTGeometry1", FatalException, description);
        abort();
      }
    } else {
      G4ExceptionDescription description;
      description << "ProcessHits: Unrecognized geometry in SCT sensitive detector. Please contact the maintainer of the SCT Detector Description.";
      G4Exception("SctSensorSD", "UnrecognizedSCTGeometry2", FatalException, description);
      abort();
    }
    if (verboseLevel>5){
      G4cout << "In the SCT Barrel" << G4endl;
      G4cout << "----- side       : " << side << G4endl;
      G4cout << "----- eta_module : " << etaMod << G4endl;
      G4cout << "----- phi_module : " << phiMod << G4endl;
      G4cout << "----- layer      : " << layerDisk << G4endl;
      G4cout << "----- barrel_ec  : " << brlEcap<< G4endl;
    }
  } else {
    // Endcap geometry from GeoModel. Flag it with a 500 or 600 in the IdTag
    // The level in the hierarchy is different for different geometries
    // 500 For pre DC3
    // 600 For DC3 and later (including SLHC)
    int sensorCopyNoTest =  sensorCopyNo/100;
    if(sensorCopyNoTest == 5 || sensorCopyNoTest == 6) {
      int signZ = (coord1z > 0) ? 1 : -1;
      brlEcap = 2 * signZ;
      side      = myTouch->GetVolume(0)->GetCopyNo() % 2;
      if (sensorCopyNoTest == 5) { // DC2 and Rome
        etaMod    = myTouch->GetVolume(3)->GetCopyNo();
        layerDisk = myTouch->GetVolume(4)->GetCopyNo();
      } else { // DC3 and later and SLHC
        etaMod    = myTouch->GetVolume(2)->GetCopyNo();
        layerDisk = myTouch->GetVolume(3)->GetCopyNo();
      }
      int phiNumber = myTouch->GetVolume(1)->GetCopyNo();
      // Number of modules in ring and module which becomes 0 is encoded in the copy number.
      // This is in order to recreate the correct id in the negative endcap.
      int maxModules = (phiNumber & 0x00ff0000) >> 16;
      int moduleZero = (phiNumber & 0xff000000) >> 24; // When phiMod = moduleZero -> 0 after inverting
      phiMod = phiNumber & 0x0000ffff;
      // If -ve endcap then invert numbering
      // maxModules is non-zero, but check for maxModules != 0  safeguards
      // against divide by zero (in case we create a geometry without the encoding).
      if (maxModules && signZ < 0) phiMod = (maxModules + moduleZero - phiMod) % maxModules;
 
      // Some printout
      if (verboseLevel>5){
        G4cout << "In the SCT EndCap" << G4endl;
        G4cout << "----- side       : " << side << G4endl;
        G4cout << "----- phi_module : " << phiMod<< G4endl;
        G4cout << "----- eta_module : " << etaMod << G4endl;
        G4cout << "----- layerDisk  : " << layerDisk << G4endl;
        G4cout << "----- barrel_ec  : " << brlEcap << G4endl;
      }
    } else {
      // Do not expect other numbers. Need to fix SCT_GeoModel or SCT_SLHC_GeoModel if this occurs.
      G4ExceptionDescription description;
      description << "ProcessHits: Unrecognized geometry in SCT sensitive detector. Please contact the maintainer of the SCT Detector Description.";
      G4Exception("SctSensorSD", "UnrecognizedSCTGeometry3", FatalException, description);
      abort();
    }
  }
  return;
}

Initialize()

void SctSensorSD::Initialize ( G4HCofThisEvent *  )

finaloverride

Definition at line 39 of file SctSensorSD.cxx.

{
  if (!m_HitColl.isValid()) m_HitColl = std::make_unique<SiHitCollection>();
}

ProcessHits()

G4bool SctSensorSD::ProcessHits ( G4Step *  aStep, G4TouchableHistory *    )

override

Definition at line 46 of file SctSensorSD.cxx.

{
  double edep = aStep->GetTotalEnergyDeposit();
  if(edep==0.) {
    if(aStep->GetTrack()->GetDefinition()!=G4Geantino::GeantinoDefinition() &&
       aStep->GetTrack()->GetDefinition()!=G4ChargedGeantino::ChargedGeantinoDefinition())
      return false;
  }
  edep *= CLHEP::MeV;
  //
  // Get the Touchable History:
  //
  const G4TouchableHistory *myTouch = dynamic_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
  if (not myTouch) {
    G4cout << "SctSensorSD::ProcessHits bad dynamic_cast" << G4endl;
    return false;
  }
  //
  // Get the hit coordinates. Start and End Point
  //
  G4ThreeVector coord1 = aStep->GetPreStepPoint()->GetPosition();
  G4ThreeVector coord2 = aStep->GetPostStepPoint()->GetPosition();
  //
  // Calculate the local step begin and end position.
  // From a G4 FAQ:
  // http://geant4-hn.slac.stanford.edu:5090/HyperNews/public/get/geometry/17/1.html
  //
  const G4AffineTransform transformation = myTouch->GetHistory()->GetTopTransform();
  G4ThreeVector localPosition1 = transformation.TransformPoint(coord1);
  G4ThreeVector localPosition2 = transformation.TransformPoint(coord2);
  //
  // Get it into a vector in local coords and with the right units:
  //
  HepGeom::Point3D<double> lP1,lP2;
 
  lP1[SiHit::xEta] = localPosition1[2]*CLHEP::mm;
  lP1[SiHit::xPhi] = localPosition1[1]*CLHEP::mm;
  lP1[SiHit::xDep] = localPosition1[0]*CLHEP::mm;
 
  lP2[SiHit::xEta] = localPosition2[2]*CLHEP::mm;
  lP2[SiHit::xPhi] = localPosition2[1]*CLHEP::mm;
  lP2[SiHit::xDep] = localPosition2[0]*CLHEP::mm;
 
  // Now Navigate the history to know in what detector the step is:
  // and finally set the ID of det element in which the hit is.
  //
  //G4int History;
  //
  // Is it TB, barrel or endcap?
  //
  int brlEcap = 0;
  int layerDisk = 0;
  int etaMod = 0;
  int phiMod = 0;
  int side = 0;
  this->indexMethod(myTouch, coord1.z(), brlEcap, layerDisk, etaMod, phiMod, side);
  // get the HepMcParticleLink from the TrackHelper
  TrackHelper trHelp(aStep->GetTrack());
  m_HitColl->Emplace(lP1,
                     lP2,
                     edep,
                     aStep->GetPreStepPoint()->GetGlobalTime(),//use the global time. i.e. the time from the beginning of the event
                     trHelp.GenerateParticleLink(),
                     1,brlEcap,layerDisk,etaMod,phiMod,side);
  return true;
}

Member Data Documentation

m_HitColl

SG::WriteHandle<SiHitCollection> SctSensorSD::m_HitColl

protected

Definition at line 52 of file SctSensorSD.h.

---

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

• SctSensorSD.h
• SctSensorSD.cxx