Ex03MCApplication Class Reference

Implementation of the TVirtualMCApplication. More...

#include <Ex03MCApplication.h>

Inheritance diagram for Ex03MCApplication:

Public Member Functions
	Ex03MCApplication (const char *name, const char *title)
	Ex03MCApplication ()
virtual	~Ex03MCApplication ()
void	InitMC (const char *setup)
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)
Ex03DetectorConstruction *	GetDetectorConstruction () const
Ex03CalorimeterSD *	GetCalorimeterSD () const
Ex03PrimaryGenerator *	GetPrimaryGenerator () const
void	SetOldGeometry (Bool_t oldGeometry=kTRUE)

Private Member Functions
	Ex03MCApplication (const Ex03MCApplication &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.
Ex03MCStack *	fStack
	VMC stack.
Ex03DetectorConstruction *	fDetConstruction
	Dector construction.
Ex03CalorimeterSD *	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.

Detailed Description

Implementation of the TVirtualMCApplication.

Date

06/03/2003

Author

I. Hrivnacova; IPN, Orsay

Definition at line 39 of file Ex03MCApplication.h.

Constructor & Destructor Documentation

Ex03MCApplication() [1/3]

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

Standard constructor

Parameters

name	The MC application name
title	The MC application description

Definition at line 44 of file Ex03MCApplication.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)
{
  /// Standard constructor
  /// \param name   The MC application name
  /// \param title  The MC application description
 
  cout << "--------------------------------------------------------------"
       << endl;
  cout << " VMC Example E03" << endl;
  cout << "--------------------------------------------------------------"
       << endl;
 
  // Create a user stack
  fStack = new Ex03MCStack(1000);
 
  // Create detector construction
  fDetConstruction = new Ex03DetectorConstruction();
 
  // Create a calorimeter SD
  fCalorimeterSD = new Ex03CalorimeterSD("Calorimeter", fDetConstruction);
 
  // Create a primary generator
  fPrimaryGenerator = new Ex03PrimaryGenerator(fStack);
 
  // Constant magnetic field (in kiloGauss)
  fMagField = new TGeoUniformMagField();
}

Ex03MCApplication() [2/3]

Ex03MCApplication::Ex03MCApplication

(

)

Default constructor

Definition at line 120 of file Ex03MCApplication.cxx.

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

~Ex03MCApplication()

Ex03MCApplication::~Ex03MCApplication ( )

virtual

Destructor

Definition at line 138 of file Ex03MCApplication.cxx.

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

Ex03MCApplication() [3/3]

Ex03MCApplication::Ex03MCApplication ( const Ex03MCApplication & origin )

private

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

Parameters

origin

The source MC application

Definition at line 86 of file Ex03MCApplication.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),
    fIsMaster(kFALSE)
{
  /// Copy constructor for cloning application on workers (in multithreading
  /// mode) \param origin   The source MC application
 
  // Create new user stack
  fStack = new Ex03MCStack(1000);
 
  // Create a calorimeter SD
  fCalorimeterSD =
    new Ex03CalorimeterSD(*(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()

void Ex03MCApplication::InitMC

(

const char *

setup

)

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

Parameters

setup

The name of the configuration macro

Definition at line 175 of file Ex03MCApplication.cxx.

{
  /// Initialize MC.
  /// The selection of the concrete MC is done in the macro.
  /// \param setup The name of the configuration macro
 
  fVerbose.InitMC();
 
  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
 
  gMC->SetStack(fStack);
  gMC->SetMagField(fMagField);
  gMC->Init();
  gMC->BuildPhysics();
 
  RegisterStack();
}

RunMC()

void Ex03MCApplication::RunMC

(

Int_t

nofEvents

)

Run MC.

Parameters

nofEvents

Number of events to be processed

Definition at line 214 of file Ex03MCApplication.cxx.

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

FinishRun()

void Ex03MCApplication::FinishRun

(

)

Finish MC run.

Definition at line 226 of file Ex03MCApplication.cxx.

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

ReadEvent()

void Ex03MCApplication::ReadEvent

(

Int_t

i

)

Read i -th event and prints hits.

Parameters

i	The number of event to be read

Definition at line 273 of file Ex03MCApplication.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);
  }
 
  fCalorimeterSD->Register();
  RegisterStack();
  fRootManager->ReadEvent(i);
}

CloneForWorker()

TVirtualMCApplication * Ex03MCApplication::CloneForWorker ( ) const

virtual

Definition at line 239 of file Ex03MCApplication.cxx.

{
  return new Ex03MCApplication(*this);
}

InitOnWorker()

void Ex03MCApplication::InitOnWorker ( )

virtual

Definition at line 245 of file Ex03MCApplication.cxx.

{
  // cout << "Ex03MCApplication::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 Ex03MCApplication::FinishRunOnWorker ( )

virtual

Definition at line 263 of file Ex03MCApplication.cxx.

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

ConstructGeometry()

void Ex03MCApplication::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 288 of file Ex03MCApplication.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) {
    fDetConstruction->ConstructMaterials();
    fDetConstruction->ConstructGeometry();
    // TGeoManager::Import("geometry.root");
    // gMC->SetRootGeometry();
  }
  else {
    Ex03DetectorConstructionOld detConstructionOld;
    detConstructionOld.ConstructMaterials();
    detConstructionOld.ConstructGeometry();
  }
}

InitGeometry()

void Ex03MCApplication::InitGeometry ( )

virtual

Initialize geometry

Definition at line 310 of file Ex03MCApplication.cxx.

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

AddParticles()

void Ex03MCApplication::AddParticles ( )

virtual

Example of user defined particle with user defined decay mode

Definition at line 324 of file Ex03MCApplication.cxx.

{
  /// Example of user defined particle with user defined decay mode
 
  fVerbose.AddParticles();
 
  // Define particle
  gMC->DefineParticle(1000020050, "He5", kPTHadron, 5.03427, 2.0, 0.002, "Ion",
    0.0, 0, 1, 0, 0, 0, 0, 0, 5, kFALSE);
 
  // 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
 
  gMC->SetDecayMode(1000020050, bratio, mode);
 
  // 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)
 
  gMC->SetDecayMode(kK0Short, bratio2, mode2);
}

AddIons()

void Ex03MCApplication::AddIons ( )

virtual

Example of user defined ion

Definition at line 373 of file Ex03MCApplication.cxx.

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

GeneratePrimaries()

void Ex03MCApplication::GeneratePrimaries ( )

virtual

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

Definition at line 383 of file Ex03MCApplication.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 Ex03MCApplication::BeginEvent ( )

virtual

User actions at beginning of event

Definition at line 396 of file Ex03MCApplication.cxx.

{
  /// User actions at beginning of event
 
  fVerbose.BeginEvent();
 
  // Clear TGeo tracks (if filled)
  if (TString(gMC->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 Ex03MCApplication::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 420 of file Ex03MCApplication.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 Ex03MCApplication::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 434 of file Ex03MCApplication.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 Ex03MCStack
 
      cout << "      Current track " << fStack->GetCurrentTrack()->GetName()
           << "  is a decay product of Parent ID = "
           << fStack->GetCurrentParentTrackNumber() << endl;
    }
  }
}

Stepping()

void Ex03MCApplication::Stepping ( )

virtual

User actions at each step

Definition at line 461 of file Ex03MCApplication.cxx.

{
  /// User actions at each step
 
  // Work around for Fluka VMC, which does not call
  // MCApplication::PreTrack()
  //
 
  static Int_t trackId = 0;
  if (TString(gMC->GetName()) == "TFluka" &&
      gMC->GetStack()->GetCurrentTrackNumber() != trackId) {
    fVerbose.PreTrack();
    trackId = gMC->GetStack()->GetCurrentTrackNumber();
  }
 
  fVerbose.Stepping();
 
  fCalorimeterSD->ProcessHits();
}

PostTrack()

void Ex03MCApplication::PostTrack ( )

virtual

User actions after finishing of each track

Definition at line 482 of file Ex03MCApplication.cxx.

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

FinishPrimary()

void Ex03MCApplication::FinishPrimary ( )

virtual

User actions after finishing of a primary track

Definition at line 490 of file Ex03MCApplication.cxx.

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

FinishEvent()

void Ex03MCApplication::FinishEvent ( )

virtual

User actions after finishing of an event

Definition at line 502 of file Ex03MCApplication.cxx.

{
  /// User actions after finishing of an event
 
  fVerbose.FinishEvent();
 
  // Geant3 + TGeo
  // (use TGeo functions for visualization)
  if (TString(gMC->GetName()) == "TGeant3TGeo") {
 
    // Draw volume
    gGeoManager->SetVisOption(0);
    gGeoManager->SetTopVisible();
    gGeoManager->GetTopVolume()->Draw();
 
    // Draw tracks (if filled)
    // Available when this feature is activated via
    // gMC->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 Ex03MCApplication::SetPrintModulo ( Int_t value )

inline

Set the event modulus number to be printed

Parameters

value

The new event modulus number value

Definition at line 109 of file Ex03MCApplication.h.

{
  fPrintModulo = value;
}

SetVerboseLevel()

void Ex03MCApplication::SetVerboseLevel ( Int_t verboseLevel )

inline

Set verbosity

Parameters

verboseLevel

The new verbose level value

Definition at line 116 of file Ex03MCApplication.h.

{
  fVerbose.SetLevel(verboseLevel);
}

SetControls()

void Ex03MCApplication::SetControls ( Bool_t isControls )

inline

Switch on/off special process controls

Parameters

isControls

If true, special process controls setting is activated

Definition at line 156 of file Ex03MCApplication.h.

{
  fIsControls = isControls;
}

SetField()

void Ex03MCApplication::SetField ( Double_t bz )

inline

Definition at line 123 of file Ex03MCApplication.h.

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

GetDetectorConstruction()

Ex03DetectorConstruction * Ex03MCApplication::GetDetectorConstruction ( ) const

inline

Returns

The detector construction

Definition at line 130 of file Ex03MCApplication.h.

{
  return fDetConstruction;
}

GetCalorimeterSD()

Ex03CalorimeterSD * Ex03MCApplication::GetCalorimeterSD ( ) const

inline

Returns

The calorimeter sensitive detector

Definition at line 136 of file Ex03MCApplication.h.

{
  return fCalorimeterSD;
}

GetPrimaryGenerator()

Ex03PrimaryGenerator * Ex03MCApplication::GetPrimaryGenerator ( ) const

inline

Returns

The primary generator

Definition at line 142 of file Ex03MCApplication.h.

{
  return fPrimaryGenerator;
}

SetOldGeometry()

void Ex03MCApplication::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 149 of file Ex03MCApplication.h.

{
  fOldGeometry = oldGeometry;
}

RegisterStack()

void Ex03MCApplication::RegisterStack ( ) const

private

Register stack in the Root manager.

Definition at line 160 of file Ex03MCApplication.cxx.

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

Member Data Documentation

fRootManager

TMCRootManager* Ex03MCApplication::fRootManager

mutableprivate

Root manager.

Definition at line 89 of file Ex03MCApplication.h.

fPrintModulo

Int_t Ex03MCApplication::fPrintModulo

private

The event modulus number to be printed.

Definition at line 90 of file Ex03MCApplication.h.

fEventNo

Int_t Ex03MCApplication::fEventNo

private

Event counter.

Definition at line 91 of file Ex03MCApplication.h.

fVerbose

TMCVerbose Ex03MCApplication::fVerbose

private

VMC verbose helper.

Definition at line 92 of file Ex03MCApplication.h.

fStack

Ex03MCStack* Ex03MCApplication::fStack

private

VMC stack.

Definition at line 93 of file Ex03MCApplication.h.

fDetConstruction

Ex03DetectorConstruction* Ex03MCApplication::fDetConstruction

private

Dector construction.

Definition at line 94 of file Ex03MCApplication.h.

fCalorimeterSD

Ex03CalorimeterSD* Ex03MCApplication::fCalorimeterSD

private

Calorimeter SD.

Definition at line 95 of file Ex03MCApplication.h.

fPrimaryGenerator

Ex03PrimaryGenerator* Ex03MCApplication::fPrimaryGenerator

private

Primary generator.

Definition at line 96 of file Ex03MCApplication.h.

fMagField

TGeoUniformMagField* Ex03MCApplication::fMagField

private

Magnetic field.

Definition at line 97 of file Ex03MCApplication.h.

fOldGeometry

Bool_t Ex03MCApplication::fOldGeometry

private

Option for geometry definition.

Definition at line 98 of file Ex03MCApplication.h.

fIsControls

Bool_t Ex03MCApplication::fIsControls

private

Option to activate special controls.

Definition at line 99 of file Ex03MCApplication.h.

fIsMaster

Bool_t Ex03MCApplication::fIsMaster

private

If is on master thread.

Definition at line 100 of file Ex03MCApplication.h.

---

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

• examples/E03/E03a/include/Ex03MCApplication.h
• examples/E03/E03a/src/Ex03MCApplication.cxx