VMC::Monopole::MCApplication Class Reference

Implementation of the TVirtualMCApplication. More...

#include <MCApplication.h>

Inheritance diagram for VMC::Monopole::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 ()
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	SetBinSize (Double_t binSize)
DetectorConstruction *	GetDetectorConstruction () const

Static Public Member Functions
static MCApplication *	Instance ()

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

Private Attributes
TMCRootManager *	fRootManager
	Root manager.
Ex03MCStack *	fStack
	The VMC stack.
DetectorConstruction *	fDetConstruction
	Dector construction.
TGeoUniformMagField *	fMagField
	Magnetic field.
Double_t	fBinSize
	Edep histogram bin size.
Double_t	fOffsetX
	The Edep histogram offset.
Double_t	fProjRange
	Projected range.
Double_t	fProjRange2
	Projected range square.
Int_t	fImedAl
	The Aluminium medium Id.
Int_t	fNofEvents
	Number of events.
Bool_t	fIsMaster
	If is on master thread.

Detailed Description

Implementation of the TVirtualMCApplication.

Date

05/04/2002

Author

I. Hrivnacova; IPN, Orsay

Definition at line 40 of file MCApplication.h.

Constructor & Destructor Documentation

MCApplication() [1/3]

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

Standard constructor

Parameters

name	The MC application name
title	The MC application description

Definition at line 57 of file MCApplication.cxx.

  : TVirtualMCApplication(name, title),
    fRootManager(0),
    fStack(0),
    fDetConstruction(0),
    fMagField(0),
    fBinSize(0.),
    fOffsetX(0.),
    fProjRange(0.),
    fProjRange2(0.),
    fImedAl(0),
    fNofEvents(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();
 
  // Constant magnetic field (in kiloGauss)
  fMagField = new TGeoUniformMagField(0, 0, 2);
}

MCApplication() [2/3]

VMC::Monopole::MCApplication::MCApplication

(

)

Default constructor

Definition at line 108 of file MCApplication.cxx.

  : TVirtualMCApplication(),
    fRootManager(0),
    fStack(0),
    fDetConstruction(0),
    fMagField(0),
    fBinSize(0.5),
    fImedAl(0),
    fNofEvents(0),
    fIsMaster(kTRUE)
{
  /// Default constructor
}

~MCApplication()

VMC::Monopole::MCApplication::~MCApplication ( )

virtual

Destructor

Definition at line 123 of file MCApplication.cxx.

{
  /// Destructor
 
  delete fStack;
  if (fIsMaster) delete fDetConstruction;
  delete fMagField;
  delete gMC;
}

MCApplication() [3/3]

VMC::Monopole::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 86 of file MCApplication.cxx.

  : TVirtualMCApplication(origin.GetName(), origin.GetTitle()),
    fRootManager(0),
    fStack(0),
    fDetConstruction(origin.fDetConstruction),
    fMagField(0),
    fImedAl(origin.fImedAl),
    fNofEvents(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);
 
  // Constant magnetic field (in kiloGauss)
  fMagField = new TGeoUniformMagField(origin.fMagField->GetFieldValue()[0],
    origin.fMagField->GetFieldValue()[1], origin.fMagField->GetFieldValue()[2]);
}

Member Function Documentation

Instance()

MCApplication * VMC::Monopole::MCApplication::Instance ( )

inlinestatic

Returns

The MC application instance

Definition at line 98 of file MCApplication.h.

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

InitMC()

void VMC::Monopole::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 153 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
 
  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();
 
  // Notify detector construction about initialization
  fDetConstruction->SetGeometryInitialized(true);
 
  // initialize projected range, tallies, Ebeam, and book histograms
  // Double_t maxStepSize = 0.5;
  fProjRange = 0.;
  fProjRange2 = 0.;
 
  Double_t length = fDetConstruction->GetAbsorberSizeX() * 10.; // in mm
  Double_t binSize = fBinSize * 10.;                            // in mm
  Int_t numBins = G4lrint(length / binSize);
 
  // Define offset for filling histogram
  fOffsetX = -0.5 * length; // in mm
 
  // cout << "numBins, length, offsetX "
  //      << numBins << ", " << length << ", " << fOffsetX << endl;
 
  // Create histograms
  // fHistograms.push_back(new TH1D("Edep_x", "Edep (MeV/mm) along absorber
  // (mm)", numBins, 0, length));
  fHistograms.push_back(
    new TH1D("h1", "Edep (MeV/mm) along absorber (mm)", numBins, 0, length));
  // fHistograms.push_back(new TH1D("h2", "DEDX (MeV/mm) of proton", 100,
  // -3., 7.)); fHistograms.push_back(new TH1D("h3", "DEDX (MeV/mm) of
  // monopole", 100, -3., 7.)); fHistograms.push_back(new TH1D("h4", "Range(mm)
  // of proton", 100, -3., 7.)); fHistograms.push_back(new TH1D("h5", "Range(mm)
  // of monopole", 100, -3., 7.)); Control histograms fHistograms.push_back(new
  // TH1D("x", "x", 100, -10, 10)); fHistograms.push_back(new TH1D("Edep (MeV)",
  // "Edep (Mev)", 100, 0, 100));
 
  RegisterStack();
}

RunMC()

void VMC::Monopole::MCApplication::RunMC

(

Int_t

nofEvents

)

Run MC.

Parameters

nofEvents

Number of events to be processed

Definition at line 221 of file MCApplication.cxx.

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

FinishRun()

void VMC::Monopole::MCApplication::FinishRun

(

)

Finish MC run.

Definition at line 232 of file MCApplication.cxx.

{
  /// Finish MC run.
 
  // Get monopole mass from TDatabasePDG
  TDatabasePDG* pdgDB = TDatabasePDG::Instance();
  TParticlePDG* monopole = pdgDB->GetParticle(60000000);
  Double_t mass = 0;
  if (monopole) {
    mass = monopole->Mass();
  }
 
  // run conditions
  cout << "\n The run consists of " << fNofEvents << " monopole "
       << " of " << mass << " GeV "
       << " through " << fDetConstruction->GetAbsorberSizeX() << " cm of "
       << fDetConstruction->GetAbsorberMaterial()
       // << " (density: "  << G4BestUnit(density,"Volumic Mass") << ")"
       << endl;
 
  // compute projected range and straggling
 
  fProjRange /= fNofEvents;
  fProjRange2 /= fNofEvents;
  Double_t rms = fProjRange2 - fProjRange * fProjRange;
  if (rms > 0.) {
    rms = std::sqrt(rms);
  }
  else {
    rms = 0.;
  }
 
  cout.precision(5);
  cout << "\n projected Range= " << fProjRange << " cm "
       << "   rms= " << rms << " cm " << endl;
 
  /*
    G4double ekin[100], dedxproton[100], dedxmp[100];
    G4EmCalculator calc;
    calc.SetVerbose(0);
    G4int i;
    for(i = 0; i < 100; ++i) {
      ekin[i] = std::pow(10., 0.1*G4double(i)) * keV;
      dedxproton[i] =
        calc.ComputeElectronicDEDX(ekin[i], "proton", matName);
      dedxmp[i] =
        calc.ComputeElectronicDEDX(ekin[i], "monopole", matName);
    }
 
    if(GetVerbose() > 0){
      G4cout << "### Stopping Powers" << G4endl;
      for(i=0; i<100; i++) {
        G4cout << " E(MeV)= " << ekin[i] << "  dedxp= " << dedxproton[i]
               << " dedxmp= " << dedxmp[i]
               << G4endl;
      }
    }
    G4cout << "### End of stopping power table" << G4endl;
 
    if(GetVerbose() > 0){
      G4cout << "Range table for " << matName << G4endl;
    }
 
    for(i=0; i<100; ++i) {
      G4double e = std::log10(ekin[i] / MeV) + 0.05;
      fHisto->Fill(1, e, dedxproton[i]);
      fHisto->Fill(2, e, dedxmp[i]);
      fHisto->Fill(3, e,
    std::log10(calc.GetRange(ekin[i],"proton",matName)/mm)); fHisto->Fill(4, e,
    std::log10(calc.GetRange(ekin[i],"monopole",matName)/mm));
    }
  */
  // normalize histogram
  // G4double fac = (mm/MeV) / (nEvents * fBinLength);
  Double_t fac = 1. / (fNofEvents * fBinSize * 10.);
  fHistograms[0]->Scale(fac);
 
  if (fRootManager) {
    fRootManager->WriteAll();
    fRootManager->Close();
  }
}

ConstructGeometry()

void VMC::Monopole::MCApplication::ConstructGeometry ( )

virtual

Construct geometry using detector contruction class. The detector contruction class is using TGeo functions.

Definition at line 350 of file MCApplication.cxx.

{
  /// Construct geometry using detector contruction class.
  /// The detector contruction class is using TGeo functions.
 
  fDetConstruction->ConstructGeometry();
}

InitGeometry()

void VMC::Monopole::MCApplication::InitGeometry ( )

virtual

Initialize geometry.

Definition at line 359 of file MCApplication.cxx.

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

GeneratePrimaries()

void VMC::Monopole::MCApplication::GeneratePrimaries ( )

virtual

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

Definition at line 367 of file MCApplication.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;
 
  // Monopole
  Int_t pdg = 60000000;
 
  // Polarization
  Double_t polx = 0.;
  Double_t poly = 0.;
  Double_t polz = 0.;
 
  // Position
  Double_t vx = -0.5 * (fDetConstruction->GetWorldSizeX()) + 1e-10; // + 1*um
  Double_t vy = 0.;
  Double_t vz = 0.;
  Double_t tof = 0.08;
 
  // Energy (in GeV)
  Double_t kinEnergy = 100.;
  Double_t mass = 100.;
  Double_t e = mass + kinEnergy;
 
  // Particle momentum
  Double_t px, py, pz;
  px = sqrt(e * e - mass * mass);
  py = 0.;
  pz = 0.;
 
  // 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 VMC::Monopole::MCApplication::BeginEvent ( )

virtual

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

Definition at line 408 of file MCApplication.cxx.

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

BeginPrimary()

void VMC::Monopole::MCApplication::BeginPrimary ( )

virtual

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

Definition at line 415 of file MCApplication.cxx.

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

PreTrack()

void VMC::Monopole::MCApplication::PreTrack ( )

virtual

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

Definition at line 422 of file MCApplication.cxx.

{
  /// User actions at beginning of each track.
  /// Print info message.
 
  // Check the monopole properties set via configuration macro
  // cout << endl;
  // std::cout << "&&&&& Track mass: " << gMC->TrackMass() << std::endl;
  // std::cout << "&&&&& Track charge: " << gMC->TrackCharge() << std::endl;
}

Stepping()

void VMC::Monopole::MCApplication::Stepping ( )

virtual

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

Definition at line 434 of file MCApplication.cxx.

{
  /// User actions at each step.
  /// Print track position, the current volume and current medium names.
 
  static TLorentzVector prevPosition;
 
  TLorentzVector position;
  gMC->TrackPosition(position);
 
  Double_t edep = gMC->Edep();
  if (edep <= 0.) {
    prevPosition = position;
    return;
  }
 
  edep *= 1.e03; // convert in MeV
 
  // Bragg curve
  Double_t x = position.X() * 10.;
  Double_t dx = (position.X() - prevPosition.X()) * 10.;
  // cout << "prevPosition "
  //      << prevPosition.X() << ", "  << prevPosition.Y() << ", "  <<
  //      prevPosition.Z() << endl;
  // cout << "position "
  //      << position.X() << ", "  << position.Y() << ", "  << position.Z() <<
  //      endl;
  // cout << "dx= " << dx << endl;
 
  Double_t random = gRandom->Uniform(0, 1);
  x += dx * random - fOffsetX;
 
  fHistograms[0]->Fill(x, edep);
  // fHistograms[1]->Fill(x);
  // fHistograms[2]->Fill(edep);
 
  // keep position from previous step
  prevPosition = position;
}

PostTrack()

void VMC::Monopole::MCApplication::PostTrack ( )

virtual

User actions after finishing of each track

Definition at line 475 of file MCApplication.cxx.

{
  /// User actions after finishing of each track
 
  // Skip secondary tracks
  if (fStack->GetCurrentParentTrackNumber() != -1) return;
 
  TLorentzVector position;
  gMC->TrackPosition(position);
 
  Double_t x = position.X() - fOffsetX / 10.; // in cm
  fProjRange += x;
  fProjRange2 += x * x;
}

FinishPrimary()

void VMC::Monopole::MCApplication::FinishPrimary ( )

virtual

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

Definition at line 491 of file MCApplication.cxx.

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

FinishEvent()

void VMC::Monopole::MCApplication::FinishEvent ( )

virtual

User actions after finishing of an event

Definition at line 498 of file MCApplication.cxx.

{
  /// User actions after finishing of an event
 
  // fRootManager->Fill();
  fStack->Reset();
}

SetBinSize()

void VMC::Monopole::MCApplication::SetBinSize

(

Double_t

binSize

)

Set Edep histogram bin size (in cm)

Definition at line 507 of file MCApplication.cxx.

{
  /// Set Edep histogram bin size (in cm)
 
  if (fDetConstruction->GetGeometryInitialized()) {
    cerr << "Geometry alredy initialized: cannot set Edep histogram bin size"
         << endl;
    return;
  }
 
  fBinSize = binSize;
}

GetDetectorConstruction()

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

inline

Returns

the detector construction

Definition at line 104 of file MCApplication.h.

{
  return fDetConstruction;
}

RegisterStack()

void VMC::Monopole::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);
  }
}

Member Data Documentation

fRootManager

TMCRootManager* VMC::Monopole::MCApplication::fRootManager

mutableprivate

Root manager.

Definition at line 80 of file MCApplication.h.

fStack

Ex03MCStack* VMC::Monopole::MCApplication::fStack

private

The VMC stack.

Definition at line 81 of file MCApplication.h.

fDetConstruction

DetectorConstruction* VMC::Monopole::MCApplication::fDetConstruction

private

Dector construction.

Definition at line 82 of file MCApplication.h.

fMagField

TGeoUniformMagField* VMC::Monopole::MCApplication::fMagField

private

Magnetic field.

Definition at line 83 of file MCApplication.h.

fBinSize

Double_t VMC::Monopole::MCApplication::fBinSize

private

Edep histogram bin size.

Definition at line 84 of file MCApplication.h.

fOffsetX

Double_t VMC::Monopole::MCApplication::fOffsetX

private

The Edep histogram offset.

Definition at line 85 of file MCApplication.h.

fProjRange

Double_t VMC::Monopole::MCApplication::fProjRange

private

Projected range.

Definition at line 86 of file MCApplication.h.

fProjRange2

Double_t VMC::Monopole::MCApplication::fProjRange2

private

Projected range square.

Definition at line 87 of file MCApplication.h.

fImedAl

Int_t VMC::Monopole::MCApplication::fImedAl

private

The Aluminium medium Id.

Definition at line 88 of file MCApplication.h.

fNofEvents

Int_t VMC::Monopole::MCApplication::fNofEvents

private

Number of events.

Definition at line 89 of file MCApplication.h.

fIsMaster

Bool_t VMC::Monopole::MCApplication::fIsMaster

private

If is on master thread.

Definition at line 90 of file MCApplication.h.

---

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

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