Ex02MCApplication Class Reference

Implementation of the TVirtualMCApplication. More...

#include <Ex02MCApplication.h>

Inheritance diagram for Ex02MCApplication:

Public Member Functions
	Ex02MCApplication (const char *name, const char *title)
	Ex02MCApplication ()
virtual	~Ex02MCApplication ()
void	InitMC (const char *setup)
void	RunMC (Int_t nofEvents)
void	FinishRun ()
virtual TVirtualMCApplication *	CloneForWorker () const
virtual void	InitOnWorker ()
virtual void	FinishRunOnWorker ()
virtual void	ConstructGeometry ()
virtual void	InitGeometry ()
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	ReadEvent (Int_t i)
void	SetOldGeometry (Bool_t oldGeometry=kTRUE)

Private Member Functions
	Ex02MCApplication (const Ex02MCApplication &origin)
void	RegisterStack () const

Private Attributes
TMCRootManager *	fRootManager
	Root manager.
Ex02MCStack *	fStack
	VMC stack.
Ex02DetectorConstruction	fDetConstruction
	Dector construction.
Ex02TrackerSD *	fTrackerSD
	Tracker SD.
TVirtualMagField *	fMagField
	Magnetic field.
Bool_t	fOldGeometry
	Option for geometry definition.

Detailed Description

Implementation of the TVirtualMCApplication.

Date

21/04/2002

Author

I. Hrivnacova; IPN, Orsay

Definition at line 37 of file Ex02MCApplication.h.

Constructor & Destructor Documentation

Ex02MCApplication() [1/3]

Ex02MCApplication::Ex02MCApplication	(	const char *	name,
		const char *	title )

Standard constructor

Parameters

name	The MC application name
title	The MC application description

Definition at line 39 of file Ex02MCApplication.cxx.

  : TVirtualMCApplication(name, title),
    fRootManager(0),
    fStack(0),
    fDetConstruction(),
    fTrackerSD(0),
    fMagField(0),
    fOldGeometry(kFALSE)
{
  /// Standard constructor
  /// \param name   The MC application name
  /// \param title  The MC application description
 
  // cout << "Ex02MCApplication::Ex02MCApplication " << this << endl;
 
  // Create application data
 
  // Create SD
  fTrackerSD = new Ex02TrackerSD("Tracker Chamber");
  // Create a user stack
  fStack = new Ex02MCStack(100);
  // Constant magnetic field (in kiloGauss)
  fMagField = new Ex02MagField(20., 0., 0.);
  // It si also possible to use TGeoUniformMagField class:
  // fMagField = new TGeoUniformMagField(20., 0., 0.);
}

Ex02MCApplication() [2/3]

Ex02MCApplication::Ex02MCApplication

(

)

Default constructor

Definition at line 94 of file Ex02MCApplication.cxx.

  : TVirtualMCApplication(),
    fRootManager(0),
    fStack(0),
    fDetConstruction(),
    fTrackerSD(),
    fMagField(0),
    fOldGeometry(kFALSE)
{
  /// Default constructor
}

~Ex02MCApplication()

Ex02MCApplication::~Ex02MCApplication ( )

virtual

Destructor

Definition at line 107 of file Ex02MCApplication.cxx.

{
  /// Destructor
 
  // cout << "Ex02MCApplication::~Ex02MCApplication " << this << endl;
 
  delete fRootManager;
  delete fStack;
  delete fTrackerSD;
  delete fMagField;
  delete gMC;
 
  // cout << "Done Ex02MCApplication::~Ex02MCApplication " << this << endl;
}

Ex02MCApplication() [3/3]

Ex02MCApplication::Ex02MCApplication ( const Ex02MCApplication & origin )

private

Copy constructor (for clonig on worker thread in MT mode).

Parameters

origin

The source object (on master).

Definition at line 67 of file Ex02MCApplication.cxx.

  : TVirtualMCApplication(origin.GetName(), origin.GetTitle()),
    fRootManager(0),
    fStack(0),
    fDetConstruction(origin.fDetConstruction),
    fTrackerSD(0),
    fMagField(0),
    fOldGeometry(kFALSE)
{
  /// Copy constructor (for clonig on worker thread in MT mode).
  /// \param origin  The source object (on master).
 
  // cout << "Ex02MCApplication::Ex02MCApplication " << this << endl;
 
  // Create application data
 
  // Create SD
  fTrackerSD = new Ex02TrackerSD(*(origin.fTrackerSD));
  // Create a user stack
  fStack = new Ex02MCStack(100);
  // Constant magnetic field (in kiloGauss)
  fMagField = new Ex02MagField(20., 0., 0.);
  // It si also possible to use TGeoUniformMagField class:
  // fMagField = new TGeoUniformMagField(20., 0., 0.);
}

Member Function Documentation

InitMC()

void Ex02MCApplication::InitMC

(

const char *

setup

)

Initialize MC from Config.C macro The selection of the concrete MC is done in the macro.

Parameters

setup

The name of the configuration macro

Definition at line 141 of file Ex02MCApplication.cxx.

{
  /// Initialize MC from Config.C macro
  /// The selection of the concrete MC is done in the macro.
  /// \param setup The name of the configuration macro
 
  if (TString(setup) != "") {
    gROOT->LoadMacro(setup);
    gInterpreter->ProcessLine("Config()");
    if (!gMC) {
      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 (!gMC->IsMT()) {
    fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite);
    // fRootManager->SetDebug(true);
  }
#else
  // Create Root manager
  fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite);
  // fRootManager->SetDebug(true);
#endif
 
  // Set data to MC
  gMC->SetStack(fStack);
  gMC->SetMagField(fMagField);
 
  // Init MC
  gMC->Init();
  gMC->BuildPhysics();
 
  RegisterStack();
}

RunMC()

void Ex02MCApplication::RunMC

(

Int_t

nofEvents

)

Run MC.

Parameters

nofEvents

Number of events to be processed

Definition at line 181 of file Ex02MCApplication.cxx.

{
  /// Run MC.
  /// \param nofEvents Number of events to be processed
 
  gMC->ProcessRun(nofEvents);
  FinishRun();
}

FinishRun()

void Ex02MCApplication::FinishRun

(

)

Finish MC run.

Definition at line 191 of file Ex02MCApplication.cxx.

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

CloneForWorker()

TVirtualMCApplication * Ex02MCApplication::CloneForWorker ( ) const

virtual

Definition at line 204 of file Ex02MCApplication.cxx.

{
  return new Ex02MCApplication(*this);
}

InitOnWorker()

void Ex02MCApplication::InitOnWorker ( )

virtual

Definition at line 210 of file Ex02MCApplication.cxx.

{
  // cout << "Ex02MCApplication::InitForWorker " << this << endl;
 
  // Create Root manager
  Int_t threadRank = 1;
        // The real thread rank will be set in MCRootManager
  fRootManager = new TMCRootManager(GetName(), TMCRootManager::kWrite, threadRank);
  // fRootManager->SetDebug(true);
 
  // Set data to MC
  gMC->SetStack(fStack);
  gMC->SetMagField(fMagField);
 
  RegisterStack();
}

FinishRunOnWorker()

void Ex02MCApplication::FinishRunOnWorker ( )

virtual

Definition at line 228 of file Ex02MCApplication.cxx.

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

ConstructGeometry()

void Ex02MCApplication::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 238 of file Ex02MCApplication.cxx.

{
  /// Construct geometry using detector contruction class.
  /// The detector contruction class is using TGeo functions or
  /// TVirtualMC functions (if oldGeometry is selected)
 
  // Cannot use Root geometry if not supported with
  // selected MC
  if (!fOldGeometry && !gMC->IsRootGeometrySupported()) {
    cerr << "Selected MC does not support TGeo geometry" << endl;
    cerr << "Exiting program" << endl;
    exit(1);
  }
 
  if (!fOldGeometry) {
    cout << "Geometry will be defined via TGeo" << endl;
    fDetConstruction.ConstructMaterials();
    fDetConstruction.ConstructGeometry();
  }
  else {
    cout << "Geometry will be defined via VMC" << endl;
    Ex02DetectorConstructionOld detConstructionOld;
    detConstructionOld.ConstructMaterials();
    detConstructionOld.ConstructGeometry();
  }
}

InitGeometry()

void Ex02MCApplication::InitGeometry ( )

virtual

Initialize geometry

Definition at line 266 of file Ex02MCApplication.cxx.

{
  /// Initialize geometry
 
  fTrackerSD->Initialize();
}

GeneratePrimaries()

void Ex02MCApplication::GeneratePrimaries ( )

virtual

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

Definition at line 274 of file Ex02MCApplication.cxx.

{
  /// Fill the user stack (derived from TVirtualMCStack) with primary particles.
 
  // Track ID (filled by stack)
  Int_t ntr;
 
  // Option: to be tracked
  Int_t toBeDone = 1;
 
  // Particle type
  // Int_t pdg  = 0;    // geantino
  Int_t pdg = kProton;
 
  // Polarization
  Double_t polx = 0.;
  Double_t poly = 0.;
  Double_t polz = 0.;
 
  // Position
  Double_t vx = 0.;
  Double_t vy = 0.;
  Double_t vz = -0.5 * (fDetConstruction.GetWorldFullLength());
  Double_t tof = 0.;
 
  // Energy
  Double_t kinEnergy = 3.0;
  Double_t mass = 0.9382723;
  Double_t e = mass + kinEnergy;
 
  // Momentum
  Double_t px, py, pz;
  px = 0.;
  py = 0.;
  pz = sqrt(e * e - mass * mass);
 
  // Add particle to stack
  fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof, polx,
    poly, polz, kPPrimary, ntr, 1., 0);
}

BeginEvent()

void Ex02MCApplication::BeginEvent ( )

virtual

User actions at beginning of event. Nothing to be done this example

Definition at line 316 of file Ex02MCApplication.cxx.

{
  /// User actions at beginning of event.
  /// Nothing to be done this example
 
  fStack->SetObjectNumber();
}

BeginPrimary()

void Ex02MCApplication::BeginPrimary ( )

virtual

User actions at beginning of a primary track. Nothing to be done this example

Definition at line 325 of file Ex02MCApplication.cxx.

{
  /// User actions at beginning of a primary track.
  /// Nothing to be done this example
}

PreTrack()

void Ex02MCApplication::PreTrack ( )

virtual

User actions at beginning of each track. Nothing to be done this example

Definition at line 332 of file Ex02MCApplication.cxx.

{
  /// User actions at beginning of each track.
  /// Nothing to be done this example
}

Stepping()

void Ex02MCApplication::Stepping ( )

virtual

User actions at each step

Definition at line 339 of file Ex02MCApplication.cxx.

{
  /// User actions at each step
 
  fTrackerSD->ProcessHits();
}

PostTrack()

void Ex02MCApplication::PostTrack ( )

virtual

User actions at each step. Nothing to be done this example

Definition at line 347 of file Ex02MCApplication.cxx.

{
  /// User actions at each step.
  /// Nothing to be done this example
}

FinishPrimary()

void Ex02MCApplication::FinishPrimary ( )

virtual

User actions after finishing of a primary track. Nothing to be done this example

Definition at line 354 of file Ex02MCApplication.cxx.

{
  /// User actions after finishing of a primary track.
  /// Nothing to be done this example
}

FinishEvent()

void Ex02MCApplication::FinishEvent ( )

virtual

User actions after finishing of an event Nothing to be done this example

Definition at line 361 of file Ex02MCApplication.cxx.

{
  /// User actions after finishing of an event
  /// Nothing to be done this example
 
  // Geant4 own visualization is activated via G4 macro (g4config.in)
 
  // TGeo visualization
  if (gGeoManager && gGeoManager->GetListOfTracks() &&
      gGeoManager->GetTrack(0) &&
      ((TVirtualGeoTrack*)gGeoManager->GetTrack(0))->HasPoints()) {
 
    gGeoManager->SetVisOption(0);
    gGeoManager->SetTopVisible();
    gGeoManager->DrawTracks("/*"); // this means all tracks
  }
 
  fRootManager->Fill();
 
  fTrackerSD->EndOfEvent();
 
  fStack->Print();
  fStack->Reset();
}

ReadEvent()

void Ex02MCApplication::ReadEvent

(

Int_t

i

)

Read i -th event and prints hits.

Parameters

i	The number of event to be read

Definition at line 387 of file Ex02MCApplication.cxx.

{
  /// Read \em i -th event and prints hits.
  /// \param i The number of event to be read
 
  if (!fRootManager) {
    fRootManager = new TMCRootManager(GetName(), TMCRootManager::kRead);
  }
 
  fTrackerSD->Register();
  RegisterStack();
  fRootManager->ReadEvent(i);
 
  fStack->Print();
  fTrackerSD->Print();
}

SetOldGeometry()

void Ex02MCApplication::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 89 of file Ex02MCApplication.h.

{
  fOldGeometry = oldGeometry;
}

RegisterStack()

void Ex02MCApplication::RegisterStack ( ) const

private

Register stack in the Root manager.

Definition at line 126 of file Ex02MCApplication.cxx.

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

Member Data Documentation

fRootManager

TMCRootManager* Ex02MCApplication::fRootManager

mutableprivate

Root manager.

Definition at line 75 of file Ex02MCApplication.h.

fStack

Ex02MCStack* Ex02MCApplication::fStack

private

VMC stack.

Definition at line 76 of file Ex02MCApplication.h.

fDetConstruction

Ex02DetectorConstruction Ex02MCApplication::fDetConstruction

private

Dector construction.

Definition at line 77 of file Ex02MCApplication.h.

fTrackerSD

Ex02TrackerSD* Ex02MCApplication::fTrackerSD

private

Tracker SD.

Definition at line 78 of file Ex02MCApplication.h.

fMagField

TVirtualMagField* Ex02MCApplication::fMagField

private

Magnetic field.

Definition at line 79 of file Ex02MCApplication.h.

fOldGeometry

Bool_t Ex02MCApplication::fOldGeometry

private

Option for geometry definition.

Definition at line 80 of file Ex02MCApplication.h.

---

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

• examples/E02/include/Ex02MCApplication.h
• examples/E02/src/Ex02MCApplication.cxx