d5/ddc/GFMaterialEffects_8h_source.html

/* Copyright 2008-2009, Technische Universitaet Muenchen,

   Authors: Christian Hoeppner & Sebastian Neubert


   This file is part of GENFIT.


   GENFIT is free software: you can redistribute it and/or modify

   it under the terms of the GNU Lesser General Public License as published

   by the Free Software Foundation, either version 3 of the License, or

   (at your option) any later version.


   GENFIT is distributed in the hope that it will be useful,

   but WITHOUT ANY WARRANTY; without even the implied warranty of

   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

   GNU Lesser General Public License for more details.


   You should have received a copy of the GNU Lesser General Public License

   along with GENFIT.  If not, see <http://www.gnu.org/licenses/>.

*/


#ifndef GFMATERIALEFFECTS_H

#define GFMATERIALEFFECTS_H


#include<iostream>

#include "RKTools.h"

#include "TObject.h"

#include <vector>

#include "TVector3.h"

#include "TMatrixT.h"

#include "GFPointPath.h"


class GFMaterialEffects : public TObject {

private:

  GFMaterialEffects();

  virtual ~GFMaterialEffects();


  static GFMaterialEffects* finstance;


public:

  static GFMaterialEffects* getInstance();

  static void destruct();


  void setNoEffects(bool opt = true) {fNoEffects = opt;}


  void setEnergyLossBetheBloch(bool opt = true) {fEnergyLossBetheBloch = opt; fNoEffects = false;}

  void setNoiseBetheBloch(bool opt = true) {fNoiseBetheBloch = opt; fNoEffects = false;}

  void setNoiseCoulomb(bool opt = true) {fNoiseCoulomb = opt; fNoEffects = false;}

  void setEnergyLossBrems(bool opt = true) {fEnergyLossBrems = opt; fNoEffects = false;}

  void setNoiseBrems(bool opt = true) {fNoiseBrems = opt; fNoEffects = false;}


  void setMscModel(const std::string& modelName);


  double effects(const std::vector<GFPointPath>& points,

                 const double& mom,

                 const int& pdg,

                 double& xx0,

                 const bool& doNoise = false,

                 double* noise7x7 = NULL,

                 const double* jacobian7x7 = NULL,

                 const TVector3* directionBefore = NULL,

                 const TVector3* directionAfter = NULL);


  double stepper(const double& maxStep, // maximum step. unsigned!

                 const double& maxAngleStep, // maximum step due to curvature. unsigned!

                 const double& posx,

                 const double& posy,

                 const double& posz,

                 const double& dirx,

                 const double& diry,

                 const double& dirz,

                 const double& mom, // momentum

                 double& relMomLoss, // relative momloss for the step will be added

                 const int& pdg);


private:

  void getParticleParameters(double mom);


  void getMaterialParameters(TGeoMaterial* mat);


  double energyLossBetheBloch(const double& mom);


  void noiseBetheBloch(const double& mom,

                       double* noise) const;


  void noiseCoulomb(const double& mom,

                    double* noise,

                    const double* jacobian,

                    const TVector3* directionBefore,

                    const TVector3* directionAfter) const;


  double energyLossBrems(const double& mom) const;


  void noiseBrems(const double& mom,

                  double* noise) const;


  bool fNoEffects;


  bool fEnergyLossBetheBloch;

  bool fNoiseBetheBloch;

  bool fNoiseCoulomb;

  bool fEnergyLossBrems;

  bool fNoiseBrems;


  const double me; // electron mass (GeV)


  double fstep; // stepsize


  // cached values for energy loss and noise calculations

  double fbeta;

  double fdedx;

  double fgamma;

  double fgammaSquare;


  double fmatDensity;

  double fmatZ;

  double fmatA;

  double fradiationLength;

  double fmEE; // mean excitation energy


  int fpdg;

  double fcharge;

  double fmass;


  int fMscModelCode;


public:

  ClassDef(GFMaterialEffects, 2);


};


#endif