#include <QweakSimPMTOnly_PMTSD.hh>

Inherits G4VSensitiveDetector.

Collaboration diagram for QweakSimPMTOnly_PMTSD:

Public Member Functions
	QweakSimPMTOnly_PMTSD (G4String name, QweakSimUserInformation *userInfo)

	~QweakSimPMTOnly_PMTSD ()

void	Initialize (G4HCofThisEvent *HCE)

virtual G4bool	ProcessHits (G4Step aStep, G4TouchableHistory )

G4bool	ProcessHits_constStep (const G4Step , G4TouchableHistory )

void	EndOfEvent (G4HCofThisEvent *HCE)

QweakSimTrajectory *	GetParentTrajectory (G4int parentID)

Private Attributes
QweakSimUserInformation *	myUserInfo

QweakSimPMTOnly_PMTHitsCollection *	PMTOnly_PMTHitsCollection

G4int	PMTOnlyPMT_CollectionID

Detailed Description

Definition at line 20 of file QweakSimPMTOnly_PMTSD.hh.

Constructor & Destructor Documentation

QweakSimPMTOnly_PMTSD::QweakSimPMTOnly_PMTSD	(	G4String	name,
		QweakSimUserInformation *	userInfo
	)

Definition at line 20 of file QweakSimPMTOnly_PMTSD.cc.

References myUserInfo, and PMTOnlyPMT_CollectionID.

 :G4VSensitiveDetector(name)
 {
   collectionName.insert("PMTHitCollection"); 
   PMTOnlyPMT_CollectionID = -1;
   myUserInfo = userInfo; 
 }

QweakSimPMTOnly_PMTSD::~QweakSimPMTOnly_PMTSD ( )

Definition at line 29 of file QweakSimPMTOnly_PMTSD.cc.

30 {

31 }

Member Function Documentation

void QweakSimPMTOnly_PMTSD::EndOfEvent ( G4HCofThisEvent * HCE )

Definition at line 71 of file QweakSimPMTOnly_PMTSD.cc.

72 {

73 }

QweakSimTrajectory * QweakSimPMTOnly_PMTSD::GetParentTrajectory ( G4int parentID )

Definition at line 76 of file QweakSimPMTOnly_PMTSD.cc.

References QweakSimTrajectory::GetCharge(), and QweakSimTrajectory::GetTrackID().

 {
   G4TrajectoryContainer* trajectoryContainer = 
    G4RunManager::GetRunManager()->GetCurrentEvent()->GetTrajectoryContainer();
 
   if(trajectoryContainer==0) {
       G4cout << " Could not find the TrajectoryContainer" << G4endl;
       return 0;
   }
 
   G4int i = 0;
   G4int n_entries = trajectoryContainer->entries();
   while( i < n_entries) {
     QweakSimTrajectory* trj = (QweakSimTrajectory*)((*trajectoryContainer)[i]);
     G4cout << " Current charge of possible parent track = " <<trj->GetCharge() << "--- Parent ID = "
      <<  trj->GetTrackID() << G4endl;
     if(trj->GetTrackID()==parentID) return trj;
     i++;
   }
 
   return 0;
 }

Here is the call graph for this function:

void QweakSimPMTOnly_PMTSD::Initialize ( G4HCofThisEvent * HCE )

Definition at line 34 of file QweakSimPMTOnly_PMTSD.cc.

References PMTOnly_PMTHitsCollection, and PMTOnlyPMT_CollectionID.

 {
    PMTOnly_PMTHitsCollection = new QweakSimPMTOnly_PMTHitsCollection(SensitiveDetectorName,collectionName[0]);
    
    if(PMTOnlyPMT_CollectionID<0)
      { PMTOnlyPMT_CollectionID = G4SDManager::GetSDMpointer()->GetCollectionID(PMTOnly_PMTHitsCollection); }
    HCE->AddHitsCollection(PMTOnlyPMT_CollectionID , PMTOnly_PMTHitsCollection);   
 }

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

virtual

Definition at line 43 of file QweakSimPMTOnly_PMTSD.cc.

                                                                       {
   return false;
 }

G4bool QweakSimPMTOnly_PMTSD::ProcessHits_constStep	(	const G4Step *	aStep,
		G4TouchableHistory *
	)

Definition at line 48 of file QweakSimPMTOnly_PMTSD.cc.

References PMTOnly_PMTHitsCollection, QweakSimPMTOnly_PMTHit::StoreHitID(), and QweakSimPMTOnly_PMTHit::StorePhotonEnergy().

Referenced by QweakSimSteppingAction::UserSteppingAction().

 {
   // Make sure that this is a 
   if (aStep->GetTrack()->GetDefinition() != G4OpticalPhoton::OpticalPhotonDefinition()) 
     return false;
 
   G4StepPoint*        preStepPoint  = aStep->GetPreStepPoint();
   G4StepPoint*        postStepPoint = aStep->GetPostStepPoint();
 
   if (postStepPoint->GetStepStatus() != fGeomBoundary) return false;
   if (preStepPoint->GetStepStatus() != fGeomBoundary) return false; // Entering Volume
 
   G4double currentPhotonEnergy = aStep->GetTrack()->GetTotalEnergy();
   QweakSimPMTOnly_PMTHit* aHit = new QweakSimPMTOnly_PMTHit();
 
   aHit->StorePhotonEnergy(currentPhotonEnergy);
   G4int hitCount = PMTOnly_PMTHitsCollection->insert(aHit); 
   aHit->StoreHitID(hitCount);
 
   return true;
 }