BesMagneticFieldMessenger.cc

Go to the documentation of this file.

#include "BesMagneticFieldMessenger.hh"
#include "BesMagneticField.hh"
 
#include "G4UIcmdWithADouble.hh"
#include "G4UIcmdWithADoubleAndUnit.hh"
#include "G4UIcmdWithAString.hh"
#include "G4UIcmdWithAnInteger.hh"
#include "G4UIcmdWithoutParameter.hh"
#include "G4UIdirectory.hh"
 
// The standard constructor
BesMagneticFieldMessenger::BesMagneticFieldMessenger( BesMagneticField* pEMfieldSetup )
    : fEMfieldSetup( pEMfieldSetup ) {
  F01detDir = new G4UIdirectory( "/field/" );
  F01detDir->SetGuidance( "Bes MagneticField tracking control." );
 
  StepperCmd = new G4UIcmdWithAnInteger( "/field/SetStepperType", this );
  StepperCmd->SetGuidance( "Select stepper type for magnetic field" );
  StepperCmd->SetParameterName( "choice", true );
  StepperCmd->SetDefaultValue( 4 );
  StepperCmd->SetRange( "choice>=0&&choice<=9" );
 
  MinStepCmd = new G4UIcmdWithADoubleAndUnit( "/field/SetMinStep", this );
  MinStepCmd->SetGuidance( "Please change it carefully,it is better keep it const" );
  MinStepCmd->SetParameterName( "MinStep", true );
  MinStepCmd->SetDefaultUnit( "mm" );
  MinStepCmd->SetDefaultValue( 0.01 );
  MinStepCmd->SetRange( "MinStep>=0&&MinStep<=1" );
 
  DeltaOneStepCmd = new G4UIcmdWithADoubleAndUnit( "/field/SetDeltaOneStep", this );
  DeltaOneStepCmd->SetGuidance( "Set DeltaOneStep for tracking" );
  DeltaOneStepCmd->SetParameterName( "DeltaOneStep", true );
  DeltaOneStepCmd->SetDefaultUnit( "mm" );
  DeltaOneStepCmd->SetDefaultValue( 0.01 );
  DeltaOneStepCmd->SetRange( "DeltaOneStep>=0" );
 
  DeltaIntersectionCmd = new G4UIcmdWithADoubleAndUnit( "/field/SetDeltaIntersection", this );
  DeltaIntersectionCmd->SetGuidance( "Set DeltaIntersection for tracking" );
  DeltaIntersectionCmd->SetParameterName( "DeltaIntersection", true );
  DeltaIntersectionCmd->SetDefaultUnit( "mm" );
  DeltaIntersectionCmd->SetDefaultValue( 0.001 );
  DeltaIntersectionCmd->SetRange( "DeltaIntersection>=0" );
 
  MinimumEpsilonStepCmd = new G4UIcmdWithADouble( "/field/SetMinimumEpsilonStep", this );
  MinimumEpsilonStepCmd->SetGuidance( "Set Minimum Epsilon Step for tracking" );
  MinimumEpsilonStepCmd->SetParameterName( "MinimumEpsilonStep", true );
  MinimumEpsilonStepCmd->SetDefaultValue( 5.0e-5 );
  MinimumEpsilonStepCmd->SetRange( "MinimumEpsilonStep>=0" );
 
  MaximumEpsilonStepCmd = new G4UIcmdWithADouble( "/field/SetMaximumEpsilonStep", this );
  MaximumEpsilonStepCmd->SetGuidance( "Set Maximum Epsilon Step for tracking" );
  MaximumEpsilonStepCmd->SetParameterName( "MaximumEpsilonStep", true );
  MaximumEpsilonStepCmd->SetDefaultValue( 1.0e-3 );
  MaximumEpsilonStepCmd->SetRange( "MaximumEpsilonStep>=0" );
}
 
BesMagneticFieldMessenger::~BesMagneticFieldMessenger() {
  delete MaximumEpsilonStepCmd;
  delete MinimumEpsilonStepCmd;
  delete DeltaIntersectionCmd;
  delete DeltaOneStepCmd;
  delete MinStepCmd;
  delete StepperCmd;
  delete F01detDir;
}
 
void BesMagneticFieldMessenger::SetNewValue( G4UIcommand* command, G4String newvalue ) {
  // Please don't change the order of code ,it will make some parameters useless if you do
  // that.
  if ( command == StepperCmd )
  {
    fEMfieldSetup->SetStepperType( StepperCmd->GetNewIntValue( newvalue ) );
    fEMfieldSetup->CreateStepperAndChordFinder();
  }
  if ( command == MinStepCmd )
  {
    fEMfieldSetup->SetMinStep( MinStepCmd->GetNewDoubleValue( newvalue ) );
    fEMfieldSetup->CreateStepperAndChordFinder();
  }
  if ( command == DeltaOneStepCmd )
  { fEMfieldSetup->SetDeltaOneStep( DeltaOneStepCmd->GetNewDoubleValue( newvalue ) ); }
  if ( command == DeltaIntersectionCmd )
  {
    fEMfieldSetup->SetDeltaIntersection( DeltaIntersectionCmd->GetNewDoubleValue( newvalue ) );
  }
  if ( command == MinimumEpsilonStepCmd )
  {
    fEMfieldSetup->SetMinimumEpsilonStep(
        MinimumEpsilonStepCmd->GetNewDoubleValue( newvalue ) );
  }
  if ( command == MaximumEpsilonStepCmd )
  {
    fEMfieldSetup->SetMaximumEpsilonStep(
        MaximumEpsilonStepCmd->GetNewDoubleValue( newvalue ) );
  }
}