TMBMetEx.cpp

Go to the documentation of this file.

//                                                                      //
// TMBTree MetEx class                                                    //
//                                                                      //

#include <iostream>
#include <string>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>

#include <cmath>

#include "tmb_tree/TMBMetEx.hpp"

using namespace std;

ClassImp(TMBMetEx);


//_____________________________________________________________________________

TMBMetEx::TMBMetEx()

{}

TMBMetEx::~TMBMetEx() 
{}


void TMBMetEx::getRings(const float *&RingEMx,
                         const float *&RingEMy,
                         const float *&RingHDx,
                         const float *&RingHDy) const
{
  RingEMx = &_RingEMx[0];
  RingEMy = &_RingEMy[0];
  RingHDx = &_RingHDx[0];
  RingHDy = &_RingHDy[0];
}

void TMBMetEx::setRings( const float *&aRingEMx,
                         const float *&aRingEMy,
                         const float *&aRingHDx,
                         const float *&aRingHDy)
{
  for(int iring=0; iring < (int )metid::NRING; iring++){
    _RingEMx[iring] = aRingEMx[iring];
    _RingEMy[iring] = aRingEMy[iring];
    _RingHDx[iring] = aRingHDx[iring];
    _RingHDy[iring] = aRingHDy[iring];
  }
}

void TMBMetEx::print(std::ostream &os) const
{
  print_METweta (os);
  print_METT   (os);
  print_METTM  (os);
  print_METTAS (os);
  print_METTBS (os);
  print_METTAN (os);
  print_METTBN (os);
  print_METCM  (os);
  print_METCAS (os);
  print_METCBS (os);
  print_METCAN (os);
  print_METCBN (os);
  print_METICD (os);
  print_METNADA(os);
  print_METMUON(os);
  print_METDM  (os);
  print_METEM  (os);
  print_METMG  (os);
  print_METFH  (os);
  print_METCH  (os);
  print_METED  (os);
  print_METNG  (os);
  print_MET42  (os);
  print_rings  (os);
}


void TMBMetEx::print_rings(std::ostream& os) const
{
  os << "MET Rings :      EMx,         EMy,         HDx,          HDy" <<endl;
  for (int iring=0; iring != (int)metid::NRING; iring++){
    os << " RingEMx[" << setw(2) << iring << "] = " 
       << setw(10) << _RingEMx[iring] << " , "
       << setw(10) << _RingEMy[iring] << " , "
       << setw(10) << _RingHDx[iring] << " , "
       << setw(10) << _RingHDy[iring] << endl;
  }  
}

void  TMBMetEx::getVisMu(float &visMuX, float &visMuY, float &visMuPt) const
{
  visMuX  = -1. * _metMUON.getmex();;
  visMuY  = -1. * _metMUON.getmey();;
  visMuPt = -1. * _metMUON.getmet();;
}


void TMBMetEx::getMETT   (float &SETT,     float &METTx,    float &METTy,    float &METT)    const
{
  SETT  = _metT.getset();
  METTx = _metT.getmex();
  METTy = _metT.getmey();
  METT  = _metT.getmet();
}
void TMBMetEx::getMETTM  (float &SETTM,    float &METTMx,   float &METTMy,   float &METTM)   const
{
  SETTM  = _metTM.getset();
  METTMx = _metTM.getmex();
  METTMy = _metTM.getmey();
  METTM  = _metTM.getmet();
}
void TMBMetEx::getVETTAS (float &SETTAS,   float &VETTASx,  float &VETTASy,  float &VETTAS)  const
{
  SETTAS  = _metTAS.getset();
  VETTASx = -1. * _metTAS.getmex();
  VETTASy = -1. * _metTAS.getmey();
  VETTAS  = _metTAS.getmet();
}
void TMBMetEx::getVETTBS (float &SETTBS,   float &VETTBSx,  float &VETTBSy,  float &VETTBS)  const
{
  SETTBS  = _metTBS.getset();
  VETTBSx = -1. * _metTBS.getmex();
  VETTBSy = -1. * _metTBS.getmey();
  VETTBS  = _metTBS.getmet();
}
void TMBMetEx::getVETTAN (float &SETTAN,   float &VETTANx,  float &VETTANy,  float &VETTAN)  const
{
  SETTAN  = _metTAN.getset();
  VETTANx = -1. * _metTAN.getmex();
  VETTANy = -1. * _metTAN.getmey();
  VETTAN  = _metTAN.getmet();
}
void TMBMetEx::getVETTBN (float &SETTBN,   float &VETTBNx,  float &VETTBNy,  float &VETTBN)  const
{
  SETTBN  = _metTBN.getset();
  VETTBNx = -1. * _metTBN.getmex();
  VETTBNy = -1. * _metTBN.getmey();
  VETTBN  = _metTBN.getmet();
}

void TMBMetEx::getMETCM  (float &SETCM,    float &METCMx,   float &METCMy,   float &METCM)   const
{
  SETCM  = _metCM.getset();
  METCMx = _metCM.getmex();
  METCMy = _metCM.getmey();
  METCM  = _metCM.getmet();
}
void TMBMetEx::getVETCAS (float &SETCAS,   float &VETCASx,  float &VETCASy,  float &VETCAS)  const
{
  SETCAS  = _metCAS.getset();
  VETCASx = -1. * _metCAS.getmex();
  VETCASy = -1. * _metCAS.getmey();
  VETCAS  = _metCAS.getmet();
}
void TMBMetEx::getVETCBS (float &SETCBS,   float &VETCBSx,  float &VETCBSy,  float &VETCBS)  const
{
  SETCBS  = _metCBS.getset();
  VETCBSx = -1. * _metCBS.getmex();
  VETCBSy = -1. * _metCBS.getmey();
  VETCBS  = _metCBS.getmet();
}
void TMBMetEx::getVETCAN (float &SETCAN,   float &VETCANx,  float &VETCANy,  float &VETCAN)  const
{
  SETCAN  = _metCAN.getset();
  VETCANx = -1. * _metCAN.getmex();
  VETCANy = -1. * _metCAN.getmey();
  VETCAN  = _metCAN.getmet();
}
void TMBMetEx::getVETCBN (float &SETCBN,   float &VETCBNx,  float &VETCBNy,  float &VETCBN)  const
{
  SETCBN  = _metCBN.getset();
  VETCBNx = -1. * _metCBN.getmex();
  VETCBNy = -1. * _metCBN.getmey();
  VETCBN  = _metCBN.getmet();
}
void TMBMetEx::getVETICD (float &SETICD,   float &VETICDx,  float &VETICDy,  float &VETICD)  const
{
  SETICD  = _metICD.getset();
  VETICDx = -1. * _metICD.getmex();
  VETICDy = -1. * _metICD.getmey();
  VETICD  = _metICD.getmet();
}
void TMBMetEx::getVETNADA(float &SETNADA,  float &VETNADAx, float &VETNADAy, float &VETNADA) const
{
  SETNADA  = _metNADA.getset();
  VETNADAx = -1. * _metNADA.getmex();
  VETNADAy = -1. * _metNADA.getmey();
  VETNADA  = _metNADA.getmet();
}
void TMBMetEx::getVETMUON(float &VETMUONx, float &VETMUONy, float &VETMUONz, float &VETMUON) const
{
  VETMUONx = -1. * _metMUON.getmex();
  VETMUONy = -1. * _metMUON.getmey();
  VETMUONz = -1. * _metMUON.getmez();
  VETMUON  = _metMUON.getmet();
}

void TMBMetEx::getVETCM  (float &SETCM,    float &VETCMx,   float &VETCMy,   float &METCM)   const
{
  SETCM  = _metCM.getset();
  VETCMx = -1. * _metCM.getmex();
  VETCMy = -1. * _metCM.getmey();
  METCM  = _metCM.getmet();
}

void TMBMetEx::getVETDM  (float &SETDM,    float &VETDMx,   float &VETDMy,   float &METDM)   const
{
  SETDM  = _metDM.getset();
  VETDMx = -1. * _metDM.getmex();
  VETDMy = -1. * _metDM.getmey();
  METDM  = _metDM.getmet();
}
void TMBMetEx::getMETDM  (float &SETDM,    float &METDMx,   float &METDMy,   float &METDM)   const
{
  SETDM  = _metDM.getset();
  METDMx = _metDM.getmex();
  METDMy = _metDM.getmey();
  METDM  = _metDM.getmet();
}
void TMBMetEx::getVETEM  (float &SETEM,    float &VETEMx,   float &VETEMy,   float &METEM)   const
{
  SETEM  = _metEM.getset();
  VETEMx = -1. * _metEM.getmex();
  VETEMy = -1. * _metEM.getmey();
  METEM  = _metEM.getmet();
}
void TMBMetEx::getVETMG  (float &SETMG,    float &VETMGx,   float &VETMGy,   float &METMG)   const
{
  SETMG  = _metMG.getset();
  VETMGx = -1. * _metMG.getmex();
  VETMGy = -1. * _metMG.getmey();
  METMG  = _metMG.getmet();
}
void TMBMetEx::getVETFH  (float &SETFH,    float &VETFHx,   float &VETFHy,   float &METFH)   const
{
  SETFH  = _metFH.getset();
  VETFHx = -1. * _metFH.getmex();
  VETFHy = -1. * _metFH.getmey();
  METFH  = _metFH.getmet();
}
void TMBMetEx::getVETCH  (float &SETCH,    float &VETCHx,   float &VETCHy,   float &METCH)   const
{
  SETCH  = _metCH.getset();
  VETCHx = -1. * _metCH.getmex();
  VETCHy = -1. * _metCH.getmey();
  METCH  = _metCH.getmet();
}
void TMBMetEx::getVETED  (float &EDGE,     float &SETED,    float &VETEDx,   float &VETEDy, float &METED) const
{
  EDGE   = 0.;
  SETED  = _metED.getset();
  VETEDx = -1. * _metED.getmex();
  VETEDy = -1. * _metED.getmey();
  METED  = _metED.getmet();
}
void TMBMetEx::getVETNG  (float &SETNG,    float &VETNGx,   float &VETNGy,   float &METNG)   const
{
  SETNG  = _metNG.getset();
  VETNGx = -1. * _metNG.getmex();
  VETNGy = -1. * _metNG.getmey();
  METNG  = _metNG.getmet();
}
void TMBMetEx::getVET42  (float &SET42,    float &VET42x,   float &VET42y,   float &MET42)   const
{
  SET42  = _metT42.getset();
  VET42x = -1. * _metT42.getmex();
  VET42y = -1. * _metT42.getmey();
  MET42  = _metT42.getmet();
}


//______________________________________________________________________________

---