MuonWallSD Class Reference

#include <MuonWallSD.h>

Inheritance diagram for MuonWallSD:

Collaboration diagram for MuonWallSD:

Public Types
using	HitVectorBuilder = TileHitVectorCellBuilder<NCells>

Public Member Functions
	MuonWallSD (const std::string &name, const std::string &hitCollectionName, int verbose)
	~MuonWallSD ()=default
void	Initialize (G4HCofThisEvent *) override final
G4bool	ProcessHits (G4Step *, G4TouchableHistory *) override final

Static Public Attributes
static constexpr int	NCells = 18

Private Member Functions
HitVectorBuilder *	GetHitCollection ()

Private Attributes
const TileTBID *	m_tileTBID {}
Identifier	m_id [NCells]
const std::string	m_hitCollectionName
HitVectorBuilder *	m_hitCollection {}

Static Private Attributes
static const int	s_nCellMu = 14
static const int	s_nCellS = 4

Detailed Description

Definition at line 31 of file MuonWallSD.h.

Member Typedef Documentation

HitVectorBuilder

using MuonWallSD::HitVectorBuilder = TileHitVectorCellBuilder<NCells>

Definition at line 38 of file MuonWallSD.h.

Constructor & Destructor Documentation

MuonWallSD()

MuonWallSD::MuonWallSD	(	const std::string &	name,
		const std::string &	hitCollectionName,
		int	verbose )

Definition at line 29 of file MuonWallSD.cxx.

    : G4VSensitiveDetector(name),
    m_hitCollectionName(hitCollectionName)
{
  verboseLevel = std::max(verboseLevel, verbose);
 
  SmartIF<StoreGateSvc> detStore{Gaudi::svcLocator()->service("DetectorStore")};
  if ( !detStore ) {
    G4ExceptionDescription description;
    description << "Constructor: DetectorStoreSvc not found!";
    G4Exception("MuonWallSD", "NoDetStore", FatalException, description);
    abort();
  } else if (verboseLevel >= 5) {
    G4cout << "DetectorStoreSvc initialized" << G4endl;
  }
 
  if (detStore->retrieve(m_tileTBID).isFailure()) {
    G4ExceptionDescription description;
    description << "Constructor: No TileTBID helper!";
    G4Exception("MuonWallSD", "NoTileTBIDHelper", FatalException, description);
    abort();
  } else if (verboseLevel >= 5) {
    G4cout << "TileTBID helper retrieved" << G4endl;
  }
 
  int type=TileTBID::ADC_TYPE, module=TileTBID::BACK_WALL;
  for (int channel=0; channel<s_nCellMu; ++channel) {
    m_id[channel] = m_tileTBID->channel_id(type,module,channel);
  }
 
  module=TileTBID::CRACK_WALL;
  for (int channel=0; channel<s_nCellS; ++channel) {
    m_id[channel+s_nCellMu] = m_tileTBID->channel_id(type,module,channel);
  }
}

~MuonWallSD()

MuonWallSD::~MuonWallSD ( )

default

Member Function Documentation

GetHitCollection()

MuonWallSD::HitVectorBuilder * MuonWallSD::GetHitCollection ( )

private

Definition at line 65 of file MuonWallSD.cxx.

{
  auto* eventManager = G4EventManager::GetEventManager();
  if (!eventManager) {
    return nullptr;
  }
 
  auto* eventInfo = dynamic_cast<AtlasG4EventUserInfo*>(eventManager->GetUserInformation());
  if (!eventInfo) {
    return nullptr;
  }
 
  auto hitCollections = eventInfo->GetHitCollectionMap();
  return hitCollections ? hitCollections->Find<HitVectorBuilder>(m_hitCollectionName) : nullptr;
}

Initialize()

void MuonWallSD::Initialize ( G4HCofThisEvent * )

finaloverride

Definition at line 81 of file MuonWallSD.cxx.

                                              {
  if (verboseLevel >= 5) {
    G4cout << "MuonWallSD::Initialize()" << G4endl;
  }
 
  m_hitCollection = GetHitCollection();
}

ProcessHits()

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

finaloverride

Definition at line 89 of file MuonWallSD.cxx.

                                                                   {
  if (verboseLevel >= 10) {
    G4cout << "MuonWallSD::ProcessHits" << G4endl;
  }
 
  const G4TouchableHistory* theTouchable = static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
  const G4VPhysicalVolume* physVol = theTouchable->GetVolume();
  const G4LogicalVolume* logiVol = physVol->GetLogicalVolume();
  const G4String nameLogiVol = logiVol->GetName();
  const G4int nScinti = physVol->GetCopyNo();
 
  const G4double edep = aStep->GetTotalEnergyDeposit() * aStep->GetTrack()->GetWeight();
  G4double stepl = 0.;
 
  if (aStep->GetTrack()->GetDefinition()->GetPDGCharge() != 0.){ // FIXME not-equal check on double
 
    stepl = aStep->GetStepLength();
  }
 
  if ((edep == 0.) && (stepl == 0.)) { //FIXME equality check on double
 
    return false;
  }
 
  int ind;
 
  if(nameLogiVol.find("MuScintillatorLayer") !=G4String::npos) {
    // for muon wall, nScinti-1 is the correct indice.
    ind = nScinti-1;
  } else if(nameLogiVol.find("S1") !=G4String::npos) {
    ind = s_nCellMu + 0;
  } else if(nameLogiVol.find("S2") !=G4String::npos) {
    ind = s_nCellMu + 1;
  } else if(nameLogiVol.find("S3") !=G4String::npos) {
    ind = s_nCellMu + 2;
  } else {
    ind = s_nCellMu + 3;
  }
 
  if (verboseLevel >= 10) {
    HitVectorBuilder* hitCollection = m_hitCollection ? m_hitCollection : GetHitCollection();
    G4cout << ((hitCollection && hitCollection->HasHit(ind))?"Additional hit in ":"First hit in ")
           << ((ind<s_nCellMu)?"MuonWall ":"beam counter S")
           << ((ind<s_nCellMu)?nScinti:(ind-s_nCellMu+1))
           << " time=" << aStep->GetPostStepPoint()->GetGlobalTime()
           << " ene=" << edep << G4endl;
  }
 
  HitVectorBuilder* hitCollection = m_hitCollection ? m_hitCollection : GetHitCollection();
  if (!hitCollection) {
    if (verboseLevel >= 5) {
      G4cout << "MuonWallSD::ProcessHits WARNING hit collection is not available" << G4endl;
    }
    return false;
  }
  m_hitCollection = hitCollection;
  hitCollection->AddHit(ind, m_id[ind], edep);
 
  return true;
}

Member Data Documentation

m_hitCollection

HitVectorBuilder* MuonWallSD::m_hitCollection {}

private

Definition at line 53 of file MuonWallSD.h.

{};

m_hitCollectionName

const std::string MuonWallSD::m_hitCollectionName

private

Definition at line 52 of file MuonWallSD.h.

m_id

Identifier MuonWallSD::m_id[NCells]

private

Definition at line 51 of file MuonWallSD.h.

m_tileTBID

const TileTBID* MuonWallSD::m_tileTBID {}

private

Definition at line 44 of file MuonWallSD.h.

{};

NCells

int MuonWallSD::NCells = 18

staticconstexpr

Definition at line 37 of file MuonWallSD.h.

s_nCellMu

const int MuonWallSD::s_nCellMu = 14

staticprivate

Definition at line 46 of file MuonWallSD.h.

s_nCellS

const int MuonWallSD::s_nCellS = 4

staticprivate

Definition at line 47 of file MuonWallSD.h.

---

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

• MuonWallSD.h
• MuonWallSD.cxx