SusyNtuple/Lepton.h

// Dear emacs, this is -*- c++ -*-
#ifndef SUSYNTUPLE_LEPTON_H
#define SUSYNTUPLE_LEPTON_H

#include "SusyNtuple/Particle.h"

#include "TBits.h"

namespace Susy
{
class Lepton : public Particle
{
public:
    Lepton() :
        trigBits(m_nTriggerBits)
        {
            clear();
        }
    virtual ~Lepton(){};
    Lepton(const Lepton &);
    Lepton& operator=(const Lepton &);

    // public member vars
    int q;                    
    bool isBaseline; 
    bool isSignal; 

    //Iso var -- same fer ele & mu
    //Etcone to be drop later in run-2
    float etconetopo20;       
    float etconetopo30;       
    float ptcone20;           
    float ptcone30;           
    float ptvarcone20;        
    float ptvarcone30;        

    // IsolationSelectionTool flags
    bool isoGradientLoose;    
    bool isoGradient;         
    bool isoLooseTrackOnly;   
    bool isoLoose;            
    bool isoTight;            


    float d0;                 
    float errD0;              
    float d0sigBSCorr;        
    float z0;                 
    float errZ0;              

    unsigned int mcType;      
    unsigned int mcOrigin;    

    bool matched2TruthLepton; 
    int truthType;            

    float effSF;              
    //float errEffSF;           ///< Uncertainty on the efficiency scale factor (for electron from LH)

    TBits   trigBits;         
    static const size_t m_nTriggerBits=64;


    float d0Sig() const {
      return d0/errD0;
    }
    float z0SinTheta() const {
      return z0 * sin(Theta());
    }

    // Polymorphism, baby!!
    virtual bool isEle() const { return false; }
    virtual bool isMu()  const { return false; }
    void setState(int sys){ resetTLV(); }

    virtual void print() const {};

    void clear(){
        q = 0;
        etconetopo20 = etconetopo30
                     = ptcone20    = ptcone30 
                     = ptvarcone20 = ptvarcone30 = 0;
        isoGradientLoose = isoGradient = isoLooseTrackOnly = isoLoose = isoTight = false;
        d0 = errD0 = d0sigBSCorr = z0 = errZ0 = 0;
        mcType = mcOrigin = 0;
        matched2TruthLepton = false;
        truthType = -1;
        effSF = 1;
        //errEffSF = 0;
        isBaseline = isSignal = false;
        Particle::clear();
    }

    ClassDef(Lepton, 20);
};
} // Susy
#endif