HadrontherapyDetectorSD Class Reference

#include <HadrontherapyDetectorSD.hh>

Inheritance diagram for HadrontherapyDetectorSD:

Public Member Functions
	HadrontherapyDetectorSD (G4String name)
	~HadrontherapyDetectorSD ()
void	Initialize (G4HCofThisEvent *)
G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROhist)
void	EndOfEvent (G4HCofThisEvent *HCE)
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
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

Data Fields
std::ofstream	ofs

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 45 of file HadrontherapyDetectorSD.hh.

Constructor & Destructor Documentation

HadrontherapyDetectorSD::HadrontherapyDetectorSD

(

G4String

name

)

Definition at line 49 of file HadrontherapyDetectorSD.cc.

References G4VSensitiveDetector::collectionName, and G4CollectionNameVector::insert().

                                                             :
    G4VSensitiveDetector(name)
{ 
    G4String HCname;
    collectionName.insert(HCname="HadrontherapyDetectorHitsCollection");
    HitsCollection = NULL; 
    sensitiveDetectorName = name;

}

HadrontherapyDetectorSD::~HadrontherapyDetectorSD

(

)

Definition at line 60 of file HadrontherapyDetectorSD.cc.

{ 
}

Member Function Documentation

void HadrontherapyDetectorSD::EndOfEvent ( G4HCofThisEvent *  HCE )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 232 of file HadrontherapyDetectorSD.cc.

References G4VSensitiveDetector::GetCollectionID().

{
 
    static G4int HCID = -1;
    if(HCID < 0)
    { 
    HCID = GetCollectionID(0); 
    }

    HCE -> AddHitsCollection(HCID,HitsCollection);
}

void HadrontherapyDetectorSD::Initialize ( G4HCofThisEvent *  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 65 of file HadrontherapyDetectorSD.cc.

References G4VSensitiveDetector::collectionName.

{
  
    HitsCollection = new HadrontherapyDetectorHitsCollection(sensitiveDetectorName,
        collectionName[0]);
}

G4bool HadrontherapyDetectorSD::ProcessHits ( G4Step *  aStep, G4TouchableHistory *  ROhist  )

virtual

Implements G4VSensitiveDetector.

Definition at line 73 of file HadrontherapyDetectorSD.cc.

References HadrontherapyLet::GetInstance(), HadrontherapyAnalysisManager::GetInstance(), HadrontherapyMatrix::GetInstance(), G4VSensitiveDetector::GetName(), G4ParticleDefinition::GetPDGEncoding(), G4Step::GetPreStepPoint(), G4VTouchable::GetReplicaNumber(), G4StepPoint::GetTouchable(), python.hepunit::keV, and python.hepunit::MeV.

{

 
  if (aStep -> GetPreStepPoint() -> GetPhysicalVolume() -> GetName() != "RODetectorZDivisionPhys") return false;
    

    // Get kinetic energy
    G4Track * theTrack = aStep  ->  GetTrack();
    G4double kineticEnergy =  theTrack -> GetKineticEnergy();  

    G4ParticleDefinition *particleDef = theTrack -> GetDefinition();
    //Get particle name  
    G4String particleName =  particleDef -> GetParticleName();  

    G4int pdg = particleDef ->GetPDGEncoding();

    // Get unique track_id (in an event)
    G4int trackID = theTrack -> GetTrackID();

    G4double energyDeposit = aStep -> GetTotalEnergyDeposit();

    G4double DX = aStep -> GetStepLength();
    G4int Z = particleDef-> GetAtomicNumber();
    G4int A = particleDef-> GetAtomicMass();

  
 // Read voxel indexes: i is the x index, k is the z index
 const G4VTouchable* touchable = aStep->GetPreStepPoint()->GetTouchable();
 G4int k  = touchable->GetReplicaNumber(0);
 G4int i  = touchable->GetReplicaNumber(2);
 G4int j  = touchable->GetReplicaNumber(1);

    

#ifdef G4ANALYSIS_USE_ROOT
    HadrontherapyAnalysisManager* analysis = HadrontherapyAnalysisManager::GetInstance();
#endif

    HadrontherapyMatrix* matrix = HadrontherapyMatrix::GetInstance();
    HadrontherapyLet* let = HadrontherapyLet::GetInstance();
    
    
  //  ******************** let ***************************
     if (let)
    {
     if ( !(Z==0 && A==1) ) // All but not neutrons 
    {
      if( energyDeposit>0. && DX >0. )
        {
          if (pdg !=22) // not gamma
        {
          let -> FillEnergySpectrum(trackID, particleDef,energyDeposit, DX, i, j, k);
        }
          else if (kineticEnergy > 50.*keV) // gamma cut
        {
            let -> FillEnergySpectrum(trackID, particleDef,energyDeposit, DX, i , j, k);
        }

        }
    }
    }
    
    
    
    //  ******************** let ***************************
  
    
    if (matrix)
    {

    // Increment Fluences & accumulate energy spectra
    // Hit voxels are marked with track_id throught hitTrack matrix
    G4int* hitTrack = matrix -> GetHitTrack(i,j,k); // hitTrack MUST BE cleared at every eventAction!
    if ( *hitTrack != trackID )
    {
                 *hitTrack = trackID;
        /*
         * Fill FLUENCE data for every single nuclide 
         * Exclude e-, neutrons, gamma, ...
         */
        if ( Z >= 1) matrix -> Fill(trackID, particleDef, i, j, k, 0, true);
            

#ifdef G4ANALYSIS_USE_ROOT
/*
        // Fragments kinetic energy (ntuple)
        if (trackID !=1 && Z>=1) 
        {
        // First step kinetic energy for every fragment 
         analysis -> FillKineticFragmentTuple(i, j, k, A, Z, kineticEnergy/MeV);
        }    
        // Kinetic energy spectra for primary particles 

        if ( trackID == 1 && i == 0) 
        {
        // First step kinetic energy for primaries only
        analysis -> FillKineticEnergyPrimaryNTuple(i, j, k, kineticEnergy/MeV);
        }
*/
#endif
    }    

    if(energyDeposit != 0)
    {
/*
 *  This method will fill a dose matrix for every single nuclide. 
 *  A method of the HadrontherapyMatrix class (StoreDoseFluenceAscii())
 *  is called automatically at the end of main (or via the macro command /analysis/writeDoseFile.
 *  It permits to store all dose/fluence data into a single plane ASCII file. 
*/      
        // if (A==1 && Z==1) // primary and sec. protons 
        if ( Z>=1 )    //  exclude e-, neutrons, gamma, ...
            matrix -> Fill(trackID, particleDef, i, j, k, energyDeposit);
        /*
         * Create a hit with the information of position is in the detector     
         */
        HadrontherapyDetectorHit* detectorHit = new HadrontherapyDetectorHit();       
        detectorHit -> SetEdepAndPosition(i, j, k, energyDeposit); 
        HitsCollection -> insert(detectorHit);
    }
    }

#ifdef G4ANALYSIS_USE_ROOT
    if(energyDeposit != 0)
    {  
    if(trackID != 1)
    {
        if (particleName == "proton")
        analysis -> SecondaryProtonEnergyDeposit(i, energyDeposit/MeV);

        else if (particleName == "neutron")
        analysis -> SecondaryNeutronEnergyDeposit(i, energyDeposit/MeV);

        else if (particleName == "alpha")
        analysis -> SecondaryAlphaEnergyDeposit(i, energyDeposit/MeV);

        else if (particleName == "gamma")
        analysis -> SecondaryGammaEnergyDeposit(i, energyDeposit/MeV);

        else if (particleName == "e-")
        analysis -> SecondaryElectronEnergyDeposit(i, energyDeposit/MeV);

        else if (particleName == "triton")
        analysis -> SecondaryTritonEnergyDeposit(i, energyDeposit/MeV);

        else if (particleName == "deuteron")
        analysis -> SecondaryDeuteronEnergyDeposit(i, energyDeposit/MeV);

        else if (particleName == "pi+" || particleName == "pi-" ||  particleName == "pi0")
        analysis -> SecondaryPionEnergyDeposit(i, energyDeposit/MeV);       
    }
    }
#endif

    return true;
}

Field Documentation

std::ofstream HadrontherapyDetectorSD::ofs

Definition at line 51 of file HadrontherapyDetectorSD.hh.

---

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

• HadrontherapyDetectorSD.hh
• HadrontherapyDetectorSD.cc