EnvelopeSlice.h

Go to the documentation of this file.

#ifndef _ENVELOPE_SLICE_H
#define _ENVELOPE_SLICE_H

#include <cmath>

enum EnvelopeSliceParameter {
    SLI_z    = 0,
    SLI_beta = 1,
    SLI_x    = 2,
    SLI_px   = 3,
    SLI_y    = 4,
    SLI_py   = 5,
    SLI_x0   = 6,
    SLI_px0  = 7,
    SLI_y0   = 8,
    SLI_py0  = 9,
    SLNPAR
};

class EnvelopeSlice {

public:

    EnvelopeSlice();

    virtual ~EnvelopeSlice() {}

    double p[SLNPAR];

    double p_old[SLNPAR];

    inline double computeGamma() {
        return sqrt(1.0 / (1.0 - p[SLI_beta] * p[SLI_beta]));
    }

    // backup the slice
    void backup();

    // restore a slice
    void restore();

    // check if a slice has changed
    int check();

    inline int is_valid() { return isValid_m; }

    bool hasSliceBeenEmitted_m;

private:
    int isValid_m;
    int backupedValid_m;

};

#endif