RecoCollectionFP420.cc

Go to the documentation of this file.

// File: RecoCollectionFP420.cc
// Date: 11.2007
// Description: RecoCollectionFP420 for FP420
// Modifications:
#include "DataFormats/FP420Cluster/interface/RecoCollectionFP420.h"
#include <algorithm>
#include <iostream>
using namespace std;
//#define mydigidebug

void RecoCollectionFP420::put(RecoCollectionFP420::Range input, unsigned int stationID) {
  // put in RecHits of stationID

  // store size of vector before put
  RecoCollectionFP420::IndexRange inputRange;

  // put in RecoFP420s from input
  bool first = true;

  // fill input in temporary vector for sorting
  std::vector<RecoFP420> temporary;
  RecoCollectionFP420::ContainerIterator sort_begin = input.first;
  RecoCollectionFP420::ContainerIterator sort_end = input.second;
#ifdef mydigidebug
  std::cout << "   !!!!!!!!!!!!!!!!    RecoCollectionFP420:: !!!!  put !!!!           start " << std::endl;
#endif
  for (; sort_begin != sort_end; ++sort_begin) {
#ifdef mydigidebug
    std::cout << "put: temporary.push_back " << std::endl;
#endif
    temporary.push_back(*sort_begin);
  }
  std::sort(temporary.begin(), temporary.end());

  // iterators over input
  RecoCollectionFP420::ContainerIterator begin = temporary.begin();
  RecoCollectionFP420::ContainerIterator end = temporary.end();
  for (; begin != end; ++begin) {
    container_.push_back(*begin);
    if (first) {
      inputRange.first = container_.size() - 1;
      //     inputRange.first = container_.size();
      first = false;
    }
  }

  // since we start from 0, then the last element will be size-1
  if (!container_.empty()) {
    inputRange.second = container_.size() - 1;
  } else {
    inputRange.first = container_.size();
    inputRange.second = container_.size();
  }

  //  inputRange.second = container_.size()-1;
  //  //inputRange.second = container_.size();

#ifdef mydigidebug
  std::cout << "put: container_.size() = " << container_.size() << std::endl;
  std::cout << "put:  inputRange.first = " << inputRange.first << std::endl;
  std::cout << "put:  inputRange.second = " << inputRange.second << std::endl;
#endif

  // fill map
  map_[stationID] = inputRange;
}

const RecoCollectionFP420::Range RecoCollectionFP420::get(unsigned int stationID) const {
  // get RecHits of stationID

#ifdef mydigidebug
  std::cout << "RecoCollectionFP420::get:stationID= " << stationID << std::endl;
#endif
  RecoCollectionFP420::RegistryIterator returnIndex = map_.find(stationID);
  RecoCollectionFP420::IndexRange returnIndexRange = returnIndex->second;
#ifdef mydigidebug
  std::cout << "RecoCollectionFP420::get1: returnIndexRange.first= " << returnIndexRange.first << std::endl;
  std::cout << "RecoCollectionFP420::get1: returnIndexRange.second= " << returnIndexRange.second << std::endl;
#endif

  RecoCollectionFP420::Range returnRange;
  returnRange.first = container_.begin() + returnIndexRange.first;
  if (returnIndexRange.second != 0) {
    returnRange.second = container_.begin() + returnIndexRange.second + 1;
  } else {
    returnRange.second = container_.begin() + returnIndexRange.second;
  }
#ifdef mydigidebug
  std::cout << "RecoCollectionFP420::get2: container_.size() = " << container_.size() << std::endl;
  std::cout << "RecoCollectionFP420::get2: returnIndexRange.first= " << returnIndexRange.first << std::endl;
  std::cout << "RecoCollectionFP420::get2: returnIndexRange.second= " << returnIndexRange.second << std::endl;
#endif
  return returnRange;
}

void RecoCollectionFP420::clear() { container_.clear(); }
void RecoCollectionFP420::putclear(RecoCollectionFP420::Range input, unsigned int stationID) {
  RecoCollectionFP420::IndexRange inputRange;

  std::vector<RecoFP420> temporary;
  RecoCollectionFP420::ContainerIterator sort_begin = input.first;
  RecoCollectionFP420::ContainerIterator sort_end = input.second;
  for (; sort_begin != sort_end; ++sort_begin) {
    temporary.push_back(*sort_begin);
  }
  std::sort(temporary.begin(), temporary.end());

  //    temporary.clear();
  RecoCollectionFP420::ContainerIterator begin = temporary.begin();
  RecoCollectionFP420::ContainerIterator end = temporary.end();
  for (; begin != end; ++begin) {
    container_.push_back(*begin);
  }
  //container_.clear();
  inputRange.first = container_.size() - container_.size();
  inputRange.second = container_.size() - container_.size();

#ifdef mydigidebug
  std::cout << "putclear: container_.size() = " << container_.size() << std::endl;
  std::cout << "putclear:  inputRange.first = " << inputRange.first << std::endl;
  std::cout << "putclear:  inputRange.second = " << inputRange.second << std::endl;
#endif

  map_[stationID] = inputRange;
}

const std::vector<unsigned int> RecoCollectionFP420::stationIDs() const {
  // returns vector of stationIDs in map

#ifdef mydigidebug
  std::cout << "RecoCollectionFP420::stationIDs:start " << std::endl;
#endif
  RecoCollectionFP420::RegistryIterator begin = map_.begin();
  RecoCollectionFP420::RegistryIterator end = map_.end();

  std::vector<unsigned int> output;

  for (; begin != end; ++begin) {
    output.push_back(begin->first);
  }

  return output;
}