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