BCMSensorSD Class Reference

#include <BCMSensorSD.h>

Inheritance diagram for BCMSensorSD:

Collaboration diagram for BCMSensorSD:

Public Member Functions
	BCMSensorSD (const std::string &name, const std::string &hitCollectionName)
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.

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

Private Attributes
std::string	m_HitCollName
	Name of the hit collection.
SiHitCollection *	m_HitColl {nullptr}
	Pointer to the hit collection.
AtlasG4EventUserInfo *	m_g4UserEventInfo {nullptr}

Detailed Description

Definition at line 24 of file BCMSensorSD.h.

Constructor & Destructor Documentation

BCMSensorSD()

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

Definition at line 31 of file BCMSensorSD.cxx.

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

Member Function Documentation

AddHit()

template<class... Args>

void BCMSensorSD::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 41 of file BCMSensorSD.h.

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

FRIEND_TEST() [1/2]

BCMSensorSD::FRIEND_TEST ( BCMSensorSDtest , AddHit  )

private

FRIEND_TEST() [2/2]

BCMSensorSD::FRIEND_TEST ( BCMSensorSDtest , ProcessHits  )

private

Initialize()

void BCMSensorSD::Initialize ( G4HCofThisEvent * )

finaloverride

Definition at line 39 of file BCMSensorSD.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 BCMSensorSD::ProcessHits ( G4Step * aStep, G4TouchableHistory *  )

finaloverride

Definition at line 50 of file BCMSensorSD.cxx.

{
  G4double edep = aStep->GetTotalEnergyDeposit();
  edep *= CLHEP::MeV;
 
  //if there is no energy deposition skip everything
  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 << "BCMSensorSD::ProcessHits bad dynamic_cast" << G4endl;
    return false;
  }
  
  int BEcopyNo =  myTouch->GetVolume()->GetCopyNo();
 
  // 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;
 
  //BCM hit stuff
  if(BEcopyNo == 11950 || BEcopyNo == 11951)
    {
      TrackHelper trHelp(aStep->GetTrack());
      //primary or not
      int primaren = 0;
      if(trHelp.IsPrimary())
        primaren = 1;
      else if(trHelp.IsRegeneratedPrimary())
        primaren = 2;
      else if(trHelp.IsSecondary())
        primaren = 3;
      else if(trHelp.IsRegisteredSecondary())
        primaren = 4;
 
      int produced_in_diamond = 0;
      if(aStep->GetTrack()->GetLogicalVolumeAtVertex()->GetName() == "Pixel::bcmDiamondLog")
        produced_in_diamond = 1;
      else if(aStep->GetTrack()->GetLogicalVolumeAtVertex()->GetName() == "Pixel::bcmModLog")
        produced_in_diamond = 2;
      else if(aStep->GetTrack()->GetLogicalVolumeAtVertex()->GetName() == "Pixel::bcmWallLog")
        produced_in_diamond = 3;
      // Temporary solution to get EventContext from a Geant4 thread 
      m_HitColl->Emplace(lP1, lP2, edep, aStep->GetPreStepPoint()->GetGlobalTime(),
                         trHelp.GenerateParticleLink(m_g4UserEventInfo ? m_g4UserEventInfo->GetEventStore() : nullptr),
                         0, 0, myTouch->GetVolume(1)->GetCopyNo()-951, BEcopyNo - 11950, primaren, produced_in_diamond);
    }
  return true;
}

Member Data Documentation

m_g4UserEventInfo

AtlasG4EventUserInfo* BCMSensorSD::m_g4UserEventInfo {nullptr}

private

Definition at line 45 of file BCMSensorSD.h.

{nullptr};

m_HitColl

SiHitCollection* BCMSensorSD::m_HitColl {nullptr}

private

Pointer to the hit collection.

Definition at line 44 of file BCMSensorSD.h.

{nullptr}; 

m_HitCollName

std::string BCMSensorSD::m_HitCollName

private

Name of the hit collection.

Definition at line 43 of file BCMSensorSD.h.

---

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

• BCMSensorSD.h
• BCMSensorSD.cxx