SctSensorGmxSD Class Reference

#include <SctSensorGmxSD.h>

Inheritance diagram for SctSensorGmxSD:

Collaboration diagram for SctSensorGmxSD:

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

Protected Attributes
std::string	m_HitCollName
SiHitCollection *	m_HitColl {nullptr}
AtlasG4EventUserInfo *	m_g4UserEventInfo {nullptr}

Private Member Functions
	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)

Private Attributes
GeoModelIO::ReadGeoModel *	m_sqlreader {nullptr}

Detailed Description

Definition at line 24 of file SctSensorGmxSD.h.

Constructor & Destructor Documentation

SctSensorGmxSD()

SctSensorGmxSD::SctSensorGmxSD	(	const std::string &	name,
		const std::string &	hitCollectionName,
		GeoModelIO::ReadGeoModel *	sqlreader = nullptr )

Definition at line 33 of file SctSensorGmxSD.cxx.

  : SctSensorSD( name, hitCollectionName) {
      m_sqlreader = sqlreader;
   }

~SctSensorGmxSD()

SctSensorGmxSD::~SctSensorGmxSD

(

)

Definition at line 38 of file SctSensorGmxSD.cxx.

{ }

Member Function Documentation

AddHit()

template<class... Args>

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

inlineinherited

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 45 of file SctSensorSD.h.

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

FRIEND_TEST() [1/3]

SctSensorSD::FRIEND_TEST ( SctSensorSDtest , AddHit  )

privateinherited

FRIEND_TEST() [2/3]

SctSensorSD::FRIEND_TEST ( SctSensorSDtest , indexMethod  )

privateinherited

FRIEND_TEST() [3/3]

SctSensorSD::FRIEND_TEST ( SctSensorSDtest , ProcessHits  )

privateinherited

indexMethod()

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

privateinherited

Definition at line 119 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 * )

finaloverrideinherited

Definition at line 39 of file SctSensorSD.cxx.

{
  // ISF calls G4SDManager::PrepareNewEvent() before the Geant4 event loop starts...
  if(auto* eventManger = G4EventManager::GetEventManager()){
    if(auto* eventInfo = static_cast<AtlasG4EventUserInfo*>(eventManger->GetUserInformation())){
      m_HitColl = eventInfo->GetHitCollectionMap()->Find<SiHitCollection>(m_HitCollName);
      m_g4UserEventInfo = eventInfo;
    }
  }
}

ProcessHits()

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

finaloverride

Definition at line 40 of file SctSensorGmxSD.cxx.

                                                                        {
  double edep = aStep->GetTotalEnergyDeposit();
  if (edep == 0.) {
    if (aStep->GetTrack()->GetDefinition() != G4Geantino::GeantinoDefinition() &&
       aStep->GetTrack()->GetDefinition() != G4ChargedGeantino::ChargedGeantinoDefinition())
      return false;
  }
  edep *= CLHEP::MeV;
 
  const G4TouchableHistory *myTouch = dynamic_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
  if (not myTouch) {
    G4cout << "SctSensorGmxSD::ProcessHits bad dynamic_cast" << G4endl;
    return false;
  }
  //
  // Get the hit start and end point local coordinates
  //
  G4ThreeVector coord1 = aStep->GetPreStepPoint()->GetPosition();
  G4ThreeVector coord2 = aStep->GetPostStepPoint()->GetPosition();
  const G4AffineTransform transformation = myTouch->GetHistory()->GetTopTransform();
  G4ThreeVector localPosition1 = transformation.TransformPoint(coord1);
  G4ThreeVector localPosition2 = transformation.TransformPoint(coord2);
  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;
 
  
  // get the HepMcParticleLink from the TrackHelper
  TrackHelper trHelp(aStep->GetTrack());
  auto particleLink = trHelp.GenerateParticleLink(m_g4UserEventInfo ? m_g4UserEventInfo->GetEventStore() : nullptr);
 
  if(m_sqlreader){
    //if sqlite inputs, Identifier indices come from PhysVol Name  
    std::string physVolName = myTouch->GetVolume(0)->GetName();
  
    int hitIdOfWafer = SiHitIdHelper::GetHelper()->buildHitIdFromStringITk(1,physVolName);                                                   
 
    m_HitColl->Emplace(lP1,
                    lP2,
                    edep,
                    aStep->GetPreStepPoint()->GetGlobalTime(),//use the global time. i.e. the time from the beginning of the event
                    std::move(particleLink),
                    hitIdOfWafer);
    return true;
  }
 
  // if not from SQLite, we assume that the Identifier has already been written in as the copy number 
  // (it should hsave done if GeoModel building ran within Athena)
  //
  //    Get the indexes of which detector the hit is in
  //
  const int id = myTouch->GetVolume(0)->GetCopyNo();
 
  m_HitColl->Emplace(lP1, lP2, edep, aStep->GetPreStepPoint()->GetGlobalTime(), std::move(particleLink), id);
 
  return true;
}

Member Data Documentation

m_g4UserEventInfo

AtlasG4EventUserInfo* SctSensorSD::m_g4UserEventInfo {nullptr}

protectedinherited

Definition at line 53 of file SctSensorSD.h.

{nullptr};

m_HitColl

SiHitCollection* SctSensorSD::m_HitColl {nullptr}

protectedinherited

Definition at line 52 of file SctSensorSD.h.

{nullptr};

m_HitCollName

std::string SctSensorSD::m_HitCollName

protectedinherited

Definition at line 51 of file SctSensorSD.h.

m_sqlreader

GeoModelIO::ReadGeoModel* SctSensorGmxSD::m_sqlreader {nullptr}

private

Definition at line 37 of file SctSensorGmxSD.h.

{nullptr};

---

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

• SctSensorGmxSD.h
• SctSensorGmxSD.cxx