DummyChargeFlipProcess Class Reference

#include <DummyChargeFlipProcess.h>

Inheritance diagram for DummyChargeFlipProcess:

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

Private Attributes
G4ParticleChange *	fPartChange

Detailed Description

Definition at line 8 of file DummyChargeFlipProcess.h.

Constructor & Destructor Documentation

DummyChargeFlipProcess()

DummyChargeFlipProcess::DummyChargeFlipProcess

(

const G4String &

processName = "Dummy"

)

Definition at line 13 of file DummyChargeFlipProcess.cc.

References fPartChange.

                                                                    : G4HadronicProcess(pname, fHadronic) {
  AddDataSet(new G4NeutronElasticXS());
  fPartChange = new G4ParticleChange();
}

~DummyChargeFlipProcess()

DummyChargeFlipProcess::~DummyChargeFlipProcess ( )

override

Definition at line 18 of file DummyChargeFlipProcess.cc.

References fPartChange.

{ delete fPartChange; }

Member Function Documentation

IsApplicable()

G4bool DummyChargeFlipProcess::IsApplicable ( const G4ParticleDefinition &  aParticleType )

override

Definition at line 20 of file DummyChargeFlipProcess.cc.

                                                                                     {
  return (aParticleType.GetParticleType() == "rhadron");
}

PostStepDoIt()

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

override

Definition at line 24 of file DummyChargeFlipProcess.cc.

References fPartChange.

                                                                                            {
  fPartChange->Initialize(aTrack);
  const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();

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

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

  G4int numberOfSecondaries = 1;
  fPartChange->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;
  if (randomParticle < 0.333)
    newType = particleTable->FindParticle(1009213);
  else if (randomParticle > 0.667)
    newType = particleTable->FindParticle(-1009213);
  else
    newType = particleTable->FindParticle(1009113);

  //G4cout << "RHADRON: New charge = " << newType->GetPDGCharge() << G4endl;

  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
  fPartChange->AddSecondary(secondary);

  // Kill the parent particle
  fPartChange->ProposeTrackStatus(fStopAndKill);
  fPartChange->ProposeLocalEnergyDeposit(energyDeposit);
  fPartChange->ProposeGlobalTime(finalGlobalTime);

  ClearNumberOfInteractionLengthLeft();

  return fPartChange;
}

Member Data Documentation

fPartChange

G4ParticleChange* DummyChargeFlipProcess::fPartChange

private

Definition at line 19 of file DummyChargeFlipProcess.h.

Referenced by DummyChargeFlipProcess(), PostStepDoIt(), and ~DummyChargeFlipProcess().