SensitiveDetector Class Reference

#include <SensitiveDetector.h>

Inheritance diagram for SensitiveDetector:

Public Types
enum	coordinates { WorldCoordinates, LocalCoordinates }

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

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, SensitiveDetectorCatalog &  clg, edm::ParameterSet const &  p  )

explicit

Definition at line 14 of file SensitiveDetector.cc.

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

SensitiveDetector::~SensitiveDetector ( )

virtual

Definition at line 20 of file SensitiveDetector.cc.

{}

Member Function Documentation

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

virtual

Definition at line 30 of file SensitiveDetector.cc.

References findQualityFiles::v.

Referenced by BscSD::BscSD(), CaloSD::CaloSD(), FiberSD::FiberSD(), FP420SD::FP420SD(), HFChamberSD::HFChamberSD(), HFWedgeSD::HFWedgeSD(), MuonSensitiveDetector::MuonSensitiveDetector(), PLTSensitiveDetector::PLTSensitiveDetector(), 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, TkAccumulatingSensitiveDetector, PLTSensitiveDetector, FiberSD, HFWedgeSD, HFChamberSD, and CaloTrkProcessing.

Local3DPoint SensitiveDetector::ConvertToLocal3DPoint

(

const G4ThreeVector &

point

)

Definition at line 72 of file SensitiveDetector.cc.

Referenced by PLTSensitiveDetector::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 )

virtual

Reimplemented in FP420SD, BscSD, CaloSD, MuonSensitiveDetector, TotemSD, PLTSensitiveDetector, TkAccumulatingSensitiveDetector, CaloTrkProcessing, FiberSD, HFWedgeSD, and HFChamberSD.

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 PLTSensitiveDetector::createHit(), TkAccumulatingSensitiveDetector::createHit(), MuonSensitiveDetector::createHit(), PLTSensitiveDetector::updateHit(), TkAccumulatingSensitiveDetector::updateHit(), and MuonSensitiveDetector::updateHit().

{
    currentStep = s;
    G4StepPoint * postStepPoint = currentStep->GetPostStepPoint();
    G4StepPoint * preStepPoint  = currentStep->GetPreStepPoint();
    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, and TkAccumulatingSensitiveDetector.

Definition at line 41 of file SensitiveDetector.h.

References nameOfSD(), and groupFilesInBlocks::temp.

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

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

void SensitiveDetector::Initialize ( G4HCofThisEvent *  eventHC )

virtual

Reimplemented in FP420SD, BscSD, CaloSD, TotemSD, FiberSD, CaloTrkProcessing, HFWedgeSD, and HFChamberSD.

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 PLTSensitiveDetector::closeHit(), TkAccumulatingSensitiveDetector::closeHit(), PLTSensitiveDetector::createHit(), TkAccumulatingSensitiveDetector::createHit(), MuonSensitiveDetector::createHit(), and MuonSensitiveDetector::ProcessHits().

{
    currentStep = s;
    G4StepPoint * preStepPoint = currentStep->GetPreStepPoint();
    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().

{ return name; }

void SensitiveDetector::NaNTrap

(

G4Step *

step

)

Definition at line 77 of file SensitiveDetector.cc.

References gather_cfg::cout, edm::isNotFinite(), NULL, and reco::return().

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

{

    if ( aStep == NULL ) return ;
    
    G4Track* CurrentTrk = aStep->GetTrack() ;
    G4ThreeVector CurrentPos = CurrentTrk->GetPosition() ;
    G4ThreeVector CurrentMom = CurrentTrk->GetMomentum() ;
    G4VPhysicalVolume* pCurrentVol = CurrentTrk->GetVolume() ;
    G4String NameOfVol ;
    if ( pCurrentVol != NULL )
    {
       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 )
    {
       // std::cout << " NaN detected in volume " << NameOfVol << std::endl ;
       throw SimG4Exception( "SimG4CoreSensitiveDetector: Corrupted Event - NaN detected (position)" );
    }

    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 )
    {
       std::cout << " NaN detected in volume " << NameOfVol << std::endl ;
       throw SimG4Exception( "SimG4CoreSensitiveDetector: Corrupted Event - NaN detected (3-momentum)" ) ;
    }

   return;

}

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

pure virtual

Implemented in FP420SD, BscSD, CaloSD, MuonSensitiveDetector, TotemSD, PLTSensitiveDetector, TkAccumulatingSensitiveDetector, HCalSD, FiberSD, CaloTrkProcessing, HFWedgeSD, HFChamberSD, DreamSD, and ZdcSD.

void SensitiveDetector::Register

(

)

Definition at line 24 of file SensitiveDetector.cc.

Referenced by BscSD::BscSD(), CaloSD::CaloSD(), FiberSD::FiberSD(), FP420SD::FP420SD(), HFChamberSD::HFChamberSD(), HFWedgeSD::HFWedgeSD(), MuonSensitiveDetector::MuonSensitiveDetector(), PLTSensitiveDetector::PLTSensitiveDetector(), 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, FiberSD, HFWedgeSD, PLTSensitiveDetector, TkAccumulatingSensitiveDetector, HCalSD, CastorSD, CaloTrkProcessing, HFChamberSD, HcalTB02SD, ECalSD, EcalTBH4BeamSD, HcalTB06BeamSD, DreamSD, and ZdcSD.

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(), dataset.Dataset::__getFileInfoList(), cuy.divideElement::__init__(), cuy.plotElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), cuy.graphElement::__init__(), dirstructure.Comparison::__make_image(), dirstructure.Directory::__repr__(), dqm_interfaces.DirID::__repr__(), dirstructure.Comparison::__repr__(), dirstructure.Directory::calcStats(), CaloTrkProcessing::CaloTrkProcessing(), utils.StatisticalTest::get_status(), CastorSD::getEnergyDeposit(), HcalTB06BeamSD::getNames(), HCalSD::HCalSD(), DreamSD::initMap(), HcalTB02SD::initMap(), ECalSD::initMap(), nameOfSD(), HCalSD::plotProfile(), dirstructure.Directory::print_report(), HcalTB06BeamSD::setDetUnitId(), and Vispa.Views.PropertyView.Property::valueChanged().