GammaRayTelTrackerSD Class Reference

#include <GammaRayTelTrackerSD.hh>

Inheritance diagram for GammaRayTelTrackerSD:

Public Member Functions
	GammaRayTelTrackerSD (G4String)
	~GammaRayTelTrackerSD ()
void	Initialize (G4HCofThisEvent *)
G4bool	ProcessHits (G4Step *astep, G4TouchableHistory *ROHist)
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 GammaRayTelTrackerSD.hh.

Constructor & Destructor Documentation

GammaRayTelTrackerSD::GammaRayTelTrackerSD

(

G4String

name

)

Definition at line 54 of file GammaRayTelTrackerSD.cc.

References G4VSensitiveDetector::collectionName, GammaRayTelDetectorConstruction::GetNbOfTKRLayers(), GammaRayTelDetectorConstruction::GetNbOfTKRStrips(), GammaRayTelDetectorConstruction::GetNbOfTKRTiles(), G4RunManager::GetRunManager(), G4RunManager::GetUserDetectorConstruction(), and G4CollectionNameVector::insert().

                                                       :G4VSensitiveDetector(name)
{
 G4RunManager* runManager = G4RunManager::GetRunManager();
  Detector =
    (GammaRayTelDetectorConstruction*)(runManager->GetUserDetectorConstruction());
  
  G4int NbOfTKRTiles  =  Detector->GetNbOfTKRTiles();
  NbOfTKRStrips  = Detector->GetNbOfTKRStrips();
  NbOfTKRLayers  = Detector->GetNbOfTKRLayers();  
  NbOfTKRStrips = NbOfTKRStrips*NbOfTKRTiles;  
  
  NbOfTKRChannels = NbOfTKRStrips* NbOfTKRTiles * NbOfTKRLayers;
  
  ThitXID = new G4int[NbOfTKRChannels];
  ThitYID = new G4int[NbOfTKRChannels];
  collectionName.insert("TrackerCollection");
}

GammaRayTelTrackerSD::~GammaRayTelTrackerSD

(

)

Definition at line 74 of file GammaRayTelTrackerSD.cc.

{
  delete [] ThitXID;
  delete [] ThitYID;
}

Member Function Documentation

void GammaRayTelTrackerSD::clear ( void  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 219 of file GammaRayTelTrackerSD.cc.

{} 

void GammaRayTelTrackerSD::DrawAll ( void  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 224 of file GammaRayTelTrackerSD.cc.

{} 

void GammaRayTelTrackerSD::EndOfEvent ( G4HCofThisEvent *  HCE )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 200 of file GammaRayTelTrackerSD.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,TrackerCollection);


  for (G4int i=0;i<NbOfTKRChannels;i++) 
    {
      ThitXID[i] = -1;
      ThitYID[i] = -1;
    };
}

void GammaRayTelTrackerSD::Initialize ( G4HCofThisEvent *  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 82 of file GammaRayTelTrackerSD.cc.

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

{
  TrackerCollection = new GammaRayTelTrackerHitsCollection
    (SensitiveDetectorName,collectionName[0]);

 for (G4int i=0;i<NbOfTKRChannels;i++)
   {
     ThitXID[i] = -1;
     ThitYID[i] = -1;
   };
}

void GammaRayTelTrackerSD::PrintAll ( void  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 229 of file GammaRayTelTrackerSD.cc.

{} 

G4bool GammaRayTelTrackerSD::ProcessHits ( G4Step *  astep, G4TouchableHistory *  ROHist  )

virtual

Implements G4VSensitiveDetector.

Definition at line 96 of file GammaRayTelTrackerSD.cc.

References GammaRayTelTrackerHit::AddSil(), G4VPhysicalVolume::GetCopyNo(), G4VPhysicalVolume::GetName(), GammaRayTelDetectorConstruction::GetNbOfTKRStrips(), GammaRayTelDetectorConstruction::GetNbOfTKRTiles(), G4StepPoint::GetPosition(), G4Step::GetPreStepPoint(), G4Step::GetTotalEnergyDeposit(), G4StepPoint::GetTouchable(), G4TouchableHistory::GetVolume(), G4THitsCollection< T >::insert(), python.hepunit::keV, GammaRayTelTrackerHit::SetNSilPlane(), GammaRayTelTrackerHit::SetNStrip(), GammaRayTelTrackerHit::SetPlaneType(), and GammaRayTelTrackerHit::SetPos().

{ 
   
  G4double edep = aStep->GetTotalEnergyDeposit();
  if ((edep/keV == 0.)) return false;      
  
  G4int StripTotal = Detector->GetNbOfTKRStrips();
  G4int TileTotal  = Detector->GetNbOfTKRTiles();  

  // This TouchableHistory is used to obtain the physical volume
  // of the hit
  G4TouchableHistory* theTouchable
    = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable());
  
  //G4VPhysicalVolume* phys_tile = theTouchable->GetVolume();  
  
  G4VPhysicalVolume* plane = theTouchable->GetVolume(1);  
  
  G4int PlaneNumber = 0;
  PlaneNumber=plane->GetCopyNo();
  G4String PlaneName = plane->GetName();

  // The RO History is used to obtain the real strip
  // of the hit

  G4int StripNumber = 0;
  G4VPhysicalVolume* strip = 0;
  strip = ROhist->GetVolume();

  //  if (strip){
    G4String StripName = strip->GetName();
    StripNumber= strip->GetCopyNo();  
    //  }
  //  ROhist->MoveUpHistory();
  G4VPhysicalVolume* tile = ROhist->GetVolume(1); 
  G4int TileNumber = tile->GetCopyNo();  
  G4String TileName = tile->GetName();   
  
  G4int NTile = (TileNumber%TileTotal);  
  G4int j=0;
  
  G4int NChannel = 0;
  
  for (j=0;j<TileTotal;j++)
    {
      if(NTile==j) StripNumber += StripTotal*NTile;
    }  
  
  NChannel = PlaneNumber*TileTotal*StripTotal + StripNumber;

  /*  G4cout << NChannel << " Channel Number" << G4endl;
      G4cout << " Plane Number = " << PlaneNumber << " " << PlaneName 
      << G4endl;
      G4cout << StripName << " " << StripNumber << G4endl;      */ 
  
  if (PlaneName == "TKRDetectorX" )
    // The hit is on an X silicon plane
    {
      // This is a new hit
      if (ThitXID[NChannel]==-1)
    {       
      GammaRayTelTrackerHit* TrackerHit = new GammaRayTelTrackerHit;
      TrackerHit->SetPlaneType(1);
      TrackerHit->AddSil(edep);
      TrackerHit->SetPos(aStep->GetPreStepPoint()->GetPosition());
      TrackerHit->SetNSilPlane(PlaneNumber);
      TrackerHit->SetNStrip(StripNumber);
      ThitXID[NChannel] = 
        TrackerCollection->insert(TrackerHit) -1;
    }
      else // This is not new
    {
      (*TrackerCollection)[ThitXID[NChannel]]->AddSil(edep);
          // G4cout << "X" << PlaneNumber << " " << StripNumber << G4endl;
    }
    }
  
  if (PlaneName == "TKRDetectorY")
    // The hit is on an Y silicon plane    
    {   
      // This is a new hit
      if (ThitYID[NChannel]==-1)
    {       
      GammaRayTelTrackerHit* TrackerHit = new GammaRayTelTrackerHit;
      TrackerHit->SetPlaneType(0);
      TrackerHit->AddSil(edep);
      TrackerHit->SetPos(aStep->GetPreStepPoint()->GetPosition());
      TrackerHit->SetNSilPlane(PlaneNumber);
      TrackerHit->SetNStrip(StripNumber);
      ThitYID[NChannel] = 
        TrackerCollection->insert(TrackerHit)-1;
    }
      else // This is not new
    {
      (*TrackerCollection)[ThitYID[NChannel]]->AddSil(edep);
          // G4cout << "Y" << PlaneNumber << " " << StripNumber << G4endl;
    }
    }
  
  return true;
}

---

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

• GammaRayTelTrackerSD.hh
• GammaRayTelTrackerSD.cc