G4DNAMolecularDissociation Class Reference

#include <G4DNAMolecularDissociation.hh>

Inheritance diagram for G4DNAMolecularDissociation:

Collaboration diagram for G4DNAMolecularDissociation:

Public Member Functions
	G4DNAMolecularDissociation (const G4String &processName="DNAMolecularDecay", G4ProcessType type=fDecay)
virtual	~G4DNAMolecularDissociation ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)
G4VParticleChange *	AtRestDoIt (const G4Track &track, const G4Step &step)
G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &step)
void	SetVerbose (G4int)
void	SetDecayDisplacer (const G4ParticleDefinition *, G4VMolecularDecayDisplacer *)
G4VMolecularDecayDisplacer *	GetDecayDisplacer (const G4ParticleDefinition *)

Protected Member Functions
virtual G4VParticleChange *	DecayIt (const G4Track &, const G4Step &)
virtual G4double	GetMeanLifeTime (const G4Track &, G4ForceCondition *)
virtual G4double	GetMeanFreePath (const G4Track &, G4double, G4ForceCondition *)

Private Types
typedef std::map< const G4ParticleDefinition *, G4VMolecularDecayDisplacer * >	DecayDisplacementMap

Private Member Functions
	G4DNAMolecularDissociation ()
	G4DNAMolecularDissociation (const G4DNAMolecularDissociation &right)
G4DNAMolecularDissociation &	operator= (const G4DNAMolecularDissociation &right)

Private Attributes
G4bool	fDecayAtFixedTime
DecayDisplacementMap	fDecayDisplacementMap
G4int	fVerbose

Detailed Description

G4DNAMolecularDissociation should be called only for molecules. It will dissociate the molecules using the decay associated to this molecule and if a displacement scheme has been registered, it will place the products to the expected position.

Definition at line 66 of file G4DNAMolecularDissociation.hh.

Member Typedef Documentation

DecayDisplacementMap

typedef std::map<const G4ParticleDefinition*, G4VMolecularDecayDisplacer*> G4DNAMolecularDissociation::DecayDisplacementMap

private

Definition at line 127 of file G4DNAMolecularDissociation.hh.

Constructor & Destructor Documentation

G4DNAMolecularDissociation() [1/3]

G4DNAMolecularDissociation::G4DNAMolecularDissociation	(	const G4String &	processName = "DNAMolecularDecay",
		G4ProcessType	type = fDecay
	)

Definition at line 52 of file G4DNAMolecularDissociation.cc.

                                               :
    G4VITRestDiscreteProcess(processName, type)
{
  // set Process Sub Type
  SetProcessSubType(59); // DNA sub-type
  enableAlongStepDoIt = false;
  enablePostStepDoIt = true;
  enableAtRestDoIt = true;

  fVerbose = 0;

#ifdef G4VERBOSE
  if(verboseLevel > 1)
  {
    G4cout << "G4MolecularDissociationProcess constructor " << "  Name:"
           << processName << G4endl;
  }
#endif

  pParticleChange = &aParticleChange;

  fDecayAtFixedTime = true;
  fProposesTimeStep = true;
}

~G4DNAMolecularDissociation()

G4DNAMolecularDissociation::~G4DNAMolecularDissociation ( )

virtual

Definition at line 80 of file G4DNAMolecularDissociation.cc.

{
  DecayDisplacementMap::iterator it = fDecayDisplacementMap.begin();

  for(; it != fDecayDisplacementMap.end(); it++)
  {
    if(it->second)
    {
      delete it->second;
      it->second = 0;
    }
  }
  fDecayDisplacementMap.clear();
}

Here is the call graph for this function:

G4DNAMolecularDissociation() [2/3]

G4DNAMolecularDissociation::G4DNAMolecularDissociation ( )

private

Here is the caller graph for this function:

G4DNAMolecularDissociation() [3/3]

G4DNAMolecularDissociation::G4DNAMolecularDissociation ( const G4DNAMolecularDissociation &  right )

private

Definition at line 98 of file G4DNAMolecularDissociation.cc.

                                                                    :
    G4VITRestDiscreteProcess(right)
{
  fDecayAtFixedTime = right.fDecayAtFixedTime;
  fDecayDisplacementMap = right.fDecayDisplacementMap;
  fVerbose = right.fVerbose;
}

Here is the call graph for this function:

Member Function Documentation

AtRestDoIt()

G4VParticleChange * G4DNAMolecularDissociation::AtRestDoIt ( const G4Track &  track, const G4Step &  step  )

inline

Definition at line 151 of file G4DNAMolecularDissociation.hh.

{
    ClearNumberOfInteractionLengthLeft();
    ClearInteractionTimeLeft();
    return DecayIt(track, step);
}

Here is the call graph for this function:

Here is the caller graph for this function:

AtRestGetPhysicalInteractionLength()

G4double G4DNAMolecularDissociation::AtRestGetPhysicalInteractionLength ( const G4Track &  track, G4ForceCondition *  condition  )

inline

Definition at line 138 of file G4DNAMolecularDissociation.hh.

{
    if(fDecayAtFixedTime)
    {
        return GetMeanLifeTime(track, condition);
    }

    return G4VITRestDiscreteProcess::AtRestGetPhysicalInteractionLength(track, condition);
//    return G4VITRestProcess::AtRestGetPhysicalInteractionLength(track, condition);
}

Here is the call graph for this function:

DecayIt()

G4VParticleChange * G4DNAMolecularDissociation::DecayIt ( const G4Track &  track, const G4Step &    )

protectedvirtual

Definition at line 141 of file G4DNAMolecularDissociation.cc.

{
  // DEBUG
  // G4cout << "Is calling G4MolecularDecayProcess::DecayIt" << G4endl;

  aParticleChange.Initialize(track);
  const G4Molecule * theMotherMolecule = GetMolecule(track);
  const G4MoleculeDefinition* moleculeDefinition = theMotherMolecule
      ->GetDefinition();

  // DEBUG
  // G4cout <<"Calling G4MolecularDecayProcess::DecayIt"<<G4endl;
  // G4cout << "The mother molecule state : " << G4endl;
  // theMotherMolecule -> PrintState();

  if(moleculeDefinition->GetDecayTable())
  {
    const vector<const G4MolecularDissociationChannel*>* DecayVector =
        (theMotherMolecule->GetDecayChannel());

    if(DecayVector == 0)
    {
      G4ExceptionDescription exceptionDescription;
      theMotherMolecule->PrintState();
      exceptionDescription << "No decay channel was found for the molecule : "
                           << theMotherMolecule->GetName() << G4endl;
      G4Exception("G4DNAMolecularDissociation::DecayIt",
                  "G4DNAMolecularDissociation::NoDecayChannel",
                  FatalException,
                  exceptionDescription);
      return &aParticleChange;
    }

    G4int DecayVectorSize = DecayVector->size();
    // DEBUG
    // G4cout<< "Number of decay channels : " << DecayVectorSize <<G4endl;
    G4double RdmValue = G4UniformRand();

    const G4MolecularDissociationChannel* decayChannel(0);
    G4int i = 0;
    do
    {
      decayChannel = (*DecayVector)[i];
      if(RdmValue < decayChannel->GetProbability()) break;
      RdmValue -= decayChannel->GetProbability();
      i++;
    }
    while(i < DecayVectorSize);

    // DEBUG
    // G4cout << "Selected Decay channel : "
    //        << decayChannel->GetName() << G4endl;

    G4double decayEnergy = decayChannel->GetEnergy();
    G4int nbProducts = decayChannel->GetNbProducts();

    if(decayEnergy)
    {
      // DEBUG
      // G4cout << "Deposit energy :"
      //        << decayChannel->GetEnergy()/eV << " eV" << G4endl;

      aParticleChange.ProposeLocalEnergyDeposit(decayChannel->GetEnergy());
    }

    if(nbProducts)
    {

      // DEBUG
      // G4cout << "Number of products :" << nbProducts << G4endl;

      vector<G4ThreeVector> ProductsDisplacement(nbProducts);
      G4ThreeVector theMotherMoleculeDisplacement;

      DecayDisplacementMap::iterator it =
          fDecayDisplacementMap.find(moleculeDefinition);

      if(it != fDecayDisplacementMap.end())
      {
        G4VMolecularDecayDisplacer* displacer = it->second;
        ProductsDisplacement = displacer->GetProductsDisplacement(decayChannel);
        theMotherMoleculeDisplacement =
            displacer->GetMotherMoleculeDisplacement(decayChannel);
      }
      else
      {
        G4ExceptionDescription errMsg;
        errMsg << "No G4MolecularDecayProcess::theDecayDisplacementMap["
               << theMotherMolecule->GetName() + "]";
        G4Exception("G4MolecularDecayProcess::DecayIt",
                    "DNAMolecularDecay001",
                    FatalErrorInArgument,
                    errMsg);
      }

      aParticleChange.SetNumberOfSecondaries(nbProducts);

//      G4cout << " nbProducts = " << nbProducts << G4endl;
//      theMotherMolecule->PrintState();

#ifdef G4VERBOSE
      if(fVerbose)
      {
        G4cout << "Decay Process : " << theMotherMolecule->GetName()
               << " (trackID :" << track.GetTrackID() << ") "
               << decayChannel->GetName() << G4endl;
      }
#endif

      G4ITNavigator* navigator =
          G4ITTransportationManager::GetTransportationManager()
              ->GetNavigatorForTracking();

      for(G4int j = 0; j < nbProducts; j++)
      {
        G4Molecule* product = new G4Molecule(decayChannel->GetProduct(j));

        G4ThreeVector displacement = theMotherMoleculeDisplacement
            + ProductsDisplacement[j];
        double mag_displacement = displacement.mag();
        G4ThreeVector displacement_direction = displacement / mag_displacement;

        double prNewSafety = DBL_MAX;

        navigator->CheckNextStep(track.GetPosition(),
                                 displacement_direction,
                                 mag_displacement,
                                 prNewSafety); // returns a value

        if(prNewSafety < mag_displacement) mag_displacement = prNewSafety;

//        const G4AffineTransform& transform = navigator
//            ->GetGlobalToLocalTransform();
//
//        G4ThreeVector localPoint =
//            transform.TransformPoint(track.GetPosition());
//
//        if(track.GetTouchable()->GetSolid()->Inside(localPoint) != EInside::kInside)
//        {
//          G4Exception("G4DNAMolecularDissociation::DecayIt",
//                      "OUTSIDE_OF_MOTHER_VOLUME",
//                      FatalException,
//                      "Product has been placed outside of the volume "
//                      "containing the mother molecule");
//        }

        G4Track* secondary =
            product->BuildTrack(track.GetGlobalTime(),
                                track.GetPosition() + displacement_direction
                                    * mag_displacement);

        secondary->SetTrackStatus(fAlive);
#ifdef G4VERBOSE
        if(fVerbose)
        {
          G4cout << "Product : " << product->GetName() << G4endl;
        }
#endif
        // add the secondary track in the List
        aParticleChange.G4VParticleChange::AddSecondary(secondary);
      }
#ifdef G4VERBOSE
      if(fVerbose) G4cout << "-------------" << G4endl;
#endif
    }
    //DEBUG
    else if(fVerbose && decayEnergy)
    {
      G4cout << "No products for this channel" << G4endl;
      G4cout<<"-------------"<<G4endl;
    }
    /*
     else if(!decayEnergy && !nbProducts)
     {
     G4ExceptionDescription errMsg;
     errMsg << "There is no products and no energy specified in the molecular "
     "dissociation channel";
     G4Exception("G4MolecularDissociationProcess::DecayIt",
     "DNAMolecularDissociation002",
     FatalErrorInArgument,
     errMsg);
     }
     */
  }

  aParticleChange.ProposeTrackStatus(fStopAndKill);

  return &aParticleChange;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetDecayDisplacer()

G4VMolecularDecayDisplacer * G4DNAMolecularDissociation::GetDecayDisplacer

(

const G4ParticleDefinition *

molDef

)

Definition at line 345 of file G4DNAMolecularDissociation.cc.

{
  return fDecayDisplacementMap[molDef];
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetMeanFreePath()

virtual G4double G4DNAMolecularDissociation::GetMeanFreePath ( const G4Track &  , G4double  , G4ForceCondition *    )

inlineprotectedvirtual

Definition at line 112 of file G4DNAMolecularDissociation.hh.

    {
      return 0;
    }

Here is the call graph for this function:

GetMeanLifeTime()

G4double G4DNAMolecularDissociation::GetMeanLifeTime ( const G4Track &  track, G4ForceCondition *    )

protectedvirtual

Definition at line 132 of file G4DNAMolecularDissociation.cc.

{
  G4double output = GetMolecule(track)->GetDecayTime() - track.GetProperTime();
  return (output > 0 ? output : 0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

IsApplicable()

G4bool G4DNAMolecularDissociation::IsApplicable ( const G4ParticleDefinition &  aParticleType )

virtual

Definition at line 109 of file G4DNAMolecularDissociation.cc.

{
  if(aParticleType.GetParticleType() == "Molecule")
  {
#ifdef G4VERBOSE

    if(fVerbose > 1)
    {
      G4cout << "G4MolecularDissociation::IsApplicable(";
      G4cout << aParticleType.GetParticleName() << ",";
      G4cout << aParticleType.GetParticleType() << ")" << G4endl;
    }
#endif
    return (true);
  }
  else
  {
    return false;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

G4DNAMolecularDissociation& G4DNAMolecularDissociation::operator= ( const G4DNAMolecularDissociation &  right )

private

Here is the caller graph for this function:

PostStepDoIt()

G4VParticleChange* G4DNAMolecularDissociation::PostStepDoIt ( const G4Track &  track, const G4Step &  step  )

inline

Definition at line 93 of file G4DNAMolecularDissociation.hh.

    {
      return AtRestDoIt(track, step);
    }

Here is the call graph for this function:

PostStepGetPhysicalInteractionLength()

G4double G4DNAMolecularDissociation::PostStepGetPhysicalInteractionLength ( const G4Track &  track, G4double  previousStepSize, G4ForceCondition *  condition  )

virtual

Definition at line 353 of file G4DNAMolecularDissociation.cc.

{
  return 0; //-1*(-log(1-G4UniformRand())*100*1e-15*s);
}

Here is the caller graph for this function:

SetDecayDisplacer()

void G4DNAMolecularDissociation::SetDecayDisplacer	(	const G4ParticleDefinition *	molDef,
		G4VMolecularDecayDisplacer *	aDisplacer
	)

Definition at line 335 of file G4DNAMolecularDissociation.cc.

{
  fDecayDisplacementMap[molDef] = aDisplacer;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetVerbose()

void G4DNAMolecularDissociation::SetVerbose ( G4int  verbose )

inline

Definition at line 133 of file G4DNAMolecularDissociation.hh.

{
    fVerbose = verbose ;
}

Here is the caller graph for this function:

Member Data Documentation

fDecayAtFixedTime

G4bool G4DNAMolecularDissociation::fDecayAtFixedTime

private

Definition at line 125 of file G4DNAMolecularDissociation.hh.

fDecayDisplacementMap

DecayDisplacementMap G4DNAMolecularDissociation::fDecayDisplacementMap

private

Definition at line 128 of file G4DNAMolecularDissociation.hh.

fVerbose

G4int G4DNAMolecularDissociation::fVerbose

private

Definition at line 130 of file G4DNAMolecularDissociation.hh.

---

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

• G4DNAMolecularDissociation.hh
• G4DNAMolecularDissociation.cc