GENFIT
Rev:NoNumberAvailable
|
Stepper and energy loss/noise matrix calculation. More...
#include <MaterialEffects.h>
Public Member Functions | |
void | init (AbsMaterialInterface *matIfc) |
set the material interface here. Material interface classes must be derived from AbsMaterialInterface. More... | |
bool | isInitialized () |
void | setNoEffects (bool opt=true) |
void | setEnergyLossBetheBloch (bool opt=true) |
void | setNoiseBetheBloch (bool opt=true) |
void | setNoiseCoulomb (bool opt=true) |
void | setEnergyLossBrems (bool opt=true) |
void | setNoiseBrems (bool opt=true) |
void | ignoreBoundariesBetweenEqualMaterials (bool opt=true) |
void | setMscModel (const std::string &modelName) |
Select the multiple scattering model that will be used during track fit. More... | |
double | effects (const std::vector< RKStep > &steps, int materialsFXStart, int materialsFXStop, const double &mom, const int &pdg, M7x7 *noise=nullptr) |
Calculates energy loss in the traveled path, optional calculation of noise matrix. More... | |
void | stepper (const RKTrackRep *rep, M1x7 &state7, const double &mom, double &relMomLoss, const int &pdg, MaterialProperties ¤tMaterial, StepLimits &limits, bool varField=true) |
Returns maximum length so that a specified momentum loss will not be exceeded. More... | |
void | setDebugLvl (unsigned int lvl=1) |
void | drawdEdx (int pdg=11) |
Static Public Member Functions | |
static MaterialEffects * | getInstance () |
static void | destruct () |
Private Member Functions | |
MaterialEffects () | |
virtual | ~MaterialEffects () |
void | getParticleParameters () |
sets charge_, mass_ More... | |
void | getMomGammaBeta (double Energy, double &mom, double &gammaSquare, double &gamma, double &betaSquare) const |
double | momentumLoss (double stepSign, double mom, bool linear) |
Returns momentum loss. More... | |
double | dEdx (double Energy) const |
Calculate dEdx for a given energy. More... | |
double | dEdxBetheBloch (double betaSquare, double gamma, double gammasquare) const |
Uses Bethe Bloch formula to calculate dEdx. More... | |
void | noiseBetheBloch (M7x7 &noise, double mom, double betaSquare, double gamma, double gammaSquare) const |
calculation of energy loss straggeling More... | |
void | noiseCoulomb (M7x7 &noise, const M1x3 &direction, double momSquare, double betaSquare) const |
calculation of multiple scattering More... | |
double | dEdxBrems (double mom) const |
Returns dEdx. More... | |
void | noiseBrems (M7x7 &noise, double momSquare, double betaSquare) const |
calculation of energy loss straggeling More... | |
Private Attributes | |
bool | noEffects_ |
bool | energyLossBetheBloch_ |
bool | noiseBetheBloch_ |
bool | noiseCoulomb_ |
bool | energyLossBrems_ |
bool | noiseBrems_ |
bool | ignoreBoundariesBetweenEqualMaterials_ |
const double | me_ |
double | stepSize_ |
double | dEdx_ |
double | E_ |
double | matDensity_ |
double | matZ_ |
double | matA_ |
double | radiationLength_ |
double | mEE_ |
int | pdg_ |
int | charge_ |
double | mass_ |
int | mscModelCode_ |
AbsMaterialInterface * | materialInterface_ |
depending on this number a specific msc model is chosen in the noiseCoulomb function. More... | |
unsigned int | debugLvl_ |
Static Private Attributes | |
static MaterialEffects * | instance_ = nullptr |
Stepper and energy loss/noise matrix calculation.
It provides functionality to limit the stepsize of an extrapolation in order not to exceed a specified maximum momentum loss. After propagation, the energy loss for the given length and (optionally) the noise matrix can be calculated. You have to set which energy-loss and noise mechanisms you want to use. At the moment, per default all energy loss and noise options are ON.
Definition at line 50 of file MaterialEffects.h.
|
private |
Definition at line 42 of file MaterialEffects.cc.
Referenced by getInstance().
|
privatevirtual |
Definition at line 66 of file MaterialEffects.cc.
References materialInterface_.
|
private |
Calculate dEdx for a given energy.
Definition at line 437 of file MaterialEffects.cc.
References dEdxBetheBloch(), dEdxBrems(), energyLossBetheBloch_, energyLossBrems_, and getMomGammaBeta().
Referenced by drawdEdx(), and momentumLoss().
|
private |
Uses Bethe Bloch formula to calculate dEdx.
Definition at line 454 of file MaterialEffects.cc.
References charge_, mass_, matA_, matDensity_, matZ_, me_, mEE_, and genfit::Exception::setFatal().
Referenced by dEdx().
|
private |
Returns dEdx.
Can be called with any pdg, but only calculates dEdx for electrons and positrons (otherwise returns 0). Uses a gaussian approximation (Bethe-Heitler formula with Migdal corrections). For positrons, dEdx is weighed with a correction factor.
Definition at line 609 of file MaterialEffects.cc.
References matA_, matDensity_, matZ_, me_, and pdg_.
Referenced by dEdx().
|
static |
Definition at line 77 of file MaterialEffects.cc.
References instance_.
void genfit::MaterialEffects::drawdEdx | ( | int | pdg = 11 | ) |
Definition at line 802 of file MaterialEffects.cc.
References dEdx(), energyLossBetheBloch_, energyLossBrems_, genfit::AbsMaterialInterface::getMaterialParameters(), getParticleParameters(), genfit::AbsMaterialInterface::initTrack(), mass_, matA_, matDensity_, materialInterface_, matZ_, mEE_, pdg_, radiationLength_, and stepSize_.
double genfit::MaterialEffects::effects | ( | const std::vector< RKStep > & | steps, |
int | materialsFXStart, | ||
int | materialsFXStop, | ||
const double & | mom, | ||
const int & | pdg, | ||
M7x7 * | noise = nullptr |
||
) |
Calculates energy loss in the traveled path, optional calculation of noise matrix.
Definition at line 112 of file MaterialEffects.cc.
References debugLvl_, genfit::debugOut, E_, energyLossBetheBloch_, energyLossBrems_, getMomGammaBeta(), getParticleParameters(), matA_, matDensity_, materialInterface_, matZ_, mEE_, momentumLoss(), noEffects_, noiseBetheBloch(), noiseBetheBloch_, noiseBrems(), noiseBrems_, noiseCoulomb(), noiseCoulomb_, pdg_, radiationLength_, genfit::Exception::setFatal(), and stepSize_.
Referenced by genfit::RKTrackRep::Extrap().
|
static |
Definition at line 71 of file MaterialEffects.cc.
References instance_, and MaterialEffects().
Referenced by genfit::RKTrackRep::estimateStep(), and genfit::RKTrackRep::Extrap().
|
private |
Definition at line 358 of file MaterialEffects.cc.
References mass_, and genfit::Exception::setFatal().
|
private |
sets charge_, mass_
Definition at line 350 of file MaterialEffects.cc.
References charge_, mass_, and pdg_.
Referenced by drawdEdx(), effects(), and stepper().
|
inline |
Definition at line 76 of file MaterialEffects.h.
References ignoreBoundariesBetweenEqualMaterials_.
void genfit::MaterialEffects::init | ( | AbsMaterialInterface * | matIfc | ) |
set the material interface here. Material interface classes must be derived from AbsMaterialInterface.
Definition at line 85 of file MaterialEffects.cc.
References materialInterface_.
|
inline |
Definition at line 67 of file MaterialEffects.h.
References materialInterface_.
|
private |
|
private |
calculation of energy loss straggeling
For the energy loss straggeling, different formulas are used for different regions:
Needs dEdx_, which is calculated in momentumLoss, so it has to be called afterwards!
Definition at line 478 of file MaterialEffects.cc.
References charge_, dEdx_, mass_, matA_, matDensity_, matZ_, me_, and stepSize_.
Referenced by effects().
|
private |
calculation of energy loss straggeling
Can be called with any pdg, but only calculates straggeling for electrons and positrons.
Definition at line 775 of file MaterialEffects.cc.
References charge_, pdg_, radiationLength_, and stepSize_.
Referenced by effects().
|
private |
calculation of multiple scattering
This function first calcuates a MSC variance based on the current material and step length 2 different formulas for the MSC variance are implemeted. One can select the formula via "setMscModel". With the MSC variance and the current direction of the track a full 7D noise matrix is calculated. This noise matrix is the additional noise at the end of fStep in the 7D globa cooridnate system taking even the (co)variances of the position coordinates into account.
Definition at line 538 of file MaterialEffects.cc.
References genfit::RKMatrix< nRows, nCols >::begin(), charge_, genfit::RKMatrix< nRows, nCols >::end(), matZ_, mscModelCode_, radiationLength_, and stepSize_.
Referenced by effects().
void genfit::MaterialEffects::setDebugLvl | ( | unsigned int | lvl = 1 | ) |
Definition at line 795 of file MaterialEffects.cc.
References debugLvl_, materialInterface_, and genfit::AbsMaterialInterface::setDebugLvl().
|
inline |
Definition at line 71 of file MaterialEffects.h.
References energyLossBetheBloch_, and noEffects_.
|
inline |
Definition at line 74 of file MaterialEffects.h.
References energyLossBrems_, and noEffects_.
void genfit::MaterialEffects::setMscModel | ( | const std::string & | modelName | ) |
Select the multiple scattering model that will be used during track fit.
At the moment two model are available GEANE and Highland. GEANE is the model was was present in Genfit first. Note that using this function has no effect if setNoiseCoulomb(false) is set.
Definition at line 96 of file MaterialEffects.cc.
References genfit::errorOut, mscModelCode_, genfit::Exception::setFatal(), and genfit::Exception::what().
|
inline |
Definition at line 69 of file MaterialEffects.h.
References noEffects_.
|
inline |
Definition at line 72 of file MaterialEffects.h.
References noEffects_, and noiseBetheBloch_.
|
inline |
Definition at line 75 of file MaterialEffects.h.
References noEffects_, and noiseBrems_.
|
inline |
Definition at line 73 of file MaterialEffects.h.
References noEffects_, and noiseCoulomb_.
void genfit::MaterialEffects::stepper | ( | const RKTrackRep * | rep, |
M1x7 & | state7, | ||
const double & | mom, | ||
double & | relMomLoss, | ||
const int & | pdg, | ||
MaterialProperties & | currentMaterial, | ||
StepLimits & | limits, | ||
bool | varField = true |
||
) |
Returns maximum length so that a specified momentum loss will not be exceeded.
The stepper returns the maximum length that the particle may travel, so that a specified relative momentum loss will not be exceeded, or the next material boundary is reached. The material crossed are stored together with their stepsizes.
Definition at line 206 of file MaterialEffects.cc.
References debugLvl_, genfit::debugOut, genfit::AbsMaterialInterface::findNextBoundary(), genfit::StepLimits::getLowestLimitSignedVal(), genfit::StepLimits::getLowestLimitVal(), genfit::AbsMaterialInterface::getMaterialParameters(), getParticleParameters(), genfit::StepLimits::getStepSign(), ignoreBoundariesBetweenEqualMaterials_, genfit::AbsMaterialInterface::initTrack(), matA_, matDensity_, materialInterface_, matZ_, mEE_, momentumLoss(), noEffects_, pdg_, genfit::MaterialProperties::Print(), radiationLength_, genfit::RKTrackRep::RKPropagate(), genfit::Exception::setFatal(), genfit::StepLimits::setLimit(), genfit::MaterialProperties::setMaterialProperties(), stepSize_, genfit::stp_boundary, and genfit::stp_momLoss.
Referenced by genfit::RKTrackRep::estimateStep().
|
private |
Definition at line 194 of file MaterialEffects.h.
Referenced by dEdxBetheBloch(), getParticleParameters(), noiseBetheBloch(), noiseBrems(), and noiseCoulomb().
|
private |
Definition at line 201 of file MaterialEffects.h.
Referenced by effects(), momentumLoss(), setDebugLvl(), and stepper().
|
private |
Definition at line 185 of file MaterialEffects.h.
Referenced by momentumLoss(), and noiseBetheBloch().
|
private |
Definition at line 186 of file MaterialEffects.h.
Referenced by effects(), and momentumLoss().
|
private |
Definition at line 172 of file MaterialEffects.h.
Referenced by dEdx(), drawdEdx(), effects(), and setEnergyLossBetheBloch().
|
private |
Definition at line 175 of file MaterialEffects.h.
Referenced by dEdx(), drawdEdx(), effects(), and setEnergyLossBrems().
|
private |
Definition at line 178 of file MaterialEffects.h.
Referenced by ignoreBoundariesBetweenEqualMaterials(), and stepper().
|
staticprivate |
Definition at line 57 of file MaterialEffects.h.
Referenced by destruct(), and getInstance().
|
private |
Definition at line 195 of file MaterialEffects.h.
Referenced by dEdxBetheBloch(), drawdEdx(), getMomGammaBeta(), getParticleParameters(), momentumLoss(), and noiseBetheBloch().
|
private |
Definition at line 189 of file MaterialEffects.h.
Referenced by dEdxBetheBloch(), dEdxBrems(), drawdEdx(), effects(), noiseBetheBloch(), and stepper().
|
private |
Definition at line 187 of file MaterialEffects.h.
Referenced by dEdxBetheBloch(), dEdxBrems(), drawdEdx(), effects(), noiseBetheBloch(), and stepper().
|
private |
depending on this number a specific msc model is chosen in the noiseCoulomb function.
Definition at line 199 of file MaterialEffects.h.
Referenced by drawdEdx(), effects(), init(), isInitialized(), setDebugLvl(), stepper(), and ~MaterialEffects().
|
private |
Definition at line 188 of file MaterialEffects.h.
Referenced by dEdxBetheBloch(), dEdxBrems(), drawdEdx(), effects(), noiseBetheBloch(), noiseCoulomb(), and stepper().
|
private |
Definition at line 180 of file MaterialEffects.h.
Referenced by dEdxBetheBloch(), dEdxBrems(), and noiseBetheBloch().
|
private |
Definition at line 191 of file MaterialEffects.h.
Referenced by dEdxBetheBloch(), drawdEdx(), effects(), and stepper().
|
private |
Definition at line 197 of file MaterialEffects.h.
Referenced by noiseCoulomb(), and setMscModel().
|
private |
Definition at line 170 of file MaterialEffects.h.
Referenced by effects(), setEnergyLossBetheBloch(), setEnergyLossBrems(), setNoEffects(), setNoiseBetheBloch(), setNoiseBrems(), setNoiseCoulomb(), and stepper().
|
private |
Definition at line 173 of file MaterialEffects.h.
Referenced by effects(), and setNoiseBetheBloch().
|
private |
Definition at line 176 of file MaterialEffects.h.
Referenced by effects(), and setNoiseBrems().
|
private |
Definition at line 174 of file MaterialEffects.h.
Referenced by effects(), and setNoiseCoulomb().
|
private |
Definition at line 193 of file MaterialEffects.h.
Referenced by dEdxBrems(), drawdEdx(), effects(), getParticleParameters(), noiseBrems(), and stepper().
|
private |
Definition at line 190 of file MaterialEffects.h.
Referenced by drawdEdx(), effects(), noiseBrems(), noiseCoulomb(), and stepper().
|
private |
Definition at line 182 of file MaterialEffects.h.
Referenced by drawdEdx(), effects(), momentumLoss(), noiseBetheBloch(), noiseBrems(), noiseCoulomb(), and stepper().