SusyNtuple/Tau.h

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

#include "SusyNtuple/Particle.h"
#include "SusyNtuple/SusyNtSys.h"

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

    int q;                    
    int author;               
    int nTrack;               
    float eleBDT;             
    float jetBDT;             
    int jetBDTSigLoose;       
    int jetBDTSigMedium;      
    int jetBDTSigTight;       
    int eleBDTLoose;          
    int eleBDTMedium;         
    int eleBDTTight;          
    int muonVeto;             

    bool trueTau;             

    int truthType;            

//obsolete -remove 
    bool matched2TruthLepton; 
    int detailedTruthType;    

    // Taus have a lot of SFs
    float looseEffSF;         
    float errLooseEffSF;      
    float mediumEffSF;        
    float errMediumEffSF;     
    float tightEffSF;         
    float errTightEffSF;      
    float looseEVetoSF;       
    float errLooseEVetoSF;    
    float mediumEVetoSF;      
    float errMediumEVetoSF;   
    float tightEVetoSF;       
    float errTightEVetoSF;    
 
    // Systematic factors
    float sme_total_up;       
    float sme_total_dn;       

    long long trigFlags;      


    bool matchTrig(long long mask) const {
      return (trigFlags & mask) == mask;
    }

    void setState(int sys);

    void print() const;

    void clear(){
        q = author = nTrack = 0;
        eleBDT = jetBDT = 0;
        jetBDTSigLoose = jetBDTSigMedium = jetBDTSigTight = 0;
        eleBDTLoose = eleBDTMedium = eleBDTTight = 0;
        muonVeto = 0;
        trueTau = false;
        matched2TruthLepton = false;
        truthType = detailedTruthType = -1;
        looseEffSF = mediumEffSF = tightEffSF = 1;
        errLooseEffSF = errMediumEffSF = errTightEffSF = 0;
        looseEVetoSF = mediumEVetoSF = tightEVetoSF = 1;
        errLooseEVetoSF = errMediumEVetoSF = errTightEVetoSF = 0;
        sme_total_up = sme_total_dn = 0;
        trigFlags = 0;
        Particle::clear();
    }

    ClassDef(Tau, 7);
  };
} // Susy
#endif