Ex06MCApplication Class Reference

Implementation of the TVirtualMCApplication. More...

#include <Ex06MCApplication.h>

Inheritance diagram for Ex06MCApplication:

Public Member Functions
	Ex06MCApplication (const char *name, const char *title)
	Ex06MCApplication ()
virtual	~Ex06MCApplication ()
void	InitMC (const char *setup)
void	RunMC (Int_t nofEvents)
virtual TVirtualMCApplication *	CloneForWorker () const
virtual void	InitOnWorker ()
virtual void	Merge (TVirtualMCApplication *localMCApplication)
virtual void	ConstructGeometry ()
virtual void	ConstructOpGeometry ()
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	FinishRun ()
void	SetVerboseLevel (Int_t verboseLevel)
Ex06PrimaryGenerator *	GetPrimaryGenerator () const
void	SetOldGeometry (Bool_t oldGeometry=kTRUE)
void	SetTestStackPopper (Bool_t option=kFALSE)

Static Public Member Functions
static Ex06MCApplication *	Instance ()

Private Member Functions
	Ex06MCApplication (const Ex06MCApplication &origin)
void	GenerateFeedback (Int_t nofPhotons, TLorentzVector momentum)

Private Attributes
Int_t	fGammaCounter
	Optical photons counter.
Int_t	fFeedbackCounter
	Feedback photons counter.
Int_t	fRunGammaCounter
	Optical photons counter2.
Int_t	fRunFeedbackCounter
	Feedback photons counter2.
TMCVerbose	fVerbose
	VMC verbose helper.
Ex03MCStack *	fStack
	VMC stack.
TVirtualMagField *	fMagField
	The magnetic field.
Ex06DetectorConstruction *	fDetConstruction
	Dector construction.
Ex06PrimaryGenerator *	fPrimaryGenerator
	Primary generator.
Bool_t	fOldGeometry
	Option for geometry definition.
Bool_t	fTestStackPopper
	Option for stack popper test.
Bool_t	fIsMaster
	If is on master thread.

Detailed Description

Implementation of the TVirtualMCApplication.

Date

16/05/2005

Author

I. Hrivnacova; IPN, Orsay

Definition at line 37 of file Ex06MCApplication.h.

Constructor & Destructor Documentation

Ex06MCApplication() [1/3]

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

Standard constructor

Parameters

name	The MC application name
title	The MC application description

Definition at line 40 of file Ex06MCApplication.cxx.

  : TVirtualMCApplication(name, title),
    fGammaCounter(0),
    fFeedbackCounter(0),
    fRunGammaCounter(0),
    fRunFeedbackCounter(0),
    fVerbose(0),
    fStack(0),
    fMagField(0),
    fDetConstruction(0),
    fPrimaryGenerator(0),
    fOldGeometry(kFALSE),
    fTestStackPopper(kFALSE),
    fIsMaster(kTRUE)
{
  /// Standard constructor
  /// \param name   The MC application name
  /// \param title  The MC application description
 
  // Create a user stack
  fStack = new Ex03MCStack(1000);
 
  // create magnetic field (with zero value)
  fMagField = new TGeoUniformMagField();
 
  // Create detector construction
  fDetConstruction = new Ex06DetectorConstruction();
 
  // Create a primary generator
  fPrimaryGenerator = new Ex06PrimaryGenerator(fStack);
}

Ex06MCApplication() [2/3]

Ex06MCApplication::Ex06MCApplication

(

)

Default constructor

Definition at line 104 of file Ex06MCApplication.cxx.

  : TVirtualMCApplication(),
    fGammaCounter(0),
    fFeedbackCounter(0),
    fRunGammaCounter(0),
    fRunFeedbackCounter(0),
    fVerbose(0),
    fStack(0),
    fMagField(0),
    fDetConstruction(0),
    fPrimaryGenerator(0),
    fOldGeometry(kFALSE),
    fTestStackPopper(kFALSE),
    fIsMaster(kTRUE)
{
  /// Default constructor
}

~Ex06MCApplication()

Ex06MCApplication::~Ex06MCApplication ( )

virtual

Destructor

Definition at line 123 of file Ex06MCApplication.cxx.

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

Ex06MCApplication() [3/3]

Ex06MCApplication::Ex06MCApplication ( const Ex06MCApplication & origin )

private

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

Parameters

origin

The source object (on master).

Definition at line 73 of file Ex06MCApplication.cxx.

  : TVirtualMCApplication(origin.GetName(), origin.GetTitle()),
    fGammaCounter(0),
    fFeedbackCounter(0),
    fRunGammaCounter(0),
    fRunFeedbackCounter(0),
    fVerbose(origin.fVerbose),
    fStack(0),
    fMagField(0),
    fDetConstruction(origin.fDetConstruction),
    fPrimaryGenerator(0),
    fOldGeometry(origin.fOldGeometry),
    fTestStackPopper(origin.fTestStackPopper),
    fIsMaster(kFALSE)
{
  /// Copy constructor (for clonig on worker thread in MT mode).
  /// \param origin  The source object (on master).
 
  // Create a user stack
  fStack = new Ex03MCStack(1000);
 
  // create magnetic field (with zero value)
  // TODO: check copying
  fMagField = new TGeoUniformMagField();
 
  // Create a primary generator
  fPrimaryGenerator =
    new Ex06PrimaryGenerator(*(origin.fPrimaryGenerator), fStack);
}

Member Function Documentation

Instance()

Ex06MCApplication * Ex06MCApplication::Instance ( )

inlinestatic

Returns

The singleton instance

Definition at line 105 of file Ex06MCApplication.h.

{
  return (Ex06MCApplication*)(TVirtualMCApplication::Instance());
}

InitMC()

void Ex06MCApplication::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 143 of file Ex06MCApplication.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.)");
    }
  }
 
  gMC->SetStack(fStack);
  gMC->SetMagField(fMagField);
  gMC->Init();
  gMC->BuildPhysics();
}

RunMC()

void Ex06MCApplication::RunMC

(

Int_t

nofEvents

)

Run MC.

Parameters

nofEvents

Number of events to be processed

Definition at line 167 of file Ex06MCApplication.cxx.

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

CloneForWorker()

TVirtualMCApplication * Ex06MCApplication::CloneForWorker ( ) const

virtual

Definition at line 180 of file Ex06MCApplication.cxx.

{
  return new Ex06MCApplication(*this);
}

InitOnWorker()

void Ex06MCApplication::InitOnWorker ( )

virtual

Definition at line 186 of file Ex06MCApplication.cxx.

{
  // cout << "Ex06MCApplication::InitForWorker " << this << endl;
 
  // Set data to MC
  gMC->SetStack(fStack);
  gMC->SetMagField(fMagField);
}

Merge()

void Ex06MCApplication::Merge ( TVirtualMCApplication * localMCApplication )

virtual

Definition at line 196 of file Ex06MCApplication.cxx.

{
  // cout << "Ex06MCApplication::Merge " << this << endl;
 
  Ex06MCApplication* ex06LocalMCApplication =
    static_cast<Ex06MCApplication*>(localMCApplication);
 
  fRunGammaCounter += ex06LocalMCApplication->fRunGammaCounter;
  fRunFeedbackCounter += ex06LocalMCApplication->fRunFeedbackCounter;
}

ConstructGeometry()

void Ex06MCApplication::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 208 of file Ex06MCApplication.cxx.

{
  /// Construct geometry using detector contruction class.
  /// The detector contruction class is using TGeo functions or
  /// TVirtualMC functions (if oldGeometry is selected)
 
  fVerbose.ConstructGeometry();
 
  // 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;
    Ex06DetectorConstructionOld detConstructionOld;
    detConstructionOld.ConstructMaterials();
    detConstructionOld.ConstructGeometry();
  }
}

ConstructOpGeometry()

void Ex06MCApplication::ConstructOpGeometry ( )

virtual

Define material optical properties

Definition at line 238 of file Ex06MCApplication.cxx.

{
  /// Define material optical properties
 
  fVerbose.ConstructGeometry();
 
  fDetConstruction->ConstructOpGeometry();
}

InitGeometry()

void Ex06MCApplication::InitGeometry ( )

virtual

Initialize geometry

Definition at line 248 of file Ex06MCApplication.cxx.

{
  /// Initialize geometry
 
  fVerbose.InitGeometry();
}

GeneratePrimaries()

void Ex06MCApplication::GeneratePrimaries ( )

virtual

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

Definition at line 256 of file Ex06MCApplication.cxx.

{
  /// Fill the user stack (derived from TVirtualMCStack) with primary particles.
 
  fVerbose.GeneratePrimaries();
 
  fPrimaryGenerator->GeneratePrimaries();
}

BeginEvent()

void Ex06MCApplication::BeginEvent ( )

virtual

User actions at beginning of event

Definition at line 266 of file Ex06MCApplication.cxx.

{
  /// User actions at beginning of event
 
  fVerbose.BeginEvent();
 
  fGammaCounter = 0;
  fFeedbackCounter = 0;
}

BeginPrimary()

void Ex06MCApplication::BeginPrimary ( )

virtual

User actions at beginning of a primary track

Definition at line 277 of file Ex06MCApplication.cxx.

{
  /// User actions at beginning of a primary track
 
  fVerbose.BeginPrimary();
}

PreTrack()

void Ex06MCApplication::PreTrack ( )

virtual

User actions at beginning of each track

Definition at line 285 of file Ex06MCApplication.cxx.

{
  /// User actions at beginning of each track
 
  fVerbose.PreTrack();
 
  if (gMC->TrackPid() == 50000050) {
    fGammaCounter++;
    fRunGammaCounter++;
  }
  if (gMC->TrackPid() == 50000051) {
    fFeedbackCounter++;
    fRunFeedbackCounter++;
  }
}

Stepping()

void Ex06MCApplication::Stepping ( )

virtual

User actions at each step

Definition at line 302 of file Ex06MCApplication.cxx.

{
  /// User actions at each step
 
  fVerbose.Stepping();
 
  // 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();
    if (gMC->TrackPid() == 50000050) fGammaCounter++;
    if (gMC->TrackPid() == 50000051) fFeedbackCounter++;
  }
 
  if (!fTestStackPopper) return;
 
  // Stack popper test
  // Add 1 feedback photon (50000051, fixed momentum) when a charged track
  // enters TANK and 3 feedback photons (momentum in opposite direction to
  // parent photon)
  //  when a photon is stopped in TANK
 
  // Charged particles entering in TANK
  TString volName = gMC->CurrentVolName();
  if (volName != "TANK") return;
 
  if ((gMC->TrackCharge() != 0.) && (gMC->IsTrackEntering())) {
    // 1 keV
    Double_t energy = 1e-06;
    // TLorentzVector momentum(energy, 0., 0., energy);
    // workaround for a problem in Geant4 11 in G4OpBoundaryProcess
    TLorentzVector momentum(
      energy, 0.1 * energy, 0.1 * energy, 0.98994949 * energy);
    GenerateFeedback(1, momentum);
  }
 
  if ((gMC->TrackPid() == 50000050) && (gMC->Edep() > 0.)) {
    TLorentzVector momentum;
    gMC->TrackMomentum(momentum);
    momentum = -1. * momentum;
    GenerateFeedback(3, momentum);
  }
}

PostTrack()

void Ex06MCApplication::PostTrack ( )

virtual

User actions after finishing of each track

Definition at line 378 of file Ex06MCApplication.cxx.

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

FinishPrimary()

void Ex06MCApplication::FinishPrimary ( )

virtual

User actions after finishing of a primary track

Definition at line 386 of file Ex06MCApplication.cxx.

{
  /// User actions after finishing of a primary track
 
  fVerbose.FinishPrimary();
}

FinishEvent()

void Ex06MCApplication::FinishEvent ( )

virtual

User actions after finishing of an event

Definition at line 394 of file Ex06MCApplication.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
    }
  }
 
  cout << "Number of optical photons produced in this event : " << fGammaCounter
       << endl;
 
  if (fTestStackPopper) {
    cout << "Number of feedback photons produced in this event : "
         << fFeedbackCounter << endl;
  }
  fStack->Reset();
}

FinishRun()

void Ex06MCApplication::FinishRun

(

)

User actions after finishing of a run

Definition at line 430 of file Ex06MCApplication.cxx.

{
  /// User actions after finishing of a run
 
  fVerbose.FinishRun();
 
  cout << "Number of optical photons produced in this run : "
       << fRunGammaCounter << endl;
 
  if (fTestStackPopper) {
    cout << "Number of feedback photons produced in this run : "
         << fRunFeedbackCounter << endl;
  }
}

SetVerboseLevel()

void Ex06MCApplication::SetVerboseLevel ( Int_t verboseLevel )

inline

Set verbosity

Parameters

verboseLevel

The new verbose level value

Definition at line 112 of file Ex06MCApplication.h.

{
  fVerbose.SetLevel(verboseLevel);
}

GetPrimaryGenerator()

Ex06PrimaryGenerator * Ex06MCApplication::GetPrimaryGenerator ( ) const

inline

Returns

The primary generator

Definition at line 118 of file Ex06MCApplication.h.

{
  return fPrimaryGenerator;
}

SetOldGeometry()

void Ex06MCApplication::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 125 of file Ex06MCApplication.h.

{
  fOldGeometry = oldGeometry;
}

SetTestStackPopper()

void Ex06MCApplication::SetTestStackPopper ( Bool_t option = kFALSE )

inline

Switch on/off the stack popper test

Parameters

option

stack popper test activation option

Definition at line 132 of file Ex06MCApplication.h.

{
  fTestStackPopper = option;
}

GenerateFeedback()

void Ex06MCApplication::GenerateFeedback ( Int_t nofPhotons, TLorentzVector momentum )

private

Generate FeedBack photons for the current particle.

Definition at line 352 of file Ex06MCApplication.cxx.

{
  /// Generate FeedBack photons for the current particle.
 
  for (Int_t i = 0; i < nofPhotons; ++i) {
    // same position as the parent track
    TLorentzVector position;
    gMC->TrackPosition(position);
    // Feedback photon
    Int_t pdgEncoding = 50000051;
    // Fixed polarization
    Double_t polarization[3];
    polarization[0] = 0.3;
    polarization[1] = 0.4;
    polarization[2] = 0.866025403784438597;
 
    Int_t ntrack;
    gMC->GetStack()->PushTrack(1, gMC->GetStack()->GetCurrentTrackNumber(),
      pdgEncoding, momentum.X(), momentum.Y(), momentum.Z(), momentum.T(),
      position.X(), position.Y(), position.Z(), position.T(), polarization[0],
      polarization[1], polarization[2], kPFeedBackPhoton, ntrack, 1.0, 0);
  }
}

Member Data Documentation

fGammaCounter

Int_t Ex06MCApplication::fGammaCounter

private

Optical photons counter.

Definition at line 86 of file Ex06MCApplication.h.

fFeedbackCounter

Int_t Ex06MCApplication::fFeedbackCounter

private

Feedback photons counter.

Definition at line 87 of file Ex06MCApplication.h.

fRunGammaCounter

Int_t Ex06MCApplication::fRunGammaCounter

private

Optical photons counter2.

Definition at line 88 of file Ex06MCApplication.h.

fRunFeedbackCounter

Int_t Ex06MCApplication::fRunFeedbackCounter

private

Feedback photons counter2.

Definition at line 89 of file Ex06MCApplication.h.

fVerbose

TMCVerbose Ex06MCApplication::fVerbose

private

VMC verbose helper.

Definition at line 90 of file Ex06MCApplication.h.

fStack

Ex03MCStack* Ex06MCApplication::fStack

private

VMC stack.

Definition at line 91 of file Ex06MCApplication.h.

fMagField

TVirtualMagField* Ex06MCApplication::fMagField

private

The magnetic field.

Definition at line 92 of file Ex06MCApplication.h.

fDetConstruction

Ex06DetectorConstruction* Ex06MCApplication::fDetConstruction

private

Dector construction.

Definition at line 93 of file Ex06MCApplication.h.

fPrimaryGenerator

Ex06PrimaryGenerator* Ex06MCApplication::fPrimaryGenerator

private

Primary generator.

Definition at line 94 of file Ex06MCApplication.h.

fOldGeometry

Bool_t Ex06MCApplication::fOldGeometry

private

Option for geometry definition.

Definition at line 95 of file Ex06MCApplication.h.

fTestStackPopper

Bool_t Ex06MCApplication::fTestStackPopper

private

Option for stack popper test.

Definition at line 96 of file Ex06MCApplication.h.

fIsMaster

Bool_t Ex06MCApplication::fIsMaster

private

If is on master thread.

Definition at line 97 of file Ex06MCApplication.h.

---

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

• examples/E06/include/Ex06MCApplication.h
• examples/E06/src/Ex06MCApplication.cxx