CSCGangedWireGrouping.cc

Go to the documentation of this file.

#include <Geometry/CSCGeometry/src/CSCGangedWireGrouping.h>

#include <FWCore/MessageLogger/interface/MessageLogger.h>

#include <iostream>
#include <iomanip>
#include <algorithm>
#include <cstddef>

CSCGangedWireGrouping::CSCGangedWireGrouping(const Container& consecutiveGroups,
                                             const Container& wiresInConsecutiveGroups,
                                             int numberOfGroups)
    : theNumberOfGroups(numberOfGroups) {
  int countGroups = 0;       // counter for no. of wire groups
  int firstWire = 1;         // (virtual) wire number which is 1st in a group
  int countConsecutive = 0;  // counter for the sections in DDD

  for (int igs : consecutiveGroups) {
    if (igs != 0) {
      // igs is number of consecutive groups each containing
      // an identical number of wires
      countGroups += igs;
      for (int ic = 0; ic != igs; ++ic) {
        theFirstWireOfEachWireGroup.emplace_back(firstWire);
        int wiresInGroup = wiresInConsecutiveGroups[countConsecutive];
        theNumberOfWiresPerWireGroup.emplace_back(wiresInGroup);
        firstWire += wiresInGroup;  // ready for next group
      }
    } else {
      // zero means a gap - just add in the no. of virtual wires in gap
      firstWire += wiresInConsecutiveGroups[countConsecutive];
    }
    ++countConsecutive;  // ready for next set of consecutive groups
  }

  theNumberOfWires = firstWire - 1;  // at end of loop we're on 1st for next

  if (countGroups != numberOfGroups) {
    edm::LogError("CSC") << "CSCGangedWireGrouping: ERROR in parsing wire info from DDD..."
                         << "\n";
    edm::LogError("CSC") << "groups expected = " << numberOfGroups << " groups seen = " << countGroups << "\n";
    edm::LogError("CSC") << "Please report this error to Tim.Cox@cern.ch"
                         << "\n";
  }

  LogTrace("CSCWireGeometry|CSC") << "CSCGangedWireGrouping constructor complete,"
                                  << " wire group list follows... ";
  LogTrace("CSCWireGeometry|CSC") << "Size of group buffers = " << theFirstWireOfEachWireGroup.size() << " and "
                                  << theNumberOfWiresPerWireGroup.size();
  LogTrace("CSCWireGeometry|CSC") << " wg#    1st wire  #wires";
  for (size_t i = 0; i != theFirstWireOfEachWireGroup.size(); ++i) {
    LogTrace("CSCWireGeometry|CSC") << std::setw(4) << i + 1 << std::setw(12) << theFirstWireOfEachWireGroup[i]
                                    << std::setw(8) << theNumberOfWiresPerWireGroup[i];
  }
  LogTrace("CSCWireGeometry|CSC") << "Total no. of wires = " << theNumberOfWires;
}

int CSCGangedWireGrouping::numberOfWiresPerGroup(int wireGroup) const {
  if (wireGroup > 0 && wireGroup <= theNumberOfGroups)
    return theNumberOfWiresPerWireGroup[wireGroup - 1];
  else
    return 0;
}

int CSCGangedWireGrouping::wireGroup(int wire) const {
  // Return wire group number (start at 1) containing (virtual) 'wire'
  // Return 0 if 'wire' is not in a wire group; this includes
  // wires out outside the chamber, and lying in 'dead' regions

  int wireG = 0;

  // upper_bound works on a sorted range and points to first element
  // _succeeding_ supplied value.

  CIterator it = upper_bound(theFirstWireOfEachWireGroup.begin(), theFirstWireOfEachWireGroup.end(), wire);

  // We are now pointing to the wire group _after_ the required one
  // (unless we are at begin() or end() when we just return wireG=0)
  ptrdiff_t pd = it - theFirstWireOfEachWireGroup.begin();

  if (pd > 0) {
    size_t id = --pd;  //@@ Need to step back one. IS THIS SANE CODE?
    int wiresInGroup = theNumberOfWiresPerWireGroup[id];
    int firstWire = theFirstWireOfEachWireGroup[id];

    // Require wire not past end of group (may be in a dead region, or
    // bigger than total wires in chamber)
    if (wire < (firstWire + wiresInGroup))
      wireG = ++id;
  }
  return wireG;
}

float CSCGangedWireGrouping::middleWireOfGroup(int wireGroup) const {
  // Return exact wire number for group with odd number of wires
  // Return half-integer wire number for group with even number of wires
  // i.e. first + 0.5 * float(#) - 0.5
  // Return 0 if out of range

  float middleWire = 0.;
  if (wireGroup > 0 && wireGroup <= theNumberOfGroups) {
    middleWire = theFirstWireOfEachWireGroup[wireGroup - 1] + 0.5 * theNumberOfWiresPerWireGroup[wireGroup - 1] - 0.5;
  }
  return middleWire;
}