utilities.hh

Go to the documentation of this file.

// Common utilities

#ifndef H_UTILITIES
#define H_UTILITIES

#include <vector>
#include <string>
#include <iostream>
#include <stdexcept>
#include "DataFormats/PatCandidates/interface/Electron.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <fastjet/PseudoJet.hh>
#include "TString.h"

typedef std::vector<const reco::Candidate*> vCands;
typedef reco::Candidate::LorentzVector LVector;

#define ERROR(x) do{throw cms::Exception("babymaker") << "In " << __FILE__ << " at line " << __LINE__ << " (in function " << __func__ << "): " << x << std::endl;}while(false)
#define DBG(x) do{edm::LogError("babymaker") << "In " << __FILE__ << " at line " << __LINE__ << " (in function " << __func__ << "): " << x << std::endl;}while(false)

namespace utilities{

  enum SysEnum{kSysJER, kSysJECUp, kSysJECDn, kSysLast};
  bool contains(const std::string &s, const std::string &pat);
  float sumMass(const LVector &a, const LVector &b);
  float sumPt(const LVector &a, const LVector &b);
  float dPhi(float phi1, float phi2);
  float dR(float phi1, float phi2, float eta1, float eta2);
  bool greaterPt(const reco::Candidate *a, const reco::Candidate *b);
  bool greaterM(const fastjet::PseudoJet &a, const fastjet::PseudoJet &b);
  float getMT(float pt1, float phi1, float pt2, float phi2);
  float getMT(float m1, float pt1, float phi1,
          float m2, float pt2, float phi2);
  float getMT2(float pt1, float phi1, float pt2, float phi2, float met, float met_phi);
  float getMT2(float m1, float pt1, float phi1,
           float m2, float pt2, float phi2,
           float met, float met_phi);
  std::string execute(const std::string &cmd);
  TString roundNumber(double num, int decimals, double denom=1.);
  TString addCommas(double num);

}

#endif