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"

#include "TDatabasePDG.h"
#include "TGeoMaterial.h"
#include "TGeoManager.h"

#include "larreco/Genfit/GFAbsEnergyLoss.h"
#include "larreco/Genfit/GFGeoMatManager.h"

  
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