#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 Attributes
SG::WriteHandle< SiHitCollection >	m_HitColl
GeoModelIO::ReadGeoModel *	m_sqlreader

Detailed Description

Definition at line 33 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 35 of file HGTDSensorGmxSD.cxx.

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

◆ ~HGTDSensorGmxSD()

virtual HGTDSensorGmxSD::~HGTDSensorGmxSD ( )

inlinevirtual

Definition at line 42 of file HGTDSensorGmxSD.h.

42{}

Member Function Documentation

◆ Initialize()

void HGTDSensorGmxSD::Initialize ( G4HCofThisEvent * )

finaloverride

Definition at line 44 of file HGTDSensorGmxSD.cxx.

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

◆ ProcessHits()

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

finaloverride

Definition at line 49 of file HGTDSensorGmxSD.cxx.

{
    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

SG::WriteHandle<SiHitCollection> HGTDSensorGmxSD::m_HitColl

private

Definition at line 52 of file HGTDSensorGmxSD.h.

◆ m_sqlreader

GeoModelIO::ReadGeoModel* HGTDSensorGmxSD::m_sqlreader