Implementation of the TVirtualMCApplication. More...

#include <Ex03cMCApplication.h>

Inheritance diagram for Ex03cMCApplication:

Public Member Functions
	Ex03cMCApplication (const char name, const char title, Bool_t isMulti=kFALSE, Bool_t splitSimulation=kFALSE)

	Ex03cMCApplication ()

virtual	~Ex03cMCApplication ()

void	InitMC (const char *setup)

void	InitMC (std::initializer_list< const char * > setupMacros)

void	InitMC ()

void	RunMC (Int_t nofEvents)

void	FinishRun ()

void	ReadEvent (Int_t i)

virtual TVirtualMCApplication *	CloneForWorker () const

virtual void	InitOnWorker ()

virtual void	FinishRunOnWorker ()

virtual void	ConstructGeometry ()

virtual void	InitGeometry ()

virtual void	AddParticles ()

virtual void	AddIons ()

virtual void	GeneratePrimaries ()

virtual void	BeginEvent ()

virtual void	BeginPrimary ()

virtual void	PreTrack ()

virtual void	Stepping ()

virtual void	PostTrack ()

virtual void	FinishPrimary ()

virtual void	FinishEvent ()

void	SetPrintModulo (Int_t value)

void	SetVerboseLevel (Int_t verboseLevel)

void	SetControls (Bool_t isConstrols)

void	SetField (Double_t bz)

void	SetDebug (Int_t debug)
	Set debug option for multiple run.

Ex03cDetectorConstruction *	GetDetectorConstruction () const

Ex03cCalorimeterSD *	GetCalorimeterSD () const

Ex03PrimaryGenerator *	GetPrimaryGenerator () const

void	SetOldGeometry (Bool_t oldGeometry=kTRUE)

Private Member Functions
	Ex03cMCApplication (const Ex03cMCApplication &origin)

void	RegisterStack () const

Private Attributes
TMCRootManager *	fRootManager
	Root manager.

Int_t	fPrintModulo
	The event modulus number to be printed.

Int_t	fEventNo
	Event counter.

TMCVerbose	fVerbose
	VMC verbose helper.

Ex03cMCStack *	fStack
	VMC stack.

Ex03cDetectorConstruction *	fDetConstruction
	Dector construction.

Ex03cCalorimeterSD *	fCalorimeterSD
	Calorimeter SD.

Ex03PrimaryGenerator *	fPrimaryGenerator
	Primary generator.

TGeoUniformMagField *	fMagField
	Magnetic field.

Bool_t	fOldGeometry
	Option for geometry definition.

Bool_t	fIsControls
	Option to activate special controls.

Bool_t	fIsMaster
	If is on master thread.

Bool_t	fIsMultiRun
	Flag if having multiple engines.

Bool_t	fSplitSimulation
	Split geometry given user criteria.

Int_t	fG3Id
	engine ID of Geant3

Int_t	fG4Id
	engine ID of Geant4

Int_t	fDebug
	debug option for multiple run

Detailed Description

Implementation of the TVirtualMCApplication.

A variant of the Ex03MCApplication class updated for multiple engine runs.

Date: 21/08/2019

Author: Benedikt Volkel, CERN

Definition at line 45 of file Ex03cMCApplication.h.

Constructor & Destructor Documentation

◆ Ex03cMCApplication() [1/3]

Ex03cMCApplication::Ex03cMCApplication	(	const char *	name,
		const char *	title,
		Bool_t	isMulti = kFALSE,
		Bool_t	splitSimulation = kFALSE )

Standard constructor

Parameters

name	The MC application name
title	The MC application description

Definition at line 48 of file Ex03cMCApplication.cxx.

  : TVirtualMCApplication(name, title),
    fRootManager(0),
    fPrintModulo(1),
    fEventNo(0),
    fVerbose(0),
    fStack(0),
    fDetConstruction(0),
    fCalorimeterSD(0),
    fPrimaryGenerator(0),
    fMagField(0),
    fOldGeometry(kFALSE),
    fIsControls(kFALSE),
    fIsMaster(kTRUE),
    fIsMultiRun(isMulti),
    fSplitSimulation(splitSimulation),
    fG3Id(-1),
    fG4Id(-1),
    fDebug(0)
{
  /// Standard constructor
  /// \param name   The MC application name
  /// \param title  The MC application description
 
  if (splitSimulation && !isMulti) {
    Fatal("Ex03cMCApplication",
      "Cannot split simulation between engines without \"isMulti\" being "
      "switched on");
  }
 
  cout << "--------------------------------------------------------------"
       << endl;
  cout << " VMC Example E03c: new version with multiple engines" << endl;
  cout << "--------------------------------------------------------------"
       << endl;
 
  // Create a user stack
  fStack = new Ex03cMCStack(1000);
 
  // Create detector construction
  fDetConstruction = new Ex03cDetectorConstruction();
 
  // Create a calorimeter SD
  fCalorimeterSD = new Ex03cCalorimeterSD("Calorimeter", fDetConstruction);
 
  // Create a primary generator
  fPrimaryGenerator = new Ex03PrimaryGenerator(fStack);
 
  // Constant magnetic field (in kiloGauss)
  fMagField = new TGeoUniformMagField();
 
  if (isMulti) {
    RequestMCManager();
    fMCManager->SetUserStack(fStack);
  }
}

◆ Ex03cMCApplication() [2/3]

Ex03cMCApplication::Ex03cMCApplication ( )

Default constructor

Definition at line 147 of file Ex03cMCApplication.cxx.

  : TVirtualMCApplication(),
    fRootManager(0),
    fPrintModulo(1),
    fEventNo(0),
    fStack(0),
    fDetConstruction(0),
    fCalorimeterSD(0),
    fPrimaryGenerator(0),
    fMagField(0),
    fOldGeometry(kFALSE),
    fIsControls(kFALSE),
    fIsMaster(kTRUE),
    fIsMultiRun(kFALSE),
    fSplitSimulation(kFALSE),
    fG3Id(-1),
    fG4Id(-1),
    fDebug(0)
{
  /// Default constructor
}

◆ ~Ex03cMCApplication()

Ex03cMCApplication::~Ex03cMCApplication ( )

virtual

Destructor

Definition at line 170 of file Ex03cMCApplication.cxx.

{
  /// Destructor
 
  // cout << "Ex03cMCApplication::~Ex03cMCApplication " << this << endl;
 
  delete fRootManager;
  delete fStack;
  if (fIsMaster) delete fDetConstruction;
  delete fCalorimeterSD;
  delete fPrimaryGenerator;
  delete fMagField;
  if (!fIsMultiRun) {
    delete fMC;
  }
 
  // cout << "Done Ex03cMCApplication::~Ex03cMCApplication " << this << endl;
}

◆ Ex03cMCApplication() [3/3]

Ex03cMCApplication::Ex03cMCApplication ( const Ex03cMCApplication & origin )

private

Copy constructor for cloning application on workers (in multithreading mode)

Parameters

origin The source MC application

Definition at line 107 of file Ex03cMCApplication.cxx.

  : TVirtualMCApplication(origin.GetName(), origin.GetTitle()),
    fRootManager(0),
    fPrintModulo(origin.fPrintModulo),
    fEventNo(0),
    fVerbose(origin.fVerbose),
    fStack(0),
    fDetConstruction(origin.fDetConstruction),
    fCalorimeterSD(0),
    fPrimaryGenerator(0),
    fMagField(0),
    fOldGeometry(origin.fOldGeometry),
    fIsControls(origin.fIsControls),
    fIsMaster(kFALSE),
    fIsMultiRun(origin.fIsMultiRun),
    fSplitSimulation(origin.fSplitSimulation),
    fG3Id(origin.fG3Id),
    fG4Id(origin.fG4Id),
    fDebug(origin.fDebug)
{
  /// Copy constructor for cloning application on workers (in multithreading
  /// mode) \param origin   The source MC application
 
  // Create new user stack
  fStack = new Ex03cMCStack(1000);
 
  // Create a calorimeter SD
  fCalorimeterSD =
    new Ex03cCalorimeterSD(*(origin.fCalorimeterSD), fDetConstruction);
 
  // Create a primary generator
  fPrimaryGenerator =
    new Ex03PrimaryGenerator(*(origin.fPrimaryGenerator), fStack);
 
  // Constant magnetic field (in kiloGauss)
  fMagField = new TGeoUniformMagField(origin.fMagField->GetFieldValue()[0],
    origin.fMagField->GetFieldValue()[1], origin.fMagField->GetFieldValue()[2]);
}

Member Function Documentation

◆ InitMC() [1/3]

void Ex03cMCApplication::InitMC ( const char * setup )

Initialize MC. The selection (and creation) of the concrete MC is done in the macro. If no macro is provided, the MC must be created before calling this function.

Parameters

setup The name of the configuration macro.

Definition at line 209 of file Ex03cMCApplication.cxx.

{
  /// Initialize MC.
  /// The selection (and creation) of the concrete MC is done in the macro.
  /// If no macro is provided, the MC must be created before calling this
  /// function. \param setup The name of the configuration macro.
 
  fVerbose.InitMC();
 
  if (TString(setup) != "") {
    gROOT->LoadMacro(setup);
    gInterpreter->ProcessLine("Config()");
    if (!fMC) {
      Fatal(
        "InitMC", "Processing Config() has failed. (No MC is instantiated.)");
    }
  }
 
// MT support available from root v 5.34/18
#if ROOT_VERSION_CODE >= 336402
  // Create Root manager
  if (!fMC->IsMT()) {
    fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite);
    // fRootManager->SetDebug(true);
  }
#else
  // Create Root manager
  fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite);
// fRootManager->SetDebug(true);
#endif
 
  fMC->SetStack(fStack);
  fMC->SetRootGeometry();
  fMC->SetMagField(fMagField);
  fMC->Init();
  fMC->BuildPhysics();
 
  RegisterStack();
 
  if (fDebug > 0) {
    Info("InitMC", "Single run initialised");
  }
}

◆ InitMC() [2/3]

void Ex03cMCApplication::InitMC ( std::initializer_list< const char * > setupMacros )

Initialize MC. New function for initialization with multiple engines.

Parameters

setupMacros The names of the configuration macros.

Definition at line 254 of file Ex03cMCApplication.cxx.

{
  /// Initialize MC.
  /// New function for initialization with multiple engines.
  /// \param setupMacros The names of the configuration macros.
 
  fVerbose.InitMC();
 
  if (!fIsMultiRun) {
    Fatal("InitMC",
      "Initialisation of multiple engines not supported in single run");
  }
 
  if (setupMacros.size() > 0) {
    for (const char* setup : setupMacros) {
      gROOT->LoadMacro(setup);
      gInterpreter->ProcessLine("Config()");
      if (!fMC) {
        Fatal(
          "InitMC", "Processing Config() has failed. (No MC is instantiated.)");
      }
      // gInterpreter->UnloadAllSharedLibraryMaps();
      gInterpreter->UnloadFile(setup);
    }
  }
 
  InitMC();
}

◆ InitMC() [3/3]

void Ex03cMCApplication::InitMC ( )

Initialize MC. New function for initialization with multiple engines. The MCs must be created before calling this function.

Definition at line 284 of file Ex03cMCApplication.cxx.

{
  /// Initialize MC.
  /// New function for initialization with multiple engines.
  /// The MCs must be created before calling this function.
 
  if (!fIsMultiRun) {
    Fatal("InitMC",
      "Initialisation of multiple engines not supported in single run");
  }
 
// MT support available from root v 5.34/18
#if ROOT_VERSION_CODE >= 336402
  // Create Root manager
  if (!fMC->IsMT()) {
    fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite);
    // fRootManager->SetDebug(true);
  }
#else
  // Create Root manager
  fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite);
// fRootManager->SetDebug(true);
#endif
 
  fMCManager->Init([this](TVirtualMC* mc) {
    mc->SetRootGeometry();
    mc->SetMagField(this->fMagField);
    mc->Init();
    mc->BuildPhysics();
  });
  RegisterStack();
 
  fG3Id = fMCManager->GetEngineId("TGeant3TGeo");
  fG4Id = fMCManager->GetEngineId("TGeant4");
  if (fDebug > 0) {
    Info("InitMC", "Multi run initialised");
    std::cout << "Engine IDs\n"
              << "TGeant3TGeo: " << fG3Id << "\n"
              << "TGeant4: " << fG4Id << std::endl;
  }
}

◆ RunMC()

void Ex03cMCApplication::RunMC ( Int_t nofEvents )

Run MC.

Parameters

nofEvents Number of events to be processed

Definition at line 327 of file Ex03cMCApplication.cxx.

{
  /// Run MC.
  /// \param nofEvents Number of events to be processed
 
  fVerbose.RunMC(nofEvents);
 
  // Prepare a timer
  TTimer timer;
 
  if (!fIsMultiRun) {
    if (fDebug > 0) {
      Info("RunMC", "Start single run");
      std::cout << "Simulation entirely done with engine " << fMC->GetName()
                << std::endl;
      timer.Start();
    }
    fMC->ProcessRun(nofEvents);
  }
  else {
    if (fDebug > 0) {
      Info("RunMC", "Start multi run");
      if (fSplitSimulation) {
        std::cout << "GAPX simulated with engine "
                  << fMCManager->GetEngine(0)->GetName() << "\n"
                  << "ABSO simulated with engine "
                  << fMCManager->GetEngine(1)->GetName() << std::endl;
      }
      else {
        std::cout << "Simulation entirely done with engine "
                  << fMCManager->GetCurrentEngine()->GetName() << std::endl;
      }
      timer.Start();
    }
    fMCManager->Run(nofEvents);
  }
  if (fDebug > 0) {
    timer.Stop();
    Info("RunMC", "Transport finished.");
    timer.Print();
  }
  FinishRun();
}

◆ FinishRun()

void Ex03cMCApplication::FinishRun ( )

Finish MC run.

Definition at line 372 of file Ex03cMCApplication.cxx.

{
  /// Finish MC run.
 
  fVerbose.FinishRun();
  // cout << "Ex03cMCApplication::FinishRun: " << endl;
  if (fRootManager) {
    fRootManager->WriteAll();
    fRootManager->Close();
  }
}

◆ ReadEvent()

void Ex03cMCApplication::ReadEvent ( Int_t i )

Read i -th event and prints hits.

Parameters

i	The number of event to be read

Definition at line 417 of file Ex03cMCApplication.cxx.

{
  /// Read \em i -th event and prints hits.
  /// \param i The number of event to be read
 
  fCalorimeterSD->Register();
  RegisterStack();
  fRootManager->ReadEvent(i);
}

◆ CloneForWorker()

TVirtualMCApplication * Ex03cMCApplication::CloneForWorker ( ) const

virtual

Definition at line 385 of file Ex03cMCApplication.cxx.

{
  return new Ex03cMCApplication(*this);
}

◆ InitOnWorker()

void Ex03cMCApplication::InitOnWorker ( )

virtual

Definition at line 391 of file Ex03cMCApplication.cxx.

{
  // cout << "Ex03cMCApplication::InitForWorker " << this << endl;
 
  // Create Root manager
  fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite);
  // fRootManager->SetDebug(true);
 
  // Set data to MC
  fMC->SetStack(fStack);
  fMC->SetMagField(fMagField);
 
  RegisterStack();
}

◆ FinishRunOnWorker()

void Ex03cMCApplication::FinishRunOnWorker ( )

virtual

Definition at line 407 of file Ex03cMCApplication.cxx.

{
  // cout << "Ex03cMCApplication::FinishWorkerRun: " << endl;
  if (fRootManager) {
    fRootManager->WriteAll();
    fRootManager->Close();
  }
}

◆ ConstructGeometry()

void Ex03cMCApplication::ConstructGeometry ( )

virtual

Construct geometry using detector contruction class. The detector contruction class is using TGeo functions or TVirtualMC functions (if oldGeometry is selected)

Definition at line 428 of file Ex03cMCApplication.cxx.

{
  /// Construct geometry using detector contruction class.
  /// The detector contruction class is using TGeo functions or
  /// TVirtualMC functions (if oldGeometry is selected)
 
  fVerbose.ConstructGeometry();
 
  if (fOldGeometry) {
    Fatal("ConstructGeometry", "Cannot run with old geomtry");
  }
  fDetConstruction->ConstructMaterials();
  fDetConstruction->ConstructGeometry();
  // TGeoManager::Import("geometry.root");
  // fMC->SetRootGeometry();
}

◆ InitGeometry()

void Ex03cMCApplication::InitGeometry ( )

virtual

Initialize geometry

Definition at line 446 of file Ex03cMCApplication.cxx.

{
  /// Initialize geometry
 
  fVerbose.InitGeometry();
 
  fDetConstruction->SetCuts();
 
  if (fIsControls) fDetConstruction->SetControls();
 
  fCalorimeterSD->Initialize();
}

◆ AddParticles()

void Ex03cMCApplication::AddParticles ( )

virtual

Example of user defined particle with user defined decay mode

Definition at line 460 of file Ex03cMCApplication.cxx.

{
  /// Example of user defined particle with user defined decay mode
 
  fVerbose.AddParticles();
 
  // Define the 2 body  phase space decay  for He5
  Int_t mode[6][3];
  Float_t bratio[6];
 
  for (Int_t kz = 0; kz < 6; kz++) {
    bratio[kz] = 0.;
    mode[kz][0] = 0;
    mode[kz][1] = 0;
    mode[kz][2] = 0;
  }
  bratio[0] = 100.;
  mode[0][0] = kNeutron;   // neutron (2112)
  mode[0][1] = 1000020040; // alpha
 
  // Overwrite a decay mode already defined in MCs
  // Kaon Short: 310 normally decays in two modes
  // pi+, pi-  68.61 %
  // pi0, pi0  31.39 %
  // and we force only the mode pi0, pi0
 
  Int_t mode2[6][3];
  Float_t bratio2[6];
 
  for (Int_t kz = 0; kz < 6; kz++) {
    bratio2[kz] = 0.;
    mode2[kz][0] = 0;
    mode2[kz][1] = 0;
    mode2[kz][2] = 0;
  }
  bratio2[0] = 100.;
  mode2[0][0] = kPi0; // pi0 (111)
  mode2[0][1] = kPi0; // pi0 (111)
 
  // Define particle
  fMC->DefineParticle(1000020050, "He5", kPTHadron, 5.03427, 2.0, 0.002, "Ion",
    0.0, 0, 1, 0, 0, 0, 0, 0, 5, kFALSE);
  fMC->SetDecayMode(1000020050, bratio, mode);
  fMC->SetDecayMode(kK0Short, bratio2, mode2);
}

◆ AddIons()

void Ex03cMCApplication::AddIons ( )

virtual

Example of user defined ion

Definition at line 507 of file Ex03cMCApplication.cxx.

{
  /// Example of user defined ion
 
  fVerbose.AddIons();
  fMC->DefineIon("MyIon", 34, 70, 12, 0.);
}

◆ GeneratePrimaries()

void Ex03cMCApplication::GeneratePrimaries ( )

virtual

Fill the user stack (derived from TVirtualMCStack) with primary particles.

Definition at line 516 of file Ex03cMCApplication.cxx.

{
  /// Fill the user stack (derived from TVirtualMCStack) with primary particles.
 
  fVerbose.GeneratePrimaries();
 
  TVector3 origin(fDetConstruction->GetWorldSizeX(),
    fDetConstruction->GetCalorSizeYZ(), fDetConstruction->GetCalorSizeYZ());
 
  fPrimaryGenerator->GeneratePrimaries(origin);
}

◆ BeginEvent()

void Ex03cMCApplication::BeginEvent ( )

virtual

User actions at beginning of event

Definition at line 529 of file Ex03cMCApplication.cxx.

{
  /// User actions at beginning of event
 
  fVerbose.BeginEvent();
 
  // Clear TGeo tracks (if filled)
  if (!fIsMultiRun && TString(fMC->GetName()) == "TGeant3TGeo" &&
      gGeoManager->GetListOfTracks() && gGeoManager->GetTrack(0) &&
      ((TVirtualGeoTrack*)gGeoManager->GetTrack(0))->HasPoints()) {
 
    gGeoManager->ClearTracks();
    // if (gPad) gPad->Clear();
  }
 
  fEventNo++;
  if (fEventNo % fPrintModulo == 0) {
    cout << "\n---> Begin of event: " << fEventNo << endl;
    // ??? How to do this in VMC
    // HepRandom::showEngineStatus();
  }
}

◆ BeginPrimary()

void Ex03cMCApplication::BeginPrimary ( )

virtual

User actions at beginning of a primary track. If test for user defined decay is activated, the primary track ID is printed on the screen.

Definition at line 553 of file Ex03cMCApplication.cxx.

{
  /// User actions at beginning of a primary track.
  /// If test for user defined decay is activated,
  /// the primary track ID is printed on the screen.
 
  fVerbose.BeginPrimary();
 
  if (fPrimaryGenerator->GetUserDecay()) {
    cout << "   Primary track ID = " << fStack->GetCurrentTrackNumber() << endl;
  }
}

◆ PreTrack()

void Ex03cMCApplication::PreTrack ( )

virtual

User actions at beginning of each track If test for user defined decay is activated, the decay products of the primary track (K0Short) are printed on the screen.

Definition at line 567 of file Ex03cMCApplication.cxx.

{
  /// User actions at beginning of each track
  /// If test for user defined decay is activated,
  /// the decay products of the primary track (K0Short)
  /// are printed on the screen.
 
  fVerbose.PreTrack();
 
  // print info about K0Short decay products
  if (fPrimaryGenerator->GetUserDecay()) {
    Int_t parentID = fStack->GetCurrentParentTrackNumber();
 
    if (parentID >= 0 &&
        fStack->GetParticle(parentID)->GetPdgCode() == kK0Short &&
        fStack->GetCurrentTrack()->GetUniqueID() == kPDecay) {
      // The production process is saved as TParticle unique ID
      // via Ex03cMCStack
 
      cout << "      Current track " << fStack->GetCurrentTrack()->GetName()
           << "  is a decay product of Parent ID = "
           << fStack->GetCurrentParentTrackNumber() << endl;
    }
  }
}

◆ Stepping()

void Ex03cMCApplication::Stepping ( )

virtual

User actions at each step

Definition at line 594 of file Ex03cMCApplication.cxx.

{
  /// User actions at each step
 
  fVerbose.Stepping();
 
  fCalorimeterSD->ProcessHits();
 
  TLorentzVector pos;
  TLorentzVector mom;
 
  fMC->TrackPosition(pos);
  fMC->TrackMomentum(mom);
 
  if (fDebug > 1) {
    std::cout << "Current engine " << fMC->GetName() << "\n"
              << "Track ID=" << fStack->GetCurrentTrackNumber() << "\n"
              << "Momentum\n"
              << "E=" << mom.T() << ", Px=" << mom.X() << ", Py=" << mom.Y()
              << ", Pz=" << mom.Z() << std::endl;
  }
 
  // Now transfer track
  if (fSplitSimulation) {
    Int_t targetId = -1;
    if (fMC->GetId() == 0 && strcmp(fMC->CurrentVolName(), "ABSO") == 0) {
      targetId = 1;
    }
    else if (fMC->GetId() == 1 && strcmp(fMC->CurrentVolName(), "GAPX") == 0) {
      targetId = 0;
    }
    if (targetId > -1) {
      if (fDebug > 1) {
        Info("Stepping", "Transfer track");
      }
      fMCManager->TransferTrack(targetId);
    }
  }
}

◆ PostTrack()

void Ex03cMCApplication::PostTrack ( )

virtual

User actions after finishing of each track

Definition at line 635 of file Ex03cMCApplication.cxx.

{
  /// User actions after finishing of each track
  fVerbose.PostTrack();
}

◆ FinishPrimary()

void Ex03cMCApplication::FinishPrimary ( )

virtual

User actions after finishing of a primary track

Definition at line 642 of file Ex03cMCApplication.cxx.

{
  /// User actions after finishing of a primary track
 
  fVerbose.FinishPrimary();
 
  if (fPrimaryGenerator->GetUserDecay()) {
    cout << endl;
  }
}

◆ FinishEvent()

void Ex03cMCApplication::FinishEvent ( )

virtual

User actions after finishing of an event

Definition at line 654 of file Ex03cMCApplication.cxx.

{
  /// User actions after finishing of an event
 
  fVerbose.FinishEvent();
 
  // Geant3 + TGeo
  // (use TGeo functions for visualization)
  if (!fIsMultiRun && TString(fMC->GetName()) == "TGeant3TGeo") {
 
    // Draw volume
    gGeoManager->SetVisOption(0);
    gGeoManager->SetTopVisible();
    gGeoManager->GetTopVolume()->Draw();
 
    // Draw tracks (if filled)
    // Available when this feature is activated via
    // fMC->SetCollectTracks(kTRUE);
    if (gGeoManager->GetListOfTracks() && gGeoManager->GetTrack(0) &&
        ((TVirtualGeoTrack*)gGeoManager->GetTrack(0))->HasPoints()) {
 
      gGeoManager->DrawTracks("/*"); // this means all tracks
    }
  }
 
  fRootManager->Fill();
 
  if (fEventNo % fPrintModulo == 0) fCalorimeterSD->PrintTotal();
 
  fCalorimeterSD->EndOfEvent();
 
  fStack->Reset();
}

◆ SetPrintModulo()

void Ex03cMCApplication::SetPrintModulo ( Int_t value )

inline

Set the event modulus number to be printed

Parameters

value The new event modulus number value

Definition at line 127 of file Ex03cMCApplication.h.

{
  fPrintModulo = value;
}

◆ SetVerboseLevel()

void Ex03cMCApplication::SetVerboseLevel ( Int_t verboseLevel )

inline

Set verbosity

Parameters

verboseLevel The new verbose level value

Definition at line 134 of file Ex03cMCApplication.h.

{
  fVerbose.SetLevel(verboseLevel);
}

◆ SetControls()

void Ex03cMCApplication::SetControls ( Bool_t isControls )

inline

Switch on/off special process controls

Parameters

isControls If true, special process controls setting is activated

Definition at line 177 of file Ex03cMCApplication.h.

{
  fIsControls = isControls;
}

◆ SetField()

void Ex03cMCApplication::SetField ( Double_t bz )

inline

Set magnetic field

Parameters

bz	The new field value in z

Definition at line 141 of file Ex03cMCApplication.h.

{
  fMagField->SetFieldValue(0., 0., bz);
}

◆ SetDebug()

void Ex03cMCApplication::SetDebug ( Int_t debug )

inline

Set debug option for multiple run.

Definition at line 147 of file Ex03cMCApplication.h.

147{ fDebug = debug; }

◆ GetDetectorConstruction()

Ex03cDetectorConstruction * Ex03cMCApplication::GetDetectorConstruction ( ) const

inline

Returns: The detector construction

Definition at line 151 of file Ex03cMCApplication.h.

{
  return fDetConstruction;
}

◆ GetCalorimeterSD()

Ex03cCalorimeterSD * Ex03cMCApplication::GetCalorimeterSD ( ) const

inline

Returns: The calorimeter sensitive detector

Definition at line 157 of file Ex03cMCApplication.h.

{
  return fCalorimeterSD;
}

◆ GetPrimaryGenerator()

Ex03PrimaryGenerator * Ex03cMCApplication::GetPrimaryGenerator ( ) const

inline

Returns: The primary generator

Definition at line 163 of file Ex03cMCApplication.h.

{
  return fPrimaryGenerator;
}

◆ SetOldGeometry()

void Ex03cMCApplication::SetOldGeometry ( Bool_t oldGeometry = kTRUE )

inline

Switch on/off the old geometry definition (via VMC functions)

Parameters

oldGeometry If true, geometry definition via VMC functions

Definition at line 170 of file Ex03cMCApplication.h.

{
  fOldGeometry = oldGeometry;
}

◆ RegisterStack()

void Ex03cMCApplication::RegisterStack ( ) const

private

Register stack in the Root manager.

Definition at line 194 of file Ex03cMCApplication.cxx.

{
  /// Register stack in the Root manager.
 
  if (fRootManager) {
    // cout << "Ex03cMCApplication::RegisterStack: " << endl;
    fRootManager->Register("stack", "Ex03cMCStack", &fStack);
  }
}

Member Data Documentation

◆ fRootManager

TMCRootManager* Ex03cMCApplication::fRootManager

mutableprivate

Root manager.

Definition at line 99 of file Ex03cMCApplication.h.

◆ fPrintModulo

Int_t Ex03cMCApplication::fPrintModulo

private

The event modulus number to be printed.

Definition at line 100 of file Ex03cMCApplication.h.

◆ fEventNo

Int_t Ex03cMCApplication::fEventNo

private

Event counter.

Definition at line 101 of file Ex03cMCApplication.h.

◆ fVerbose

TMCVerbose Ex03cMCApplication::fVerbose

private

VMC verbose helper.

Definition at line 102 of file Ex03cMCApplication.h.

◆ fStack

Ex03cMCStack* Ex03cMCApplication::fStack

private

VMC stack.

Definition at line 103 of file Ex03cMCApplication.h.

◆ fDetConstruction

Ex03cDetectorConstruction* Ex03cMCApplication::fDetConstruction

private

Dector construction.

Definition at line 104 of file Ex03cMCApplication.h.

◆ fCalorimeterSD

Ex03cCalorimeterSD* Ex03cMCApplication::fCalorimeterSD

private

Calorimeter SD.

Definition at line 105 of file Ex03cMCApplication.h.

◆ fPrimaryGenerator

Ex03PrimaryGenerator* Ex03cMCApplication::fPrimaryGenerator

private

Primary generator.

Definition at line 106 of file Ex03cMCApplication.h.

◆ fMagField

TGeoUniformMagField* Ex03cMCApplication::fMagField

private

Magnetic field.

Definition at line 107 of file Ex03cMCApplication.h.

◆ fOldGeometry

Bool_t Ex03cMCApplication::fOldGeometry

private

Option for geometry definition.

Definition at line 108 of file Ex03cMCApplication.h.

◆ fIsControls

Bool_t Ex03cMCApplication::fIsControls

private

Option to activate special controls.

Definition at line 109 of file Ex03cMCApplication.h.

◆ fIsMaster

Bool_t Ex03cMCApplication::fIsMaster

private

If is on master thread.

Definition at line 110 of file Ex03cMCApplication.h.

◆ fIsMultiRun

Bool_t Ex03cMCApplication::fIsMultiRun

private

Flag if having multiple engines.

Definition at line 111 of file Ex03cMCApplication.h.

◆ fSplitSimulation

Bool_t Ex03cMCApplication::fSplitSimulation

private

Split geometry given user criteria.

Definition at line 112 of file Ex03cMCApplication.h.

◆ fG3Id

Int_t Ex03cMCApplication::fG3Id

private

engine ID of Geant3

Definition at line 113 of file Ex03cMCApplication.h.

◆ fG4Id

Int_t Ex03cMCApplication::fG4Id

private

engine ID of Geant4

Definition at line 114 of file Ex03cMCApplication.h.

◆ fDebug

Int_t Ex03cMCApplication::fDebug

private

debug option for multiple run

Definition at line 115 of file Ex03cMCApplication.h.

The documentation for this class was generated from the following files:

examples/E03/E03c/include/Ex03cMCApplication.h
examples/E03/E03c/src/Ex03cMCApplication.cxx

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Ex03cMCApplication() [1/3]

◆ Ex03cMCApplication() [2/3]

◆ ~Ex03cMCApplication()

◆ Ex03cMCApplication() [3/3]

Member Function Documentation

◆ InitMC() [1/3]

◆ InitMC() [2/3]

◆ InitMC() [3/3]

◆ RunMC()

◆ FinishRun()

◆ ReadEvent()

◆ CloneForWorker()

◆ InitOnWorker()

◆ FinishRunOnWorker()

◆ ConstructGeometry()

◆ InitGeometry()

◆ AddParticles()

◆ AddIons()

◆ GeneratePrimaries()

◆ BeginEvent()

◆ BeginPrimary()

◆ PreTrack()

◆ Stepping()

◆ PostTrack()

◆ FinishPrimary()

◆ FinishEvent()

◆ SetPrintModulo()

◆ SetVerboseLevel()

◆ SetControls()

◆ SetField()

◆ SetDebug()

◆ GetDetectorConstruction()

◆ GetCalorimeterSD()

◆ GetPrimaryGenerator()

◆ SetOldGeometry()

◆ RegisterStack()

Member Data Documentation

◆ fRootManager

◆ fPrintModulo

◆ fEventNo

◆ fVerbose

◆ fStack

◆ fDetConstruction

◆ fCalorimeterSD

◆ fPrimaryGenerator

◆ fMagField

◆ fOldGeometry

◆ fIsControls

◆ fIsMaster

◆ fIsMultiRun

◆ fSplitSimulation

◆ fG3Id

◆ fG4Id

◆ fDebug