Ex01MCApplication Class Reference

Implementation of the TVirtualMCApplication. More...

#include <Ex01MCApplication.h>

Inheritance diagram for Ex01MCApplication:

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

Static Public Member Functions
static Ex01MCApplication *	Instance ()

Private Member Functions
void	ConstructMaterials ()
void	ConstructVolumes ()

Private Attributes
TVirtualMCStack *	fStack
	The VMC stack.
TVirtualMagField *	fMagField
	The magnetic field.
Int_t	fImedAr
	The Argon gas medium Id.
Int_t	fImedAl
	The Aluminium medium Id.
Int_t	fImedPb
	The Lead medium Id.
Bool_t	fOldGeometry
	Option for geometry definition.

Detailed Description

Implementation of the TVirtualMCApplication.

Date

05/04/2002

Author

I. Hrivnacova; IPN, Orsay

Definition at line 31 of file Ex01MCApplication.h.

Constructor & Destructor Documentation

Ex01MCApplication() [1/2]

Ex01MCApplication::Ex01MCApplication	(	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 Ex01MCApplication.cxx.

  : TVirtualMCApplication(name, title),
    fStack(0),
    fMagField(0),
    fImedAr(0),
    fImedAl(0),
    fImedPb(0),
    fOldGeometry(kFALSE)
{
  /// Standard constructor
  /// \param name   The MC application name
  /// \param title  The MC application description
 
  // create a user stack
  fStack = new Ex01MCStack(100);
 
  // create magnetic field (with zero value)
  fMagField = new TGeoUniformMagField();
}

Ex01MCApplication() [2/2]

Ex01MCApplication::Ex01MCApplication

(

)

Default constructor

Definition at line 61 of file Ex01MCApplication.cxx.

  : TVirtualMCApplication(),
    fStack(0),
    fMagField(0),
    fImedAr(0),
    fImedAl(0),
    fImedPb(0),
    fOldGeometry(kFALSE)
{
  /// Default constructor
}

~Ex01MCApplication()

Ex01MCApplication::~Ex01MCApplication ( )

virtual

Destructor

Definition at line 74 of file Ex01MCApplication.cxx.

{
  /// Destructor
 
  delete fStack;
  delete fMagField;
  delete gMC;
}

Member Function Documentation

Instance()

Ex01MCApplication * Ex01MCApplication::Instance ( )

inlinestatic

Returns

The MC application instance

Definition at line 81 of file Ex01MCApplication.h.

{
  /// \return The MC application instance
  return (Ex01MCApplication*)(TVirtualMCApplication::Instance());
}

InitMC()

void Ex01MCApplication::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 222 of file Ex01MCApplication.cxx.

{
  /// Initialize MC.
  /// 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.)");
    }
  }
 
  gMC->SetStack(fStack);
  gMC->SetMagField(fMagField);
  gMC->Init();
  gMC->BuildPhysics();
}

RunMC()

void Ex01MCApplication::RunMC

(

Int_t

nofEvents

)

Run MC.

Parameters

nofEvents

Number of events to be processed

Definition at line 244 of file Ex01MCApplication.cxx.

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

FinishRun()

void Ex01MCApplication::FinishRun

(

)

Finish MC run.

Definition at line 254 of file Ex01MCApplication.cxx.

{
  /// Finish MC run.
}

CloneForWorker()

TVirtualMCApplication * Ex01MCApplication::CloneForWorker ( ) const

virtual

Definition at line 260 of file Ex01MCApplication.cxx.

{
  return new Ex01MCApplication(GetName(), GetTitle());
}

InitOnWorker()

void Ex01MCApplication::InitOnWorker ( )

virtual

Definition at line 266 of file Ex01MCApplication.cxx.

{
  gMC->SetStack(fStack);
  gMC->SetMagField(fMagField);
}

ConstructGeometry()

void Ex01MCApplication::ConstructGeometry ( )

virtual

Construct geometry using TGeo functions or TVirtualMC functions (if oldGeometry is selected)

Definition at line 273 of file Ex01MCApplication.cxx.

{
  /// Construct geometry 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;
    ConstructMaterials();
    ConstructVolumes();
  }
  else {
    cout << "Geometry will be defined via VMC" << endl;
    Ex01DetectorConstructionOld detConstructionOld;
    detConstructionOld.ConstructMaterials();
    detConstructionOld.ConstructVolumes();
  }
}

InitGeometry()

void Ex01MCApplication::InitGeometry ( )

virtual

Initialize geometry.

Definition at line 300 of file Ex01MCApplication.cxx.

{
  /// Initialize geometry.
 
  fImedAr = gMC->MediumId("ArgonGas");
  fImedAl = gMC->MediumId("Aluminium");
  fImedPb = gMC->MediumId("Lead");
}

GeneratePrimaries()

void Ex01MCApplication::GeneratePrimaries ( )

virtual

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

Definition at line 310 of file Ex01MCApplication.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;
 
  // Geantino
  Int_t pdg = 0;
 
  // Polarization
  Double_t polx = 0.;
  Double_t poly = 0.;
  Double_t polz = 0.;
 
  // Position
  Double_t vx = -200.;
  Double_t vy = 0.;
  Double_t vz = 0.;
  Double_t tof = 0.;
 
  // Momentum
  Double_t px, py, pz, e;
  px = 1.;
  py = 0.;
  pz = 0.;
  e = 1.;
 
  // Add particle to stack
  fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof, polx,
    poly, polz, kPPrimary, ntr, 1., 0);
 
  // Change direction and add particle to stack
  px = 1.;
  py = 0.1;
  pz = 0.;
  fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof, polx,
    poly, polz, kPPrimary, ntr, 1., 0);
 
  // Change direction and add particle to stack
  px = 1.;
  py = 0.;
  pz = 0.1;
  fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof, polx,
    poly, polz, kPPrimary, ntr, 1., 0);
}

BeginEvent()

void Ex01MCApplication::BeginEvent ( )

virtual

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

Definition at line 361 of file Ex01MCApplication.cxx.

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

BeginPrimary()

void Ex01MCApplication::BeginPrimary ( )

virtual

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

Definition at line 368 of file Ex01MCApplication.cxx.

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

PreTrack()

void Ex01MCApplication::PreTrack ( )

virtual

User actions at beginning of each track. Print info message.

Definition at line 375 of file Ex01MCApplication.cxx.

{
  /// User actions at beginning of each track.
  /// Print info message.
 
  cout << endl;
  cout << "Starting new track" << endl;
}

Stepping()

void Ex01MCApplication::Stepping ( )

virtual

User actions at each step. Print track position, the current volume and current medium names.

Definition at line 385 of file Ex01MCApplication.cxx.

{
  /// User actions at each step.
  /// Print track position, the current volume and current medium names.
 
  TLorentzVector position;
  gMC->TrackPosition(position);
 
  cout << "Track " << position.X() << " " << position.Y() << " " << position.Z()
       << "  in " << gMC->CurrentVolName() << "  ";
 
  if (gMC->CurrentMedium() == fImedAr) cout << "ArgonGas";
  if (gMC->CurrentMedium() == fImedAl) cout << "Aluminium";
  if (gMC->CurrentMedium() == fImedPb) cout << "Lead";
 
  cout << endl;
 
  // // Test other TVirtualMC::TrackPosition() functions
 
  // Double_t dx, dy, dz;
  // gMC->TrackPosition(dx, dy, dz);
 
  // Float_t x, y, z;
  // gMC->TrackPosition(x, y, z);
 
  // cout << "Track position (double): "  << dx << " " << dy << " " << dz
  //      << "  (float): "   << x << " " << y << " " << z << endl;
}

PostTrack()

void Ex01MCApplication::PostTrack ( )

virtual

User actions after finishing of each track Nothing to be done this example

Definition at line 415 of file Ex01MCApplication.cxx.

{
  /// User actions after finishing of each track
  /// Nothing to be done this example
}

FinishPrimary()

void Ex01MCApplication::FinishPrimary ( )

virtual

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

Definition at line 422 of file Ex01MCApplication.cxx.

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

FinishEvent()

void Ex01MCApplication::FinishEvent ( )

virtual

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

Definition at line 429 of file Ex01MCApplication.cxx.

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

SetOldGeometry()

void Ex01MCApplication::SetOldGeometry ( Bool_t oldGeometry = kTRUE )

inline

Select old geometry definition (via TVirtualMC)

Definition at line 87 of file Ex01MCApplication.h.

{
  /// Select old geometry definition (via TVirtualMC)
  fOldGeometry = oldGeometry;
}

TestVMCGeometryGetters()

void Ex01MCApplication::TestVMCGeometryGetters

(

)

Test (new) TVirtualMC functions: GetTransform(), GetShape(), GetMaterial(), GetMedium()

Definition at line 436 of file Ex01MCApplication.cxx.

{
  /// Test (new) TVirtualMC functions:
  /// GetTransform(), GetShape(), GetMaterial(), GetMedium()
 
  // Get transformation of 10th layer
  //
  TString volPath = "/EXPH_1/CALB_1/LAYB_9";
  TGeoHMatrix matrix;
  Bool_t result = gMC->GetTransformation(volPath, matrix);
  if (result) {
    cout << "Transformation for " << volPath.Data() << ": " << endl;
    matrix.Print();
  }
  else {
    cerr << "Volume path " << volPath.Data() << " not found" << endl;
  }
  cout << endl;
 
  volPath = "/EXPH_1/CALB_1/LAYB_100";
  result = gMC->GetTransformation(volPath, matrix);
  if (result) {
    cout << "Transformation for " << volPath.Data() << ": " << endl;
    matrix.Print();
  }
  else {
    cerr << "Volume path " << volPath.Data() << " not found" << endl;
  }
  cout << endl;
 
  volPath = "/EXPH_1/CALB_1/LAYB_9";
  result = gMC->GetTransformation(volPath, matrix);
  if (result) {
    cout << "Transformation for " << volPath.Data() << ": " << endl;
    matrix.Print();
  }
  else {
    cerr << "Volume path " << volPath.Data() << " not found" << endl;
  }
  cout << endl;
 
  // Get shape
  //
  volPath = "/EXPH_1/CALB_1/LAYB_9";
  TString shapeType;
  TArrayD par;
  result = gMC->GetShape(volPath, shapeType, par);
  if (result) {
    cout << "Shape for " << volPath.Data() << ": " << endl;
    cout << shapeType.Data() << "  parameters: ";
    for (Int_t ipar = 0; ipar < par.GetSize(); ipar++)
      cout << par.At(ipar) << ",  ";
    cout << endl;
  }
  else {
    cerr << "Volume path " << volPath.Data() << " not found" << endl;
  }
  cout << endl;
 
  // Get material by material ID
  //
  TString matName;
  Int_t imat = 2;
  Double_t a, z, density, radl, inter;
  TArrayD mpar;
#if ROOT_VERSION_CODE >= ROOT_VERSION(5, 30, 0)
  result = gMC->GetMaterial(imat, matName, a, z, density, radl, inter, mpar);
  if (result) {
    cout << "Material with ID " << imat << ": " << endl;
    cout << matName.Data() << "  Aeff = " << a << "  Zeff = " << z
         << "  density = " << density << "  radl = " << radl
         << "  inter = " << inter << endl;
    if (mpar.GetSize() > 0) {
      cout << " User defined parameters: ";
      for (Int_t ipar = 0; ipar < par.GetSize(); ipar++)
        cout << mpar.At(ipar) << ",  ";
      cout << endl;
    }
  }
  else {
    cerr << "Material with ID " << imat << " not found" << endl;
  }
  cout << endl;
#endif
 
  // Get material by volume name
  //
  TString volName = "LAYB";
  mpar.Set(0);
  result =
    gMC->GetMaterial(volName, matName, imat, a, z, density, radl, inter, mpar);
  if (result) {
    cout << "Material for " << volName.Data() << " volume: " << endl;
    cout << matName.Data() << "  " << imat << "  Aeff = " << a
         << "  Zeff = " << z << "  density = " << density << "  radl = " << radl
         << "  inter = " << inter << endl;
    if (mpar.GetSize() > 0) {
      cout << " User defined parameters: ";
      for (Int_t ipar = 0; ipar < par.GetSize(); ipar++)
        cout << mpar.At(ipar) << ",  ";
      cout << endl;
    }
  }
  else {
    cerr << "Volume " << volName.Data() << " not found" << endl;
  }
  cout << endl;
 
  // Get medium
  //
  TString medName;
  Int_t imed, nmat, isvol, ifield;
  Double_t fieldm, tmaxfd, stemax, deemax, epsil, stmin;
  result = gMC->GetMedium(volName, medName, imed, nmat, isvol, ifield, fieldm,
    tmaxfd, stemax, deemax, epsil, stmin, mpar);
  if (result) {
    cout << "Medium for " << volName.Data() << " volume: " << endl;
    cout << medName.Data() << "  " << imed << "  nmat = " << nmat
         << "  isvol = " << isvol << "  ifield = " << ifield
         << "  fieldm = " << fieldm << "  tmaxfd = " << tmaxfd
         << "  stemax = " << stemax << "  deemax = " << deemax
         << "  epsil = " << epsil << "  stmin = " << stmin << endl;
    if (mpar.GetSize() > 0) {
      cout << " User defined parameters: ";
      for (Int_t ipar = 0; ipar < par.GetSize(); ipar++)
        cout << mpar.At(ipar) << ",  ";
      cout << endl;
    }
  }
  else {
    cerr << "Volume " << volName.Data() << " not found" << endl;
  }
  cout << endl;
 
  // Test getters non-existing position/volume
  //
 
  // Transformation
  volPath = "/EXPH_1/CALB_1/LAYB_100";
  result = gMC->GetTransformation(volPath, matrix);
  cout << "GetTransformation: Volume path " << volPath.Data();
  if (!result)
    cout << " not found" << endl;
  else
    cout << " found" << endl;
 
  // Shape
  result = gMC->GetShape(volPath, shapeType, par);
  cout << "GetShape: Volume path " << volPath.Data();
  if (!result)
    cout << " not found" << endl;
  else
    cout << " found" << endl;
 
  // Material
  volName = "XYZ";
  result =
    gMC->GetMaterial(volName, matName, imat, a, z, density, radl, inter, mpar);
  cout << "GetMaterial: Volume name " << volName.Data();
  if (!result)
    cout << " not found" << endl;
  else
    cout << " found" << endl;
 
  // Medium
  result = gMC->GetMedium(volName, medName, imed, nmat, isvol, ifield, fieldm,
    tmaxfd, stemax, deemax, epsil, stmin, mpar);
  cout << "GetMedium: Volume name " << volName.Data();
  if (!result)
    cout << " not found" << endl;
  else
    cout << " found" << endl;
}

ConstructMaterials()

void Ex01MCApplication::ConstructMaterials ( )

private

Construct materials using TGeo modeller

Definition at line 88 of file Ex01MCApplication.cxx.

{
  /// Construct materials using TGeo modeller
 
  // Create Root geometry manager
  new TGeoManager("E01_geometry", "E01 VMC example geometry");
 
  Double_t a;       // Mass of a mole in g/mole
  Double_t z;       // Atomic number
  Double_t density; // Material density in g/cm3
 
  a = 39.95;
  z = 18.;
  density = 1.782e-03;
  TGeoMaterial* matAr = new TGeoMaterial("ArgonGas", a, z, density);
 
  a = 26.98;
  z = 13.;
  density = 2.7;
  TGeoMaterial* matAl = new TGeoMaterial("Aluminium", a, z, density);
 
  a = 207.19;
  z = 82.;
  density = 11.35;
  TGeoMaterial* matLead = new TGeoMaterial("Lead", a, z, density);
 
  /*
    // Set material IDs
    // This step is needed, only if user wants to use the material Ids
    // in his application. Be aware that the material Ids vary
    // with each concrete MC.
    // It is recommended to use Media Ids instead, which values
    // set by user are preserved in all MCs
    Int_t imat = 0;
    matAr->SetUniqueID(imat++);
    matAl->SetUniqueID(imat++);
    matLead->SetUniqueID(imat++);
  */
 
  //
  // Tracking medias
  //
 
  Double_t param[20];
  param[0] = 0;     // isvol  - Not used
  param[1] = 2;     // ifield - User defined magnetic field
  param[2] = 10.;   // fieldm - Maximum field value (in kiloGauss)
  param[3] = -20.;  // tmaxfd - Maximum angle due to field deflection
  param[4] = -0.01; // stemax - Maximum displacement for multiple scat
  param[5] = -.3;   // deemax - Maximum fractional energy loss, DLS
  param[6] = .001;  // epsil - Tracking precision
  param[7] = -.8;   // stmin
  for (Int_t i = 8; i < 20; ++i) param[i] = 0.;
 
  fImedAr = 1;
  new TGeoMedium("ArgonGas", fImedAr, matAr, param);
 
  fImedAl = 2;
  new TGeoMedium("Aluminium", fImedAl, matAl, param);
 
  fImedPb = 3;
  new TGeoMedium("Lead", fImedPb, matLead, param);
}

ConstructVolumes()

void Ex01MCApplication::ConstructVolumes ( )

private

Contruct volumes using TGeo modeller

Definition at line 153 of file Ex01MCApplication.cxx.

{
  /// Contruct volumes using TGeo modeller
 
  //------------------------------ experimental hall (world volume)
  //------------------------------ beam line along x axis
 
  Double_t* ubuf = 0;
 
  Double_t expHall[3];
  expHall[0] = 300.;
  expHall[1] = 100.;
  expHall[2] = 100.;
  TGeoVolume* top = gGeoManager->Volume("EXPH", "BOX", fImedAr, expHall, 3);
  gGeoManager->SetTopVolume(top);
 
  //------------------------------ a tracker tube
 
  Double_t trackerTube[3];
  trackerTube[0] = 0.;
  trackerTube[1] = 60.;
  trackerTube[2] = 50.;
  gGeoManager->Volume("TRTU", "TUBE", fImedAl, trackerTube, 3);
 
  Double_t posX = -100.;
  Double_t posY = 0.;
  Double_t posZ = 0.;
  gGeoManager->Node("TRTU", 1, "EXPH", posX, posY, posZ, 0, kTRUE, ubuf);
 
  //------------------------------ a calorimeter block
 
  Double_t calBox[3];
  calBox[0] = 100.;
  calBox[1] = 50.;
  calBox[2] = 50.;
  gGeoManager->Volume("CALB", "BOX", fImedPb, calBox, 3);
 
  posX = 100.;
  posY = 0.;
  posZ = 0.;
  gGeoManager->Node("CALB", 1, "EXPH", posX, posY, posZ, 0, kTRUE, ubuf);
 
  //------------------------------ calorimeter layers
 
  Double_t layerBox[3];
  layerBox[0] = 1.;
  layerBox[1] = 40.;
  layerBox[2] = 40.;
  gGeoManager->Volume("LAYB", "BOX", fImedAl, layerBox, 3);
 
  for (Int_t i = 0; i < 19; i++) {
    posX = (i - 9) * 10.;
    posY = 0.;
    posZ = 0.;
    gGeoManager->Node("LAYB", i, "CALB", posX, posY, posZ, 0, kTRUE, ubuf);
  }
 
  // close geometry
  gGeoManager->CloseGeometry();
 
  // notify VMC about Root geometry
  gMC->SetRootGeometry();
}

Member Data Documentation

fStack

TVirtualMCStack* Ex01MCApplication::fStack

private

The VMC stack.

Definition at line 69 of file Ex01MCApplication.h.

fMagField

TVirtualMagField* Ex01MCApplication::fMagField

private

The magnetic field.

Definition at line 70 of file Ex01MCApplication.h.

fImedAr

Int_t Ex01MCApplication::fImedAr

private

The Argon gas medium Id.

Definition at line 71 of file Ex01MCApplication.h.

fImedAl

Int_t Ex01MCApplication::fImedAl

private

The Aluminium medium Id.

Definition at line 72 of file Ex01MCApplication.h.

fImedPb

Int_t Ex01MCApplication::fImedPb

private

The Lead medium Id.

Definition at line 73 of file Ex01MCApplication.h.

fOldGeometry

Bool_t Ex01MCApplication::fOldGeometry

private

Option for geometry definition.

Definition at line 74 of file Ex01MCApplication.h.

---

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

• examples/E01/include/Ex01MCApplication.h
• examples/E01/src/Ex01MCApplication.cxx