#include <QweakSimCerenkov_PMTSD.hh>

Inherits G4VSensitiveDetector.

Collaboration diagram for QweakSimCerenkov_PMTSD:

Public Member Functions
	QweakSimCerenkov_PMTSD (G4String name, QweakSimUserInformation *myUserInfo)

virtual	~QweakSimCerenkov_PMTSD ()

void	Initialize (G4HCofThisEvent *HCE)

G4bool	ProcessHits (G4Step , G4TouchableHistory )

G4bool	ProcessHits_constStep (const G4Step aStep, G4TouchableHistory ROhist)

void	EndOfEvent (G4HCofThisEvent *HCE)

QweakSimTrajectory *	GetParentTrajectory (G4int parentID)

Private Attributes
QweakSimUserInformation *	myUserInfo

QweakSimCerenkovDetector_PMTHitsCollection *	CerenkovDetector_PMTHitsCollection

G4int	CerenkovDetectorPMT_CollectionID

Detailed Description

Definition at line 15 of file QweakSimCerenkov_PMTSD.hh.

Constructor & Destructor Documentation

QweakSimCerenkov_PMTSD::QweakSimCerenkov_PMTSD	(	G4String	name,
		QweakSimUserInformation *	myUserInfo
	)

Definition at line 17 of file QweakSimCerenkov_PMTSD.cc.

References CerenkovDetectorPMT_CollectionID, and myUserInfo.

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

QweakSimCerenkov_PMTSD::~QweakSimCerenkov_PMTSD ( )

virtual

Definition at line 26 of file QweakSimCerenkov_PMTSD.cc.

 {
 
 }

Member Function Documentation

void QweakSimCerenkov_PMTSD::EndOfEvent ( G4HCofThisEvent * HCE )

Definition at line 84 of file QweakSimCerenkov_PMTSD.cc.

 {
 
 }

QweakSimTrajectory * QweakSimCerenkov_PMTSD::GetParentTrajectory ( G4int parentID )

Definition at line 90 of file QweakSimCerenkov_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;
   }
 //  TrajectoryVector* vect = container->GetVector();
 //
 //    (G4VTrajectory**) tr = vect->begin();
 //
 //   while(tr != vect->end())
 //   {
 //     QweakSimTrajectory* tr1 = (QweakSimTrajectory*)(*tr);
 //     if(tr1->GetTrackID()==parentID) return tr1;
 //     tr++;
 //   }
 
 G4int i = 0;
 
 G4int n_trajectories = trajectoryContainer->entries();
 
 while(i != n_trajectories)
 {
 
     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 QweakSimCerenkov_PMTSD::Initialize ( G4HCofThisEvent * HCE )

Definition at line 32 of file QweakSimCerenkov_PMTSD.cc.

References CerenkovDetector_PMTHitsCollection, and CerenkovDetectorPMT_CollectionID.

 {
    CerenkovDetector_PMTHitsCollection = new QweakSimCerenkovDetector_PMTHitsCollection(SensitiveDetectorName,collectionName[0]);
    
    if(CerenkovDetectorPMT_CollectionID<0)
      { CerenkovDetectorPMT_CollectionID = G4SDManager::GetSDMpointer()->GetCollectionID(CerenkovDetector_PMTHitsCollection); }
    HCE->AddHitsCollection(CerenkovDetectorPMT_CollectionID , CerenkovDetector_PMTHitsCollection);   
 }

G4bool QweakSimCerenkov_PMTSD::ProcessHits	(	G4Step *	,
		G4TouchableHistory *
	)

inline

Definition at line 23 of file QweakSimCerenkov_PMTSD.hh.

                                                      {
     // ProcessHits_constStep(aStep, ROhist);
     return true;
   };

G4bool QweakSimCerenkov_PMTSD::ProcessHits_constStep	(	const G4Step *	aStep,
		G4TouchableHistory *	ROhist
	)

Definition at line 42 of file QweakSimCerenkov_PMTSD.cc.

References CerenkovDetector_PMTHitsCollection, myUserInfo, QweakSimUserInformation::SetCurrentPMTHit(), QweakSimCerenkov_PMTHit::StoreDetectorID(), QweakSimCerenkov_PMTHit::StoreHitID(), QweakSimCerenkov_PMTHit::StoreHitTime(), QweakSimCerenkov_PMTHit::StorePhotonEnergy(), and QweakSimCerenkov_PMTHit::StorePMTID().

Referenced by QweakSimSteppingAction::UserSteppingAction().

 {
 
   if (aStep->GetTrack()->GetDefinition() != G4OpticalPhoton::OpticalPhotonDefinition()) 
     return false;
 
   G4StepPoint*        preStepPoint  = aStep->GetPreStepPoint();
   G4StepPoint*        postStepPoint = aStep->GetPostStepPoint();
 //   G4VPhysicalVolume*  thePrePV      = preStepPoint->GetPhysicalVolume();
 //   G4VPhysicalVolume*  thePostPV     = postStepPoint->GetPhysicalVolume();
 
   G4TouchableHistory* theTouchable  = (G4TouchableHistory*)(preStepPoint->GetTouchable());
 //   G4VPhysicalVolume*  physVol       = theTouchable->GetVolume();
 
   if (postStepPoint->GetStepStatus() != fGeomBoundary) return false;
 
   //if( physVol->GetName().compare("Cathode_Physical")!=0 ) return false;
 
   if (preStepPoint->GetStepStatus() != fGeomBoundary) return false; // Entering Geometry
 
   G4double currentPhotonEnergy = aStep->GetTrack()->GetTotalEnergy();
   G4double currentHitTime      = aStep->GetTrack()->GetGlobalTime();
   G4int MotherCopyNo      = theTouchable->GetVolume(1)->GetCopyNo();
   QweakSimCerenkov_PMTHit* aHit = new QweakSimCerenkov_PMTHit();
   G4int MotherReplicaNo2  = theTouchable->GetReplicaNumber(3);        // Several MotherVolumes
 
 //  std::cout<<"store hit info: detector "<<MotherReplicaNo2<<"  PMTID "<<MotherCopyNo<<std::endl;
   aHit->StoreDetectorID(MotherReplicaNo2); // which octant (number to be converted) 
   aHit->StorePMTID(MotherCopyNo);          // left or right pmt
 
   aHit->StorePhotonEnergy(currentPhotonEnergy);
   aHit->StoreHitTime(currentHitTime);
   G4int hitCount = CerenkovDetector_PMTHitsCollection->insert(aHit); 
   aHit->StoreHitID(hitCount);
   //aHit->SetHitValid(False);
   myUserInfo->SetCurrentPMTHit(aHit,MotherCopyNo);
 
   return true;
 }