#include <GammaRayTelAnticoincidenceSD.hh>

Inheritance diagram for GammaRayTelAnticoincidenceSD:

Public Member Functions
	GammaRayTelAnticoincidenceSD (G4String)

	~GammaRayTelAnticoincidenceSD ()

void	Initialize (G4HCofThisEvent *)

G4bool	ProcessHits (G4Step astep, G4TouchableHistory )

void	EndOfEvent (G4HCofThisEvent *)

void	clear ()

void	DrawAll ()

void	PrintAll ()

Public Member Functions inherited from G4VSensitiveDetector
	G4VSensitiveDetector (G4String name)

	G4VSensitiveDetector (const G4VSensitiveDetector &right)

virtual	~G4VSensitiveDetector ()

const G4VSensitiveDetector &	operator= (const G4VSensitiveDetector &right)

G4int	operator== (const G4VSensitiveDetector &right) const

G4int	operator!= (const G4VSensitiveDetector &right) const

G4bool	Hit (G4Step *aStep)

void	SetROgeometry (G4VReadOutGeometry *value)

void	SetFilter (G4VSDFilter *value)

G4int	GetNumberOfCollections () const

G4String	GetCollectionName (G4int id) const

void	SetVerboseLevel (G4int vl)

void	Activate (G4bool activeFlag)

G4bool	isActive () const

G4String	GetName () const

G4String	GetPathName () const

G4String	GetFullPathName () const

G4VReadOutGeometry *	GetROgeometry () const

G4VSDFilter *	GetFilter () const

virtual G4VSensitiveDetector *	Clone () const

Additional Inherited Members
Protected Member Functions inherited from G4VSensitiveDetector
virtual G4int	GetCollectionID (G4int i)

Protected Attributes inherited from G4VSensitiveDetector
G4CollectionNameVector	collectionName

G4String	SensitiveDetectorName

G4String	thePathName

G4String	fullPathName

G4int	verboseLevel

G4bool	active

G4VReadOutGeometry *	ROgeometry

G4VSDFilter *	filter

Detailed Description

Definition at line 51 of file GammaRayTelAnticoincidenceSD.hh.

Constructor & Destructor Documentation

GammaRayTelAnticoincidenceSD::GammaRayTelAnticoincidenceSD ( G4String name )

Definition at line 52 of file GammaRayTelAnticoincidenceSD.cc.

References G4VSensitiveDetector::collectionName, G4cout, G4endl, GammaRayTelDetectorConstruction::GetNbOfACDLateralTiles(), GammaRayTelDetectorConstruction::GetNbOfACDTopTiles(), G4RunManager::GetRunManager(), G4RunManager::GetUserDetectorConstruction(), and G4CollectionNameVector::insert().

 :G4VSensitiveDetector(name)
 {
   G4RunManager* runManager = G4RunManager::GetRunManager();
   Detector =
     (GammaRayTelDetectorConstruction*)(runManager->GetUserDetectorConstruction());
   
   NbOfACDLateralTiles  =  Detector->GetNbOfACDLateralTiles();
   NbOfACDTopTiles  =  Detector->GetNbOfACDTopTiles(); 
 
   G4cout <<  NbOfACDLateralTiles << " LAT " << G4endl;
   G4cout <<  NbOfACDTopTiles << " TOP " << G4endl;
   
   HitLateralID = new G4int[NbOfACDLateralTiles];
   HitTopID = new G4int[NbOfACDTopTiles];
   collectionName.insert("AnticoincidenceCollection");
 }

GammaRayTelAnticoincidenceSD::~GammaRayTelAnticoincidenceSD ( )

Definition at line 72 of file GammaRayTelAnticoincidenceSD.cc.

 {
   delete [] HitLateralID;
   delete [] HitTopID;
 }

Member Function Documentation

void GammaRayTelAnticoincidenceSD::clear ( void )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 210 of file GammaRayTelAnticoincidenceSD.cc.

211 {}

void GammaRayTelAnticoincidenceSD::DrawAll ( void )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 215 of file GammaRayTelAnticoincidenceSD.cc.

216 {}

void GammaRayTelAnticoincidenceSD::EndOfEvent ( G4HCofThisEvent * HCE )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 186 of file GammaRayTelAnticoincidenceSD.cc.

References G4HCofThisEvent::AddHitsCollection(), G4VSensitiveDetector::collectionName, G4SDManager::GetCollectionID(), and G4SDManager::GetSDMpointer().

 {
   static G4int HCID = -1;
   if(HCID<0)
     { 
       HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
     }
   HCE->AddHitsCollection(HCID,AnticoincidenceCollection);
 
   for (G4int i=0;i<NbOfACDLateralTiles;i++)
     {
       HitLateralID[i] = -1;
     };
   
   for (G4int j=0;j<NbOfACDTopTiles;j++) 
     {
       
       HitTopID[j] = -1;
     };
 
 }

void GammaRayTelAnticoincidenceSD::Initialize ( G4HCofThisEvent * )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 80 of file GammaRayTelAnticoincidenceSD.cc.

References G4VSensitiveDetector::collectionName, and G4VSensitiveDetector::SensitiveDetectorName.

 {
   AnticoincidenceCollection = new GammaRayTelAnticoincidenceHitsCollection
     (SensitiveDetectorName,collectionName[0]);
   for (G4int i=0;i<NbOfACDLateralTiles;i++)
     {
       HitLateralID[i] = -1;
     };
   
   for (G4int j=0;j<NbOfACDTopTiles;j++) 
     {
       
       HitTopID[j] = -1;
     };
 }

void GammaRayTelAnticoincidenceSD::PrintAll ( void )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 220 of file GammaRayTelAnticoincidenceSD.cc.

221 {}

G4bool GammaRayTelAnticoincidenceSD::ProcessHits	(	G4Step *	astep,
		G4TouchableHistory *
	)

virtual

Implements G4VSensitiveDetector.

Definition at line 98 of file GammaRayTelAnticoincidenceSD.cc.

References GammaRayTelAnticoincidenceHit::AddEnergy(), G4VPhysicalVolume::GetCopyNo(), G4VPhysicalVolume::GetName(), G4StepPoint::GetPosition(), G4Step::GetPreStepPoint(), G4Step::GetTotalEnergyDeposit(), G4StepPoint::GetTouchable(), G4TouchableHistory::GetVolume(), G4THitsCollection< T >::insert(), python.hepunit::keV, GammaRayTelAnticoincidenceHit::SetACDTileNumber(), GammaRayTelAnticoincidenceHit::SetACDType(), and GammaRayTelAnticoincidenceHit::SetPos().

 { 
   
   G4double edep = aStep->GetTotalEnergyDeposit();
   if ((edep/keV == 0.)) return false;      
   
   // This TouchableHistory is used to obtain the physical volume
   // of the hit
   G4TouchableHistory* theTouchable
     = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable());
   
   G4VPhysicalVolume* acd_tile = theTouchable->GetVolume();  
 
   G4int ACDTileNumber=acd_tile->GetCopyNo();
   G4String ACDTileName = acd_tile->GetName();
   
   //  G4cout << ACDTileName << " " << edep/keV << G4endl;
 
   if (ACDTileName == "ACT" )
     // The hit is on an top ACD tile (ACDType 0)
     {
       // This is a new hit
       if (HitTopID[ACDTileNumber]==-1)
     {       
       GammaRayTelAnticoincidenceHit* 
         AnticoincidenceHit = new GammaRayTelAnticoincidenceHit;
       AnticoincidenceHit->SetACDType(0);
       AnticoincidenceHit->AddEnergy(edep);
       AnticoincidenceHit->SetPos(aStep->GetPreStepPoint()->GetPosition());
       AnticoincidenceHit->SetACDTileNumber(ACDTileNumber);
       HitTopID[ACDTileNumber] = 
         AnticoincidenceCollection->insert(AnticoincidenceHit) -1;
     }
       else // This is not new
     {
       (*AnticoincidenceCollection)
         [HitTopID[ACDTileNumber]]->AddEnergy(edep);
     }
     }
  
   if (ACDTileName == "ACL1")
     // The hit is on an lateral (left-right) ACD tile (ACDType 1)
     {   
       // This is a new hit
       if (HitLateralID[ACDTileNumber]==-1)
     {       
       GammaRayTelAnticoincidenceHit* 
         AnticoincidenceHit = new GammaRayTelAnticoincidenceHit;
       AnticoincidenceHit->SetACDType(1);
       AnticoincidenceHit->AddEnergy(edep);
       AnticoincidenceHit->SetPos(aStep->GetPreStepPoint()->GetPosition());
       AnticoincidenceHit->SetACDTileNumber(ACDTileNumber);
       HitLateralID[ACDTileNumber] = 
         AnticoincidenceCollection->insert(AnticoincidenceHit) -1;
     }
       else // This is not new
     {
       (*AnticoincidenceCollection)
         [HitLateralID[ACDTileNumber]]->AddEnergy(edep);
     }
     }
 
    if (ACDTileName == "ACL2")
     // The hit is on an lateral (rear - front) ACD tile (ACDType 2)
     {   
       // This is a new hit
       if (HitLateralID[ACDTileNumber]==-1)
     {       
       GammaRayTelAnticoincidenceHit* 
         AnticoincidenceHit = new GammaRayTelAnticoincidenceHit;
       AnticoincidenceHit->SetACDType(2);
       AnticoincidenceHit->AddEnergy(edep);
       AnticoincidenceHit->SetPos(aStep->GetPreStepPoint()->GetPosition());
       AnticoincidenceHit->SetACDTileNumber(ACDTileNumber);
       HitLateralID[ACDTileNumber] = 
         AnticoincidenceCollection->insert(AnticoincidenceHit) -1;
     }
       else // This is not new
     {
       (*AnticoincidenceCollection)
         [HitLateralID[ACDTileNumber]]->AddEnergy(edep);
     }
     }
   return true;
 }

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation