cmc C++ API Reference¶

struct beta_dist¶

#include <I3CascadeSimulation.h>

Local definition of beta distribution.

The struct defines a pdf function, ie. the beta distribtuion not normed, though. It provides a rvs method to draw random samples from the beta distribution using gsl_rng

The parameters used for the beta distribibution (alpha and beta) are defined.

Functions are passed to I3MetropolisHastings and need to be static.

Public Functions

SET_LOGGER ("I3CascadeSimulation")

Public Static Functions

static inline double pdf(const double &x)¶

Beta Function $ x^{\alpha-1} (1-x)^{\beta-1}$

Parameters:: x –

static inline double rvs(const gsl_rng *random)¶

Returns a random variate drawn from the beta distribution

Uses the gsl_ran_beta method

Parameters:: random – gsl_rng instance

Public Static Attributes

static double alpha¶

static double beta¶

class I3CascadeDevelopment¶

#include <I3CascadeDevelopment.h>

This abstract class defines the interface of a cascade development description.

Cascade longitudinal base class for parametrization and simulation classes

Rcs: I3CascadeDevelopment.h 32874 2007-06-04 18:22:36Z bvoigt

An implementation of this class would define a Simulate method that calculates the length and the longitudinal energy deposit profile of a cascade of given energy.

Version
Rcs: 32874
Date
Rcs: 2007-06-04 12:22:36 -0600 (Mon, 04 Jun 2007)
Author: Bernhard Voigt bernhard.voigt@desy.de Last changed by:
Rcs: bvoigt

After this method has been invoked, the length and longitudinal profile are accessible via the corresponding methods of this class.

Implementations can be either a pure parametrization or a simulation of the shower development.

Author: Bernhard Voigt

Subclassed by I3CascadeParametrization, I3CascadeSimulation

Public Functions

inline virtual ~I3CascadeDevelopment()¶: Virtual destructor

virtual void Simulate(I3Particle::ParticleType type, double energy) = 0¶

Simulates the cascade development

This method must be called in order to get the length and/or longitudinal profile.

Parameters:

type – EM or hadronic
energy – energy of the shower

double GetLength()¶

Returns the length of a shower units of radiation length

I3CascadeDevelopment::Simulate must be called previously in order to get the shower length

Returns:: the length of the shower in units of radiation length

double GetEnergyDeposit(int a, int b)¶

Returns the energy loss of the shower between a and b (units are radiation length)

I3CascadeDevelopment::Simulate must be called beforehand.

Parameters:

a – start position (units are radiation length)
b – end position (units are radiation length)

Returns:

fractional energy loss between a and b (excluding b)

inline std::vector<double> GetEnergyLossProfile()¶: Returns the energy loss profile filled in the simulate method.

Public Static Attributes

static const double RADIATION_LENGTH = 0.358 / 0.9216¶

radiation length in ice in meter

Implementation I3CascadeDevelopment class.

$Id$

Version: $Revision$
Date: $LastChangedDate$
Author: Bernhard Voigt bernhard.voigt@desy.de Last changed by: $LastChangedBy$

Protected Attributes

std::vector<double> energyLossProfile_¶: pointer to the vector holding the energy loss profile at every radiation length filled in Simulate method.

Private Functions

SET_LOGGER ("I3CascadeDevelopment")

class I3CascadeMCService : public I3PropagatorService ¶

#include <I3CascadeMCService.h>

Cascade Monte Carlo Service.

This module provides the functionality to replace cascade-like in the MCTree with with a list of sub-cascades to simulate the cascade longitudinal development.

Public Functions

inline ~I3CascadeMCService()¶

I3CascadeMCService(I3RandomServicePtr r)¶: Constructor

virtual std::vector<I3Particle> Propagate(I3Particle &p, DiagnosticMapPtr frame, I3FramePtr)¶

virtual void SetRandomNumberGenerator(I3RandomServicePtr random)¶

void Simulate(I3Particle &c, std::vector<I3Particle> &d)¶

void SetEnergyThresholdMuons(double t)¶

void SetMaxMuons(int i)¶

void SetThresholdSplit(double t)¶

void SetSegmentMaxEnergy(double t)¶

void SetEnergyThresholdSimulation(double t)¶

void SetStepWidth(int w)¶

Private Functions

SET_LOGGER ("I3CascadeMCService")

Private Members

I3RandomServicePtr random_¶: random service provided by the frame

double energyThresholdMuons_¶: hadronic particles with an energy exceeding this threshold will generate muons.

int maxMuons_¶: maximal number of generated muons.

double energyThresholdSplit_¶: particles with an energy exceeding this threshold will be split

double energyThresholdSimulation_¶: particles with an energy exceeding this threshold will be split and the the longitudinal development is simulated

int stepWidth_¶: distance between individual sub-cascades in units of radiation length

boost::shared_ptr<I3CascadeSplit> splitter_¶: cascade splitting class

boost::shared_ptr<I3CascadeMuonSplit> muonSplitter_¶: cascade moun splitting class

class I3CascadeMuonSplit : public I3CascadeMCCommon::Sampler ¶

#include <I3CascadeMuonSplit.h>

This class defines methods to generate muons from a cascade.

Generation of muons for a realistic treatment of the longitudinal development of cascasdes

Rcs: I3CascadeSplit.h 33096 2007-06-11 15:09:40Z bvoigt

There is only one method of interest to the user which is the split method. It takes a I3Particle of a cascade typ, generates a set of muons, reduces the energy of the cascade and returns the corresponding I3Particles for the muons.

Version
Rcs: 33096
Date
Rcs: 2007-06-11 17:09:40 +0200 (Mon, 11 Jun 2007)
Author: Sebastian Panknin panknin@physik.hu-berlin.de Last changed by:
Rcs: bvoigt

Author: Sebastian Panknin

Public Functions

I3CascadeMuonSplit(I3RandomServicePtr random)¶

Creates a I3CascadeMuonSplit object

Parameters:: random – service pointer

virtual ~I3CascadeMuonSplit()¶: Default destructor

std::vector<I3Particle> GenerateMuons(I3Particle &cascade)¶

Generate Muons for a given I3Particle.

For a (hadronic) cascade-like particle muons according to the cascade energy are generated. The number and energies of the muons is taken from a parameterization. The muons have the same origian and the same direction as the orgriginal cascade. The are positivly charged (I3Particle::MuPlus). The energy of the cascade is reduced by the energy of the muons.

Parameters:: cascade – I3Particle
Returns:: vector of I3Particles

void SetEnergyThresholdMuons(double threshold)¶: set the energyThresholdMuons

void SetMaxMuons(int max)¶: set the maxMuons number

Public Static Attributes

static const double DEFAULT_ENERGY_THRESHOLD = 1 * I3Units::GeV ¶

static const int DEFAULT_MAX_MUONS = 10¶

Private Functions

SET_LOGGER ("I3CascadeMuonSplit")

I3CascadeMuonSplit(const I3CascadeMuonSplit &cascadeMuonSplit)¶

Private Members

double energyThresholdMuons_¶: energy treshold for muons to consider

int maxMuons_¶: maximal number of muons to consider

class I3CascadeParametrization : public I3CascadeDevelopment, public I3CascadeMCCommon::Sampler ¶

#include <I3CascadeParametrization.h>

This class implements the longitudinal shower profile of electromagnetic cascades for energies below the LPM Threshold (roughly 1PeV in Ice)

Cascade longitudinal energy profile parametrization

Rcs: I3CascadeParametrization.h 32874 2007-06-04 18:22:36Z bvoigt

The parametrization is take from the Particle Data Booklet, Section 26.5 The parameters a,b of the dE/dt function described their have been determined by Christoher Wiebusch using GEANT 3. They are taken from his thesis.

Version
Rcs: 32874
Date
Rcs: 2007-06-04 12:22:36 -0600 (Mon, 04 Jun 2007)
Author: Bernhard Voigt bernhard.voigt@desy.de Last changed by:
Rcs: bvoigt

Author: Bernhard Voigt

Public Functions

I3CascadeParametrization()¶: Default constructor

inline ~I3CascadeParametrization()¶: Default destructor

virtual void Simulate(I3Particle::ParticleType type, double energy)¶

Evaluates the longitudinal shower profile parametrization for an EM shower of given energy

The parametrization is taken from PDG Booklet Section (Particle Passage through matter)

$\frac{dE}{dt} = E\ b\ \frac{(bt)^{a-1}\ e^{-bt}}{\Gamma(a)}$

The parameters used are a, b which have been determined by Christopher Wiebusch and quoted in his thesis.

The energy deposit is calculated in steps of one radiation length. It is internally stored and can be accessed using the I3CascadeDevelopment::GetEnergyDeposit method.

Parameters:

type – EM or hadronic
energy – energy of the shower

double dE_dt(double energy, double t)¶

Calculates the energy loss function dE/dt

Parameters:

energy – of the shower
t – is the depth in radiation length units

inline void SetThreshold(double threshold)¶

Sets the threshold for the parametrization

Calculation of energy deposit is quit, when the deposit in one radiation length step is less than this threshold, default is 1 TeV

Parameters:: threshold –

inline double GetThreshold()¶: Returns the threshold of energy deposit calculation

Public Static Attributes

static const double DEFAULT_THRESHOLD = 1 * I3Units::TeV ¶

threshold for energy deposit calculation

Implementation I3CascadeParametrization class.

$Id$

Version: $Revision$
Date: $LastChangedDate$
Author: Bernhard Voigt bernhard.voigt@desy.de Last changed by: $LastChangedBy$

Private Functions

SET_LOGGER ("I3CascadeParametrization")

Private Members

double threshold_¶

Threshold for calculation energy deposit.

Calculation of energy deposit is quit, when the deposit in one radiation length step is less than this threshold

Private Static Attributes

static const double A_0¶

parameter for ‘a’ parameter calculation, no units.

From Wiebusch thesis

static const double A_1¶

static const double B¶

paramter for dE/dt function

From Wiebusch thesis

class I3CascadeSimulation : public I3CascadeDevelopment, public I3CascadeMCCommon::Sampler ¶

#include <I3CascadeSimulation.h>

This class implements a Cascade development simulation including the LPM effect.

Description of the simulation can be found here: http://wiki.icecube.wisc.edu/index.php/Cascade_Simulation

Cross section formulars used are defined in I3CascadeSimulationCrossSections

Author: Bernhard Voigt

Public Functions

I3CascadeSimulation(I3RandomServicePtr random)¶

Constructor.

Constructor with random service given

same as above, except random is not created

Parameters:: random – random number service from icetray frame

~I3CascadeSimulation()¶

Destructor.

Frees the gsl objects

Destructor - frees all gsl struct that have been initialized

virtual void SetRandomNumberGenerator(I3RandomServicePtr)¶

virtual void Simulate(I3Particle::ParticleType type, double energy)¶

Simulates a cascade of the given energy

Parameters:

type – EM or hadronic
energy –

inline void SetThreshold(double thresholdEnergy)¶

Sets the energy threshold down to which particles are tracked (Default is 50 TeV)

This shouldn’t be lower than 1 GeV, since this is the limit for the Bremsstrahlung interaction

Parameters:: thresholdEnergy –

inline double GetThreshold()¶

Returns the threshold energy down to which particles are tracked

Returns:: threshold energy

Public Static Attributes

static const double DEFAULT_THRESHOLD = 50 * I3Units::TeV ¶: default energy threshold down to which particles are tracked

Private Functions

SET_LOGGER ("I3CascadeSimulation")

I3CascadeSimulation(const I3CascadeSimulation&)¶

double PairProductionMeanFreePath(double energy)¶

Cacluates the pair productio mean free path for the given energy

Parameters:: energy –

double BremsstrahlungMeanFreePath(double energy)¶

Cacluates the bremsstrahlung radiation length (cut off energy is defined as 1 GeV)

Parameters:: energy –

double SampleBremsstrahlungCrossSection(double energy)¶

Samples the bremsstrahlung cross section to get a fractional energy of a secondary particle produced in an interaction with a incident particle of the given energy

cmc C++ API Reference¶

Table of Contents

Previous topic

Next topic

This Page