GSNovaModel.h

Go to the documentation of this file.

#ifndef SNOVAMODEL__H
#define SNOVAMODEL__H

// C++ includes
#include <list>

// GENIE includes
#include "GENIE/Framework/ParticleData/PDGCodeList.h"

// NOvASoft includes
#include "GenieSNova/src/GenParticle.h"
#include "GenieSNova/src/GeneratorBase.h"
#include "GenieSNova/src/GFluxGenerator.h"

namespace genie {
    namespace supernova {
        class GSNovaModel: public GeneratorBase {

            public:
                GSNovaModel(std::string filename);
                void ReadModel(std::string filename);

                /// Get model features
                const genie::PDGCodeList& PdgList();

                /// Populate model with generators
                void AddGenerator(GFluxGenerator);

                /// Generate next particle
                void GenerateNext(GenParticle&);

                /// Check if there are no active generators left
                bool End();

                /// Reactivate all generators
                void Reload();
                void SetEmin(double Emin);
                void SetProbScale(double Scale);

                /// Access to generators
                const std::list<GFluxGenerator> ActiveGenerators();

            protected:
                std::list<GFluxGenerator>::iterator SelectGenerator();
                genie::PDGCodeList fPdgList; /// list of particles

                /// Flux containers
                std::list<GFluxGenerator> fGeneratorsActive;
                std::list<GFluxGenerator> fGeneratorsDead;

        };

    } // end of namespace supernova
} // end of namespace genie
#endif