GFMaterialEffects.h

Go to the documentation of this file.

/* 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 "TObject.h"
#include <vector>
#include "TVector3.h"
#include "TMatrixT.h"

class GFMaterialEffects : public TObject {
 private:
  GFMaterialEffects();
  virtual ~GFMaterialEffects();

  static GFMaterialEffects *finstance;

 public:
  static GFMaterialEffects *getInstance();
  static void destruct();

  void setEnergyLossBetheBloch(bool opt = true) { fEnergyLossBetheBloch = opt; }
  void setNoiseBetheBloch(bool opt = true) { fNoiseBetheBloch = opt; }
  void setNoiseCoulomb(bool opt = true) { fNoiseCoulomb = opt; }
  void setEnergyLossBrems(bool opt = true) { fEnergyLossBrems = opt; }
  void setNoiseBrems(bool opt = true) { fNoiseBrems = opt; }

  double effects(const std::vector<TVector3> &points, const std::vector<double> &pointPaths, const double &mom, const int &pdg, const bool &doNoise = false,
                 TMatrixT<double> *noise = nullptr, const TMatrixT<double> *jacobian = nullptr, const TVector3 *directionBefore = nullptr,
                 const TVector3 *directionAfter = nullptr);


  double stepper(const double &maxDist, const double &posx, const double &posy, const double &posz, const double &dirx, const double &diry, const double &dirz, const double &mom,
                 const int &pdg);

  double stepper(const double &maxDist, const TVector3 &pos, const TVector3 &dir, const double &mom, const int &pdg)
  {
    return stepper(maxDist, pos.X(), pos.Y(), pos.Z(), dir.X(), dir.Y(), dir.Z(), mom, pdg);
  }

 private:
  void getParameters();

  void calcBeta(double mom);


  double energyLossBetheBloch(const double &mom);


  void noiseBetheBloch(const double &mom, TMatrixT<double> *noise) const;


  void noiseCoulomb(const double &mom, TMatrixT<double> *noise, const TMatrixT<double> *jacobian, const TVector3 *directionBefore, const TVector3 *directionAfter) const;


  double energyLossBrems(const double &mom) const;


  void noiseBrems(const double &mom, TMatrixT<double> *noise) const;

  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;

 public:
  ClassDef(GFMaterialEffects, 1)
};

#endif