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

Private Attributes
SG::WriteHandle< SiHitCollection >	m_HitColl
GeoModelIO::ReadGeoModel *	m_sqlreader

Detailed Description

Definition at line 28 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 37 of file PixelSensorGmxSD.cxx.

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

~PixelSensorGmxSD()

virtual PixelSensorGmxSD::~PixelSensorGmxSD ( )

inlinevirtual

Definition at line 38 of file PixelSensorGmxSD.h.

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

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 49 of file PixelSensorGmxSD.h.

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

FRIEND_TEST() [1/3]

PixelSensorGmxSD::FRIEND_TEST ( PixelSensorGmxSDtest  , AddHit    )

private

FRIEND_TEST() [2/3]

PixelSensorGmxSD::FRIEND_TEST ( PixelSensorGmxSDtest  , Initialize    )

private

FRIEND_TEST() [3/3]

PixelSensorGmxSD::FRIEND_TEST ( PixelSensorGmxSDtest  , ProcessHits    )

private

Initialize()

void PixelSensorGmxSD::Initialize ( G4HCofThisEvent *  )

finaloverridevirtual

Definition at line 45 of file PixelSensorGmxSD.cxx.

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

ProcessHits()

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

finaloverridevirtual

Definition at line 51 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());
 
    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
                     trHelp.GenerateParticleLink(),
                     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(),
                     trHelp.GenerateParticleLink(),
                     id);
  return true; 
}

Member Data Documentation

m_HitColl

SG::WriteHandle<SiHitCollection> PixelSensorGmxSD::m_HitColl

private

Definition at line 53 of file PixelSensorGmxSD.h.

m_sqlreader

GeoModelIO::ReadGeoModel* PixelSensorGmxSD::m_sqlreader

private

Definition at line 54 of file PixelSensorGmxSD.h.

---

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

• PixelSensorGmxSD.h
• PixelSensorGmxSD.cxx