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