logintpack Namespace Reference

Functions
int16_t	pack16log (double x, double lmin, double lmax, uint16_t base=32768)
int16_t	pack16logCeil (double x, double lmin, double lmax, uint16_t base=32768)
int16_t	pack16logClosed (double x, double lmin, double lmax, uint16_t base=32768)
int8_t	pack8log (double x, double lmin, double lmax, uint8_t base=128)
int8_t	pack8logCeil (double x, double lmin, double lmax, uint8_t base=128)
int8_t	pack8logClosed (double x, double lmin, double lmax, uint8_t base=128)
double	unpack16log (int16_t i, double lmin, double lmax, uint16_t base=32768)
double	unpack16logClosed (int16_t i, double lmin, double lmax, uint16_t base=32768)
	reverse of pack8logClosed More...
double	unpack8log (int8_t i, double lmin, double lmax, uint8_t base=128)
double	unpack8logClosed (int8_t i, double lmin, double lmax, uint8_t base=128)
	reverse of pack8logClosed More...

Variables
constexpr int8_t	smallestNegative = -1
constexpr int8_t	smallestPositive = 0

Function Documentation

int16_t logintpack::pack16log ( double  x, double  lmin, double  lmax, uint16_t  base = 32768  )

inline

Definition at line 26 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, delta, JetChargeProducer_cfi::exp, checklumidiff::l, cmsBatch::log, and alignCSCRings::r.

Referenced by CompressionElement::pack().

    {
            if(base>32768) base=32768;
                const float delta=(log(1.+exp((lmax-lmin)/base))-log(2.))*base/(lmax-lmin);
        const double l = std::log(std::abs(x));
        const double centered = (l-lmin)/(lmax-lmin)*base;
        int16_t  r=std::floor(centered);
        if(centered-r>delta) r+=1;
        if(centered >= base-1) r=base-1;
        if(centered < 0) r=0;
        if(x<0) r = r==0 ? -1 : -r;
        return r;
    }

int16_t logintpack::pack16logCeil ( double  x, double  lmin, double  lmax, uint16_t  base = 32768  )

inline

Definition at line 13 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, checklumidiff::l, cmsBatch::log, and alignCSCRings::r.

        {
                if(base>32768) base=32768;
                const double l = std::log(std::abs(x));
                const double centered = (l-lmin)/(lmax-lmin)*base;
                int16_t  r=std::ceil(centered);
                if(centered >= base-1) r=base-1;
                if(centered < 0) r=0;
                if(x<0) r = r==0 ? -1 : -r;
                return r;
        }

int16_t logintpack::pack16logClosed ( double  x, double  lmin, double  lmax, uint16_t  base = 32768  )

inline

pack a value x distributed in [-1,1], with guarantee that -1 and 1 are preserved exactly in packing and unpacking. tries to keep the best precision for x close to the endpoints, sacrifying that in the middle

Definition at line 42 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, checklumidiff::l, cmsBatch::log, and alignCSCRings::r.

    {
            if(base>32768) base=32768;
        const double l = std::log(std::abs(x));
        const double centered = (l-lmin)/(lmax-lmin)*(base-1);
        int16_t  r=round(centered);
        if(centered >= base-1) r=base-1;
        if(centered < 0) r=0;
        if(x<0) r = r==0 ? -1 : -r;
        return r;
    }

int8_t logintpack::pack8log ( double  x, double  lmin, double  lmax, uint8_t  base = 128  )

inline

Definition at line 91 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, checklumidiff::l, cmsBatch::log, and alignCSCRings::r.

    {
            if(base>128) base=128;
        const double l = std::log(std::abs(x));
        const double centered = (l-lmin)/(lmax-lmin)*base;
        int8_t  r=centered;
        if(centered >= base-1) r=base-1;
        if(centered < 0) r=0;
        if(x<0) r = r==0 ? -1 : -r;
        return r;
    }

int8_t logintpack::pack8logCeil ( double  x, double  lmin, double  lmax, uint8_t  base = 128  )

inline

Definition at line 78 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, checklumidiff::l, cmsBatch::log, and alignCSCRings::r.

        {
                if(base>128) base=128;
                const double l = std::log(std::abs(x));
                const double centered = (l-lmin)/(lmax-lmin)*base;
                int8_t  r=std::ceil(centered);
                if(centered >= base-1) r=base-1;
                if(centered < 0) r=0;
                if(x<0) r = r==0 ? -1 : -r;
                return r;
        }

int8_t logintpack::pack8logClosed ( double  x, double  lmin, double  lmax, uint8_t  base = 128  )

inline

pack a value x distributed in [-1,1], with guarantee that -1 and 1 are preserved exactly in packing and unpacking. tries to keep the best precision for x close to the endpoints, sacrifying that in the middle

Definition at line 106 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, checklumidiff::l, cmsBatch::log, and alignCSCRings::r.

Referenced by pat::PackedCandidate::setPuppiWeight().

    {
            if(base>128) base=128;
        const double l = std::log(std::abs(x));
        const double centered = (l-lmin)/(lmax-lmin)*(base-1);
        int8_t  r=round(centered);
        if(centered >= base-1) r=base-1;
        if(centered < 0) r=0;
        if(x<0) r = r==0 ? -1 : -r;
        return r;
    }

double logintpack::unpack16log ( int16_t  i, double  lmin, double  lmax, uint16_t  base = 32768  )

inline

Definition at line 55 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, JetChargeProducer_cfi::exp, checklumidiff::l, and heppy_batch::val.

Referenced by CompressionElement::unpack().

    {
            if(base>32768) base=32768;
            const double basef=base;
        const double l=lmin+std::abs(i)/basef*(lmax-lmin);
        const double val=std::exp(l);
        if(i<0) return -val; else return val;
    }

double logintpack::unpack16logClosed ( int16_t  i, double  lmin, double  lmax, uint16_t  base = 32768  )

inline

reverse of pack8logClosed

Definition at line 65 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, JetChargeProducer_cfi::exp, checklumidiff::l, and heppy_batch::val.

    {
            if(base>32768) base=32768;
            const double basef=base-1;
        double l=lmin+std::abs(i)/basef*(lmax-lmin);
        if (std::abs(i) == base-1) l = lmax;
        const double val=std::exp(l);
        if(i<0) return -val; else return val;
    }

double logintpack::unpack8log ( int8_t  i, double  lmin, double  lmax, uint8_t  base = 128  )

inline

Definition at line 119 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, JetChargeProducer_cfi::exp, checklumidiff::l, and heppy_batch::val.

    {
            if(base>128) base=128;
            const double basef=base;
        const double l=lmin+std::abs(i)/basef*(lmax-lmin);
        const double val=std::exp(l);
        if(i<0) return -val; else return val;
    }

double logintpack::unpack8logClosed ( int8_t  i, double  lmin, double  lmax, uint8_t  base = 128  )

inline

reverse of pack8logClosed

Definition at line 130 of file liblogintpack.h.

References funct::abs(), runEdmFileComparison::base, JetChargeProducer_cfi::exp, checklumidiff::l, and heppy_batch::val.

Referenced by pat::PackedCandidate::puppiWeight(), and pat::PackedCandidate::puppiWeightNoLep().

    {
            if(base>128) base=128;
            const double basef=base-1;
        double l=lmin+std::abs(i)/basef*(lmax-lmin);
        if (std::abs(i) == base-1) l = lmax;
        const double val=std::exp(l);
        if(i<0) return -val; else return val;
    }

Variable Documentation

constexpr int8_t logintpack::smallestNegative = -1

Definition at line 12 of file liblogintpack.h.

constexpr int8_t logintpack::smallestPositive = 0

Definition at line 9 of file liblogintpack.h.