CorsikaShowerFileParticleIterator.h

Go to the documentation of this file.

#ifndef _io_CorsikaShowerFileParticleIterator_h_
#define _io_CorsikaShowerFileParticleIterator_h_

#include <utl/VShowerFileParticleIterator.h>
#include <io/RawCorsikaFile.h>
#include <io/CorsikaBlock.h>
#include <io/CorsikaUtilities.h>
#include <utl/Particle.h>
#include <utl/CoordinateSystemPtr.h>
#include <utl/Validated.h>


namespace utl {
  class Particle;
}

namespace io {

  namespace Corsika {
    class RawFile;
  }

  class CorsikaShowerFile;


  class CorsikaShowerFileParticleIterator : public utl::VShowerFileParticleIterator {

  public:
    typedef io::Corsika::RawFile::PositionType PositionType;

    // The following two constructors do exactly the same for different block types.
    // It cannot be a template because of a bug in gcc:
    // http://gcc.gnu.org/bugzilla/show_bug.cgi?id=9050
    CorsikaShowerFileParticleIterator(Corsika::RawFile& rawFile,
                                      const io::Corsika::EventHeader& header,
                                      const Corsika::AtmosphereParameters& params,
                                      const double timeOffset,
                                      const PositionType& startPosition,
                                      const unsigned int observationLevel,
                                      const bool isThinned,
                                      const double wmaxEM,
                                      const bool CERFIL = false,
                                      const double CERwlMin = 0,
                                      const double CERwlMax = 0);

    virtual ~CorsikaShowerFileParticleIterator() { }

    virtual utl::Particle* GetOneParticle(const utl::CoordinateSystemPtr& cs);
    virtual void Rewind();

  private:
    const io::Corsika::ParticleData* GetOneParticleRecord();
    const io::Corsika::ParticleDataUnthinned* GetOneParticleRecordUnthinned();

    double GetHeightAtDepth(const double depth) const;
    double GetDepthAtHeight(const double h) const;

    utl::Particle::Source GetEMParticleSource(const float wt, const short int gen) const;

    Corsika::RawFile& fRawFile;
    const Corsika::AtmosphereParameters fAtmPars;
    const double fTimeOffset;
    const PositionType fStartPosition;
    const double fVersion;

    const unsigned int fObservationLevel;

    PositionType fCurrentPosition;
    unsigned int fParticleInBlock = 0;
    bool fIteratorValid = false;

    utl::Particle fCurrentParticle;
    io::Corsika::Block fCurrentBlock;
    io::Corsika::BlockUnthinned fCurrentBlockUnthinned;

    const bool fIsThinned;
    const double fWMaxEM;
    const bool fCERFIL;
    const double fMCERFI;
    const double fWlLowerLimit;
    const double fWlUpperLimit;
    const bool fCurvedObsLevel;
    const bool fCurved;
    const double fObsLevelHeight;
    bool fBlockBufferValid = false;

    // This is to handle Corsika's history option.
    utl::Validated<utl::Point> fProductionPoint;
    utl::Validated<utl::Particle> fParent;
    utl::Validated<utl::Particle> fGrandParent;

  };

}


#endif