ml2.cc File Reference

#include "ML2Main.hh"
#include "ML2PrimaryGenerationAction.hh"
#include "ML2WorldConstruction.hh"
#include "ML2PhysicsList.hh"
#include "ML2SteppingAction.hh"
#include "ML2EventAction.hh"
#include "ML2TrackingAction.hh"
#include "G4UImanager.hh"
#include "G4Timer.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 82 of file ml2.cc.

References G4UImanager::ApplyCommand(), SGeneralData::bCompareExp, G4RunManager::BeamOn(), CML2CInputData::bOnlyVisio, CML2WorldConstruction::checkVolumeOverlap(), CML2WorldConstruction::create(), CML2PrimaryGenerationAction::design(), SGeneralData::fileExperimentalData, SGeneralData::fileExperimentalDataOut, SInputData::generalData, CML2AcceleratorConstruction::getAcceleratorIsoCentre(), CML2CInputData::getbOnlyVisio(), CML2WorldConstruction::getCML2AcceleratorConstruction(), CML2PhantomConstruction::GetInstance(), CML2PrimaryGenerationAction::GetInstance(), CML2WorldConstruction::GetInstance(), CML2PrimaryGenerationAction::getNrecycling(), G4UImanager::GetUIpointer(), G4Timer::GetUserElapsed(), G4RunManager::Initialize(), CML2PrimaryGenerationAction::inizialize(), CML2CInputData::inputData, SGeneralData::maxNumberOfEvents, SGeneralData::nBeam, CML2WorldConstruction::newGeometry(), SGeneralData::nMaxLoop, SInputData::primaryParticleData, SGeneralData::saving_in_Selected_Voxels_every_events, SGeneralData::seed, CML2Convergence::setNewGeometry(), CLHEP::HepRandom::setTheEngine(), CLHEP::HepRandom::setTheSeed(), G4RunManager::SetUserAction(), G4RunManager::SetUserInitialization(), G4Timer::Start(), SGeneralData::StartFileInputData, G4Timer::Stop(), and CML2Convergence::stopRun().

{
    // instantiate the world class
    CML2WorldConstruction *myWorld=CML2WorldConstruction::GetInstance(); 

    G4RunManager *runManager=new G4RunManager();
    ML2PhysicsList *physics=new ML2PhysicsList();
    runManager->SetUserInitialization(physics);
    
    // build the primary generator
    CML2PrimaryGenerationAction *gun;
    gun = CML2PrimaryGenerationAction::GetInstance();

    // build the main messenger class for the input data
    CML2CInputData *myInputData;
    myInputData=new CML2CInputData();

    // initialize the primary generator variables
    gun->inizialize(&myInputData->inputData.primaryParticleData);

    // according to the number of the launching line
    if (argc==1)
    {
        myInputData->inputData.generalData.seed=1;
        myInputData->inputData.generalData.StartFileInputData="ml2.mac";
    }
    if (argc==2)
    {
        myInputData->inputData.generalData.seed=1;
        myInputData->inputData.generalData.StartFileInputData=(G4String)argv[1];
    }
    if (argc==3)
    {
        sscanf(argv[2],"%d", &myInputData->inputData.generalData.seed);
        myInputData->inputData.generalData.StartFileInputData=(G4String)argv[1];
    }

    // read the main mac file and execute the commands
    G4UImanager* UImanager = G4UImanager::GetUIpointer();
    G4String command = "/control/execute ";
    UImanager->ApplyCommand(command+myInputData->inputData.generalData.StartFileInputData); 


    // set and initialize the random generator
    CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine);
    CLHEP::HepRandom :: setTheSeed(myInputData->inputData.generalData.seed);

    // create the world class
    if (!myWorld->create(&myInputData->inputData, myInputData->getbOnlyVisio()))
    {
        return 1; // if it fails to create the world
    }
    else
    {
        // initialize the primary generator according to the choosen particle source
        gun->design(myWorld->getCML2AcceleratorConstruction()->getAcceleratorIsoCentre());
    }

    // instantiate the convergence control class and assign it to myStepAction
    CML2Convergence *convergence=new CML2Convergence(myInputData->inputData.generalData.seed, myInputData->inputData.generalData.saving_in_Selected_Voxels_every_events, myInputData->inputData.generalData.fileExperimentalData, myInputData->inputData.generalData.fileExperimentalDataOut, myInputData->inputData.generalData.bCompareExp, myInputData->inputData.generalData.maxNumberOfEvents, gun->getNrecycling(), myInputData->inputData.generalData.nMaxLoop);


    // build the ML2RunAction to assign the single phantom name at each run 
    CML2RunAction *myRunAction=new CML2RunAction(convergence, myInputData->inputData.generalData.nBeam, myInputData->bOnlyVisio);
    
    CML2SteppingAction *myStepAction=new CML2SteppingAction(convergence);
    CML2EventAction *ML2EventAction = new CML2EventAction();

    runManager->SetUserInitialization(myWorld);
    runManager->SetUserAction(myRunAction);
    runManager->SetUserAction(gun);
    runManager->SetUserAction(myStepAction);
    runManager->SetUserAction(ML2EventAction);
    runManager->SetUserAction(new CML2TrackingAction);
    runManager->Initialize();

    // performances info 
    int nLoop=0;
    G4Timer MyFullTime;
    G4double loopElapsedTime;
    G4bool bStopRun=false;
    G4bool bNewGeometry=true;
    if (myInputData->bOnlyVisio)
    {
#ifdef G4VIS_USE
        // visualization
        myWorld->newGeometry();
        convergence->setNewGeometry();
        visio(argc, argv);
#endif
    }
    else
    {
        MyFullTime.Start();
        // compute 
        while (bNewGeometry)
        {
            bNewGeometry=myWorld->newGeometry();
            convergence->setNewGeometry();
            if (bNewGeometry)
            {
                if (CML2PhantomConstruction::GetInstance()->getPhantomName()!="Dicom1")
                {CML2WorldConstruction::GetInstance()->checkVolumeOverlap();}
                std::cout<<"################ START NEW GEOMETRY ########################"<<'\n';
                myRunAction->setActualLoop(nLoop);
                while (!bStopRun)
                {
                    runManager->BeamOn(myInputData->inputData.generalData.nBeam);
                    // check if the run has to be reapeted
                    bStopRun=convergence->stopRun(); 
                }
                nLoop=0;
            std::cout<<"################ END NEW GEOMETRY ########################"<<'\n';
            }
            bStopRun=false;
        }
        MyFullTime.Stop();
        loopElapsedTime=MyFullTime.GetUserElapsed();
        std::cout << "loop elapsed time [s] : "<< loopElapsedTime << '\n';
        std::cout <<'\n';
    }
        delete runManager;
        return 0;
}