DummyChargeFlipProcess Class Reference

#include <DummyChargeFlipProcess.h>

Inheritance diagram for DummyChargeFlipProcess:

Public Member Functions
void	BuildPhysicsTable (const G4ParticleDefinition &aParticleType)
	DummyChargeFlipProcess (const G4String &processName="Dummy")
void	DumpPhysicsTable (const G4ParticleDefinition &aParticleType)
G4bool	IsApplicable (const G4ParticleDefinition &aParticleType)
G4VParticleChange *	PostStepDoIt (const G4Track &aTrack, const G4Step &aStep)
	~DummyChargeFlipProcess ()

Private Member Functions
G4double	GetMicroscopicCrossSection (const G4DynamicParticle *aParticle, const G4Element *anElement, G4double aTemp)

Detailed Description

Definition at line 18 of file DummyChargeFlipProcess.h.

Constructor & Destructor Documentation

DummyChargeFlipProcess::DummyChargeFlipProcess

(

const G4String &

processName = "Dummy"

)

Definition at line 8 of file DummyChargeFlipProcess.cc.

                                                    : 
      G4HadronicProcess(processName)
{
  AddDataSet(new G4HadronElasticDataSet);
}

DummyChargeFlipProcess::~DummyChargeFlipProcess

(

)

Definition at line 14 of file DummyChargeFlipProcess.cc.

{
}

Member Function Documentation

void DummyChargeFlipProcess::BuildPhysicsTable

(

const G4ParticleDefinition &

aParticleType

)

Definition at line 20 of file DummyChargeFlipProcess.cc.

{
   if (!G4HadronicProcess::GetCrossSectionDataStore()) {
      return;
   }
   G4HadronicProcess::GetCrossSectionDataStore()->BuildPhysicsTable(aParticleType);
}

void DummyChargeFlipProcess::DumpPhysicsTable

(

const G4ParticleDefinition &

aParticleType

)

Definition at line 51 of file DummyChargeFlipProcess.cc.

{

}

G4double DummyChargeFlipProcess::GetMicroscopicCrossSection ( const G4DynamicParticle *  aParticle, const G4Element *  anElement, G4double  aTemp  )

private

Definition at line 30 of file DummyChargeFlipProcess.cc.

{

   return 30*millibarn*element->GetN(); 
}

G4bool DummyChargeFlipProcess::IsApplicable

(

const G4ParticleDefinition &

aParticleType

)

Definition at line 40 of file DummyChargeFlipProcess.cc.

{
   if(aParticleType.GetParticleType()  == "rhadron")
    return true;
  else
    return false;
}

G4VParticleChange * DummyChargeFlipProcess::PostStepDoIt	(	const G4Track &	aTrack,
		const G4Step &	aStep
	)

Definition at line 57 of file DummyChargeFlipProcess.cc.

References gather_cfg::cout.

{
  G4ParticleChange * pc = new G4ParticleChange();
  pc->Initialize(aTrack);
  const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
  G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();

  G4double   ParentEnergy  = aParticle->GetTotalEnergy();
  G4ThreeVector ParentDirection(aParticle->GetMomentumDirection());

  G4double energyDeposit = 0.0;
  G4double finalGlobalTime = aTrack.GetGlobalTime();

  G4int numberOfSecondaries = 1;
  pc->SetNumberOfSecondaries(numberOfSecondaries);
  const G4TouchableHandle thand = aTrack.GetTouchableHandle();

     // get current position of the track
     aTrack.GetPosition();
     // create a new track object
      G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
      float randomParticle = G4UniformRand();
      G4ParticleDefinition * newType = aParticleDef;
      if(randomParticle < 0.333)
        newType=particleTable->FindParticle(1009213);
      else if(randomParticle > 0.667)
        newType=particleTable->FindParticle(-1009213);
      else
        newType=particleTable->FindParticle(1009113);
     
      std::cout << "RHADRON: New charge = " << newType->GetPDGCharge() << std::endl;
       
     G4DynamicParticle * newP =  new G4DynamicParticle(newType,ParentDirection,ParentEnergy);
     G4Track* secondary = new G4Track( newP ,
                      finalGlobalTime ,
                       aTrack.GetPosition()
     );
     // switch on good for tracking flag
     secondary->SetGoodForTrackingFlag();
     secondary->SetTouchableHandle(thand);
     // add the secondary track in the List
     pc->AddSecondary(secondary);

  // Kill the parent particle
  pc->ProposeTrackStatus( fStopAndKill ) ;
  pc->ProposeLocalEnergyDeposit(energyDeposit); 
  pc->ProposeGlobalTime( finalGlobalTime );
  // Clear NumberOfInteractionLengthLeft
  ClearNumberOfInteractionLengthLeft();



   return pc;

}