PrimaryGenerator.cxx

Go to the documentation of this file.

//------------------------------------------------
// The Virtual Monte Carlo examples
// Copyright (C) 2007 - 2016 Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see geant4_vmc/LICENSE.
// Contact: root-vmc@cern.ch
//-------------------------------------------------
 
/// \file ExGarfield/src/PrimaryGenerator.cxx
/// \brief Implementation of the ExGarfield::PrimaryGenerator class
///
/// Garfield garfieldpp example adapted to Virtual Monte Carlo.
///
/// \date 28/10/2015
/// \author I. Hrivnacova; IPN, Orsay
 
#include <TDatabasePDG.h>
#include <TPDGCode.h>
#include <TParticlePDG.h>
#include <TRandom.h>
#include <TVector3.h>
#include <TVirtualMC.h>
#include <TVirtualMCApplication.h>
#include <TVirtualMCStack.h>
 
#include "PrimaryGenerator.h"
 
/// \cond CLASSIMP
ClassImp(VMC::ExGarfield::PrimaryGenerator)
  /// \endcond
 
  namespace VMC
{
  namespace ExGarfield
  {
 
  //_____________________________________________________________________________
  PrimaryGenerator::PrimaryGenerator(TVirtualMCStack* stack)
    : TObject(),
      fStack(stack),
      fNofPrimaries(1),
      fVertexPosition(),
      fVertexDirection(TVector3(0, 0, 1))
  {
    /// Standard constructor
    /// \param stack  The VMC stack
  }
 
  //_____________________________________________________________________________
  PrimaryGenerator::PrimaryGenerator(
    const PrimaryGenerator& origin, TVirtualMCStack* stack)
    : TObject(origin),
      fStack(stack),
      fNofPrimaries(origin.fNofPrimaries),
      fVertexPosition(origin.fVertexPosition),
      fVertexDirection(origin.fVertexDirection)
  {
    /// Copy constructor (for clonig on worker thread in MT mode).
    /// \param origin    The source object (on master).
    /// \param stack  The VMC stack
  }
 
  //_____________________________________________________________________________
  PrimaryGenerator::PrimaryGenerator()
    : TObject(),
      fStack(0),
      fNofPrimaries(0),
      fVertexPosition(),
      fVertexDirection()
  {
    /// Default constructor
  }
 
  //_____________________________________________________________________________
  PrimaryGenerator::~PrimaryGenerator()
  {
    /// Destructor
  }
 
  //
  // private methods
  //
 
  //_____________________________________________________________________________
  void PrimaryGenerator::GenerateOnePrimary(const TVector3& origin)
  {
    /// Add one primary particle (kElectron) to the user stack
    /// (derived from TVirtualMCStack).
    /// \param origin  The track position
 
    // Track ID (filled by stack)
    Int_t ntr;
 
    // Option: to be tracked
    Int_t toBeDone = 1;
 
    // PDG
    Int_t pdg = kMuonPlus;
 
    // Polarization
    Double_t polx = 0.;
    Double_t poly = 0.;
    Double_t polz = 0.;
 
    // Position
    Double_t vx = origin.X();
    Double_t vy = origin.Y();
    Double_t vz = origin.Z();
    Double_t tof = 0.;
    fVertexPosition = origin;
 
    // Energy (in GeV)
    Double_t kinEnergy = 4.;
    Double_t mass = 0.1056583715;
    Double_t e = mass + kinEnergy;
    Double_t pmag = sqrt(e * e - mass * mass);
 
    // Particle momentum
    Double_t px, py, pz;
    px = pmag * fVertexDirection.X();
    py = pmag * fVertexDirection.Y();
    pz = pmag * fVertexDirection.Z();
 
    // Randomize position
    // if (fIsRandom) {
    //   vy = origin.Y()*(gRandom->Rndm() - 0.5);
    //   vz = origin.Z()*(gRandom->Rndm() - 0.5);
    // }
 
    // Add particle to stack
    fStack->PushTrack(toBeDone, -1, pdg, px, py, pz, e, vx, vy, vz, tof, polx,
      poly, polz, kPPrimary, ntr, 1., 0);
  }
 
  //
  // public methods
  //
 
  //_____________________________________________________________________________
  void PrimaryGenerator::GeneratePrimaries(const TVector3& origin)
  {
    /// Fill the user stack (derived from TVirtualMCStack) with primary
    /// particles.
 
    for (Int_t i = 0; i < fNofPrimaries; i++) GenerateOnePrimary(origin);
  }
 
  } // namespace ExGarfield
}