RadialInterval.h

Go to the documentation of this file.

#ifndef RadialInterval_H
#define RadialInterval_H

//  Created 1/11/04. F. Beaudette (CERN)
//  This class is used to ease the lateral development with
// different spot fractions in FamosShower.

#include <vector>

class RandomEngineAndDistribution;

class RadialInterval {
public:
  RadialInterval(double RC, unsigned nSpots, double energy, const RandomEngineAndDistribution* engine);
  ~RadialInterval() { ; }

  void addInterval(double, double);
  void compute();
  inline unsigned nIntervals() const { return nInter; }
  inline double getSpotEnergy(unsigned i) const {
    //    std::cout << " getSpotEnergy " << i << " " << spotE.size() << std::endl;
    return spotE[i];
  }
  inline unsigned getNumberOfSpots(unsigned i) const {
    //    std::cout << " getNumberOfSpots " << i << " " << nspots.size() << std::endl;
    return nspots[i];
  }
  inline double getUmin(unsigned i) const {
    //    std::cout << " getUmin " << i << " " << uMin.size() << std::endl;
    return uMin[i];
  }
  inline double getUmax(unsigned i) const {
    //    std::cout << " getUmax " << i << " " << uMax.size() << std::endl;
    return uMax[i];
  }

private:
  double currentRad;
  double currentEnergyFraction;
  double currentUlim;
  double theR;
  unsigned theNumberOfSpots;
  double theSpotEnergy;
  unsigned nInter;

  std::vector<double> uMin;
  std::vector<double> uMax;
  std::vector<unsigned> nspots;
  std::vector<double> spotE;
  std::vector<double> dspotsunscaled;
  std::vector<double> spotfraction;

private:
  // Fraction of the energy in rm Moliere radius
  double energyFractionInRadius(double rm);

  // Famos Random Engine
  const RandomEngineAndDistribution* random;
};
#endif