HGTDSensorGmxSD Class Reference

#include <HGTDSensorGmxSD.h>

Inheritance diagram for HGTDSensorGmxSD:

Collaboration diagram for HGTDSensorGmxSD:

Public Member Functions
	HGTDSensorGmxSD (const std::string &name, const std::string &hitCollectionName, GeoModelIO::ReadGeoModel *sqlreader=nullptr)
virtual	~HGTDSensorGmxSD ()
G4bool	ProcessHits (G4Step *, G4TouchableHistory *) override final
void	Initialize (G4HCofThisEvent *) override final

Private Member Functions
SiHitCollection *	getHitCollection () const

Private Attributes
std::string	m_hitCollectionName
SiHitCollection *	m_hitColl {}
GeoModelIO::ReadGeoModel *	m_sqlreader

Detailed Description

Definition at line 32 of file HGTDSensorGmxSD.h.

Constructor & Destructor Documentation

HGTDSensorGmxSD()

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

Definition at line 37 of file HGTDSensorGmxSD.cxx.

    : G4VSensitiveDetector( name ),
      m_hitCollectionName( hitCollectionName ),
      m_sqlreader( sqlreader )
{
 
}

~HGTDSensorGmxSD()

virtual HGTDSensorGmxSD::~HGTDSensorGmxSD ( )

inlinevirtual

Definition at line 41 of file HGTDSensorGmxSD.h.

{}

Member Function Documentation

getHitCollection()

SiHitCollection * HGTDSensorGmxSD::getHitCollection ( ) const

private

Definition at line 146 of file HGTDSensorGmxSD.cxx.

{
    auto* eventInfo = AtlasG4EventUserInfo::GetEventUserInfo();
    if (!eventInfo) {
        return nullptr;
    }
 
    auto hitCollections = eventInfo->GetHitCollectionMap();
    return hitCollections ? hitCollections->Find<SiHitCollection>(m_hitCollectionName) : nullptr;
}

Initialize()

void HGTDSensorGmxSD::Initialize ( G4HCofThisEvent * )

finaloverride

Definition at line 46 of file HGTDSensorGmxSD.cxx.

{
    m_hitColl = getHitCollection();
}

ProcessHits()

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

finaloverride

Definition at line 51 of file HGTDSensorGmxSD.cxx.

{
    if (!m_hitColl) {
        m_hitColl = getHitCollection();
        if (!m_hitColl) {
            return false;
        }
    }
 
    if (verboseLevel>5) G4cout << "Process Hit" << G4endl;
 
    G4double edep = aStep->GetTotalEnergyDeposit();
    edep *= CLHEP::MeV;
 
    if (edep==0.) {
        if (aStep->GetTrack()->GetDefinition() != G4Geantino::GeantinoDefinition() &&
            aStep->GetTrack()->GetDefinition() != G4ChargedGeantino::ChargedGeantinoDefinition())
            return false;
    }
 
    //
    // Get the Touchable History:
    //
    const G4TouchableHistory*  myTouch = dynamic_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
    if (not myTouch){
      G4cout<<"HGTDSensorGmxSD::ProcessHits: dynamic cast failed"<<G4endl;
      return false;
    }
    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;
        }
    }
 
    //
    // Get the hit coordinates. Start and End Point
    //
    G4ThreeVector startCoord = aStep->GetPreStepPoint()->GetPosition();
    G4ThreeVector endCoord   = aStep->GetPostStepPoint()->GetPosition();
 
    // Create the SiHits
 
    const G4AffineTransform transformation = myTouch->GetHistory()->GetTopTransform();
 
    G4ThreeVector localPosition1 = transformation.TransformPoint(startCoord);
    G4ThreeVector localPosition2 = transformation.TransformPoint(endCoord);
 
    HepGeom::Point3D<double> lP1,lP2;
    //TODO need to check
    lP1[SiHit::xEta] = localPosition1[1]*CLHEP::mm; //long edge of the module
    lP1[SiHit::xPhi] = localPosition1[0]*CLHEP::mm; //short edge of the module
    lP1[SiHit::xDep] = localPosition1[2]*CLHEP::mm; //depth (z)
 
    lP2[SiHit::xEta] = localPosition2[1]*CLHEP::mm;
    lP2[SiHit::xPhi] = localPosition2[0]*CLHEP::mm;
    lP2[SiHit::xDep] = localPosition2[2]*CLHEP::mm;
 
    // get the HepMcParticleLink from the TrackHelper
    TrackHelper trHelp(aStep->GetTrack());
 
    if(m_sqlreader){
        // if sqlite inputs, Identifier indices come from PhysVol Name
        std::string physVolName = myTouch->GetVolume()->GetName();
 
        int hitIdOfWafer = SiHitIdHelper::GetHelper()->buildHitIdFromStringHGTD(2,physVolName);
 
        m_hitColl->Emplace(lP1,
                           lP2,
                           edep,
                           aStep->GetPreStepPoint()->GetGlobalTime(),
                           trHelp.GenerateParticleLink(),
                           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()->GetCopyNo();
 
    m_hitColl->Emplace(lP1,
                       lP2,
                       edep,
                       aStep->GetPreStepPoint()->GetGlobalTime(),
                       trHelp.GenerateParticleLink(),
                       id);
 
    return true;
}

Member Data Documentation

m_hitColl

SiHitCollection* HGTDSensorGmxSD::m_hitColl {}

private

Definition at line 55 of file HGTDSensorGmxSD.h.

{};

m_hitCollectionName

std::string HGTDSensorGmxSD::m_hitCollectionName

private

Definition at line 53 of file HGTDSensorGmxSD.h.

m_sqlreader

GeoModelIO::ReadGeoModel* HGTDSensorGmxSD::m_sqlreader

private

Definition at line 56 of file HGTDSensorGmxSD.h.

---

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

• HGTDSensorGmxSD.h
• HGTDSensorGmxSD.cxx