SctSensor_CTB Class Reference

#include <SctSensor_CTB.h>

Inheritance diagram for SctSensor_CTB:

Collaboration diagram for SctSensor_CTB:

Public Member Functions
	SctSensor_CTB (const std::string &name, const std::string &hitCollectionName)
	~SctSensor_CTB ()
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 (SctSensor_CTBtest, Initialize)
	FRIEND_TEST (SctSensor_CTBtest, ProcessHits)
	FRIEND_TEST (SctSensor_CTBtest, AddHit)

Private Attributes
SG::WriteHandle< SiHitCollection >	m_HitColl

Detailed Description

Definition at line 27 of file SctSensor_CTB.h.

Constructor & Destructor Documentation

SctSensor_CTB()

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

Definition at line 33 of file SctSensor_CTB.cxx.

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

~SctSensor_CTB()

SctSensor_CTB::~SctSensor_CTB ( )

inline

Definition at line 37 of file SctSensor_CTB.h.

{ /* If all goes well we do not own myHitColl at this point */ }

Member Function Documentation

AddHit()

template<class... Args>

void SctSensor_CTB::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 48 of file SctSensor_CTB.h.

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

FRIEND_TEST() [1/3]

SctSensor_CTB::FRIEND_TEST ( SctSensor_CTBtest  , AddHit    )

private

FRIEND_TEST() [2/3]

SctSensor_CTB::FRIEND_TEST ( SctSensor_CTBtest  , Initialize    )

private

FRIEND_TEST() [3/3]

SctSensor_CTB::FRIEND_TEST ( SctSensor_CTBtest  , ProcessHits    )

private

Initialize()

void SctSensor_CTB::Initialize ( G4HCofThisEvent *  )

finaloverride

Definition at line 41 of file SctSensor_CTB.cxx.

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

ProcessHits()

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

finaloverride

Definition at line 48 of file SctSensor_CTB.cxx.

{
  G4double edep = aStep->GetTotalEnergyDeposit();
  if(edep==0.){
    if(aStep->GetTrack()->GetDefinition()!=G4Geantino::GeantinoDefinition() &&
       aStep->GetTrack()->GetDefinition()!=G4ChargedGeantino::ChargedGeantinoDefinition())
      return false;
  }
  edep *= CLHEP::MeV;
  //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 << "SctSensor_CTB::ProcessHits bad dynamic_cast" << G4endl;
    return false;
  }
  //
  // 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;
 
  //
  // Get it into a vector in local coords and with the right units:
 
 
  // Now Navigate the history to know in what detector the step is:
  // and finally set the ID of det element in which the hit is.
  //
  // G4int History;
  //
  // Is it TB, barrel or endcap?
  //
  int BrlEcap = 0;
  int LayerDisk = 0;
  int etaMod = 1;
  int phiMod = 0;
  int side = 0;
  //
  // In the case of the TB the positioning integers won't be initialized
  // and the identifying integer will be zero. There is no need to do
  // anything else
 
  int BEcopyNo =  myTouch->GetVolume()->GetCopyNo();
 
#ifdef SCTSD_CTB_DEBUG
  for (int i=0;i<myTouch->GetHistoryDepth();i++){
    std::string detname=myTouch->GetVolume(i)->GetLogicalVolume()->GetName();
    int copyno=myTouch->GetVolume(i)->GetCopyNo();
    if (verboseLevel>1) G4cout << "Volume "<<detname <<" Copy Nr. "<<copyno << G4endl;
  }
#endif
 
  // Testbeam geometry from GeoModel.
  //
  int BEcopyNoTest =  BEcopyNo/100;
  if(BEcopyNoTest == 5 || BEcopyNoTest == 6 || BEcopyNo == 1000 || BEcopyNo == 2000) {
    BrlEcap=0;
    etaMod=1;
 
    if(BEcopyNoTest == 5 || BEcopyNoTest == 6 ) {
      side      = myTouch->GetVolume(0)->GetCopyNo()%2;
      phiMod=myTouch->GetVolume(1)->GetCopyNo()%2;
      LayerDisk = myTouch->GetVolume(1)->GetCopyNo()/2;
    } else if ( BEcopyNo == 1000) {
      side      = myTouch->GetVolume(1)->GetCopyNo();
      phiMod=myTouch->GetVolume(2)->GetCopyNo()%2;
      LayerDisk = myTouch->GetVolume(2)->GetCopyNo()/2;
    } // otherwise use side = phiMod = LayerDisk = 0.
    // Some printout
#ifdef SCTSD_CTB_DEBUG
    if (verboseLevel>1) G4cout << "In the SCT TestBeam" << G4endl;
    if (verboseLevel>1) G4cout << " * Module Number: " << LayerDisk << "," << phiMod << " Side: " << side << G4endl;
    if (verboseLevel>5) G4cout << " Identifier will be [2.2."
                     << BrlEcap << "."
                     << LayerDisk << "."
                     << phiMod << "."
                     << etaMod << "."
                     << side << "]" << G4endl;
#endif
  } else {
    G4ExceptionDescription description;
    description << "ProcessHits: Unexpected copy number for sensor";
    G4Exception("SctSensor_CTB", "UnexpectedCopyNumberForSctSensor", FatalException, description);
    abort();
  }
  TrackHelper trHelp(aStep->GetTrack());
  m_HitColl->Emplace(lP1,
                     lP2,
                     edep,
                     aStep->GetPreStepPoint()->GetGlobalTime(),
                     trHelp.GenerateParticleLink(),
                     1,BrlEcap,LayerDisk,etaMod,phiMod,side);
  return true;
}

Member Data Documentation

m_HitColl

SG::WriteHandle<SiHitCollection> SctSensor_CTB::m_HitColl

private

Definition at line 52 of file SctSensor_CTB.h.

---

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

• SctSensor_CTB.h
• SctSensor_CTB.cxx