GammaRayTelAnticoincidenceSD Class Reference

#include <GammaRayTelAnticoincidenceSD.hh>

Inheritance diagram for GammaRayTelAnticoincidenceSD:

Collaboration 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 ()
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.

: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");
}

Here is the call graph for this function:

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.

{} 

void GammaRayTelAnticoincidenceSD::DrawAll ( void  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 215 of file GammaRayTelAnticoincidenceSD.cc.

{} 

void GammaRayTelAnticoincidenceSD::EndOfEvent ( G4HCofThisEvent *  HCE )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 186 of file GammaRayTelAnticoincidenceSD.cc.

{
  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;
    };

}

Here is the call graph for this function:

void GammaRayTelAnticoincidenceSD::Initialize ( G4HCofThisEvent *  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 80 of file GammaRayTelAnticoincidenceSD.cc.

{
  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.

{} 

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

virtual

Implements G4VSensitiveDetector.

Definition at line 98 of file GammaRayTelAnticoincidenceSD.cc.

{ 
  
  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;
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/examples/advanced/gammaray_telescope/include/GammaRayTelAnticoincidenceSD.hh
• source/geant4.10.03.p03/examples/advanced/gammaray_telescope/src/GammaRayTelAnticoincidenceSD.cc