#include <G4ReactionTableMessenger.hh>

Inheritance diagram for G4ReactionTableMessenger:

Public Member Functions
	G4ReactionTableMessenger (G4DNAMolecularReactionTable *)

	~G4ReactionTableMessenger () 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
G4DNAMolecularReactionTable *	fpTable

std::unique_ptr< G4UIcmdWithoutParameter >	fpActivateReactionUI

G4UIcmdWithAString *	fpAddReaction

G4UIcmdWithAString *	fpNewDiffContReaction

G4UIcmdWithoutParameter *	fpPrintTable

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 45 of file G4ReactionTableMessenger.hh.

Constructor & Destructor Documentation

◆ G4ReactionTableMessenger()

G4ReactionTableMessenger::G4ReactionTableMessenger ( G4DNAMolecularReactionTable * table )

Definition at line 44 of file G4ReactionTableMessenger.cc.

  : 
   fpTable(table)
  , fpActivateReactionUI(new G4UIcmdWithoutParameter("/chem/reaction/UI", this))
{
  fpNewDiffContReaction = new G4UIcmdWithAString("/chem/reaction/new", this);
  fpAddReaction = new G4UIcmdWithAString("/chem/reaction/add", this);
  fpPrintTable = new G4UIcmdWithoutParameter("/chem/reaction/print", this);
}

◆ ~G4ReactionTableMessenger()

G4ReactionTableMessenger::~G4ReactionTableMessenger ( )

override

Definition at line 56 of file G4ReactionTableMessenger.cc.

{
  delete fpNewDiffContReaction;
  delete fpAddReaction;
  delete fpPrintTable;
}

Member Function Documentation

◆ SetNewValue()

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

overridevirtual

Reimplemented from G4UImessenger.

Definition at line 64 of file G4ReactionTableMessenger.cc.

{
  if(command == fpActivateReactionUI.get())
  {
    //assert(false);
      fpTable->Reset();//release reaction data
  }
 
  if(command == fpNewDiffContReaction)
  {
    std::istringstream iss(newValue);
 
    G4String species1;
    iss >> species1;
 
    G4String species2;
    iss >> species2;
 
    double reactionRate;
    iss >> reactionRate;
 
//    G4String reactionRateUnit;
//    iss >> reactionRateUnit;
 
    double dimensionedReactionRate = reactionRate * (1e-3 * m3 / (mole * s));
//        G4UIcmdWithADoubleAndUnit::ConvertToDimensionedDouble((reactionRate
//            + G4String(" ") + reactionRateUnit).c_str());
 
    auto reactionData =
        new G4DNAMolecularReactionData(dimensionedReactionRate,
                                       species1,
                                       species2);
 
//    G4String productionRate;
//    iss >> productionRate;
//
//    if(productionRate != "" && productionRate != "X")
//    {
//      double prodRateReal = G4UIcommand::ConvertToDouble(productionRate);
//
//      if(prodRateReal == 0 || isnan(prodRateReal))
//      {
//        G4Exception("G4ReactionTableMessenger",
//                    "WRONG_PRODUCTION_RATE",
//                    FatalException, "");
//      }
//
//      double dimensionedProductionRate = prodRateReal
//          * (1e-3 * m3 / (mole * s));
//      reactionData->SetProductionRate(dimensionedProductionRate);
//    }
 
    while(!iss.eof())
    {
      G4String product;
      iss >> product;
 
      if(!product.empty())
      {
        reactionData->AddProduct(product);
      }
      else
      {
        break;
      }
    };
 
    fpTable->SetReaction(reactionData);
  }
//  else if(command == fpNewPartDiffContReactionByRadius)
//  {
//    std::istringstream iss(newValue);
//
//    G4String species1;
//    iss >> species1;
//
//    G4String species2;
//    iss >> species2;
//
//    double reactionRate;
//    iss >> reactionRate;
//
////    G4String reactionRateUnit;
////    iss >> reactionRateUnit;
//
////    G4String reactionRateUnit;
////    iss >> reactionRateUnit;
//
//    double reactionRadius;
//    iss >> reactionRadius;
//
////    G4String reactionRadiusUnit;
////    iss >> reactionRadiusUnit;
//
//    double dimensionedReactionRate = reactionRate * (1e-3 * m3 / (mole * s));
//
////    double dimensionedReactionRate =
////        G4UIcmdWithADoubleAndUnit::ConvertToDimensionedDouble((reactionRate
////            + " " + reactionRateUnit).c_str());
//
//    double dimensionedReactionRadius = reactionRadius * nm;
////        G4UIcmdWithADoubleAndUnit::ConvertToDimensionedDouble((reactionRadius
////            + " " + reactionRadiusUnit).c_str());
//
//    G4DNAMolecularReactionData* reactionData =
//        new G4DNAMolecularReactionData(dimensionedReactionRate,
//                                       species1,
//                                       species2);
//    reactionData->SetPartiallyDiffusionControlledReaction(dimensionedReactionRate,
//                                                          dimensionedReactionRadius);
//
//    while(iss.eof() == false)
//    {
//      G4String product;
//      iss >> product;
//
//      if(product != "")
//      {
//        reactionData->AddProduct(product);
//      }
//      else
//      {
//        break;
//      }
//    }
//
//    fpTable->SetReaction(reactionData);
//  }
//  else if(command == fpNewPartDiffContReactionByReactionRate)
//  {
//    std::istringstream iss(newValue);
//
//    G4String species1;
//    iss >> species1;
//
//    G4String species2;
//    iss >> species2;
//
//    double reactionRate;
//    iss >> reactionRate;
//
//    //    G4String reactionRateUnit;
//    //    iss >> reactionRateUnit;
//
//    //    G4String reactionRateUnit;
//    //    iss >> reactionRateUnit;
//
//    double activationRate;
//    iss >> activationRate;
//
//    //    G4String reactionRadiusUnit;
//    //    iss >> reactionRadiusUnit;
//
//    double dimensionedReactionRate = reactionRate * (1e-3 * m3 / (mole * s));
//
//    double dimensionedActivationRate = activationRate
//        * (1e-3 * m3 / (mole * s));
//
//    G4DNAMolecularReactionData* reactionData =
//        new G4DNAMolecularReactionData(dimensionedReactionRate,
//                                       species1,
//                                       species2);
//    reactionData->
//      SetPartiallyDiffusionControlledReactionByActivation(dimensionedReactionRate,
//                                                          dimensionedActivationRate);
//
//    while(iss.eof() == false)
//    {
//      G4String product;
//      iss >> product;
//
//      if(product != "")
//      {
//        reactionData->AddProduct(product);
//      }
//      else
//      {
//        break;
//      }
//    }
//
//    fpTable->SetReaction(reactionData);
//  }
  else if(command == fpPrintTable)
  {
    fpTable->PrintTable();
  }
  else if(command == fpAddReaction)
  {
    std::istringstream iss(newValue);
 
    //--------------------------------------------------------------------------
    // Reactants definition
 
    G4String species1;
    iss >> species1;
 
    G4String marker;
    iss >> marker; // peut etre +, ->, |
 
    G4String species2;
 
    if(marker == "+")
    {
        iss >> species2;
        iss >> marker; // peut etre ->, |
    }
 
    //--------------------------------------------------------------------------
 
    auto  reactionData =
                new G4DNAMolecularReactionData(0,
                                               species1,
                                               species2);
    //fpTable->SetReaction(reactionData);
 
    //--------------------------------------------------------------------------
    // Add products
    if(marker == "->")
    {
      iss >> marker; // doit etre = species name
 
      while(marker!="|"
          //&& marker!=""
          && !iss.eof()
          )
      {
        //G4cout << marker << G4endl;
        if(marker == "+")
        {
          iss >> marker; // doit etre species name
          continue;
        }
        if(marker != "H2O")
        {
          reactionData->AddProduct(marker);
        }
 
        iss >> marker; // peut etre species name, +, |
      };
    }
 
//    G4cout << "out of 1st loop" << G4endl;
 
    //--------------------------------------------------------------------------
    // Add reaction rate method
    G4String rateconst_method;
    iss >> rateconst_method;
    if(rateconst_method == "Fix")
    {
      iss >> marker; // must be |
      double reactionRate;
      iss >> reactionRate;
 
      double dimensionedReactionRate = reactionRate * (1e-3 * m3 / (mole * s));
      reactionData->SetObservedReactionRateConstant(dimensionedReactionRate);
      reactionData->ComputeEffectiveRadius();
      G4String markerType;
      iss >> markerType; // must be |
      if(markerType == "|")
      {
        G4int reactionType;
        iss >> reactionType;
        if(reactionType == 1)
        {
          reactionData->SetReactionType(reactionType);
        }
      }
 
 
//      G4String productionRate;
//      iss >> productionRate;
//
//      if(productionRate != "" && productionRate != "X")
//      {
//        double prodRateReal = G4UIcommand::ConvertToDouble(productionRate);
//
//        if(prodRateReal == 0 || isnan(prodRateReal))
//        {
//          G4Exception("G4ReactionTableMessenger",
//                      "WRONG_PRODUCTION_RATE",
//                      FatalException,
//                      "");
//        }
//
//        double dimensionedProductionRate = prodRateReal
//            * (1e-3 * m3 / (mole * s));
//        reactionData->SetProductionRate(dimensionedProductionRate);
//      }
    }
    else if(rateconst_method == "Arr")
    {
      iss >> marker; // must be |
      double A0  = 0;
      double E_R = 0;
 
      iss >> A0;
      iss >> E_R;
      reactionData->SetArrehniusParameterization(A0, E_R);
    }
    else if(rateconst_method == "Pol")
    {
      iss >> marker; // must be |
      std::vector<double> P = {0, 0, 0, 0, 0};
 
      size_t i = 0;
      while(i < 4) // could be changed to 5 only if marker is used as delimiter
      {
        double tmp;
        iss >> tmp;
        P[i] = tmp;
//        G4cout << newValue << G4endl;
//        G4cout << tmp << G4endl;
//        G4cout << P[i] << G4endl;
        ++i;
      };
      reactionData->SetPolynomialParameterization(P);
    }
    else if(rateconst_method == "Scale")
    {
      iss >> marker; // must be |
      double temp_K;
      iss >> temp_K;
      double reactionRateCste;
      iss >> reactionRateCste;
      double dimensionedReactionRate = reactionRateCste * (1e-3 * m3 / (mole * s));
      reactionData->SetObservedReactionRateConstant(dimensionedReactionRate);
      reactionData->SetScaledParameterization(temp_K, dimensionedReactionRate);
    }
 
//    if(iss.eof() == false)
//    {
//      iss >> marker;
//
//      if(marker == "|")
//      {
//        G4String productionRate ;
//        iss >> productionRate;
//
////        G4cout << productionRate << G4endl;
//
//        double dimProductionRate = G4UIcommand::ConvertToDouble(productionRate)* (1e-3 * m3 / (mole * s));
//
//        G4cout << " DIM PROD RATE = " << reactionData->GetReactant1()->GetName()
//             << " + " << reactionData->GetReactant2()->GetName() << " = " << dimProductionRate << G4endl;
//
//        reactionData->SetProductionRate(dimProductionRate);
//      }
//    }
    fpTable->SetReaction(reactionData);
//    G4cout << "Reaction " << species1 << " + " << species2 << " added" << G4endl;
  }
}