doiPET.cc File Reference

#include "doiPETDetectorConstruction.hh"
#include "doiPETPhysicsList.hh"
#include "doiPETAnalysis.hh"
#include "doiPETActionInitialization.hh"
#include "Randomize.hh"
#include "G4UImanager.hh"
#include "G4VisExecutive.hh"
#include "G4UIExecutive.hh"
#include "G4SystemOfUnits.hh"
#include "G4RunManager.hh"

Go to the source code of this file.

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

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 61 of file doiPET.cc.

{
        // Choose the Random engine
        G4Random::setTheEngine(new CLHEP::RanecuEngine);

#ifdef G4MULTITHREADED

        G4MTRunManager* runManager = new G4MTRunManager;
        runManager->SetNumberOfThreads(4); // Is equal to 2 by default
#else

        G4RunManager* runManager = new G4RunManager;
#endif

        runManager->SetUserInitialization(new doiPETDetectorConstruction);

        runManager->SetUserInitialization(new doiPETPhysicsList);

        // Set user action initialization
        runManager->SetUserInitialization(new doiPETActionInitialization());

        //Initialize analysis
        doiPETAnalysis* ptrAnalysis = doiPETAnalysis::GetInstance();


        G4double act  = 1000000 * becquerel;//Activity is set via run.mac file
        ptrAnalysis->SetActivity(act);

        G4double halfLife = 109.771 * 60 * s; //Halflife of F-18 as a default
        ptrAnalysis -> SetIsotopeHalfLife(halfLife);

        //Blurring specification of the scanner. see inputParameter.txt
        ptrAnalysis -> BlurringParameters();

        //Open file to write the output of the simulation
        ptrAnalysis->Open("result"); //file extention is affixed based on the type of the output (.root for root or .data for ascii)


        //
        ptrAnalysis -> PMTPosition();
        //Read reflector pattern from the inputParameter.txt file
        ptrAnalysis->ReadReflectorPattern();

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



        // Process macro or start UI session
        if (argc!=1)   // batch mode  
        {
                G4String command = "/control/execute ";
                G4String fileName = argv[1];
                UI->ApplyCommand(command+fileName);
        }

        else           //define visualization and UI terminal for interactive mode
        { 
                G4VisManager* visManager = new G4VisExecutive;
                visManager->Initialize();

                G4UIExecutive * ui = new G4UIExecutive(argc,argv);      
                ui->SessionStart();
                delete ui;

                //
                delete visManager;  
        }

        //close the file
        ptrAnalysis->Close();
        ptrAnalysis->Delete();
        
        delete runManager;
        return 0;

}