PixelSensorGmxSD Class Reference

#include <PixelSensorGmxSD.h>

Inheritance diagram for PixelSensorGmxSD:

Collaboration diagram for PixelSensorGmxSD:

Public Member Functions
	PixelSensorGmxSD (const std::string &name, const std::string &hitCollectionName, GeoModelIO::ReadGeoModel *sqlreader=nullptr)
virtual G4bool	ProcessHits (G4Step *, G4TouchableHistory *) override final
virtual 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...

Private Member Functions
	FRIEND_TEST (PixelSensorGmxSDtest, ProcessHits)
	FRIEND_TEST (PixelSensorGmxSDtest, AddHit)

Private Attributes
std::string	m_HitCollName
SiHitCollection *	m_HitColl {nullptr}
GeoModelIO::ReadGeoModel *	m_sqlreader {nullptr}

Detailed Description

Definition at line 27 of file PixelSensorGmxSD.h.

Constructor & Destructor Documentation

PixelSensorGmxSD()

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

Definition at line 36 of file PixelSensorGmxSD.cxx.

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

Member Function Documentation

AddHit()

template<class... Args>

void PixelSensorGmxSD::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 44 of file PixelSensorGmxSD.h.

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

FRIEND_TEST() [1/2]

PixelSensorGmxSD::FRIEND_TEST ( PixelSensorGmxSDtest  , AddHit    )

private

FRIEND_TEST() [2/2]

PixelSensorGmxSD::FRIEND_TEST ( PixelSensorGmxSDtest  , ProcessHits    )

private

Initialize()

void PixelSensorGmxSD::Initialize ( G4HCofThisEvent *  )

finaloverridevirtual

Definition at line 44 of file PixelSensorGmxSD.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);
    }
  }
 
}

ProcessHits()

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

finaloverridevirtual

Definition at line 56 of file PixelSensorGmxSD.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;
  }
 
  //use the global time. i.e. the time from the beginning of the event
  //
  // Get the Touchable History:
  //
  const G4TouchableHistory *myTouch = dynamic_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
  if (not myTouch) {
    G4cout << "PixelSensorGmxSD::ProcessHits bad dynamic_cast" << 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 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);
 
  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;
 
  TrackHelper trHelp(aStep->GetTrack());
  // Temporary solution to get EventContext from a Geant4 thread 
  EventContext const* eventContext{nullptr};
  if(auto* eventInfo = AtlasG4EventUserInfo::GetEventUserInfo()){
    eventContext = &eventInfo->GetEventContext();
  }
  auto mcParticleLink = eventContext ? 
    trHelp.GenerateParticleLink(*eventContext) : 
    trHelp.GenerateParticleLink();
 
    if(m_sqlreader){
        //if sqlite inputs, Identifier indices come from PhysVol Name  
        std::string physVolName = myTouch->GetVolume(0)->GetName();
 
        int hitIdOfWafer = SiHitIdHelper::GetHelper()->buildHitIdFromStringITk(0,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(mcParticleLink),
                     hitIdOfWafer);
        return true;
        
    }
    // if not from SQLite, we assume that the Identifier has already been written in as the copy number 
  // (it should have done if GeoModel building ran within Athena)
 
  int id = myTouch->GetVolume()->GetCopyNo();
 
  m_HitColl->Emplace(lP1,
                     lP2,
                     edep,
                     aStep->GetPreStepPoint()->GetGlobalTime(),
                     std::move(mcParticleLink),
                     id);
  return true; 
}

Member Data Documentation

m_HitColl

SiHitCollection* PixelSensorGmxSD::m_HitColl {nullptr}

private

Definition at line 49 of file PixelSensorGmxSD.h.

m_HitCollName

std::string PixelSensorGmxSD::m_HitCollName

private

Definition at line 48 of file PixelSensorGmxSD.h.

m_sqlreader

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

private

Definition at line 50 of file PixelSensorGmxSD.h.

---

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

• PixelSensorGmxSD.h
• PixelSensorGmxSD.cxx