VMC::Gflash::MCApplication Class Reference

Implementation of the TVirtualMCApplication. More...

#include <MCApplication.h>

Inheritance diagram for VMC::Gflash::MCApplication:

Public Member Functions
	MCApplication (const char *name, const char *title)
	MCApplication ()
virtual	~MCApplication ()
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	GeneratePrimaries ()
virtual void	BeginEvent ()
virtual void	BeginPrimary ()
virtual void	PreTrack ()
virtual void	Stepping ()
virtual void	PostTrack ()
virtual void	FinishPrimary ()
virtual void	FinishEvent ()
void	SetVerboseLevel (Int_t verboseLevel)
DetectorConstruction *	GetDetectorConstruction () const
SensitiveDetector *	GetCalorimeterSD () const
PrimaryGenerator *	GetPrimaryGenerator () const

Private Member Functions
	MCApplication (const MCApplication &origin)
void	RegisterStack () const
void	ComputeEventStatistics () const

Private Attributes
TMCRootManager *	fRootManager
	Root manager.
Int_t	fEventNo
	Event counter.
TMCVerbose	fVerbose
	VMC verbose helper.
Ex03MCStack *	fStack
	VMC stack.
DetectorConstruction *	fDetConstruction
	Dector construction.
SensitiveDetector *	fSensitiveDetector
	Calorimeter SD.
PrimaryGenerator *	fPrimaryGenerator
	Primary generator.
Bool_t	fIsMaster
	If is on master thread.
TStopwatch *	fEventTimer
	Event timer.

Detailed Description

Implementation of the TVirtualMCApplication.

Geant4 gflash example adapted to Virtual Monte Carlo.

Date

28/10/2015

Author

I. Hrivnacova; IPN, Orsay

Definition at line 47 of file MCApplication.h.

Constructor & Destructor Documentation

MCApplication() [1/3]

VMC::Gflash::MCApplication::MCApplication	(	const char *	name,
		const char *	title )

Standard constructor

Parameters

name	The MC application name
title	The MC application description

Definition at line 48 of file MCApplication.cxx.

    : TVirtualMCApplication(name, title),
      fRootManager(0),
      fEventNo(0),
      fVerbose(0),
      fStack(0),
      fDetConstruction(0),
      fSensitiveDetector(0),
      fPrimaryGenerator(0),
      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 detector construction
    fDetConstruction = new DetectorConstruction();
 
    // Create a calorimeter SD
    fSensitiveDetector = new SensitiveDetector("Calorimeter");
 
    // Create a primary generator
    fPrimaryGenerator = new PrimaryGenerator(fStack);
  }

MCApplication() [2/3]

VMC::Gflash::MCApplication::MCApplication

(

)

Default constructor

Definition at line 103 of file MCApplication.cxx.

    : TVirtualMCApplication(),
      fRootManager(0),
      fEventNo(0),
      fStack(0),
      fDetConstruction(0),
      fSensitiveDetector(0),
      fPrimaryGenerator(0),
      fIsMaster(kTRUE)
  {
    /// Default constructor
  }

~MCApplication()

VMC::Gflash::MCApplication::~MCApplication ( )

virtual

Destructor

Definition at line 117 of file MCApplication.cxx.

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

MCApplication() [3/3]

VMC::Gflash::MCApplication::MCApplication ( const MCApplication & origin )

private

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

Parameters

origin

The source MC application

Definition at line 77 of file MCApplication.cxx.

    : TVirtualMCApplication(origin.GetName(), origin.GetTitle()),
      fRootManager(0),
      fEventNo(0),
      fVerbose(origin.fVerbose),
      fStack(0),
      fDetConstruction(origin.fDetConstruction),
      fSensitiveDetector(0),
      fPrimaryGenerator(0),
      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
    fSensitiveDetector = new SensitiveDetector(*(origin.fSensitiveDetector));
 
    // Create a primary generator
    fPrimaryGenerator =
      new PrimaryGenerator(*(origin.fPrimaryGenerator), fStack);
  }

Member Function Documentation

InitMC()

void VMC::Gflash::MCApplication::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 235 of file MCApplication.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->Init();
    gMC->BuildPhysics();
 
    RegisterStack();
  }

RunMC()

void VMC::Gflash::MCApplication::RunMC

(

Int_t

nofEvents

)

Run MC.

Parameters

nofEvents

Number of events to be processed

Definition at line 273 of file MCApplication.cxx.

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

FinishRun()

void VMC::Gflash::MCApplication::FinishRun

(

)

Finish MC run.

Definition at line 285 of file MCApplication.cxx.

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

ReadEvent()

void VMC::Gflash::MCApplication::ReadEvent

(

Int_t

i

)

Read i -th event and prints hits.

Parameters

i	The number of event to be read

Definition at line 329 of file MCApplication.cxx.

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

CloneForWorker()

TVirtualMCApplication * VMC::Gflash::MCApplication::CloneForWorker ( ) const

virtual

Definition at line 298 of file MCApplication.cxx.

  {
    return new MCApplication(*this);
  }

InitOnWorker()

void VMC::Gflash::MCApplication::InitOnWorker ( )

virtual

Definition at line 304 of file MCApplication.cxx.

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

FinishRunOnWorker()

void VMC::Gflash::MCApplication::FinishRunOnWorker ( )

virtual

Definition at line 319 of file MCApplication.cxx.

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

ConstructGeometry()

void VMC::Gflash::MCApplication::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 340 of file MCApplication.cxx.

  {
    /// Construct geometry using detector contruction class.
    /// The detector contruction class is using TGeo functions or
    /// TVirtualMC functions (if oldGeometry is selected)
 
    fVerbose.ConstructGeometry();
 
    fDetConstruction->Construct();
  }

InitGeometry()

void VMC::Gflash::MCApplication::InitGeometry ( )

virtual

Initialize geometry

Definition at line 352 of file MCApplication.cxx.

  {
    /// Initialize geometry
 
    fVerbose.InitGeometry();
 
    fSensitiveDetector->Initialize();
  }

GeneratePrimaries()

void VMC::Gflash::MCApplication::GeneratePrimaries ( )

virtual

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

Definition at line 362 of file MCApplication.cxx.

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

BeginEvent()

void VMC::Gflash::MCApplication::BeginEvent ( )

virtual

User actions at beginning of event

Definition at line 374 of file MCApplication.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++;
    cout << " Start generating event Nr " << fEventNo << endl;
 
    // start timer
    fEventTimer = new TStopwatch();
    fEventTimer->Start();
  }

BeginPrimary()

void VMC::Gflash::MCApplication::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 398 of file MCApplication.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();
  }

PreTrack()

void VMC::Gflash::MCApplication::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 408 of file MCApplication.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();
  }

Stepping()

void VMC::Gflash::MCApplication::Stepping ( )

virtual

User actions at each step

Definition at line 419 of file MCApplication.cxx.

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

PostTrack()

void VMC::Gflash::MCApplication::PostTrack ( )

virtual

User actions after finishing of each track

Definition at line 440 of file MCApplication.cxx.

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

FinishPrimary()

void VMC::Gflash::MCApplication::FinishPrimary ( )

virtual

User actions after finishing of a primary track

Definition at line 448 of file MCApplication.cxx.

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

FinishEvent()

void VMC::Gflash::MCApplication::FinishEvent ( )

virtual

User actions after finishing of an event

Definition at line 456 of file MCApplication.cxx.

  {
    /// User actions after finishing of an event
 
    // VMC
 
    fVerbose.FinishEvent();
 
    fRootManager->Fill();
    // The application code
 
    fEventTimer->Stop();
    cout << endl;
    cout << "******************************************";
    cout << endl;
    cout << "Elapsed Time: " << endl;
    fEventTimer->Print();
    cout << endl;
    cout << "******************************************" << endl;
    // TO DO
    // fDtime+=fTimerIntern.GetRealElapsed();
 
    ComputeEventStatistics();
 
    fSensitiveDetector->EndOfEvent();
 
    fStack->Reset();
  }

SetVerboseLevel()

void VMC::Gflash::MCApplication::SetVerboseLevel ( Int_t verboseLevel )

inline

Set verbosity

Parameters

verboseLevel

The new verbose level value

Definition at line 107 of file MCApplication.h.

{
  fVerbose.SetLevel(verboseLevel);
}

GetDetectorConstruction()

DetectorConstruction * VMC::Gflash::MCApplication::GetDetectorConstruction ( ) const

inline

Returns

The detector construction

Definition at line 113 of file MCApplication.h.

{
  return fDetConstruction;
}

GetCalorimeterSD()

SensitiveDetector * VMC::Gflash::MCApplication::GetCalorimeterSD ( ) const

inline

Returns

The calorimeter sensitive detector

Definition at line 119 of file MCApplication.h.

{
  return fSensitiveDetector;
}

GetPrimaryGenerator()

PrimaryGenerator * VMC::Gflash::MCApplication::GetPrimaryGenerator ( ) const

inline

Returns

The primary generator

Definition at line 125 of file MCApplication.h.

{
  return fPrimaryGenerator;
}

RegisterStack()

void VMC::Gflash::MCApplication::RegisterStack ( ) const

private

Register stack in the Root manager.

Definition at line 138 of file MCApplication.cxx.

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

ComputeEventStatistics()

void VMC::Gflash::MCApplication::ComputeEventStatistics ( ) const

private

Compute event statisics

Definition at line 149 of file MCApplication.cxx.

  {
    /// Compute event statisics
 
    cout << " ------ ExGflashEventAction::End of event nr. " << fEventNo
         << "  -----" << endl;
 
    TClonesArray* hitsCollection = fSensitiveDetector->GetHitsCollection();
 
    // Hits in sensitive Detector
    int n_hit = hitsCollection->GetEntriesFast();
    cout << "  " << n_hit << " hits are stored in HitsCollection " << endl;
 
    // Get (x,y,z) of vertex of initial particles
    TVector3 vertexPosition = fPrimaryGenerator->GetVertexPosition();
    TVector3 vertexDirection = fPrimaryGenerator->GetVertexDirection();
 
    // ExGflashEventAction: Magicnumber
    // Should be retrieved from detector construction
    Double_t energyincrystal[100];
    for (Int_t i = 0; i < 100; i++) energyincrystal[i] = 0.;
 
    // For all Hits in sensitive detector
    Double_t totE = 0;
    for (Int_t i = 0; i < n_hit; i++) {
      Hit* hit = static_cast<Hit*>(hitsCollection->At(i));
 
      Double_t estep = hit->GetEdep();
      if (estep > 0) {
        totE += estep;
        Int_t num = hit->GetCrystalNum();
        energyincrystal[num] += estep;
        // cout << " Crystal Nummer " <<  num  << endl;
 
        TVector3 hitpos = hit->GetPos();
        TVector3 l(hitpos);
        // distance from shower start
        l = l - vertexPosition;
        // projection on shower axis = longitudinal profile
        TVector3 longitudinal(l);
        // shower profiles (Radial)
        TVector3 radial(vertexPosition.Cross(l));
      }
    }
 
    // Find crystal with maximum energy
    Double_t max = 0;
    Int_t index = 0;
    for (Int_t i = 0; i < 100; i++) {
      // cout << i <<"  " << energyincrystal[i] << G4endl;
      if (max < energyincrystal[i]) {
        max = energyincrystal[i];
        index = i;
      }
    }
    // cout << " NMAX  " << index << G4endl;
 
    // 3x3 matrix of crystals around the crystal with the maximum energy deposit
    Double_t e3x3 = energyincrystal[index] + energyincrystal[index + 1] +
                    energyincrystal[index - 1] + energyincrystal[index - 10] +
                    energyincrystal[index - 9] + energyincrystal[index - 11] +
                    energyincrystal[index + 10] + energyincrystal[index + 11] +
                    energyincrystal[index + 9];
 
    // 5x5 matrix of crystals around the crystal with the maximum energy deposit
    Double_t e5x5 = energyincrystal[index] + energyincrystal[index + 1] +
                    energyincrystal[index - 1] + energyincrystal[index + 2] +
                    energyincrystal[index - 2] + energyincrystal[index - 10] +
                    energyincrystal[index - 9] + energyincrystal[index - 11] +
                    energyincrystal[index - 8] + energyincrystal[index - 12] +
                    energyincrystal[index + 10] + energyincrystal[index + 11] +
                    energyincrystal[index + 9] + energyincrystal[index + 12] +
                    energyincrystal[index + 8];
 
    cout << "   e1  " << energyincrystal[index] << "   e3x3  " << e3x3
         << "   e5x5  " << e5x5 << " GeV" << endl;
 
    cout << " Total energy deposited in the calorimeter: " << totE << " (GeV)"
         << endl;
  }

Member Data Documentation

fRootManager

TMCRootManager* VMC::Gflash::MCApplication::fRootManager

mutableprivate

Root manager.

Definition at line 90 of file MCApplication.h.

fEventNo

Int_t VMC::Gflash::MCApplication::fEventNo

private

Event counter.

Definition at line 91 of file MCApplication.h.

fVerbose

TMCVerbose VMC::Gflash::MCApplication::fVerbose

private

VMC verbose helper.

Definition at line 92 of file MCApplication.h.

fStack

Ex03MCStack* VMC::Gflash::MCApplication::fStack

private

VMC stack.

Definition at line 93 of file MCApplication.h.

fDetConstruction

DetectorConstruction* VMC::Gflash::MCApplication::fDetConstruction

private

Dector construction.

Definition at line 94 of file MCApplication.h.

fSensitiveDetector

SensitiveDetector* VMC::Gflash::MCApplication::fSensitiveDetector

private

Calorimeter SD.

Definition at line 95 of file MCApplication.h.

fPrimaryGenerator

PrimaryGenerator* VMC::Gflash::MCApplication::fPrimaryGenerator

private

Primary generator.

Definition at line 96 of file MCApplication.h.

fIsMaster

Bool_t VMC::Gflash::MCApplication::fIsMaster

private

If is on master thread.

Definition at line 97 of file MCApplication.h.

fEventTimer

TStopwatch* VMC::Gflash::MCApplication::fEventTimer

private

Event timer.

Definition at line 98 of file MCApplication.h.

---

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

• examples/Gflash/include/MCApplication.h
• examples/Gflash/src/MCApplication.cxx