ProfileSimulator.h

Go to the documentation of this file.

#ifndef _ProfileSimulatorOG_ProfileSimulator_h_
#define _ProfileSimulatorOG_ProfileSimulator_h_

#include <fwk/VModule.h>
#include <evt/GaisserHillas4Parameter.h>


namespace utl {
  class RandomEngine;
}

namespace ProfileSimulatorOG {

  class ProfileSimulator : public fwk::VModule {

  public:
    ProfileSimulator();
    virtual ~ProfileSimulator();

    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event& theEvent);
    fwk::VModule::ResultFlag Finish();

    enum ProfileType {
      eFromEnergy = 0,  // generate the profile using a
                        // parametrization of Nmax from energy
      eFromFixParameters // generates the profile from
                         // a given set of GH parameters
    };

  private:
    evt::GaisserHillas4Parameter fGH;
    //long fSeed; // unused. LN.
    double fMinEnergy;
    double fMaxEnergy;
    double fEnergyCutoff;
    int fPrimary;
    int fProfileType;
    bool fGenerateEnergy;
    double fAlphaSpectrum;
    double fXMaxSig;
    //double fThetaMin; // unused. LN.
    //double fThetaMax; // unused. LN.
    //double fNCos; // unused. LN.
    //double fProtonFraction; // unused. LN.
    double fTheta;
    double fPhi;
    bool fConstantX0;

    utl::RandomEngine* fRandomEngine = nullptr;

    double DicePowerLaw(double min, double max, double index) const;

    double RandomXOne(int primary, double energy) const;

    double RandomXZero(int primary, double energy) const;

    double RandomXMax(int primary, double energy)  const;

    double RandomLambda(int primary, double energy) const;

    double RandomNMax(int primary, double energy)  const;

    double CalculateNMax(double depth, double xZero, double xMax, double energy) const;

    double CalculatedEdX(double age, double density) const;

    REGISTER_MODULE("ProfileSimulatorOG", ProfileSimulator);

  };

}


#endif