QweakSimG4.cc

// geant4 includes
#include "G4RunManager.hh"
#include "G4UImanager.hh"
#include "G4UIterminal.hh"
#include "G4UItcsh.hh"

#include "G4Version.hh"
#if G4VERSION_NUMBER < 1000
#include "G4StepLimiterBuilder.hh"
#else
#include "G4StepLimiterPhysics.hh"
#endif

#include "G4OpticalPhysics.hh"
#include "G4PhysListFactory.hh"
#include "G4VModularPhysicsList.hh"

// user includes
#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 "G4UIExecutive.hh"

#include <time.h>

//#define G4UI_USE_QT 1
#ifdef G4UI_USE_QT
    #include "G4UIQt.hh"
#endif

#ifdef G4UI_USE_XM
    #include "G4UIXm.hh"
#endif

#ifdef G4VIS_USE
    #include "G4VisExecutive.hh"
#endif

#define USE_CUSTOM_NUCLEAR_SCATTERING 1
#if USE_CUSTOM_NUCLEAR_SCATTERING
//#include "physics_lists/constructors/electromagnetic/QweakSimEmStandardPhysics.hh"
#include "physics_lists/constructors/electromagnetic/QweakSimEmLivermorePhysics.hh"
#endif

#include <vector>

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

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

  // Verbose output class
  //G4VSteppingVerbose::SetInstance( new QweakSimSteppingVerbose() ); 
  clock_t tStart=clock();
  // 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);

  /*
    vector to pass asym information around
    0: pp = \Pi (1+A_i)
    1: pm = \Pi (1-A_i)
    2: calculation flag (can be used to stop calculation during event 
       -- not used now; if needed needs to be set in steppingAction)
    3: bitInfo for physProcesses: a*2^0+b*2^1
             a=1 for debugPrint; 
             b=1 for modifyTrajectory; 
        e.g.: 3:modifyTrajectory is on and will print for debugging 
   */
  std::vector<double> asymInfo(4,-2);
  asymInfo[3]=0;//default false for both

  // 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(0,&asymInfo));
#endif
  runManager->SetUserInitialization(physlist);

  // UserAction classes
  runManager->SetUserAction( new QweakSimPrimaryGeneratorAction(myQweakSimUserInformation, myEPEvent) );
  //runManager->SetUserAction( new QweakSimPrimaryGeneratorAction( ) );
  runManager->SetUserAction( new QweakSimSteppingAction(myQweakSimUserInformation, myEPEvent,&asymInfo));
  runManager->SetUserAction( new QweakSimStackingAction() );
  runManager->SetUserAction( new QweakSimTrackingAction(myQweakSimUserInformation) );
  runManager->SetUserAction( new QweakSimEventAction(&asymInfo,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 vis/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;
  G4cout<<" Running time[s]: "<< (double) ((clock() - tStart)/CLOCKS_PER_SEC)<<G4endl;
  return 0;
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......