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 10 of file DummyChargeFlipProcess.cc.

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

DummyChargeFlipProcess::~DummyChargeFlipProcess

(

)

Definition at line 16 of file DummyChargeFlipProcess.cc.

References BuildPhysicsTable().

{
}

Member Function Documentation

void DummyChargeFlipProcess::BuildPhysicsTable

(

const G4ParticleDefinition &

aParticleType

)

Definition at line 22 of file DummyChargeFlipProcess.cc.

References GetMicroscopicCrossSection().

Referenced by ~DummyChargeFlipProcess().

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

void DummyChargeFlipProcess::DumpPhysicsTable

(

const G4ParticleDefinition &

aParticleType

)

Definition at line 53 of file DummyChargeFlipProcess.cc.

Referenced by IsApplicable().

{

}

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

private

Definition at line 32 of file DummyChargeFlipProcess.cc.

References IsApplicable().

Referenced by BuildPhysicsTable().

{

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

G4bool DummyChargeFlipProcess::IsApplicable

(

const G4ParticleDefinition &

aParticleType

)

Definition at line 42 of file DummyChargeFlipProcess.cc.

References DumpPhysicsTable().

Referenced by GetMicroscopicCrossSection().

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

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

Definition at line 59 of file DummyChargeFlipProcess.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::G4cout.

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

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

  return pc;
}