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

Private Attributes
SG::WriteHandle< SiHitCollection >	m_HitColl

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_HitColl( hitCollectionName )
{
}

~BCMSensorSD()

BCMSensorSD::~BCMSensorSD ( )

inline

Definition at line 34 of file BCMSensorSD.h.

{ /* I don't own myHitColl if all has gone well */ }

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

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

FRIEND_TEST() [1/3]

BCMSensorSD::FRIEND_TEST ( BCMSensorSDtest  , AddHit    )

private

FRIEND_TEST() [2/3]

BCMSensorSD::FRIEND_TEST ( BCMSensorSDtest  , Initialize    )

private

FRIEND_TEST() [3/3]

BCMSensorSD::FRIEND_TEST ( BCMSensorSDtest  , ProcessHits    )

private

Initialize()

void BCMSensorSD::Initialize ( G4HCofThisEvent *  )

finaloverride

Definition at line 39 of file BCMSensorSD.cxx.

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

ProcessHits()

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

finaloverride

Definition at line 44 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;
 
      m_HitColl->Emplace(lP1, lP2, edep, aStep->GetPreStepPoint()->GetGlobalTime(), trHelp.GenerateParticleLink(),
                         0, 0, myTouch->GetVolume(1)->GetCopyNo()-951, BEcopyNo - 11950, primaren, produced_in_diamond);
    }
  return true;
}

Member Data Documentation

m_HitColl

SG::WriteHandle<SiHitCollection> BCMSensorSD::m_HitColl

private

Definition at line 49 of file BCMSensorSD.h.

---

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

• BCMSensorSD.h
• BCMSensorSD.cxx