TrackerSD Class Reference

#include <TrackerSD.hh>

Inheritance diagram for TrackerSD:

Collaboration diagram for TrackerSD:

Public Member Functions
	TrackerSD (const G4String &name, const G4String &hitsCollectionName)
virtual	~TrackerSD ()
virtual void	Initialize (G4HCofThisEvent *hitCollection)
virtual G4bool	ProcessHits (G4Step *step, G4TouchableHistory *history)
virtual void	EndOfEvent (G4HCofThisEvent *hitCollection)
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
virtual void	clear ()
virtual void	DrawAll ()
virtual void	PrintAll ()
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

Private Attributes
TrackerHitsCollection *	fHitsCollection

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

Tracker sensitive detector class

The hits are accounted in hits in ProcessHits() function which is called by Geant4 kernel at each step. A hit is created with each step with non zero energy deposit.

Definition at line 52 of file TrackerSD.hh.

Constructor & Destructor Documentation

TrackerSD()

TrackerSD::TrackerSD	(	const G4String &	name,
		const G4String &	hitsCollectionName
	)

Definition at line 47 of file TrackerSD.cc.

 : G4VSensitiveDetector(name),
   fHitsCollection(NULL)
{
  collectionName.insert(hitsCollectionName);
}

Here is the call graph for this function:

~TrackerSD()

TrackerSD::~TrackerSD ( )

virtual

Definition at line 57 of file TrackerSD.cc.

{}

Member Function Documentation

EndOfEvent()

void TrackerSD::EndOfEvent ( G4HCofThisEvent *  hitCollection )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 101 of file TrackerSD.cc.

{
  G4int nofHits = fHitsCollection->entries();
     
  if ( verboseLevel>1 ) 

            G4cout << G4endl
            << "-------->Hits Collection: in this event they are " 
            << nofHits 
            << " hits in the target volume " << G4endl;
  
  // PROCESSING OF MICRODOSIMETRY Y & Z SPECTRA
  
  // main loop

  // *************************************
  // Please select herebelow :
  // 1 - the number of sampling per incident envent
  //     variable name = nbSampling
  //     (it is set to 100 by default)
  //
  //
  // 2 - the radius of the target sphere
  //     variable name = radius
  //     (it is set to 50 nm by default)
  //
  
  G4int nbSampling = 1000;
  G4double radius = 1.50*um;

  // ************************************
    
  for ( G4int k=0; k<nbSampling; k++ ) 
  {
   
  //***************
  // FIRST SAMPLING: Y
  //***************
  
  // select random hit
  G4int randHit=0; // Runs from 0 to number of hits - 1
  randHit = static_cast<G4int>( G4UniformRand()*nofHits );
  
  //G4cout << static_cast<G4int>(2.7) 
  //<< "======> random selection of hit number randHit =" 
  //       << randHit << G4endl;
  
  // get random hit position
  G4ThreeVector hitPos =  (*fHitsCollection)[randHit]->GetPos();
  //G4cout << "======> random hit position x/nm =" << hitPos.x()/nm << G4endl; 
  //G4cout << "======> random hit position y/nm =" << hitPos.y()/nm << G4endl; 
  //G4cout << "======> random hit position z/nm =" << hitPos.z()/nm << G4endl; 
  
  // select random center of sphere within radius
  G4double chord = 4.*radius/3;
  //G4double density = 1 * g/cm3;
  //G4double mass = (4./3)*CLHEP::pi*radius*radius*radius*density;
  
  G4ThreeVector randDir = G4RandomDirection();
  G4double randRadius = G4UniformRand()*radius;
  G4ThreeVector randCenterPos = randRadius*randDir + hitPos;
  
  // search for neighbouring hits in the sphere and cumulate deposited energy 
  //  in epsilon
  G4double epsilon = 0;
  G4int nbEdep = 0;
  
  for ( G4int i=0; i<nofHits; i++ ) 
  { 
    G4ThreeVector localPos = (*fHitsCollection)[i]->GetPos();; 
    if ( 
        (localPos.x()-randCenterPos.x()) * (localPos.x()-randCenterPos.x()) +
        (localPos.y()-randCenterPos.y()) * (localPos.y()-randCenterPos.y()) +
        (localPos.z()-randCenterPos.z()) * (localPos.z()-randCenterPos.z()) 
         <= radius*radius
       ) 
       
    { 
      epsilon = epsilon + (*fHitsCollection)[i]->GetEdep() ;
      nbEdep = nbEdep+1;
    }
       
  }
  
        
  //***************

  /*
  // for testing only

  G4cout << "======> for hit number #" << randHit <<
  ", we collect " 
  << nbEdep << " energy depositions in a sphere of radius " 
  << radius/nm << " nm and mass " 
  << mass/kg << " kg for a total of " 
  << epsilon/eV << " eV or " 
  << (epsilon/joule)/(mass/kg) << " Gy" << G4endl;
  G4cout << "-" << G4endl;
  
  FILE* myFile;
  myFile=fopen("yz.txt","a");
  fprintf(myFile,"%e %e %e\n",radius/nm,(epsilon/eV)/(chord/nm),
   (epsilon/joule)/(mass/kg));
  fclose(myFile);
  */
  
  // get analysis manager
  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();

  // fill ntuple
  analysisManager->FillNtupleDColumn(0, radius/nm);
  analysisManager->FillNtupleDColumn(1, (epsilon/eV)/(chord/nm));
  analysisManager->AddNtupleRow();
  
  } // END MAIN LOOP ON SAMPLING

}

Here is the call graph for this function:

Initialize()

void TrackerSD::Initialize ( G4HCofThisEvent *  hitCollection )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 62 of file TrackerSD.cc.

{
  // Create hits collection

  fHitsCollection 
    = new TrackerHitsCollection(SensitiveDetectorName, collectionName[0]); 

  // Add this collection in hce

  G4int hcID 
    = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
  hce->AddHitsCollection( hcID, fHitsCollection ); 
}

Here is the call graph for this function:

ProcessHits()

G4bool TrackerSD::ProcessHits ( G4Step *  step, G4TouchableHistory *  history  )

virtual

Implements G4VSensitiveDetector.

Definition at line 78 of file TrackerSD.cc.

{  
  // energy deposit
  G4double edep = aStep->GetTotalEnergyDeposit();

  if (edep==0.) return false;

  TrackerHit* newHit = new TrackerHit();

  newHit->SetTrackID  (aStep->GetTrack()->GetTrackID());
  newHit->SetEdep(edep);
  newHit->SetPos (aStep->GetPostStepPoint()->GetPosition());

  fHitsCollection->insert( newHit );

  //newHit->Print();

  return true;
}

Here is the call graph for this function:

Member Data Documentation

fHitsCollection

TrackerHitsCollection* TrackerSD::fHitsCollection

private

Definition at line 65 of file TrackerSD.hh.

---

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

• TrackerSD.hh
• TrackerSD.cc