BLMSensorSD Class Reference

#include <BLMSensorSD.h>

Inheritance diagram for BLMSensorSD:

Collaboration diagram for BLMSensorSD:

Public Member Functions
	BLMSensorSD (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. More...

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

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

Detailed Description

Definition at line 23 of file BLMSensorSD.h.

Constructor & Destructor Documentation

BLMSensorSD()

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

Definition at line 29 of file BLMSensorSD.cxx.

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

Member Function Documentation

AddHit()

template<class... Args>

void BLMSensorSD::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 40 of file BLMSensorSD.h.

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

FRIEND_TEST() [1/2]

BLMSensorSD::FRIEND_TEST ( BLMSensorSDtest  , AddHit    )

private

FRIEND_TEST() [2/2]

BLMSensorSD::FRIEND_TEST ( BLMSensorSDtest  , ProcessHits    )

private

Initialize()

void BLMSensorSD::Initialize ( G4HCofThisEvent *  )

finaloverride

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

finaloverride

Definition at line 47 of file BLMSensorSD.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 << "BLMSensorSD::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;
 
  //BLM hit stuff
  if(BEcopyNo == 2009)
    {
      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::blmDiamondLog")
        produced_in_diamond = 1;
      else if(aStep->GetTrack()->GetLogicalVolumeAtVertex()->GetName() == "Pixel::blmModLog")
        produced_in_diamond = 2;
      else if(aStep->GetTrack()->GetLogicalVolumeAtVertex()->GetName() == "Pixel::blmWallLog")
        produced_in_diamond = 3;
 
      m_HitColl->Emplace(lP1, lP2, edep, aStep->GetPreStepPoint()->GetGlobalTime(), trHelp.GenerateParticleLink(),
                         0, 0, myTouch->GetVolume(1)->GetCopyNo()-222, 0, primaren, produced_in_diamond);
    }
  return true;
}

Member Data Documentation

m_HitColl

SiHitCollection* BLMSensorSD::m_HitColl {nullptr}

private

Definition at line 44 of file BLMSensorSD.h.

m_HitCollName

std::string BLMSensorSD::m_HitCollName

private

Definition at line 43 of file BLMSensorSD.h.

---

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

• BLMSensorSD.h
• BLMSensorSD.cxx