exampleB03.cc File Reference

Main program of the biasing/B03 example. More...

#include <iostream>
#include "G4Types.hh"
#include "G4RunManagerFactory.hh"
#include "G4VPhysicalVolume.hh"
#include "G4UImanager.hh"
#include "G4SystemOfUnits.hh"
#include "B03DetectorConstruction.hh"
#include "B03PhysicsList.hh"
#include "B03ActionInitialization.hh"
#include "B03ImportanceDetectorConstruction.hh"
#include "G4IStore.hh"

Go to the source code of this file.

Functions
int	main (int, char **)

Detailed Description

Main program of the biasing/B03 example.

Definition in file exampleB03.cc.

Function Documentation

main()

int main	(	int	,
		char **
	)

Definition at line 83 of file exampleB03.cc.

{
  G4int numberOfEvents = 100;
 
  G4long myseed = 345354;
 
  auto* runManager = G4RunManagerFactory::CreateRunManager();
  runManager->SetNumberOfThreads(2);
 
  G4Random::setTheSeed(myseed);
 
  // create the detector      ---------------------------
  B03DetectorConstruction *detector = new B03DetectorConstruction();
  runManager->SetUserInitialization(detector);
  //  ---------------------------------------------------
 
  // create a parallel detector
  G4String parallelName("ParallelBiasingWorld");
  B03ImportanceDetectorConstruction *pdet =
    new B03ImportanceDetectorConstruction(parallelName);
  detector->RegisterParallelWorld(pdet);
 
  // G4GeometrySampler pgs(pdet->GetWorldVolume(),"neutron");
  // B03PhysicsList* physlist = new B03PhysicsList;
  // name of first parallel world:
  B03PhysicsList* physlist = new B03PhysicsList(parallelName);
  // push parallel world to store in case of multiple worlds:
  physlist->AddParallelWorldName(parallelName);
  // physlist->AddParallelWorldName(parallelName);
  //physlist->AddParallelWorldName(sparallelName);
  // physlist->AddBiasing(&pgs,parallelName);
 
  runManager->SetUserInitialization(physlist);
 
 // Set user action classes through Worker Initialization
 //
  B03ActionInitialization* actions = new B03ActionInitialization;
  runManager->SetUserInitialization(actions);
 
//   runManager->SetUserAction(new B03PrimaryGeneratorAction);
//   //  runManager->SetUserAction(new B03PrimaryGeneratorAction(ifElectron));
//   runManager->SetUserAction(new B03RunAction);
 
  runManager->Initialize();
 
  G4VPhysicalVolume& aghostWorld = pdet->GetWorldVolumeAddress();
  G4cout << " ghost world: " << pdet->GetName() << G4endl;
 
  // create an importance
  G4IStore *aIstore = G4IStore::GetInstance(pdet->GetName());
 
  // create a geometry cell for the world volume replpicanumber is 0!
  G4GeometryCell gWorldVolumeCell(aghostWorld, 0);
  // set world volume importance to 1
  aIstore->AddImportanceGeometryCell(1, gWorldVolumeCell);
 
  // set importance values and create scorers
  G4int cell(1);
  for (cell=1; cell<=18; cell++) {
    G4GeometryCell gCell = pdet->GetGeometryCell(cell);
    G4cout << " adding cell: " << cell
           << " replica: " << gCell.GetReplicaNumber()
           << " name: " << gCell.GetPhysicalVolume().GetName() << G4endl;
    G4double imp = std::pow(2.0,cell-1);
    //x    aIstore.AddImportanceGeometryCell(imp, gCell);
    aIstore->AddImportanceGeometryCell(imp, gCell.GetPhysicalVolume(), cell);
  }
 
  // creating the geometry cell and add both to the store
  //  G4GeometryCell gCell = pdet->GetGeometryCell(18);
 
  // create importance geometry cell pair for the "rest"cell
  // with the same importance as the last concrete cell
  G4GeometryCell gCell = pdet->GetGeometryCell(19);
  //  G4double imp = std::pow(2.0,18);
  G4double imp = std::pow(2.0,17);
  aIstore->AddImportanceGeometryCell(imp, gCell.GetPhysicalVolume(), 19);
 
  //temporary fix before runManager->BeamOn works...
  G4UImanager* UImanager = G4UImanager::GetUIpointer();
  G4String command1 = "/control/cout/setCoutFile fileName";
  UImanager->ApplyCommand(command1);
 
  G4String command2 = "/run/beamOn "
                    + G4UIcommand::ConvertToString(numberOfEvents);;
  UImanager->ApplyCommand(command2);
 
  //  runManager->BeamOn(numberOfEvents);
 
  // pgs.ClearSampling();
 
  delete runManager;
 
  return 0;
}