sTGCSensitiveDetector Class Reference

#include <sTGCSensitiveDetector.h>

Inheritance diagram for sTGCSensitiveDetector:

Collaboration diagram for sTGCSensitiveDetector:

Public Member Functions
	sTGCSensitiveDetector (const std::string &name, const std::string &hitCollectionName)
	construction/destruction More...
	~sTGCSensitiveDetector ()
void	Initialize (G4HCofThisEvent *HCE) override final
	member functions More...
G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROhist) override final

Private Member Functions
	FRIEND_TEST (sTGCSensitiveDetectortest, Initialize)
	FRIEND_TEST (sTGCSensitiveDetectortest, ProcessHits)

Private Attributes
SG::WriteHandle< sTGCSimHitCollection >	m_sTGCSimHitCollection
const sTgcHitIdHelper *	m_muonHelper

Detailed Description

Definition at line 15 of file sTGCSensitiveDetector.h.

Constructor & Destructor Documentation

sTGCSensitiveDetector()

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

construction/destruction

Definition at line 18 of file sTGCSensitiveDetector.cxx.

  : G4VSensitiveDetector( name )
  , m_sTGCSimHitCollection( hitCollectionName )
{
  m_muonHelper = sTgcHitIdHelper::GetHelper();
  //m_muonHelper->PrintFields();
}

~sTGCSensitiveDetector()

sTGCSensitiveDetector::~sTGCSensitiveDetector ( )

inline

Definition at line 22 of file sTGCSensitiveDetector.h.

{}

Member Function Documentation

FRIEND_TEST() [1/2]

sTGCSensitiveDetector::FRIEND_TEST ( sTGCSensitiveDetectortest  , Initialize    )

private

FRIEND_TEST() [2/2]

sTGCSensitiveDetector::FRIEND_TEST ( sTGCSensitiveDetectortest  , ProcessHits    )

private

Initialize()

void sTGCSensitiveDetector::Initialize ( G4HCofThisEvent *  HCE )

finaloverride

member functions

Definition at line 27 of file sTGCSensitiveDetector.cxx.

{
  if (!m_sTGCSimHitCollection.isValid()) m_sTGCSimHitCollection = std::make_unique<sTGCSimHitCollection>();
}

ProcessHits()

G4bool sTGCSensitiveDetector::ProcessHits ( G4Step *  aStep, G4TouchableHistory *  ROhist  )

finaloverride

Definition at line 32 of file sTGCSensitiveDetector.cxx.

{
  G4Track* currentTrack = aStep->GetTrack();
  int charge=currentTrack->GetDefinition()->GetPDGCharge();
 
  bool geantinoHit = (currentTrack->GetDefinition()==G4Geantino::GeantinoDefinition()) ||
                     (currentTrack->GetDefinition()==G4ChargedGeantino::ChargedGeantinoDefinition());
 
  if (!charge && (!geantinoHit)) return false;
  //  G4cout << "\t\t sTGCSD: Hit in a sensitive layer!!!!! " << G4endl;
  G4StepPoint* postStep=aStep->GetPostStepPoint();
  const G4Step* post_Step=aStep->GetTrack()->GetStep();
 
  Amg::Vector3D position = Amg::Hep3VectorToEigen(postStep->GetPosition());
  
  G4StepPoint* preStep = aStep->GetPreStepPoint();
  Amg::Vector3D preposition = Amg::Hep3VectorToEigen(preStep->GetPosition());
 
  int pdgCode=currentTrack->GetDefinition()->GetPDGEncoding();
 
  float globalTime=postStep->GetGlobalTime();
 
  Amg::Vector3D direction= Amg::Hep3VectorToEigen( postStep->GetMomentumDirection() );
  float depositEnergy=post_Step->GetTotalEnergyDeposit();
 
  if (depositEnergy<0.0001 && (!geantinoHit)) return false;
 
  const G4TouchableHistory* touchHist = static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
 
  // int iDepth=touchHist->GetHistoryDepth();
  //  G4cout << "\t\t\t\t Touchable history dump "<<G4endl;
  int nLayer=touchHist->GetVolume(0)->GetCopyNo();
  std::string chName=touchHist->GetVolume(1)->GetLogicalVolume()->GetName();
  //G4cout << "sTGCSensitiveDetector name: "<<chName<<G4endl;
  std::string subType=chName.substr(chName.find('-')+1);
  //G4cout << "\t\t sType: "<<subType);
  if (subType[0]!='T'&&subType[0]!='Q' ) G4cout << " something is wrong, this is no sTGC!"<<G4endl;
  std::string temp(&subType[2]);
  std::istringstream is(temp);
  int iRing;
  is>>iRing;
  // identifiers have eta naming 0-2, eta encoded in subtype is 1-3
  iRing--;
  // double phiDiff=2*M_PI;
 
  G4ThreeVector posH=postStep->GetPosition(); //posH is equivalent to position - eigen not used to avoid additional dependence on EventPrimitives
  if (subType[1]=='L') posH.rotateZ(M_PI/8.);
  double phiHit=posH.phi();
  if(phiHit<=0) phiHit+=2.*M_PI;
  int iPhi=1+int(phiHit/(M_PI/4.));
  iPhi*=2;
  if (subType[1]=='L') iPhi-=1;
 
  int iSide=1;
  if (position.z()<0) iSide=-1;
 
  int mLayer=0;
  if (subType[1]=='S')
  {
    if (subType[3]=='C') mLayer=1;
    else if (subType[3]=='P') mLayer=2;
  }
  else if (subType[1]=='L')
  {
    if (subType[3]=='P') mLayer=1;
    else if (subType[3]=='C') mLayer=2;
  }
 
  if (mLayer != 1 && mLayer !=2) G4cout << " something is wrong - multilayer index is " << mLayer << G4endl;
 
  int sTgcId = m_muonHelper->BuildsTgcHitId(subType, iPhi, iRing, mLayer,nLayer, iSide);
  TrackHelper trHelp(aStep->GetTrack());
  m_sTGCSimHitCollection->Emplace(sTgcId,globalTime,position,pdgCode,direction,depositEnergy,
                                  trHelp.GenerateParticleLink(),
                                  preStep->GetKineticEnergy(),preposition);
 
  return true;
}

Member Data Documentation

m_muonHelper

const sTgcHitIdHelper* sTGCSensitiveDetector::m_muonHelper

private

Definition at line 31 of file sTGCSensitiveDetector.h.

m_sTGCSimHitCollection

SG::WriteHandle<sTGCSimHitCollection> sTGCSensitiveDetector::m_sTGCSimHitCollection

private

Definition at line 30 of file sTGCSensitiveDetector.h.

---

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

• sTGCSensitiveDetector.h
• sTGCSensitiveDetector.cxx