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

class TGeoMaterial;

namespace genf {

  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_t>* noise = NULL,
                   const TMatrixT<Double_t>* jacobian = NULL,
                   const TVector3* directionBefore = NULL,
                   const TVector3* directionAfter = NULL);


    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:
    //  std::vector<GFAbsEnergyLoss*> fEnergyLoss;
    //GFGeoMatManager *geoMatManager;
    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;
    double MeanExcEnergy_get(int Z);
    double MeanExcEnergy_get(TGeoMaterial*);

    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)
  };

} // end namespace

#endif