HadrontherapyIonStandard.cc

Go to the documentation of this file.

//
// ********************************************************************
// * DISCLAIMER                                                       *
// *                                                                  *
// * The following disclaimer summarizes all the specific disclaimers *
// * of contributors to this software. The specific disclaimers,which *
// * govern, are listed with their locations in:                      *
// *   http://cern.ch/geant4/license                                  *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.                                                             *
// *                                                                  *
// * This  code  implementation is the  intellectual property  of the *
// * GEANT4 collaboration.                                            *
// * By copying,  distributing  or modifying the Program (or any work *
// * based  on  the Program)  you indicate  your  acceptance of  this *
// * statement, and all its terms.                                    *
// ********************************************************************
// $Id: HadrontherapyIonStandard.cc; May 2005
// ----------------------------------------------------------------------------
//                 GEANT 4 - Hadrontherapy example
// ----------------------------------------------------------------------------
// Code developed by:
//
// G.A.P. Cirrone(a)*, F. Di Rosa(a), S. Guatelli(b), G. Russo(a)
//
// (a) Laboratori Nazionali del Sud
//     of the National Institute for Nuclear Physics, Catania, Italy
// (b) National Institute for Nuclear Physics Section of Genova, genova, Italy
//
// * cirrone@lns.infn.it
// ----------------------------------------------------------------------------

#include "HadrontherapyIonStandard.hh"
#include "G4ProcessManager.hh"
#include "G4ParticleDefinition.hh"
#include "G4ParticleTable.hh"
#include "G4hMultipleScattering.hh"
#include "G4hIonisation.hh"
#include "G4ionIonisation.hh"
#include "G4hMultipleScattering.hh"
#include "G4StepLimiter.hh"

HadrontherapyIonStandard::HadrontherapyIonStandard(const G4String& name): G4VPhysicsConstructor(name)
{ }

HadrontherapyIonStandard::~HadrontherapyIonStandard()
{ }

void HadrontherapyIonStandard::ConstructProcess()
{
  G4ParticleTable::G4PTblDicIterator* theParticleIterator = G4ParticleTable::GetParticleTable()->GetIterator();
  theParticleIterator -> reset();

  while( (*theParticleIterator)() )
    {
      G4ParticleDefinition* particle = theParticleIterator -> value();
      G4ProcessManager* manager = particle -> GetProcessManager();
      G4String particleName = particle -> GetParticleName();
      G4double charge = particle -> GetPDGCharge();

      if (particleName == "GenericIon")
   {
     G4ionIonisation* ionisation = new G4ionIonisation();
     G4VProcess*  multipleScattering = new G4hMultipleScattering();
     manager -> AddProcess(multipleScattering, -1,1,1);
     manager -> AddProcess(ionisation, -1,2,2);
     manager -> AddProcess(new G4StepLimiter(),-1,-1, 3);

   }
      //protons and generic hadrons
      if (( charge != 0. ) && particleName != "e+" && particleName != "mu+" &&
     particleName != "e-" && particleName != "mu-" && particleName !="GenericIon")
   {
     if((!particle -> IsShortLived()) &&
        (particle -> GetParticleName() != "chargedgeantino"))
       {
         G4hIonisation* ionisation = new G4hIonisation();
         G4VProcess*  multipleScattering = new G4hMultipleScattering();
         manager -> AddProcess(multipleScattering, -1,1,1);
         manager -> AddProcess(ionisation, -1,2,2);
         manager -> AddProcess(new G4StepLimiter(),-1,-1, 3);
       }
   }
    }
}