ActarSim.cc

/////////////////////////////////////////////////////////////////
//*-- AUTHOR : Hector Alvarez-Pol (hapol@fpddux.usc.es)
//             Beatriz Fernandez (bfd@ns.ph.liv.ac.uk)
//             Esther Estevez (hachebarra79@yahoo.es)
//             Elisangela Benjamim (ebenjami@usc.es)
//*-- Date: 03/2006
//*-- Last Update: 22/04/20
//*-- Copyright: GENP (Univ. Santiago de Compostela)
//
// --------------------------------------------------------------
// Comments:
//    ACTARSIM is an application for the simulation of the
//    ACTive TArget detector.
//
// --------------------------------------------------------------
/////////////////////////////////////////////////////////////////

#include <time.h>

#include "G4RunManager.hh"
#include "G4UImanager.hh"
#include "G4UIterminal.hh"
#include "G4UItcsh.hh"
#include "G4UIQt.hh"

#include "Randomize.hh"

#include "G4VisManager.hh"
#include "G4UIExecutive.hh"

#include "ActarSimDetectorConstruction.hh"
#include "ActarSimPhysicsList.hh"
#include "ActarSimPrimaryGeneratorAction.hh"
#include "ActarSimRunAction.hh"
#include "ActarSimEventAction.hh"
#include "ActarSimSteppingAction.hh"
#include "ActarSimSteppingVerbose.hh"
#include "ActarSimROOTAnalysis.hh"
#include "ActarSimVisManager.hh"

int main(int argc,char** argv)
{
  // choose the Random engine
  CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine);
  G4long seed=time(0); //returns time in seconds as an integer
  CLHEP::HepRandom::setTheSeed(seed);//changes the seed of the random engine

  //my Verbose output class
  G4VSteppingVerbose::SetInstance(new ActarSimSteppingVerbose);

  // Construct the default run manager
  G4RunManager* runManager = new G4RunManager;

  // set mandatory initialization classes
  ActarSimDetectorConstruction* detector = new ActarSimDetectorConstruction;
  runManager->SetUserInitialization(detector);
  runManager->SetUserInitialization(new ActarSimPhysicsList);

  // histogramming
  //ActarSimROOTAnalysis *analysis = new ActarSimROOTAnalysis(detector);
  ActarSimROOTAnalysis *analysis = new ActarSimROOTAnalysis();

  G4UIsession* session=0;

  if (argc==1){   // Define UI session for interactive mode.
    session = new G4UIQt(argc,argv);
  }
  // visualization manager
  G4VisManager* visManager = new ActarSimVisManager;
  visManager->Initialize();

  // set mandatory user action class
  runManager->SetUserAction(new ActarSimPrimaryGeneratorAction);
  runManager->SetUserAction(new ActarSimRunAction);
  ActarSimEventAction* eventaction = new ActarSimEventAction;
  runManager->SetUserAction(eventaction);
  runManager->SetUserAction(new ActarSimSteppingAction(detector,eventaction));

  // Initialize G4 kernel -->Make it manually to allow the definition of
  //    commands in PreInit state (for instance to define the PhysicsList)
  //runManager->Initialize();

  // get the pointer to the UI manager and set verbosities
  G4UImanager* UI = G4UImanager::GetUIpointer();

  if (session) {
    // G4UIterminal is a (dumb) terminal.
    //UI->ApplyCommand("/control/execute vis2.mac");
    UI->ApplyCommand("/control/execute gui.mac");
    session->SessionStart();
    if ( (G4UImanager::GetUIpointer()) &&
	 (session ==(G4UImanager::GetUIpointer())->GetSession()) )
      delete session;
  }
  else {
    //This part should be used for batch
    G4String command = "/control/execute ";
    G4String fileName = argv[1];
    UI->ApplyCommand(command+fileName);
  }

  // job termination
  delete analysis;
  delete visManager;
  G4cout << "RunManager deleting... "<< G4endl;
  delete runManager;

return 0;
}