MJGeneratorPrimaryMessenger.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//                                                                           //
//                         MAJORANA Simulation                               //
//                                                                           //
//      This code implementation is the intellectual property of the         //
//      MAJORANA Collaboration. It is based on Geant4, an intellectual       //
//      property of the RD44 GEANT4 collaboration.                           //
//                                                                           //
//                        *********************                              //
//                                                                           //
//    Neither the authors of this software system, nor their employing       //
//    institutes, nor the agencies providing financial support for this      //
//    work  make  any representation or  warranty, express or implied,       //
//    regarding this software system or assume any liability for its use.    //
//    By copying, distributing or modifying the Program (or any work based   //
//    on on the Program) you indicate your acceptance of this statement,     //
//    and all its terms.                                                     //
//                                                                           //
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//---------------------------------------------------------------------------//
//                                                          
// $Id: MJGeneratorPrimaryMessenger.cc,v 1.4 2004/11/15 09:53:56 xliu Exp $ 
//      
// CLASS IMPLEMENTATION:  MJGeneratorPrimaryMessenger.cc
//
//---------------------------------------------------------------------------//
// 
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//
//

#include "globals.hh"
#include "G4UIcmdWithoutParameter.hh"
#include "G4UIcmdWithAString.hh"
#include "G4UIdirectory.hh"
#include "G4PhysicalVolumeStore.hh"

#include "generators/MJGeneratorPNNL.hh"
#include "generators/MJGeneratorPrimary.hh"
#include "generators/MJGeneratorRDM.hh"
#include "generators/MJGeneratorG4Gun.hh"
#include "generators/MJGeneratorTUNLFEL.hh"
#include "generators/MaGeGeneratorDecay0.hh"

#include "io/MJLogger.hh"
//---------------------------------------------------------------------------//

#include "generators/MJGeneratorPrimaryMessenger.hh" 

//---------------------------------------------------------------------------//

MJGeneratorPrimaryMessenger::MJGeneratorPrimaryMessenger
(MJGeneratorPrimary *generator):
  fGeneratorPrimary(generator)
{
  // /MJ/generator/ directory definition
  fGeneratorDirectory = new G4UIdirectory("/MJ/generator/");
  fGeneratorDirectory->SetGuidance("Control commands for generators:");
  fGeneratorDirectory->SetGuidance("/MJ/generator/select: Select generator.");
  
  // /MJ/generator/select command
  fSelectCmd = new G4UIcmdWithAString("/MJ/generator/select", this);
  fSelectCmd->SetGuidance("Selects generator for events.");
  fSelectCmd->SetGuidance("Options are:");
  fSelectCmd->SetGuidance("PNNLiso : PNNL model for radioactive isotopes.");
  fSelectCmd->SetGuidance("RDMiso: G4 model for radioactive isotopes.");
  fSelectCmd->SetGuidance("TUNLFEL: TUNL FEL beam.");
  fSelectCmd->SetGuidance("G4gun: Standard G4 gun.");
  fSelectCmd->SetGuidance("decay0: decay0 generator.");
  fSelectCmd->SetCandidates("PNNLiso RDMiso TUNLFEL G4gun decay0");

  // /MJ/generator/name
  fNameCmd = new G4UIcmdWithoutParameter("/MJ/generator/name", this);
  fNameCmd->SetGuidance("Returns name of current generator.");

  // /MJ/generator/confine
  fConfineCmd = new G4UIcmdWithAString("/MJ/generator/confine", this);
  fConfineCmd->SetGuidance("Selects confinement for the source.");
  fConfineCmd->SetGuidance("Options are:");
  fConfineCmd->SetGuidance("noconfined : source not confined");
  fConfineCmd->SetGuidance("volume : source confined in a volume.");
  fConfineCmd->SetGuidance("surface : source confined in a surface.");
  fConfineCmd->SetCandidates("noconfined volume surface");

  // /MJ/generator/volume
  fVolumeCmd = new G4UIcmdWithAString("/MJ/generator/volume", this);
  fVolumeCmd->SetGuidance("Selects the volume where the source is confined");
  fVolumeCmd->AvailableForStates(G4State_Init,G4State_Idle);
}

//---------------------------------------------------------------------------//

MJGeneratorPrimaryMessenger::MJGeneratorPrimaryMessenger(const MJGeneratorPrimaryMessenger & other)
{;}

//---------------------------------------------------------------------------//

MJGeneratorPrimaryMessenger::~MJGeneratorPrimaryMessenger()
{
  delete fGeneratorDirectory;
  delete fNameCmd;
  delete fSelectCmd;
  delete fConfineCmd;
  delete fVolumeCmd;
}

//---------------------------------------------------------------------------//

G4String MJGeneratorPrimaryMessenger::GetCurrentValue(G4UIcommand *command)
{
  G4String cv;

  // /MJ/generator/name
  if(command == fNameCmd) {
    cv = "Current generator is : " + 
      fGeneratorPrimary->GetMJGenerator()->GetGeneratorName();
  }

  return cv;
}

//---------------------------------------------------------------------------//

void MJGeneratorPrimaryMessenger::SetNewValue(G4UIcommand *command,
                                              G4String newValues)
{
  // /MJ/generator/select
  if(command == fSelectCmd) {
    if(newValues == "PNNLiso") 
      fGeneratorPrimary->SetMJGenerator(new MJGeneratorPNNL);
    else if(newValues == "RDMiso")
      fGeneratorPrimary->SetMJGenerator(new MJGeneratorRDM);
    else if(newValues == "G4gun")
      fGeneratorPrimary->SetMJGenerator(new MJGeneratorG4Gun);
    else if(newValues == "TUNLFEL")
      fGeneratorPrimary->SetMJGenerator(new MJGeneratorTUNLFEL);
    else if(newValues == "decay0")
      fGeneratorPrimary->SetMJGenerator(new MaGeGeneratorDecay0);
    else 
      MJLog(warning) << "Unknown generator. Please try again." << endlog;
  }

  if(command == fConfineCmd) {
    if(newValues == "noconfined")
      {
        fGeneratorPrimary->SetConfinementCode(noconfined);
        MJLog(trace) << "Source not confined" << endlog;
      }
    else if(newValues == "volume")
      {
        fGeneratorPrimary->SetConfinementCode(volume);
        MJLog(trace) << "Source confined in volume" << endlog;
      }
    else if(newValues == "surface")
      {
        fGeneratorPrimary->SetConfinementCode(surface);
        MJLog(trace) << "Source confined in skin surface" << endlog;
      }
  }

  if (command == fVolumeCmd)
    {
      MJLog(trace) << " choose volume command "<<endlog;
      fGeneratorPrimary->SetVolumeName(newValues);

      G4bool ifound = false;
      G4PhysicalVolumeStore* volumeStore = G4PhysicalVolumeStore::GetInstance();
      G4int nVolumes = (G4int) volumeStore->size();
      G4String candidateList;
      for(G4int i=0;i<nVolumes;i++)
        {
          candidateList += (*volumeStore)[i]->GetName();
          candidateList += ", ";
          if ((*volumeStore)[i]->GetName() == newValues) ifound = true;
        }

      if (ifound) {
        if (fGeneratorPrimary->GetConfinementCode() != noconfined) 
          {
            fGeneratorPrimary->SetVolumeName(newValues);
            MJLog(trace) << "Source confined in " << newValues << endlog;
          }
        else 
          MJLog(warning) << "Source not confined: nothing happens " << endlog;
      }
      else
        {
          MJLog(warning) << "Volume not found " << endlog;
          MJLog(warning) << "The list of volumes is: "
                         << candidateList << endlog;
        }

    }

}

//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//

---