Par01PiModel Class Reference

#include <Par01PiModel.hh>

Inheritance diagram for Par01PiModel:

Public Member Functions
	Par01PiModel (G4Region *anEnvelope)
	~Par01PiModel ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
virtual G4bool	ModelTrigger (const G4FastTrack &)
virtual void	DoIt (const G4FastTrack &, G4FastStep &)
Public Member Functions inherited from G4VFastSimulationModel
	G4VFastSimulationModel (const G4String &aName)
	G4VFastSimulationModel (const G4String &aName, G4Envelope *, G4bool IsUnique=FALSE)
virtual	~G4VFastSimulationModel ()
virtual G4bool	AtRestModelTrigger (const G4FastTrack &)
virtual void	AtRestDoIt (const G4FastTrack &, G4FastStep &)
const G4String	GetName () const
G4bool	operator== (const G4VFastSimulationModel &) const

Detailed Description

Definition at line 38 of file Par01PiModel.hh.

Constructor & Destructor Documentation

Par01PiModel::Par01PiModel

(

G4Region *

anEnvelope

)

Definition at line 35 of file Par01PiModel.cc.

                                               :
  G4VFastSimulationModel("Par01PiModel",anEnvelope)
{;}

Par01PiModel::~Par01PiModel

(

)

Definition at line 39 of file Par01PiModel.cc.

{;}

Member Function Documentation

void Par01PiModel::DoIt ( const G4FastTrack &  fastTrack, G4FastStep &  fastStep  )

virtual

Implements G4VFastSimulationModel.

Definition at line 82 of file Par01PiModel.cc.

References G4FastStep::CreateSecondaryTrack(), G4cout, G4endl, G4Gamma::GammaDefinition(), G4FastTrack::GetEnvelopeSolid(), G4Track::GetGlobalTime(), G4FastTrack::GetPrimaryTrack(), G4FastTrack::GetPrimaryTrackLocalDirection(), G4FastTrack::GetPrimaryTrackLocalPosition(), position, G4FastStep::ProposePrimaryTrackFinalPosition(), G4FastStep::SetNumberOfSecondaryTracks(), and CLHEP::Hep3Vector::setZ().

{

  //------------------------------------------------
  // The primary track continues along its direction.
  // One secondary (a photon) is added:
  //------------------------------------------------
  G4cout << "      Pion `model' applied\n" << G4endl;

  //------------------------------
  // Primary:
  //    idem as in "DefaultModel":
  //
  //------------------------------
  G4ThreeVector position;
  G4double distance;
  distance = fastTrack.GetEnvelopeSolid()->
    DistanceToOut(fastTrack.GetPrimaryTrackLocalPosition(),
                  fastTrack.GetPrimaryTrackLocalDirection());
  position = fastTrack.GetPrimaryTrackLocalPosition() + 
    distance*fastTrack.GetPrimaryTrackLocalDirection();
  
  // -- set final position:
  fastStep.ProposePrimaryTrackFinalPosition(position);
  
  //---------------------------
  // Secondary:
  //   Adds one "secondary":
  //---------------------------
  // -- First, user has to say how many secondaries will be created:
  fastStep.SetNumberOfSecondaryTracks(1);
  
  //------------------------
  // -- Build the secondary:
  //------------------------
  // -- direction:
  G4ParticleMomentum direction(fastTrack.GetPrimaryTrackLocalDirection());
  direction.setZ(direction.z()*0.5);
  direction.setY(direction.y()+direction.z()*0.1);
  direction = direction.unit(); // necessary ?
  
  // -- dynamics (Note that many constructors exists for G4DynamicParticle
  // -- see prototype/particle+matter/particles/management/include/G4DynamicParticle.hh)
  G4DynamicParticle dynamique(G4Gamma::GammaDefinition(),
                              direction,
                              fastTrack.GetPrimaryTrack()->
                              GetKineticEnergy()/2.);
  // -- position:
  G4double Dist;
  Dist = fastTrack.GetEnvelopeSolid()->
    DistanceToOut(fastTrack.GetPrimaryTrackLocalPosition(),
                  direction);
  G4ThreeVector posi;
  posi = fastTrack.GetPrimaryTrackLocalPosition() + Dist*direction;
  
  //------------------------------------
  //-- Creation of the secondary Track:
  //------------------------------------
  fastStep.CreateSecondaryTrack(dynamique, posi, 
                                fastTrack.GetPrimaryTrack()->GetGlobalTime());
  
}

G4bool Par01PiModel::IsApplicable ( const G4ParticleDefinition &  particleType )

virtual

Implements G4VFastSimulationModel.

Definition at line 42 of file Par01PiModel.cc.

References G4PionMinus::PionMinusDefinition(), and G4PionPlus::PionPlusDefinition().

{
  return
    &particleType == G4PionMinus::PionMinusDefinition() ||
    &particleType == G4PionPlus::PionPlusDefinition();
}

G4bool Par01PiModel::ModelTrigger ( const G4FastTrack &  fastTrack )

virtual

Implements G4VFastSimulationModel.

Definition at line 49 of file Par01PiModel.cc.

References G4cout, G4endl, G4Track::GetMomentum(), G4Track::GetPosition(), G4FastTrack::GetPrimaryTrack(), G4FastTrack::GetPrimaryTrackLocalDirection(), G4FastTrack::GetPrimaryTrackLocalPosition(), and CLHEP::Hep3Vector::unit().

                                                              {
  //-------------------------------------------------------------
  // UserTrigger() method: method which has to decide if
  // the parameterisation has to be applied.
  // Here ModelTrigger() asks the user (ie you) a 0/1 answer.
  //
  // Note that quantities like the local/global position/direction etc..
  // are available at this level via the fastTrack parameter (allowing 
  // to check distance from boundaries, see below  to allow the decision)
  //--------------------------------------------------------------
  G4cout << "\nPar01PiModel::ModelTrigger() called:" << G4endl;
  G4cout <<   "--------------------------------" << G4endl;
  G4cout << "(particle is a " << fastTrack.GetPrimaryTrack()->
    GetDefinition()->GetParticleName() << " )\n" << G4endl;

  // -- Examples of available informations:

  // -- position:
  G4cout << "         Track position: " <<
    fastTrack.GetPrimaryTrack()->GetPosition()  << "(global coord.)"   <<
    fastTrack.GetPrimaryTrackLocalPosition()  << "(in envelope coord.)" 
       << G4endl;

  // -- direction:
  G4cout << "         Track direction:"       <<
    fastTrack.GetPrimaryTrack()->GetMomentum().unit() << 
    "(global coord.)"      <<
    fastTrack.GetPrimaryTrackLocalDirection() << "(in envelope coord.)" << 
    G4endl;

  return true;
}

---

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

• Par01PiModel.hh
• Par01PiModel.cc