#include <G4MoleculeTableMessenger.hh>

Inheritance diagram for G4MoleculeTableMessenger:

Public Member Functions
	G4MoleculeTableMessenger ()

	~G4MoleculeTableMessenger () override

void	SetNewValue (G4UIcommand *command, G4String newValue) override

Public Member Functions inherited from G4UImessenger
	G4UImessenger ()=default

	G4UImessenger (const G4String &path, const G4String &dsc, G4bool commandsToBeBroadcasted=true)

virtual	~G4UImessenger ()

virtual G4String	GetCurrentValue (G4UIcommand *command)

G4bool	CommandsShouldBeInMaster () const

Protected Attributes
std::unique_ptr< G4UIcmdWithoutParameter >	fpPrintTable

std::unique_ptr< G4UIcmdWithAString >	fpSpecies

Protected Attributes inherited from G4UImessenger
G4UIdirectory *	baseDir = nullptr

G4String	baseDirName = ""

G4bool	commandsShouldBeInMaster = false

Additional Inherited Members
Protected Member Functions inherited from G4UImessenger
G4String	ItoS (G4int i)

G4String	LtoS (G4long l)

G4String	DtoS (G4double a)

G4String	BtoS (G4bool b)

G4int	StoI (const G4String &s)

G4long	StoL (const G4String &s)

G4double	StoD (const G4String &s)

G4bool	StoB (G4String s)

void	AddUIcommand (G4UIcommand *newCommand)

void	CreateDirectory (const G4String &path, const G4String &dsc, G4bool commandsToBeBroadcasted=true)

template<typename T >
T *	CreateCommand (const G4String &cname, const G4String &dsc)

Detailed Description

Definition at line 41 of file G4MoleculeTableMessenger.hh.

Constructor & Destructor Documentation

◆ G4MoleculeTableMessenger()

G4MoleculeTableMessenger::G4MoleculeTableMessenger ( )

Definition at line 37 of file G4MoleculeTableMessenger.cc.

  : 
   fpPrintTable(new G4UIcmdWithoutParameter("/chem/PrintSpeciesTable", this))
  , fpSpecies(new G4UIcmdWithAString("/chem/species", this))
{}

◆ ~G4MoleculeTableMessenger()

G4MoleculeTableMessenger::~G4MoleculeTableMessenger ( )

overridedefault

Member Function Documentation

◆ SetNewValue()

void G4MoleculeTableMessenger::SetNewValue	(	G4UIcommand *	command,
		G4String	newValue )

overridevirtual

Reimplemented from G4UImessenger.

Definition at line 46 of file G4MoleculeTableMessenger.cc.

{
  if(command == fpPrintTable.get())
  {
    G4MolecularConfiguration::PrintAll();
  }
  if(command == fpSpecies.get())
  {
    std::istringstream iss(newValue);
 
    G4String speciesName;
    iss >> speciesName;
 
    G4String marker;
    iss >> marker;  // must be [
 
    if(marker != "[")
    {
      G4ExceptionDescription description;
      description << " marker : " << marker << G4endl;
      G4Exception("G4MoleculeTableMessenger::SetNewValue",
                  "FAIL_SPECIES_DEFINITION04", FatalException, description);
    }
 
    G4String speciesDefName;
    iss >> speciesDefName;
 
    iss >> marker;  // peut etre |
 
    G4int charge = 0;
    if(marker == "|")
    {
      iss >> charge;
    }
    iss >> marker;  // peut etre |
    G4double diffusion_coefficient = 0;
    if(marker == "|")
    {
      iss >> diffusion_coefficient;
    }
    iss >> marker;  // peut etre |
    G4double VanDerVaalsRadius = 0;
    if(marker == "|")
    {
      iss >> VanDerVaalsRadius;  // in nm
    }
 
    auto pConf =
      G4MolecularConfiguration::GetMolecularConfiguration(speciesName);
    if(pConf != nullptr)
    {
      pConf->UnFinalize();
      if(VanDerVaalsRadius != 0)
      {
        pConf->SetVanDerVaalsRadius(VanDerVaalsRadius * nm);
      }
      if(diffusion_coefficient != 0)
      {
        pConf->SetDiffusionCoefficient(diffusion_coefficient * (m2 / s));
      }
    }
    else
    {
      auto particleDef = dynamic_cast<G4MoleculeDefinition*>(
        G4ParticleTable::GetParticleTable()->FindParticle(speciesDefName));
      if(particleDef != nullptr)
      {
        auto molConf =
          G4MolecularConfiguration::GetOrCreateMolecularConfiguration(
            particleDef, charge);
        if(molConf == nullptr)
        {
          G4ExceptionDescription description;
          description << "This molecule has not been defined" << G4endl;
          G4Exception("G4MoleculeTableMessenger::SetNewValue",
                      "FAIL_SPECIES_DEFINITION02", FatalException, description);
        }
        molConf->UnFinalize();
        if(VanDerVaalsRadius != 0)
        {
          molConf->SetVanDerVaalsRadius(VanDerVaalsRadius * nm);
        }
        if(diffusion_coefficient != 0)
        {
          molConf->SetDiffusionCoefficient(diffusion_coefficient * (m2 / s));
        }
 
        auto usedName = molConf->GetUserID();
        if(!usedName.empty())
        {
          molConf->PrintState();
          G4ExceptionDescription description;
          description << "This molecule has been defined by the name : "
                      << usedName << " . Please, use this name." << G4endl;
          G4Exception("G4MoleculeTableMessenger::SetNewValue",
                      "FAIL_SPECIES_DEFINITION", FatalException, description);
        }
        else
        {
          molConf->SetUserID(speciesName);
        }
      }
      else
      {
        auto speciesDef =
          new G4MoleculeDefinition(speciesDefName, /*mass no needed*/ 0,
                                   /*D*/ diffusion_coefficient * (m * m / s),
                                   /*charge*/ charge, 1,
                                   /*electronL*/ 0,
                                   /*radius*/ VanDerVaalsRadius * nm);
        G4bool alreadyCreated(false);
        G4MolecularConfiguration::CreateMolecularConfiguration(
          speciesName, speciesDef, alreadyCreated);
      }
    }
  }
}

Member Data Documentation

◆ fpPrintTable

std::unique_ptr<G4UIcmdWithoutParameter> G4MoleculeTableMessenger::fpPrintTable

protected

Definition at line 49 of file G4MoleculeTableMessenger.hh.

Referenced by SetNewValue().

◆ fpSpecies

std::unique_ptr<G4UIcmdWithAString> G4MoleculeTableMessenger::fpSpecies

protected

Definition at line 50 of file G4MoleculeTableMessenger.hh.