MEtoEDM< T > Class Template Reference

#include <MEtoEDMFormat.h>

Classes
struct	MEtoEDMObject

Public Types
typedef std::vector< MEtoEDMObject >	MEtoEdmObjectVector
typedef std::vector< uint32_t >	TagList

Public Member Functions
bool	CheckBinLabels (const TAxis *a1, const TAxis *a2)
const MEtoEdmObjectVector &	getMEtoEdmObject () const
bool	mergeProduct (const MEtoEDM< T > &newMEtoEDM)
template<>
bool	mergeProduct (const MEtoEDM< double > &newMEtoEDM)
template<>
bool	mergeProduct (const MEtoEDM< int > &newMEtoEDM)
template<>
bool	mergeProduct (const MEtoEDM< long long > &newMEtoEDM)
template<>
bool	mergeProduct (const MEtoEDM< TString > &newMEtoEDM)
	MEtoEDM ()
	MEtoEDM (size_t reservedSize)
void	putMEtoEdmObject (const std::string &name, const T &object)
void	swap (MEtoEDM< T > &iOther)
virtual	~MEtoEDM ()

Private Attributes
MEtoEdmObjectVector	MEtoEdmObject

Detailed Description

template<class T>
class MEtoEDM< T >

DataFormat class to hold the information from a ME tranformed into ROOT objects as appropriate

Author

M. Strang SUNY-Buffalo

Definition at line 37 of file MEtoEDMFormat.h.

Member Typedef Documentation

MEtoEdmObjectVector

template<class T>

typedef std::vector<MEtoEDMObject> MEtoEDM< T >::MEtoEdmObjectVector

Definition at line 51 of file MEtoEDMFormat.h.

TagList

template<class T>

typedef std::vector<uint32_t> MEtoEDM< T >::TagList

Definition at line 43 of file MEtoEDMFormat.h.

Constructor & Destructor Documentation

MEtoEDM() [1/2]

template<class T>

MEtoEDM< T >::MEtoEDM ( )

inline

Definition at line 39 of file MEtoEDMFormat.h.

{}

MEtoEDM() [2/2]

template<class T>

MEtoEDM< T >::MEtoEDM ( size_t  reservedSize )

inlineexplicit

Definition at line 40 of file MEtoEDMFormat.h.

References MEtoEDM< T >::MEtoEdmObject.

{ MEtoEdmObject.reserve(reservedSize); }

~MEtoEDM()

template<class T>

virtual MEtoEDM< T >::~MEtoEDM ( )

inlinevirtual

Definition at line 41 of file MEtoEDMFormat.h.

{}

Member Function Documentation

CheckBinLabels()

template<class T>

bool MEtoEDM< T >::CheckBinLabels ( const TAxis *  a1, const TAxis *  a2  )

inline

Definition at line 62 of file MEtoEDMFormat.h.

References testProducerWithPsetDescEmpty_cfi::a2, mps_fire::i, and hbheRecHitProducerGPUTask_cff::label1.

Referenced by MEtoEDM< T >::mergeProduct().

                                                        {
    // check that axis have same labels
    THashList *l1 = (const_cast<TAxis *>(a1))->GetLabels();
    THashList *l2 = (const_cast<TAxis *>(a2))->GetLabels();

    if (!l1 && !l2)
      return true;
    if (!l1 || !l2) {
      return false;
    }
    // check now labels sizes  are the same
    if (l1->GetSize() != l2->GetSize()) {
      return false;
    }
    for (int i = 1; i <= a1->GetNbins(); ++i) {
      TString label1 = a1->GetBinLabel(i);
      TString label2 = a2->GetBinLabel(i);
      if (label1 != label2) {
        return false;
      }
    }
    return true;
  }

getMEtoEdmObject()

template<class T>

const MEtoEdmObjectVector& MEtoEDM< T >::getMEtoEdmObject ( ) const

inline

Definition at line 60 of file MEtoEDMFormat.h.

References MEtoEDM< T >::MEtoEdmObject.

Referenced by MEtoEDM< T >::mergeProduct().

{ return MEtoEdmObject; }

mergeProduct() [1/5]

template<class T>

bool MEtoEDM< T >::mergeProduct ( const MEtoEDM< T > &  newMEtoEDM )

inline

Definition at line 86 of file MEtoEDMFormat.h.

References MEtoEDM< T >::CheckBinLabels(), gather_cfg::cout, MEtoEDM< T >::getMEtoEdmObject(), mps_fire::i, dqmiolumiharvest::j, MEtoEDM< T >::MEtoEdmObject, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                  {
    const MEtoEdmObjectVector &newMEtoEDMObject = newMEtoEDM.getMEtoEdmObject();
    const size_t nObjects = newMEtoEDMObject.size();
    //  NOTE: we remember the present size since we will only add content
    //        from newMEtoEDMObject after this point
    const size_t nOldObjects = MEtoEdmObject.size();

    // if the old and new are not the same size, we want to report a problem
    if (nObjects != nOldObjects) {
      std::cout << "WARNING MEtoEDM::mergeProducts(): the lists of histograms to be merged have different sizes: new="
                << nObjects << ", old=" << nOldObjects << std::endl;
    }

    for (unsigned int i = 0; i < nObjects; ++i) {
      unsigned int j = 0;
      // see if the name is already in the old container up to the point where
      // we may have added new entries in the container
      const std::string &name = newMEtoEDMObject[i].name;
      if (i < nOldObjects && (MEtoEdmObject[i].name == name)) {
        j = i;
      } else {
        j = 0;
        while (j < nOldObjects && (MEtoEdmObject[j].name != name))
          ++j;
      }
      if (j >= nOldObjects) {
        // this value is only in the new container, not the old one
#if METOEDMFORMAT_DEBUG
        std::cout << "WARNING MEtoEDM::mergeProducts(): adding new histogram '" << name << "'" << std::endl;
#endif
        MEtoEdmObject.push_back(newMEtoEDMObject[i]);
      } else if (MEtoEdmObject[j].object.CanExtendAllAxes() && newMEtoEDMObject[i].object.CanExtendAllAxes()) {
        TList list;
        list.Add((TObject *)&newMEtoEDMObject[i].object);
        if (MEtoEdmObject[j].object.Merge(&list) == -1) {
          std::cout << "ERROR MEtoEDM::mergeProducts(): merge failed for '" << name << "'" << std::endl;
        }
      } else {
        // this value is also in the new container: add the two
        if (MEtoEdmObject[j].object.GetNbinsX() == newMEtoEDMObject[i].object.GetNbinsX() &&
            MEtoEdmObject[j].object.GetXaxis()->GetXmin() == newMEtoEDMObject[i].object.GetXaxis()->GetXmin() &&
            MEtoEdmObject[j].object.GetXaxis()->GetXmax() == newMEtoEDMObject[i].object.GetXaxis()->GetXmax() &&
            MEtoEdmObject[j].object.GetNbinsY() == newMEtoEDMObject[i].object.GetNbinsY() &&
            MEtoEdmObject[j].object.GetYaxis()->GetXmin() == newMEtoEDMObject[i].object.GetYaxis()->GetXmin() &&
            MEtoEdmObject[j].object.GetYaxis()->GetXmax() == newMEtoEDMObject[i].object.GetYaxis()->GetXmax() &&
            MEtoEdmObject[j].object.GetNbinsZ() == newMEtoEDMObject[i].object.GetNbinsZ() &&
            MEtoEdmObject[j].object.GetZaxis()->GetXmin() == newMEtoEDMObject[i].object.GetZaxis()->GetXmin() &&
            MEtoEdmObject[j].object.GetZaxis()->GetXmax() == newMEtoEDMObject[i].object.GetZaxis()->GetXmax() &&
            CheckBinLabels((TAxis *)MEtoEdmObject[j].object.GetXaxis(),
                           (TAxis *)newMEtoEDMObject[i].object.GetXaxis()) &&
            CheckBinLabels((TAxis *)MEtoEdmObject[j].object.GetYaxis(),
                           (TAxis *)newMEtoEDMObject[i].object.GetYaxis()) &&
            CheckBinLabels((TAxis *)MEtoEdmObject[j].object.GetZaxis(),
                           (TAxis *)newMEtoEDMObject[i].object.GetZaxis())) {
          MEtoEdmObject[j].object.Add(&newMEtoEDMObject[i].object);
        } else {
          std::cout
              << "ERROR MEtoEDM::mergeProducts(): found histograms with different axis limits or different labels, '"
              << name << "' not merged" << std::endl;
#if METOEDMFORMAT_DEBUG
          std::cout << MEtoEdmObject[j].name << " " << newMEtoEDMObject[i].name << std::endl;
          std::cout << MEtoEdmObject[j].object.GetNbinsX() << " " << newMEtoEDMObject[i].object.GetNbinsX()
                    << std::endl;
          std::cout << MEtoEdmObject[j].object.GetXaxis()->GetXmin() << " "
                    << newMEtoEDMObject[i].object.GetXaxis()->GetXmin() << std::endl;
          std::cout << MEtoEdmObject[j].object.GetXaxis()->GetXmax() << " "
                    << newMEtoEDMObject[i].object.GetXaxis()->GetXmax() << std::endl;
          std::cout << MEtoEdmObject[j].object.GetNbinsY() << " " << newMEtoEDMObject[i].object.GetNbinsY()
                    << std::endl;
          std::cout << MEtoEdmObject[j].object.GetYaxis()->GetXmin() << " "
                    << newMEtoEDMObject[i].object.GetYaxis()->GetXmin() << std::endl;
          std::cout << MEtoEdmObject[j].object.GetYaxis()->GetXmax() << " "
                    << newMEtoEDMObject[i].object.GetYaxis()->GetXmax() << std::endl;
          std::cout << MEtoEdmObject[j].object.GetNbinsZ() << " " << newMEtoEDMObject[i].object.GetNbinsZ()
                    << std::endl;
          std::cout << MEtoEdmObject[j].object.GetZaxis()->GetXmin() << " "
                    << newMEtoEDMObject[i].object.GetZaxis()->GetXmin() << std::endl;
          std::cout << MEtoEdmObject[j].object.GetZaxis()->GetXmax() << " "
                    << newMEtoEDMObject[i].object.GetZaxis()->GetXmax() << std::endl;
#endif
        }
      }
    }
    return true;
  }

mergeProduct() [2/5]

template<>

bool MEtoEDM< double >::mergeProduct ( const MEtoEDM< double > &  newMEtoEDM )

inline

Definition at line 180 of file MEtoEDMFormat.h.

References gather_cfg::cout, MEtoEDM< T >::getMEtoEdmObject(), mps_fire::i, dqmiolumiharvest::j, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                           {
  const MEtoEdmObjectVector &newMEtoEDMObject = newMEtoEDM.getMEtoEdmObject();
  const size_t nObjects = newMEtoEDMObject.size();
  //  NOTE: we remember the present size since we will only add content
  //        from newMEtoEDMObject after this point
  const size_t nOldObjects = MEtoEdmObject.size();

  // if the old and new are not the same size, we want to report a problem
  if (nObjects != nOldObjects) {
    std::cout << "WARNING MEtoEDM::mergeProducts(): the lists of histograms to be merged have different sizes: new="
              << nObjects << ", old=" << nOldObjects << std::endl;
  }

  for (unsigned int i = 0; i < nObjects; ++i) {
    unsigned int j = 0;
    // see if the name is already in the old container up to the point where
    // we may have added new entries in the container
    const std::string &name = newMEtoEDMObject[i].name;
    if (i < nOldObjects && (MEtoEdmObject[i].name == name)) {
      j = i;
    } else {
      j = 0;
      while (j < nOldObjects && (MEtoEdmObject[j].name != name))
        ++j;
    }
    if (j >= nOldObjects) {
      // this value is only in the new container, not the old one
#if METOEDMFORMAT_DEBUG
      std::cout << "WARNING MEtoEDM::mergeProducts(): adding new histogram '" << name << "'" << std::endl;
#endif
      MEtoEdmObject.push_back(newMEtoEDMObject[i]);
    }
  }
  return true;
}

mergeProduct() [3/5]

template<>

bool MEtoEDM< int >::mergeProduct ( const MEtoEDM< int > &  newMEtoEDM )

inline

Definition at line 217 of file MEtoEDMFormat.h.

References gather_cfg::cout, MEtoEDM< T >::getMEtoEdmObject(), mps_fire::i, dqmiolumiharvest::j, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                     {
  const MEtoEdmObjectVector &newMEtoEDMObject = newMEtoEDM.getMEtoEdmObject();
  const size_t nObjects = newMEtoEDMObject.size();
  //  NOTE: we remember the present size since we will only add content
  //        from newMEtoEDMObject after this point
  const size_t nOldObjects = MEtoEdmObject.size();

  // if the old and new are not the same size, we want to report a problem
  if (nObjects != nOldObjects) {
    std::cout << "WARNING MEtoEDM::mergeProducts(): the lists of histograms to be merged have different sizes: new="
              << nObjects << ", old=" << nOldObjects << std::endl;
  }

  for (unsigned int i = 0; i < nObjects; ++i) {
    unsigned int j = 0;
    // see if the name is already in the old container up to the point where
    // we may have added new entries in the container
    const std::string &name = newMEtoEDMObject[i].name;
    if (i < nOldObjects && (MEtoEdmObject[i].name == name)) {
      j = i;
    } else {
      j = 0;
      while (j < nOldObjects && (MEtoEdmObject[j].name != name))
        ++j;
    }
    if (j >= nOldObjects) {
      // this value is only in the new container, not the old one
#if METOEDMFORMAT_DEBUG
      std::cout << "WARNING MEtoEDM::mergeProducts(): adding new histogram '" << name << "'" << std::endl;
#endif
      MEtoEdmObject.push_back(newMEtoEDMObject[i]);
    } else {
      // this value is also in the new container: add the two
      if (MEtoEdmObject[j].name.find("EventInfo/processedEvents") != std::string::npos) {
        MEtoEdmObject[j].object += (newMEtoEDMObject[i].object);
      }
      if (MEtoEdmObject[j].name.find("EventInfo/iEvent") != std::string::npos ||
          MEtoEdmObject[j].name.find("EventInfo/iLumiSection") != std::string::npos) {
        if (MEtoEdmObject[j].object < newMEtoEDMObject[i].object) {
          MEtoEdmObject[j].object = (newMEtoEDMObject[i].object);
        }
      }
    }
  }
  return true;
}

mergeProduct() [4/5]

template<>

bool MEtoEDM< long long >::mergeProduct ( const MEtoEDM< long long > &  newMEtoEDM )

inline

Definition at line 265 of file MEtoEDMFormat.h.

References gather_cfg::cout, MEtoEDM< T >::getMEtoEdmObject(), mps_fire::i, dqmiolumiharvest::j, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
  const MEtoEdmObjectVector &newMEtoEDMObject = newMEtoEDM.getMEtoEdmObject();
  const size_t nObjects = newMEtoEDMObject.size();
  //  NOTE: we remember the present size since we will only add content
  //        from newMEtoEDMObject after this point
  const size_t nOldObjects = MEtoEdmObject.size();

  // if the old and new are not the same size, we want to report a problem
  if (nObjects != nOldObjects) {
    std::cout << "WARNING MEtoEDM::mergeProducts(): the lists of histograms to be merged have different sizes: new="
              << nObjects << ", old=" << nOldObjects << std::endl;
  }

  for (unsigned int i = 0; i < nObjects; ++i) {
    unsigned int j = 0;
    // see if the name is already in the old container up to the point where
    // we may have added new entries in the container
    const std::string &name = newMEtoEDMObject[i].name;
    if (i < nOldObjects && (MEtoEdmObject[i].name == name)) {
      j = i;
    } else {
      j = 0;
      while (j < nOldObjects && (MEtoEdmObject[j].name != name))
        ++j;
    }
    if (j >= nOldObjects) {
      // this value is only in the new container, not the old one
#if METOEDMFORMAT_DEBUG
      std::cout << "WARNING MEtoEDM::mergeProducts(): adding new histogram '" << name << "'" << std::endl;
#endif
      MEtoEdmObject.push_back(newMEtoEDMObject[i]);
    } else {
      // this value is also in the new container: add the two
      if (MEtoEdmObject[j].name.find("EventInfo/processedEvents") != std::string::npos) {
        MEtoEdmObject[j].object += (newMEtoEDMObject[i].object);
      }
      if (MEtoEdmObject[j].name.find("EventInfo/iEvent") != std::string::npos ||
          MEtoEdmObject[j].name.find("EventInfo/iLumiSection") != std::string::npos) {
        if (MEtoEdmObject[j].object < newMEtoEDMObject[i].object) {
          MEtoEdmObject[j].object = (newMEtoEDMObject[i].object);
        }
      }
    }
  }
  return true;
}

mergeProduct() [5/5]

template<>

bool MEtoEDM< TString >::mergeProduct ( const MEtoEDM< TString > &  newMEtoEDM )

inline

Definition at line 313 of file MEtoEDMFormat.h.

References gather_cfg::cout, MEtoEDM< T >::getMEtoEdmObject(), mps_fire::i, dqmiolumiharvest::j, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                             {
  const MEtoEdmObjectVector &newMEtoEDMObject = newMEtoEDM.getMEtoEdmObject();
  const size_t nObjects = newMEtoEDMObject.size();
  //  NOTE: we remember the present size since we will only add content
  //        from newMEtoEDMObject after this point
  const size_t nOldObjects = MEtoEdmObject.size();

  // if the old and new are not the same size, we want to report a problem
  if (nObjects != nOldObjects) {
    std::cout << "WARNING MEtoEDM::mergeProducts(): the lists of histograms to be merged have different sizes: new="
              << nObjects << ", old=" << nOldObjects << std::endl;
  }

  for (unsigned int i = 0; i < nObjects; ++i) {
    unsigned int j = 0;
    // see if the name is already in the old container up to the point where
    // we may have added new entries in the container
    const std::string &name = newMEtoEDMObject[i].name;
    if (i < nOldObjects && (MEtoEdmObject[i].name == name)) {
      j = i;
    } else {
      j = 0;
      while (j < nOldObjects && (MEtoEdmObject[j].name != name))
        ++j;
    }
    if (j >= nOldObjects) {
      // this value is only in the new container, not the old one
#if METOEDMFORMAT_DEBUG
      std::cout << "WARNING MEtoEDM::mergeProducts(): adding new histogram '" << name << "'" << std::endl;
#endif
      MEtoEdmObject.push_back(newMEtoEDMObject[i]);
    }
  }
  return true;
}

putMEtoEdmObject()

template<class T>

void MEtoEDM< T >::putMEtoEdmObject ( const std::string &  name, const T &  object  )

inline

Definition at line 53 of file MEtoEDMFormat.h.

References MEtoEDM< T >::MEtoEdmObject, Skims_PA_cff::name, resolutioncreator_cfi::object, and groupFilesInBlocks::temp.

                                                                {
    typename MEtoEdmObjectVector::value_type temp;
    MEtoEdmObject.push_back(temp);
    MEtoEdmObject.back().name = name;
    MEtoEdmObject.back().object = object;
  }

swap()

template<class T>

void MEtoEDM< T >::swap ( MEtoEDM< T > &  iOther )

inline

Definition at line 172 of file MEtoEDMFormat.h.

References MEtoEDM< T >::MEtoEdmObject.

{ MEtoEdmObject.swap(iOther.MEtoEdmObject); }

Member Data Documentation

MEtoEdmObject

template<class T>

MEtoEdmObjectVector MEtoEDM< T >::MEtoEdmObject

private

Definition at line 175 of file MEtoEDMFormat.h.

Referenced by MEtoEDM< T >::getMEtoEdmObject(), MEtoEDM< T >::mergeProduct(), MEtoEDM< T >::MEtoEDM(), MEtoEDM< T >::putMEtoEdmObject(), and MEtoEDM< T >::swap().