CSCRecHit2D.cc

Go to the documentation of this file.

#include <DataFormats/CSCRecHit/interface/CSCRecHit2D.h>
#include <iostream>

CSCRecHit2D::CSCRecHit2D()
    : theTpeak(-999.),
      thePositionWithinStrip(-999.),
      theErrorWithinStrip(-999.),
      theEnergyDeposit(-994.),
      theQuality(0),
      theScaledWireTime(0),
      theBadStrip(0),
      theBadWireGroup(0),
      nStrips_(0),
      nWireGroups_(0),
      nTimeBins_(0),
      theLocalPosition(0., 0.),
      theLocalError(0., 0., 0.) {
  for (unsigned int i = 0; i < MAXSTRIPS; i++)
    theStrips_[i] = 0;
  for (unsigned int i = 0; i < MAXSTRIPS; i++)
    for (unsigned int j = 0; j < MAXTIMEBINS; j++)
      theADCs_[i * MAXTIMEBINS + j] = 0;
}

CSCRecHit2D::CSCRecHit2D(const CSCDetId& id,
                         const LocalPoint& pos,
                         const LocalError& err,
                         const ChannelContainer& channels,
                         const ADCContainer& adcs,
                         const ChannelContainer& wgroups,
                         float tpeak,
                         float posInStrip,
                         float errInStrip,
                         int quality,
                         short int badStrip,
                         short int badWireGroup,
                         int scaledWireTime,
                         float energyDeposit)
    : RecHit2DLocalPos(id),
      theTpeak(tpeak),
      thePositionWithinStrip(posInStrip),
      theErrorWithinStrip(errInStrip),
      theEnergyDeposit(energyDeposit),
      theQuality(quality),
      theScaledWireTime(scaledWireTime),
      theBadStrip(badStrip),
      theBadWireGroup(badWireGroup),
      theLocalPosition(pos),
      theLocalError(err) {
  nStrips_ = channels.size();
  nWireGroups_ = wgroups.size();

  if (nStrips_ > MAXSTRIPS) {
    std::cout << "CSCRecHit2D: not enough strips in DataFormat! " << unsigned(nStrips_) << std::endl;
    nStrips_ = MAXSTRIPS;
  }

  for (unsigned int i = 0; i < MAXSTRIPS; i++)
    theStrips_[i] = 0;
  for (unsigned int i = 0; i < MAXSTRIPS; i++)
    theL1APhaseBits_[i] = 0;
  for (unsigned int i = 0; i < MAXSTRIPS; i++)
    for (unsigned int j = 0; j < MAXTIMEBINS; j++)
      theADCs_[i * MAXTIMEBINS + j] = 0;

  for (unsigned int i = 0; i < nStrips_; i++) {
    theStrips_[i] = channels[i] & 0x000000FF;
    theL1APhaseBits_[i] = channels[i] & 0x0000FF00;
  }
  if (nWireGroups_ > 0) {
    //take only the low bits
    hitWire_ = wgroups[nWireGroups_ / 2] & 0x0000FFFF;
    theWGroupsBX_ = (wgroups[nWireGroups_ / 2] >> 16) & 0x0000FFFF;
  } else {
    hitWire_ = 0;
    theWGroupsBX_ = 0;
  }
  ADCContainer tmp(adcs);  //must be a bug in RangeMap!!!???
  nTimeBins_ = tmp.size() / nStrips_;
  unsigned int k = 0;
  for (unsigned int i = 0; i < nStrips_; i++)
    for (unsigned int j = 0; j < nTimeBins_; j++) {
      theADCs_[i * MAXTIMEBINS + j] = tmp[k];
      k++;
    }
}

CSCRecHit2D::~CSCRecHit2D() {}

bool CSCRecHit2D::sharesInput(const TrackingRecHit* other, TrackingRecHit::SharedInputType what) const {
  // This is to satisfy the base class virtual function

  // @@ Cast the enum (!) But what if the TRH::SIT changes?!
  CSCRecHit2D::SharedInputType cscWhat = static_cast<CSCRecHit2D::SharedInputType>(what);
  return sharesInput(other, cscWhat);
}

bool CSCRecHit2D::sharesInput(const TrackingRecHit* other, CSCRecHit2D::SharedInputType what) const {
  // Check to see if the TrackingRecHit is actually a CSCRecHit2D.
  if (other->geographicalId().subdetId() != MuonSubdetId::CSC)
    return false;

  // Now I can static cast, because the previous guarantees that this is a CSCRecHit2D
  const CSCRecHit2D* otherRecHit = static_cast<const CSCRecHit2D*>(other);

  return sharesInput(otherRecHit, what);
}

bool CSCRecHit2D::sharesInput(const CSCRecHit2D* otherRecHit, CSCRecHit2D::SharedInputType what) const {
  // Check to see if the geographical ID of the two are the same
  if (geographicalId() != otherRecHit->geographicalId())
    return false;

  // Trivial cases
  if (nStrips() == 0 && otherRecHit->nStrips() == 0 && nWireGroups() == 0 && otherRecHit->nWireGroups() == 0)
    return true;
  if ((what == allWires || what == someWires) && nWireGroups() == 0 && otherRecHit->nWireGroups() == 0)
    return true;
  if ((what == allStrips || what == someStrips) && nStrips() == 0 && otherRecHit->nStrips() == 0)
    return true;

  // Check to see if the wire containers are the same length
  if ((what == all || what == allWires) && nWireGroups() != otherRecHit->nWireGroups())
    return false;

  // Check to see if the strip containers are the same length
  if ((what == all || what == allStrips) && nStrips() != otherRecHit->nStrips())
    return false;

  bool foundWire = false;
  // Check to see if the wires are the same
  if (what != allStrips && what != someStrips) {
    //can we do better here?
    if (hitWire() != otherRecHit->hitWire())
      return false;
  }

  // Check to see if the wires are the same
  bool foundStrip = false;
  if (what != allWires && what != someWires) {
    for (unsigned int i = 0; i < nStrips(); i++) {
      bool found = false;
      for (unsigned int j = 0; j < nStrips(); j++) {
        //a strip is a channel for all but ME1/1a chambers (where 3 ganged strips are a channel)
        if (cscDetId().channel(channels(i)) == otherRecHit->cscDetId().channel(otherRecHit->channels(j))) {
          if (what == some || what == someStrips)
            return true;
          else {
            found = true;
            foundStrip = true;
            break;
          }
        }
      }
      if ((what == all || what == allStrips) && !found)
        return false;
    }
    if (what == someStrips && !foundStrip)
      return false;
  }

  // In case we were looking for "some" and found absolutely nothing.
  if (!foundWire && !foundStrip)
    return false;

  // If we made it this far, then:
  //  1) the detector IDs are the same
  //  2) the channel containers have the same number of entries
  //  3) for each entry in my channel container, I can find the same value in the other RecHit's corresponding channel container
  // I think that means we are the same.
  return true;
}

void CSCRecHit2D::print() const {
  std::cout << "CSCRecHit in CSC Detector: " << cscDetId() << std::endl;
  std::cout << "  local x = " << localPosition().x() << " +/- " << sqrt(localPositionError().xx())
            << " y = " << localPosition().y() << " +/- " << sqrt(localPositionError().yy()) << std::endl;

  std::cout << " tpeak " << theTpeak << " psoInStrip " << thePositionWithinStrip << " errorinstrip "
            << theErrorWithinStrip << " "
            << " qual " << theQuality << " wiretime " << theScaledWireTime << " tbs " << theBadStrip << " bwg "
            << theBadWireGroup << std::endl;

  std::cout << "  Channels: ";
  for (unsigned int i = 0; i < nStrips(); i++) {
    std::cout << std::dec << channels(i) << " "
              << " ("
              << "HEX: " << std::hex << channels(i) << ")"
              << " ";
  }
  std::cout << std::endl;

  std::cout << "  L1APhase: ";
  for (int i = 0; i < (int)nStrips(); i++) {
    std::cout << "|";
    for (int k = 0; k < 8; k++) {
      std::cout << ((channelsl1a(i) >> (15 - k)) & 0x1) << " ";
    }
    std::cout << "| ";
  }
  std::cout << std::endl;

  std::cout << "nWireGroups " << (int)nWireGroups() << " central wire " << hitWire_ << std::endl;
}

std::ostream& operator<<(std::ostream& os, const CSCRecHit2D& rh) {
  os << "CSCRecHit2D: "
     << "local x: " << rh.localPosition().x() << " +/- " << sqrt(rh.localPositionError().xx())
     << " y: " << rh.localPosition().y() << " +/- " << sqrt(rh.localPositionError().yy())
     << " in strip X: " << rh.positionWithinStrip() << " +/- " << rh.errorWithinStrip() << " quality: " << rh.quality()
     << " tpeak: " << rh.tpeak() << " wireTime: " << rh.wireTime() << std::endl;
  os << "strips: ";
  for (size_t iS = 0; iS < rh.nStrips(); ++iS) {
    os << rh.channels(iS) << "  ";
  }
  int nwgs = rh.nWireGroups();
  if (nwgs == 1) {
    os << "central wire: " << rh.hitWire() << " of " << nwgs << " wiregroup" << std::endl;
  } else {
    os << "central wire: " << rh.hitWire() << " of " << nwgs << " wiregroups" << std::endl;
  }
  return os;
}