/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/FastSimulation/MaterialEffects/src/MaterialEffectsSimulator.cc

Go to the documentation of this file.

#include "FastSimulation/MaterialEffects/interface/MaterialEffectsSimulator.h"
//#include "FastSimulation/Utilities/interface/RandomEngine.h"

#include <list>

using std::list;
using std::pair;

MaterialEffectsSimulator:: MaterialEffectsSimulator(const RandomEngine* engine,
                                                    double A, double Z, 
                                                    double density, double radLen) :
  A(A), Z(Z), density(density), radLen(radLen)
{ 
  random = engine;
  _theUpdatedState.clear(); 
}

MaterialEffectsSimulator::~MaterialEffectsSimulator() {
  // Don't delete the objects contained in the list
  _theUpdatedState.clear();
}

void MaterialEffectsSimulator::updateState(ParticlePropagator & Particle,
                                           double radlen)
{

  _theUpdatedState.clear();
  theClosestChargedDaughterId = -1;

  radLengths = radlen;
  if ( radLengths > 0. ) compute(Particle);

}

XYZVector 
MaterialEffectsSimulator::orthogonal(const XYZVector& aVector) const { 

  double x = fabs(aVector.X());
  double y = fabs(aVector.Y());
  double z = fabs(aVector.Z());

  if ( x < y ) 
    return x < z ? 
      XYZVector(0.,aVector.Z(),-aVector.Y()) :
      XYZVector(aVector.Y(),-aVector.X(),0.);
  else
    return y < z ? 
      XYZVector(-aVector.Z(),0.,aVector.X()) :
      XYZVector(aVector.Y(),-aVector.X(),0.);

}