SensitiveDetector Class Reference

#include <SensitiveDetector.h>

Inheritance diagram for SensitiveDetector:

Public Types
enum	coordinates { WorldCoordinates, LocalCoordinates }

Public Member Functions
virtual void	AssignSD (const std::string &vname)
virtual void	clearHits ()=0
Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point)
void	EndOfEvent (G4HCofThisEvent *eventHC) override
Local3DPoint	FinalStepPosition (G4Step *s, coordinates)
virtual std::vector< std::string >	getNames ()
void	Initialize (G4HCofThisEvent *eventHC) override
Local3DPoint	InitialStepPosition (G4Step *s, coordinates)
std::string	nameOfSD ()
void	NaNTrap (G4Step *step)
G4bool	ProcessHits (G4Step *step, G4TouchableHistory *tHistory) override=0
void	Register ()
	SensitiveDetector (std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &, edm::ParameterSet const &p)
virtual uint32_t	setDetUnitId (G4Step *step)=0
	~SensitiveDetector () override

Private Attributes
G4Step *	currentStep
std::string	name

Detailed Description

Definition at line 22 of file SensitiveDetector.h.

Member Enumeration Documentation

enum SensitiveDetector::coordinates

Enumerator
WorldCoordinates
LocalCoordinates

Definition at line 36 of file SensitiveDetector.h.

{WorldCoordinates, LocalCoordinates};

Constructor & Destructor Documentation

SensitiveDetector::SensitiveDetector ( std::string &  iname, const DDCompactView &  cpv, const SensitiveDetectorCatalog &  clg, edm::ParameterSet const &  p  )

explicit

Definition at line 14 of file SensitiveDetector.cc.

                                                :
  G4VSensitiveDetector(iname), name(iname) {}

SensitiveDetector::~SensitiveDetector ( )

override

Definition at line 20 of file SensitiveDetector.cc.

{}

Member Function Documentation

void SensitiveDetector::AssignSD ( const std::string &  vname )

virtual

Definition at line 30 of file SensitiveDetector.cc.

References findQualityFiles::v.

Referenced by Bcm1fSD::Bcm1fSD(), BHMSD::BHMSD(), BscSD::BscSD(), CaloSD::CaloSD(), FastTimerSD::FastTimerSD(), FiberSD::FiberSD(), FP420SD::FP420SD(), HFChamberSD::HFChamberSD(), HFWedgeSD::HFWedgeSD(), MuonSensitiveDetector::MuonSensitiveDetector(), PltSD::PltSD(), TkAccumulatingSensitiveDetector::TkAccumulatingSensitiveDetector(), and TotemSD::TotemSD().

{
  G4LogicalVolumeStore * theStore = G4LogicalVolumeStore::GetInstance();
  G4LogicalVolumeStore::const_iterator it;
  for (it = theStore->begin(); it != theStore->end(); it++)
    {
      G4LogicalVolume * v = *it;
      if (vname==v->GetName()) { v->SetSensitiveDetector(this); }
    }
}

virtual void SensitiveDetector::clearHits ( )

pure virtual

Implemented in FP420SD, CaloSD, BscSD, MuonSensitiveDetector, TotemSD, FastTimerSD, BHMSD, TkAccumulatingSensitiveDetector, Bcm1fSD, PltSD, FiberSD, HFWedgeSD, HFChamberSD, and CaloTrkProcessing.

Local3DPoint SensitiveDetector::ConvertToLocal3DPoint

(

const G4ThreeVector &

point

)

Definition at line 72 of file SensitiveDetector.cc.

Referenced by Bcm1fSD::createHit(), PltSD::createHit(), TkAccumulatingSensitiveDetector::createHit(), MuonSensitiveDetector::createHit(), FinalStepPosition(), MuonSensitiveDetector::FinalStepPositionVsParent(), InitialStepPosition(), and MuonSensitiveDetector::InitialStepPositionVsParent().

{
    return Local3DPoint(p.x(),p.y(),p.z());
}

void SensitiveDetector::EndOfEvent ( G4HCofThisEvent *  eventHC )

override

Definition at line 41 of file SensitiveDetector.cc.

{}

Local3DPoint SensitiveDetector::FinalStepPosition	(	G4Step *	s,
		coordinates	c
	)

Definition at line 58 of file SensitiveDetector.cc.

References ConvertToLocal3DPoint(), currentStep, alignCSCRings::s, and WorldCoordinates.

Referenced by Bcm1fSD::createHit(), PltSD::createHit(), TkAccumulatingSensitiveDetector::createHit(), MuonSensitiveDetector::createHit(), Bcm1fSD::updateHit(), PltSD::updateHit(), TkAccumulatingSensitiveDetector::updateHit(), and MuonSensitiveDetector::updateHit().

{
    currentStep = s;
    G4StepPoint * postStepPoint = currentStep->GetPostStepPoint();
    G4StepPoint * preStepPoint  = currentStep->GetPreStepPoint();
    const G4ThreeVector& globalCoordinates = postStepPoint->GetPosition();
    if (c == WorldCoordinates) return ConvertToLocal3DPoint(globalCoordinates);
    G4TouchableHistory * theTouchable = (G4TouchableHistory *)
                                        (preStepPoint->GetTouchable());
    G4ThreeVector localCoordinates = theTouchable->GetHistory()
                  ->GetTopTransform().TransformPoint(globalCoordinates);
    return ConvertToLocal3DPoint(localCoordinates); 
}

virtual std::vector<std::string> SensitiveDetector::getNames ( )

inlinevirtual

Reimplemented in FP420SD, BscSD, MuonSensitiveDetector, BHMSD, FastTimerSD, and TkAccumulatingSensitiveDetector.

Definition at line 41 of file SensitiveDetector.h.

References nameOfSD(), NaNTrap(), and groupFilesInBlocks::temp.

Referenced by CaloTrkProcessing::CaloTrkProcessing(), CaloTrkProcessing::fillHits(), HCalSD::HCalSD(), and HcalTB06BeamSD::HcalTB06BeamSD().

  {
    std::vector<std::string> temp;
    temp.push_back(nameOfSD());
    return temp;
  }

void SensitiveDetector::Initialize ( G4HCofThisEvent *  eventHC )

override

Definition at line 22 of file SensitiveDetector.cc.

{}

Local3DPoint SensitiveDetector::InitialStepPosition	(	G4Step *	s,
		coordinates	c
	)

Definition at line 45 of file SensitiveDetector.cc.

References ConvertToLocal3DPoint(), currentStep, alignCSCRings::s, and WorldCoordinates.

Referenced by Bcm1fSD::closeHit(), PltSD::closeHit(), TkAccumulatingSensitiveDetector::closeHit(), Bcm1fSD::createHit(), PltSD::createHit(), TkAccumulatingSensitiveDetector::createHit(), MuonSensitiveDetector::createHit(), and MuonSensitiveDetector::ProcessHits().

{
    currentStep = s;
    G4StepPoint * preStepPoint = currentStep->GetPreStepPoint();
    const G4ThreeVector& globalCoordinates = preStepPoint->GetPosition();
    if (c == WorldCoordinates) return ConvertToLocal3DPoint(globalCoordinates);
    G4TouchableHistory * theTouchable=(G4TouchableHistory *)
                                      (preStepPoint->GetTouchable());
    G4ThreeVector localCoordinates = theTouchable->GetHistory()
                  ->GetTopTransform().TransformPoint(globalCoordinates);
    return ConvertToLocal3DPoint(localCoordinates); 
}

std::string SensitiveDetector::nameOfSD ( )

inline

Definition at line 40 of file SensitiveDetector.h.

References name.

Referenced by getNames(), and HGCSD::HGCSD().

{ return name; }

void SensitiveDetector::NaNTrap

(

G4Step *

step

)

Definition at line 77 of file SensitiveDetector.cc.

References edm::isNotFinite(), and reco::return().

Referenced by getNames(), ZdcSD::ProcessHits(), HGCSD::ProcessHits(), HCalSD::ProcessHits(), and CaloSD::ProcessHits().

{

    if ( aStep == nullptr ) return ;
    
    G4Track* CurrentTrk = aStep->GetTrack() ;
    const G4ThreeVector& CurrentPos = CurrentTrk->GetPosition() ;
    G4ThreeVector CurrentMom = CurrentTrk->GetMomentum() ;
    G4VPhysicalVolume* pCurrentVol = CurrentTrk->GetVolume() ;
    G4String NameOfVol ;
    if ( pCurrentVol != nullptr )
    {
       NameOfVol = pCurrentVol->GetName() ;
    }
    else
    {
       NameOfVol = "CorruptedVolumeInfo" ;
    }
    
    // for simplicity... maybe edm::isNotFinite() will work on the 
    // 3-vector directly...

    double xyz[3] ;
    xyz[0] = CurrentPos.x() ;
    xyz[1] = CurrentPos.y() ;
    xyz[2] = CurrentPos.z() ;
    
    //
    // this is another trick to check on a NaN, maybe it's even CPU-faster...
    // but ler's stick to system function edm::isNotFinite(...) for now
    //
    // if ( !(xyz[0]==xyz[0]) || !(xyz[1]==xyz[1]) || !(xyz[2]==xyz[2]) )
    if( edm::isNotFinite(xyz[0]+xyz[1]+xyz[2]) != 0 )
    {
      throw SimG4Exception( "SimG4CoreSensitiveDetector: Corrupted Event - NaN detected (position) in volume " + NameOfVol);
    }

    xyz[0] = CurrentMom.x() ;
    xyz[1] = CurrentMom.y() ;
    xyz[2] = CurrentMom.z() ;
    if ( !(xyz[0]==xyz[0]) || !(xyz[1]==xyz[1]) || !(xyz[2]==xyz[2]) ||
         edm::isNotFinite(xyz[0]) != 0 || edm::isNotFinite(xyz[1]) != 0 || 
     edm::isNotFinite(xyz[2]) != 0 )
    {
      throw SimG4Exception( "SimG4CoreSensitiveDetector: Corrupted Event - NaN detected (3-momentum) in volume " + NameOfVol);
    }

   return;

}

G4bool SensitiveDetector::ProcessHits ( G4Step *  step, G4TouchableHistory *  tHistory  )

overridepure virtual

Implemented in FP420SD, BscSD, CaloSD, MuonSensitiveDetector, TotemSD, BHMSD, FastTimerSD, HCalSD, Bcm1fSD, PltSD, TkAccumulatingSensitiveDetector, FiberSD, CaloTrkProcessing, HGCSD, HFWedgeSD, HFChamberSD, DreamSD, and ZdcSD.

void SensitiveDetector::Register

(

)

Definition at line 24 of file SensitiveDetector.cc.

Referenced by Bcm1fSD::Bcm1fSD(), BHMSD::BHMSD(), BscSD::BscSD(), CaloSD::CaloSD(), FastTimerSD::FastTimerSD(), FiberSD::FiberSD(), FP420SD::FP420SD(), HFChamberSD::HFChamberSD(), HFWedgeSD::HFWedgeSD(), MuonSensitiveDetector::MuonSensitiveDetector(), PltSD::PltSD(), TkAccumulatingSensitiveDetector::TkAccumulatingSensitiveDetector(), and TotemSD::TotemSD().

{
  G4SDManager * SDman = G4SDManager::GetSDMpointer();
  SDman->AddNewDetector(this);
}

virtual uint32_t SensitiveDetector::setDetUnitId ( G4Step *  step )

pure virtual

Implemented in FP420SD, BscSD, CaloSD, MuonSensitiveDetector, TotemSD, BHMSD, FastTimerSD, FiberSD, HCalSD, HFWedgeSD, Bcm1fSD, PltSD, TkAccumulatingSensitiveDetector, ECalSD, CastorSD, CaloTrkProcessing, HFChamberSD, HcalTB02SD, HGCSD, HcalTB06BeamSD, EcalTBH4BeamSD, HGCalTB16SD01, DreamSD, ZdcSD, and AHCalSD.

Member Data Documentation

G4Step* SensitiveDetector::currentStep

private

Definition at line 52 of file SensitiveDetector.h.

Referenced by FinalStepPosition(), and InitialStepPosition().

std::string SensitiveDetector::name

private

Definition at line 51 of file SensitiveDetector.h.

Referenced by dirstructure.Directory::__create_pie_image(), dqm_interfaces.DirID::__eq__(), dirstructure.Directory::__get_full_path(), dirstructure.Comparison::__get_img_name(), dirstructure.Comparison::__make_image(), dirstructure.Directory::__repr__(), dqm_interfaces.DirID::__repr__(), dirstructure.Comparison::__repr__(), config.CFG::__str__(), dirstructure.Directory::calcStats(), CaloTrkProcessing::CaloTrkProcessing(), validation.Sample::digest(), utils.StatisticalTest::get_status(), CastorSD::getEnergyDeposit(), HcalTB06BeamSD::getNames(), HCalSD::HCalSD(), VIDSelectorBase.VIDSelectorBase::initialize(), DreamSD::initMap(), HcalTB02SD::initMap(), ECalSD::initMap(), nameOfSD(), HCalSD::plotProfile(), dirstructure.Directory::print_report(), SensitiveCaloDetector::SensitiveCaloDetector(), SensitiveTkDetector::SensitiveTkDetector(), HGCalTB16SD01::setDetUnitId(), HcalTB06BeamSD::setDetUnitId(), and Vispa.Views.PropertyView.Property::valueChanged().