SusyNtuple/Met.h

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

#include "TObject.h"
#include "TLorentzVector.h"

namespace Susy
{

class Met : public TObject
{
public:
    Met(){ clear(); }
    virtual ~Met(){};
    Met(const Met &);
    Met& operator=(const Met &);

    float Et;
    float phi;
    float sumet;

    const TLorentzVector lv() const {
      TLorentzVector tlv;
      tlv.SetPtEtaPhiE(this->Et,0,this->phi,this->Et);
      return tlv;
    }

    // MET Composition info - do we want TLorentzVectors, TVector2, or just floats?
    // TODO: clean out the obsolete term
    float refEle_et;          
    float refEle_phi;         
    float refEle_sumet;       
    float refGamma_et;        
    float refGamma_phi;       
    float refGamma_sumet;     
    float refTau_et;          
    float refTau_phi;         
    float refTau_sumet;       
    float refJet_et;          
    float refJet_phi;         
    float refJet_sumet;       
    float softTerm_et;        
    float softTerm_phi;       
    float softTerm_sumet;     
    float refMuo_et;          
    float refMuo_phi;         
    float refMuo_sumet;       
   
    // I'm not a big fan of the way this is setup, would like to improve it someday
    int sys;                  


    void print() const;

    void clear(){
      Et = phi = sumet = 0;
      refEle_et = refGamma_et = refTau_et = refJet_et = softTerm_et = refMuo_et = 0;
      refEle_phi = refGamma_phi = refTau_phi = refJet_phi = softTerm_phi = refMuo_phi = 0;
      refEle_sumet = refGamma_sumet = refTau_sumet = refJet_sumet = softTerm_sumet = refMuo_sumet = 0;
      sys = 0;
    }

    ClassDef(Met, 5);
};
} // Susy
#endif