QweakSimG4.cc File Reference

#include "G4RunManager.hh"
#include "G4UImanager.hh"
#include "G4UIterminal.hh"
#include "G4UItcsh.hh"
#include "G4Version.hh"
#include "G4StepLimiterBuilder.hh"
#include "G4OpticalPhysics.hh"
#include "G4PhysListFactory.hh"
#include "G4VModularPhysicsList.hh"
#include "G4UIExecutive.hh"
#include "QweakSimUserInformation.hh"
#include "QweakSimDetectorConstruction.hh"
#include "QweakSimPrimaryGeneratorAction.hh"
#include "QweakSimRunAction.hh"
#include "QweakSimEventAction.hh"
#include "QweakSimStackingAction.hh"
#include "QweakSimSteppingAction.hh"
#include "QweakSimTrackingAction.hh"
#include "QweakSimSteppingVerbose.hh"
#include "QweakSimAnalysis.hh"
#include "QweakSimEPEvent.hh"

Include dependency graph for QweakSimG4.cc:

Go to the source code of this file.

Macros
#define	USE_CUSTOM_NUCLEAR_SCATTERING   0

Functions
int	main (int argc, char **argv)

Detailed Description

Revision:

1.2

Date:

2005/12/27 19:09:32

Author

Klaus Hans Grimm

Definition in file QweakSimG4.cc.

Macro Definition Documentation

#define USE_CUSTOM_NUCLEAR_SCATTERING   0

Definition at line 63 of file QweakSimG4.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 71 of file QweakSimG4.cc.

                               {

  // Verbose output class
  //G4VSteppingVerbose::SetInstance( new QweakSimSteppingVerbose() ); 

  // Run manager
  G4RunManager * runManager = new G4RunManager();

  // UserInitialization classes (mandatory)
  QweakSimUserInformation*       myQweakSimUserInformation  = new QweakSimUserInformation();
  QweakSimDetectorConstruction*  myQweakSimExperiment       = new QweakSimDetectorConstruction(myQweakSimUserInformation);
  QweakSimAnalysis*              myQweakSimAnalysis         = new QweakSimAnalysis(myQweakSimUserInformation);

//jpan@nuclear.uwinnipeg.ca
  QweakSimEPEvent*               myEPEvent                  = new QweakSimEPEvent(myQweakSimUserInformation);

  runManager->SetUserInitialization(myQweakSimExperiment);


  // Calls a reference physics list for the simulation
  G4PhysListFactory factory;
  G4VModularPhysicsList* physlist = factory.GetReferencePhysList("QGSP_BERT_LIV");
  physlist->RegisterPhysics(new G4OpticalPhysics());
  #if G4VERSION_NUMBER < 1000
  physlist->RegisterPhysics(new G4StepLimiterBuilder());
  #else
  physlist->RegisterPhysics(new G4StepLimiterPhysics());
  #endif
  // Replace the standard EM with the customized version to add Pb A_T
  #if USE_CUSTOM_NUCLEAR_SCATTERING
  physlist->ReplacePhysics(new QweakSimEmLivermorePhysics());
  #endif

  runManager->SetUserInitialization(physlist);

  // Original Qweak Physics List, uncomment to use and add correct #include lines
  //runManager->SetUserInitialization(new QweakSimPhysicsList() );

  // UserAction classes
  runManager->SetUserAction( new QweakSimPrimaryGeneratorAction(myQweakSimUserInformation, myEPEvent) );
  //runManager->SetUserAction( new QweakSimPrimaryGeneratorAction( ) );
  runManager->SetUserAction( new QweakSimSteppingAction(myQweakSimUserInformation, myEPEvent));
  runManager->SetUserAction( new QweakSimStackingAction() );
  runManager->SetUserAction( new QweakSimTrackingAction(myQweakSimUserInformation) );
  runManager->SetUserAction( new QweakSimEventAction(myQweakSimAnalysis, myQweakSimUserInformation) );
  runManager->SetUserAction( new QweakSimRunAction(myQweakSimAnalysis) );  

  runManager->StoreRandomNumberStatusToG4Event(1);

  //Initialize G4 kernel
  runManager->Initialize();
  G4UIsession* session = 0;

 if (argc==1)   // Define UI session for interactive mode.
 {

     // G4UIterminal is a (dumb) terminal.
    #if defined(G4UI_USE_QT)
     session = new G4UIQt(argc,argv);
    #elif defined(G4UI_USE_XM)
     session = new G4UIXm(argc,argv);
    #elif defined(G4UI_USE_WIN32)
     session = new G4UIWin32();
    #elif defined(G4UI_USE_TCSH)
     session = new G4UIterminal(new G4UItcsh);
    #else
      session = new G4UIterminal();
    #endif

 }


#ifdef G4VIS_USE
  // Visualization, if you choose to have it!
  // 
  // Simple graded message scheme - give first letter or a digit:
  //  0) quiet,         // Nothing is printed.
  //  1) startup,       // Startup and endup messages are printed...
  //  2) errors,        // ...and errors...
  //  3) warnings,      // ...and warnings...
  //  4) confirmations, // ...and confirming messages...
  //  5) parameters,    // ...and parameters of scenes and views...
  //  6) all            // ...and everything available.

  G4VisManager* visManager = new G4VisExecutive;
  //visManager -> SetVerboseLevel (1);
  visManager ->Initialize();
#endif


  //get the pointer to the User Interface manager 
  G4UImanager * UI = G4UImanager::GetUIpointer();  

 if (session)   // Define UI session for interactive mode.
    {
      // G4UIterminal is a (dumb) terminal.
      //UI->ApplyCommand("/control/execute myVis.mac");

#if defined(G4UI_USE_XM) || defined(G4UI_USE_WIN32) || defined(G4UI_USE_QT)
      // Customize the G4UIXm,Win32 menubar with a macro file :
      UI->ApplyCommand("/control/execute gui.mac");
#endif

      session->SessionStart();
      delete session;
    }
  else           // Batch mode
    {
#ifdef G4VIS_USE
      visManager->SetVerboseLevel("quiet");
#endif
      G4String command = "/control/execute ";
      G4String fileName = argv[1];
      UI->ApplyCommand(command+fileName);
    }



#ifdef G4VIS_USE
  delete visManager;
#endif

  delete runManager;
  delete myEPEvent;

  return 0;
}