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)
	~BLMSensorSD ()
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, Initialize)
	FRIEND_TEST (BLMSensorSDtest, ProcessHits)
	FRIEND_TEST (BLMSensorSDtest, AddHit)

Private Attributes
SG::WriteHandle< SiHitCollection >	m_HitColl

Detailed Description

Definition at line 24 of file BLMSensorSD.h.

Constructor & Destructor Documentation

BLMSensorSD()

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

Definition at line 28 of file BLMSensorSD.cxx.

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

~BLMSensorSD()

BLMSensorSD::~BLMSensorSD ( )

inline

Definition at line 34 of file BLMSensorSD.h.

{}

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 45 of file BLMSensorSD.h.

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

FRIEND_TEST() [1/3]

BLMSensorSD::FRIEND_TEST ( BLMSensorSDtest  , AddHit    )

private

FRIEND_TEST() [2/3]

BLMSensorSD::FRIEND_TEST ( BLMSensorSDtest  , Initialize    )

private

FRIEND_TEST() [3/3]

BLMSensorSD::FRIEND_TEST ( BLMSensorSDtest  , ProcessHits    )

private

Initialize()

void BLMSensorSD::Initialize ( G4HCofThisEvent *  )

finaloverride

Definition at line 36 of file BLMSensorSD.cxx.

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

ProcessHits()

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

finaloverride

Definition at line 41 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

SG::WriteHandle<SiHitCollection> BLMSensorSD::m_HitColl

private

Definition at line 49 of file BLMSensorSD.h.

---

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

• BLMSensorSD.h
• BLMSensorSD.cxx