/data/refman/pasoursint/CMSSW_4_1_8_patch13/src/FastSimulation/MaterialEffects/interface/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 MultipleScatteringSimulator;
class BremsstrahlungSimulator;
class EnergyLossSimulator;
class NuclearInteractionSimulator;
class RandomEngine;

namespace edm {
  class ParameterSet;
}

class MaterialEffects
{

 public:

  MaterialEffects(const edm::ParameterSet& matEff,
                  const RandomEngine* engine);

  ~MaterialEffects();

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

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


 private:

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

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

 private:

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

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

  // debugging
  //  double myEta;

  // The random engine
  const RandomEngine* random;

};

#endif