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 MaterialEffectsSimulator;
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;
  MaterialEffectsSimulator* NuclearInteraction;

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

  // debugging
  //  double myEta;

  bool use_hardcoded;
};

#endif