MaterialEffects.h

Go to the documentation of this file.

#ifndef MaterialEffects_h
#define MaterialEffects_h

//Framework Headers
//#include "FWCore/ParameterSet/interface/ParameterSet.h"

// Geometry Headers
#include "DataFormats/GeometryVector/interface/GlobalVector.h"

#include <vector>
class FSimEvent;
class TrackerLayer;
class ParticlePropagator;
class PairProductionSimulator;
class BremsstrahlungSimulator;
class EnergyLossSimulator;
class NuclearInteractionSimulator;
class MultipleScatteringSimulator;
class MuonBremsstrahlungSimulator;
class RandomEngineAndDistribution;

namespace edm {
  class ParameterSet;
}

class MaterialEffects
{

 public:

  MaterialEffects(const edm::ParameterSet& matEff);

  ~MaterialEffects();

  void interact(FSimEvent& simEvent,
        const TrackerLayer& layer,
        ParticlePropagator& PP,
        unsigned i,
                RandomEngineAndDistribution const*);

  void save();

  inline double thickness() const { return theThickness; }

  inline double energyLoss() const { return theEnergyLoss; }

  inline MultipleScatteringSimulator* multipleScatteringSimulator() const { 
    return MultipleScattering;
  }

  inline EnergyLossSimulator* energyLossSimulator() const { 
    return EnergyLoss;
  }

  inline MuonBremsstrahlungSimulator* muonBremsstrahlungSimulator() const {
    return MuonBremsstrahlung;
  }

 private:

  double radLengths(const TrackerLayer& layer,
            ParticlePropagator& myTrack);

  GlobalVector normalVector(const TrackerLayer& layer,
                ParticlePropagator& myTrack ) const;

 private:

  PairProductionSimulator* PairProduction;
  BremsstrahlungSimulator* Bremsstrahlung;
  MuonBremsstrahlungSimulator* MuonBremsstrahlung;
  MultipleScatteringSimulator* MultipleScattering;
  EnergyLossSimulator* EnergyLoss;
  NuclearInteractionSimulator* NuclearInteraction;

  // Cuts for material effects
  double pTmin;
  GlobalVector theNormalVector;
  double theThickness;
  double theEnergyLoss;
  double theTECFudgeFactor;

  // debugging
  //  double myEta;

  bool use_hardcoded;

};

#endif