CSCHitFromStripOnly Class Reference

#include <CSCHitFromStripOnly.h>

Public Types
typedef std::array< CSCStripData, 100 >	PulseHeightMap

Public Member Functions
	CSCHitFromStripOnly (const edm::ParameterSet &ps)
bool	ganged ()
std::vector< CSCStripHit >	runStrip (const CSCDetId &id, const CSCLayer *layer, const CSCStripDigiCollection::Range &rstripd)
void	setConditions (const CSCRecoConditions *reco)
void	setGanged (bool ig)
	~CSCHitFromStripOnly ()

Private Member Functions
void	fillPulseHeights (const CSCStripDigiCollection::Range &rstripd)
	Store SCA pulseheight information from strips in digis of one layer. More...
float	findHitOnStripPosition (const std::vector< CSCStripHitData > &data, const int &centerStrip)
	Find position of hit in strip cluster in terms of strip #. More...
void	findMaxima (const CSCDetId &id)
	Find local maxima. More...
bool	isDeadStrip (const CSCDetId &id, int centralStrip, int nstrips)
	Is the strip 'bad'? More...
bool	isNearDeadStrip (const CSCDetId &id, int centralStrip, int nstrips)
	Is either neighbour 'bad'? More...
bool	isPeakOK (int iStrip, float heightCluster)
float	makeCluster (int centerStrip)
	Make clusters using local maxima. More...
CSCStripHitData	makeStripData (int centerStrip, int offset)

Private Attributes
CSCPedestalChoice *	calcped_
int	clusterSize
float	gainWeight [80]
bool	ganged_
CSCDetId	id_
const CSCLayer *	layer_
unsigned	nstrips_
const CSCRecoConditions *	recoConditions_
std::vector< float >	strips_adc
std::vector< float >	strips_adcRaw
std::vector< int >	theClosestMaximum
std::vector< int >	theConsecutiveStrips
std::vector< int >	theMaxima
PulseHeightMap	thePulseHeightMap
std::vector< int >	theStrips
float	theThresholdForAPeak
float	theThresholdForCluster
int	tmax_cluster
bool	useCalib

Static Private Attributes
static const int	theClusterSize = 3

Detailed Description

below)

Search for strip with ADC output exceeding theThresholdForAPeak. For each of these strips, build a cluster of strip of size theClusterSize (typically 5 strips). Finally, make a Strip Hit out of these clusters by finding the center-of-mass position of the hit The DetId, strip hit position, and peaking time are stored in a CSCStripHit collection.

Be careful with the ME_1/a chambers: the 48 strips are ganged into 16 channels, Each channel contains signals from three strips, each separated by 16 strips from the next.

Definition at line 34 of file CSCHitFromStripOnly.h.

Member Typedef Documentation

PulseHeightMap

typedef std::array<CSCStripData, 100> CSCHitFromStripOnly::PulseHeightMap

Definition at line 36 of file CSCHitFromStripOnly.h.

Constructor & Destructor Documentation

CSCHitFromStripOnly()

CSCHitFromStripOnly::CSCHitFromStripOnly ( const edm::ParameterSet &  ps )

explicit

Definition at line 22 of file CSCHitFromStripOnly.cc.

References calcped_, edm::ParameterSet::getParameter(), LogTrace, theThresholdForAPeak, theThresholdForCluster, and useCalib.

    : recoConditions_(nullptr), calcped_(nullptr), ganged_(false) {
  useCalib = ps.getParameter<bool>("CSCUseCalibrations");
  bool useStaticPedestals = ps.getParameter<bool>("CSCUseStaticPedestals");
  int noOfTimeBinsForDynamicPed = ps.getParameter<int>("CSCNoOfTimeBinsForDynamicPedestal");

  theThresholdForAPeak = ps.getParameter<double>("CSCStripPeakThreshold");
  theThresholdForCluster = ps.getParameter<double>("CSCStripClusterChargeCut");

  LogTrace("CSCRecHit") << "[CSCHitFromStripOnly] CSCUseStaticPedestals = " << useStaticPedestals;
  if (!useStaticPedestals)
    LogTrace("CSCRecHit") << "[CSCHitFromStripOnly] CSCNoOfTimeBinsForDynamicPedestal = " << noOfTimeBinsForDynamicPed;

  if (useStaticPedestals) {
    calcped_ = new CSCStaticPedestal();
  } else {
    if (noOfTimeBinsForDynamicPed == 1) {
      calcped_ = new CSCDynamicPedestal1();
    } else {
      calcped_ = new CSCDynamicPedestal2();  // NORMAL DEFAULT!
    }
  }
}

~CSCHitFromStripOnly()

CSCHitFromStripOnly::~CSCHitFromStripOnly

(

)

Definition at line 46 of file CSCHitFromStripOnly.cc.

References calcped_.

{ delete calcped_; }

Member Function Documentation

fillPulseHeights()

void CSCHitFromStripOnly::fillPulseHeights ( const CSCStripDigiCollection::Range &  rstripd )

private

Store SCA pulseheight information from strips in digis of one layer.

Definition at line 350 of file CSCHitFromStripOnly.cc.

References cms::cuda::assert(), calcped_, filterCSVwithJSON::copy, gainWeight, ganged(), id_, dqmiolumiharvest::j, CSCPedestalChoice::pedestal(), recoConditions_, thePulseHeightMap, tmax, and useCalib.

Referenced by runStrip().

                                                                                     {
  // Loop over strip digis in one CSCLayer and fill PulseHeightMap with pedestal-subtracted
  // SCA pulse heights.

  for (auto& ph : thePulseHeightMap)
    ph.reset();

  // for storing sca pulseheights once they may no longer be integer (e.g. after ped subtraction)
  for (CSCStripDigiCollection::const_iterator it = rstripd.first; it != rstripd.second; ++it) {
    int thisChannel = (*it).getStrip();
    auto& stripData = thePulseHeightMap[thisChannel - 1];
    auto& scaRaw = stripData.phRaw_;
    auto& sca = stripData.ph_;

    auto const& scaOri = (*it).getADCCounts();
    assert(scaOri.size() == 8);
    // Fill sca from scaRaw, implicitly converting to float
    std::copy(scaOri.begin(), scaOri.end(), scaRaw.begin());
    std::copy(scaRaw.begin(), scaRaw.end(), sca.begin());

    //@@ Find bin with largest pulseheight (_before_ ped subtraction - shouldn't matter, right?)
    int tmax = std::max_element(sca.begin(), sca.end()) - sca.begin();  // counts from 0

    // get pedestal - calculated as appropriate - for this sca pulse
    float ped = calcped_->pedestal(sca, recoConditions_, id_, thisChannel);

    // subtract the pedestal (from BOTH sets of sca pulseheights)
    std::for_each(sca.begin(), sca.end(), CSCSubtractPedestal(ped));
    std::for_each(scaRaw.begin(), scaRaw.end(), CSCSubtractPedestal(ped));

    //@@ Max in first 3 or last time bins is unacceptable, if so set to zero (why?)
    float phmax = 0.f;
    if (tmax > 2 && tmax < 7) {
      phmax = sca[tmax];
    }
    stripData.phmax_ = phmax;
    stripData.tmax_ = tmax;

    // Fill the map, possibly apply gains from cond data, and unfold ME1A channels
    // (To apply gains use CSCStripData::op*= which scales only the non-raw sca ph's;
    // but note that both sca & scaRaw are pedestal-subtracted.)

    // From StripDigi, thisChannel labels strip channel. Values phmax, tmax, scaRaw, sca belong to thisChannel
    if (useCalib)
      stripData *= gainWeight[thisChannel - 1];

    // for ganged ME1a need to duplicate values on istrip=thisChannel to iStrip+16 and iStrip+32
    if (ganged()) {
      for (int j = 1; j < 3; ++j) {
        thePulseHeightMap[thisChannel - 1 + 16 * j] = stripData;
      }
    }
  }
}

findHitOnStripPosition()

float CSCHitFromStripOnly::findHitOnStripPosition ( const std::vector< CSCStripHitData > &  data, const int &  centerStrip  )

private

Find position of hit in strip cluster in terms of strip #.

Definition at line 546 of file CSCHitFromStripOnly.cc.

References filterCSVwithJSON::copy, data, mps_fire::i, LogTrace, digitizers_cfi::strip, strips_adc, strips_adcRaw, and w().

Referenced by makeCluster().

                                                                                                                {
  float strippos = -1.;

  if (data.empty())
    return strippos;

  // biggestStrip is strip with largest pulse height
  // Use pointer subtraction

  int biggestStrip = max_element(data.begin(), data.end()) - data.begin();
  strippos = data[biggestStrip].strip() * 1.;

  // If more than one strip:  use centroid to find center of cluster
  // but only use time bin == tmax (otherwise, bias centroid).
  float sum = 0.;
  float sum_w = 0.;

  //  std::vector<float> w(4);
  // std::vector<float> wRaw(4);

  for (size_t i = 0; i != data.size(); ++i) {
    auto const& w = data[i].ph();
    auto const& wRaw = data[i].phRaw();

    // (Require ADC to be > 0.)
    // No later studies suggest that this only do harm
    /*
    for ( size_t j = 0; j != w.size(); ++j ) {
      if ( w[j] < 0. ) w[j] = 0.001;
    }
    */

    // Fill the data members
    std::copy(w.begin(), w.end(), std::back_inserter(strips_adc));
    std::copy(wRaw.begin(), wRaw.end(), std::back_inserter(strips_adcRaw));

    if (data[i].strip() < 1) {
      LogTrace("CSCRecHit") << "[CSCHitFromStripOnly::findHitOnStripPosition] problem in indexing of strip, strip= "
                            << data[i].strip();
    }
    sum_w += w[1];
    sum += w[1] * data[i].strip();
  }

  if (sum_w > 0.)
    strippos = sum / sum_w;

  return strippos;
}

findMaxima()

void CSCHitFromStripOnly::findMaxima ( const CSCDetId &  id )

private

Find local maxima.

Definition at line 411 of file CSCHitFromStripOnly.cc.

References mps_fire::i, isDeadStrip(), isPeakOK(), dqmiolumiharvest::j, dqmdumpme::k, cmsLHEtoEOSManager::l, nstrips_, theClosestMaximum, theConsecutiveStrips, theMaxima, and thePulseHeightMap.

Referenced by runStrip().

                                                       {
  theMaxima.clear();
  theConsecutiveStrips.clear();
  theClosestMaximum.clear();
  for (size_t i = 0; i != thePulseHeightMap.size(); ++i) {
    // sum 3 strips so that hits between strips are not suppressed
    float heightCluster = 0.;

    bool maximumFound = false;
    // Left edge of chamber
    if (!isDeadStrip(id, i + 1, nstrips_)) {  // Is it i or i+1
      if (i == 0) {
        heightCluster = thePulseHeightMap[i].phmax() + thePulseHeightMap[i + 1].phmax();
        // Have found a strip Hit if...
        if (thePulseHeightMap[i].phmax() >= thePulseHeightMap[i + 1].phmax() && isPeakOK(i, heightCluster)) {
          maximumFound = true;
        }
        // Right edge of chamber
      } else if (i == thePulseHeightMap.size() - 1) {
        heightCluster = thePulseHeightMap[i - 1].phmax() + thePulseHeightMap[i].phmax();
        // Have found a strip Hit if...
        if (thePulseHeightMap[i].phmax() > thePulseHeightMap[i - 1].phmax() && isPeakOK(i, heightCluster)) {
          maximumFound = true;
        }
        // Any other strips
      } else {
        heightCluster =
            thePulseHeightMap[i - 1].phmax() + thePulseHeightMap[i].phmax() + thePulseHeightMap[i + 1].phmax();
        // Have found a strip Hit if...
        if (thePulseHeightMap[i].phmax() > thePulseHeightMap[i - 1].phmax() &&
            thePulseHeightMap[i].phmax() >= thePulseHeightMap[i + 1].phmax() && isPeakOK(i, heightCluster)) {
          maximumFound = true;
        }
      }
    }
    //---- Consecutive strips with charge (real cluster); if too wide - measurement is not accurate
    if (maximumFound) {
      int numberOfConsecutiveStrips = 1;
      float testThreshold = 10.;  //---- ADC counts;
                                  //---- this is not XTalk corrected so it is correct in first approximation only
      int j = 0;
      for (int l = 0; l < 8; ++l) {
        if (j < 0)
          edm::LogWarning("FailedStripCountingWrongConsecutiveStripNumber")
              << "This should never occur!!! Contact CSC expert!";
        ++j;
        bool signalPresent = false;
        for (int k = 0; k < 2; ++k) {
          j *= -1;  //---- check from left and right
          int anotherConsecutiveStrip = i + j;
          if (anotherConsecutiveStrip >= 0 && anotherConsecutiveStrip < int(thePulseHeightMap.size())) {
            if (thePulseHeightMap[anotherConsecutiveStrip].phmax() > testThreshold) {
              ++numberOfConsecutiveStrips;
              signalPresent = true;
            }
          }
        }
        if (!signalPresent) {
          break;
        }
      }

      bool additional_maxima_found = false;
      // search for additional maxima if:
      // - hit is closer than 3 strips from the edge
      // - enough consecutive strips with signal
      // - strip charge distribution looks abnormal

      if (i > 2 && i + 3 < thePulseHeightMap.size() && numberOfConsecutiveStrips > 3) {
        //try to look for additional maxima at the left side from the main maxima

        if (((thePulseHeightMap[i + 1].phmax() >= thePulseHeightMap[i - 1].phmax() &&
              thePulseHeightMap[i + 1].phmax() >= thePulseHeightMap[i - 2].phmax() &&
              thePulseHeightMap[i + 2].phmax() <= thePulseHeightMap[i - 2].phmax()) ||
             (thePulseHeightMap[i + 1].phmax() <= thePulseHeightMap[i - 1].phmax() &&
              thePulseHeightMap[i + 1].phmax() <= thePulseHeightMap[i - 2].phmax())) &&
            //to avoid close maxima delimitation (this is already present in the code)
            thePulseHeightMap[i - 1].phmax() >= thePulseHeightMap[i - 2].phmax() &&
            //no need in a small charge maxima (might need adjustment)
            thePulseHeightMap[i - 2].phmax() > 20) {
          additional_maxima_found = true;
          theMaxima.push_back(i - 2);  //insert left maxima first
          //insert the same number of cosecutive strips, because they belong to both maximas
          theConsecutiveStrips.push_back(numberOfConsecutiveStrips);
          theMaxima.push_back(i);  //insert main maxima
          //insert the same number of cosecutive strips, because they belong to both maximas
          theConsecutiveStrips.push_back(numberOfConsecutiveStrips);

        }  //looking for additional maxima on the left

        //try to look for additional maxima at the right side from the main maxima
        if (((thePulseHeightMap[i + 1].phmax() >= thePulseHeightMap[i - 1].phmax() &&
              thePulseHeightMap[i + 2].phmax() >= thePulseHeightMap[i - 1].phmax()) ||
             (thePulseHeightMap[i + 1].phmax() <= thePulseHeightMap[i - 1].phmax() &&
              thePulseHeightMap[i + 2].phmax() <= thePulseHeightMap[i - 1].phmax() &&
              thePulseHeightMap[i + 2].phmax() >= thePulseHeightMap[i - 2].phmax())) &&
            //to avoid close maxima delimitation (this is already present in the code)
            thePulseHeightMap[i + 1].phmax() >= thePulseHeightMap[i + 2].phmax() &&
            //no need in a small charge maxima (might need adjustment)
            thePulseHeightMap[i + 2].phmax() > 20) {
          additional_maxima_found = true;
          theMaxima.push_back(i);  //insert main maxima first
          //insert the same number of cosecutive strips, because they belong to both maximas
          theConsecutiveStrips.push_back(numberOfConsecutiveStrips);
          theMaxima.push_back(i + 2);  //insert right maxima
          //insert the same number of cosecutive strips, because they belong to both maximas
          theConsecutiveStrips.push_back(numberOfConsecutiveStrips);

        }  //looking for additional maxima on the right

        //if nothing additional found fill the maxima
        if (!additional_maxima_found) {
          theMaxima.push_back(i);
          theConsecutiveStrips.push_back(numberOfConsecutiveStrips);
        }
      } else {  //not the case for looking for the additional maxima
        theMaxima.push_back(i);
        theConsecutiveStrips.push_back(numberOfConsecutiveStrips);
      }
    }  //if maximafound
  }    //all pulses
}  //find maxima procedure

ganged()

bool CSCHitFromStripOnly::ganged ( )

inline

Definition at line 48 of file CSCHitFromStripOnly.h.

References ganged_.

Referenced by fillPulseHeights(), and runStrip().

{ return ganged_; }

isDeadStrip()

bool CSCHitFromStripOnly::isDeadStrip ( const CSCDetId &  id, int  centralStrip, int  nstrips  )

private

Is the strip 'bad'?

Definition at line 602 of file CSCHitFromStripOnly.cc.

References CSCRecoConditions::badStrip(), and recoConditions_.

Referenced by findMaxima(), and runStrip().

                                                                                       {
  return recoConditions_->badStrip(id, centralStrip, nstrips);
}

isNearDeadStrip()

bool CSCHitFromStripOnly::isNearDeadStrip ( const CSCDetId &  id, int  centralStrip, int  nstrips  )

private

Is either neighbour 'bad'?

Definition at line 596 of file CSCHitFromStripOnly.cc.

References CSCRecoConditions::nearBadStrip(), and recoConditions_.

                                                                                           {
  //@@ Tim says: not sure I understand this properly... but just moved code to CSCRecoConditions
  // where it can handle the conversion from strip to channel etc.
  return recoConditions_->nearBadStrip(id, centralStrip, nstrips);
}

isPeakOK()

bool CSCHitFromStripOnly::isPeakOK ( int  iStrip, float  heightCluster  )

private

Definition at line 534 of file CSCHitFromStripOnly.cc.

References mps_fire::i, thePulseHeightMap, theThresholdForAPeak, and theThresholdForCluster.

Referenced by findMaxima().

                                                                  {
  int i = iStrip;
  bool peakOK = (thePulseHeightMap[i].phmax() > theThresholdForAPeak && heightCluster > theThresholdForCluster &&
                 // ... and proper peak time; note that the values below are used elsewhere in this file;
                 // they should become parameters or at least constants defined in appropriate place
                 thePulseHeightMap[i].tmax() > 2 && thePulseHeightMap[i].tmax() < 7);
  return peakOK;
}

makeCluster()

float CSCHitFromStripOnly::makeCluster ( int  centerStrip )

private

Make clusters using local maxima.

Definition at line 211 of file CSCHitFromStripOnly.cc.

References clusterSize, data, findHitOnStripPosition(), mps_fire::i, LogTrace, makeStripData(), nstrips_, theClusterSize, and theStrips.

Referenced by runStrip().

                                                      {
  float strippos = -1.;
  clusterSize = theClusterSize;
  std::vector<CSCStripHitData> stripDataV;

  // We only want to use strip position in terms of strip # for the strip hit. //@@ What other choice is there?

  // If the cluster size is such that you go beyond the edge of detector, shrink cluster appropriately
  for (int i = 1; i < theClusterSize / 2 + 1; ++i) {
    if (centerStrip - i < 1 || centerStrip + i > int(nstrips_)) {
      // Shrink cluster size, but keep it an odd number of strips.
      clusterSize = 2 * i - 1;
    }
  }
  for (int i = -clusterSize / 2; i <= clusterSize / 2; ++i) {
    CSCStripHitData data = makeStripData(centerStrip, i);
    stripDataV.push_back(data);
    theStrips.push_back(centerStrip + i);
  }
  strippos = findHitOnStripPosition(stripDataV, centerStrip);

  LogTrace("CSCHitFromStripOnly") << "[CSCHitFromStripOnly::makeCluster] centerStrip= " << centerStrip
                                  << " strippos=" << strippos;

  return strippos;
}

makeStripData()

CSCStripHitData CSCHitFromStripOnly::makeStripData ( int  centerStrip, int  offset  )

private

makeStripData

Definition at line 241 of file CSCHitFromStripOnly.cc.

References gpuClustering::adc, clusterSize, filterCSVwithJSON::copy, spr::find(), mps_fire::i, dqmdumpme::k, LogTrace, SiStripPI::min, hltrates_dqm_sourceclient-live_cfg::offset, Validation_hcalonly_cfi::sign, theMaxima, thePulseHeightMap, tmax, tmax_cluster, and validateGeometry_cfg::valid.

Referenced by makeCluster().

                                                                              {
  CSCStripHitData prelimData;
  int thisStrip = centerStrip + offset;

  int tmax = thePulseHeightMap[centerStrip - 1].tmax();
  tmax_cluster = tmax;

  std::vector<float> adc(4);
  std::vector<float> adcRaw(4);

  // Fill adc & adcRaw

  int istart = tmax - 1;
  int istop = std::min(tmax + 2, 7);  // there are only time bins 0-7
  adc[3] = 0.1;                       // in case it isn't filled

  if (tmax > 2 && tmax < 7) {  // for time bins 3-6
    int ibin = thisStrip - 1;
    if (thePulseHeightMap[ibin].valid()) {
      std::copy(
          thePulseHeightMap[ibin].ph().begin() + istart, thePulseHeightMap[ibin].ph().begin() + istop + 1, adc.begin());

      std::copy(thePulseHeightMap[ibin].phRaw().begin() + istart,
                thePulseHeightMap[ibin].phRaw().begin() + istop + 1,
                adcRaw.begin());
    }
  } else {
    adc[0] = 0.1;
    adc[1] = 0.1;
    adc[2] = 0.1;
    adc[3] = 0.1;
    adcRaw = adc;
    LogTrace("CSCRecHit") << "[CSCHitFromStripOnly::makeStripData] Tmax out of range: contact CSC expert!";
  }

  if (offset == 0) {
    prelimData = CSCStripHitData(thisStrip, tmax_cluster, adcRaw, adc);
  } else {
    int sign = offset > 0 ? 1 : -1;
    // If there's another maximum that would like to use part of this cluster,
    // it gets shared in proportion to the height of the maxima
    for (int i = 1; i <= clusterSize / 2; ++i) {
      // Find the direction of the offset
      int testStrip = thisStrip + sign * i;
      std::vector<int>::iterator otherMax = find(theMaxima.begin(), theMaxima.end(), testStrip - 1);

      // No other maxima found, so just store
      if (otherMax == theMaxima.end()) {
        prelimData = CSCStripHitData(thisStrip, tmax_cluster, adcRaw, adc);
      } else {
        // Another maximum found - share
        std::vector<float> adc1(4);
        std::vector<float> adcRaw1(4);
        std::vector<float> adc2(4);
        std::vector<float> adcRaw2(4);
        // In case we only copy (below) into 3 of the 4 bins i.e. when istart=5, istop=7
        adc1[3] = 0.1;
        adc2[3] = 0.1;
        adcRaw1[3] = 0.1;
        adcRaw2[3] = 0.1;

        // Fill adcN with content of time bins tmax-1 to tmax+2 (if it exists!)
        if (tmax > 2 && tmax < 7) {  // for time bin tmax from 3-6
          int ibin = testStrip - 1;
          int jbin = centerStrip - 1;
          if (thePulseHeightMap[ibin].valid()) {
            std::copy(thePulseHeightMap[ibin].ph().begin() + istart,
                      thePulseHeightMap[ibin].ph().begin() + istop + 1,
                      adc1.begin());
            std::copy(thePulseHeightMap[ibin].phRaw().begin() + istart,
                      thePulseHeightMap[ibin].phRaw().begin() + istop + 1,
                      adcRaw1.begin());
          }

          if (thePulseHeightMap[jbin].valid()) {
            std::copy(thePulseHeightMap[jbin].ph().begin() + istart,
                      thePulseHeightMap[jbin].ph().begin() + istop + 1,
                      adc2.begin());

            std::copy(thePulseHeightMap[jbin].phRaw().begin() + istart,
                      thePulseHeightMap[jbin].phRaw().begin() + istop + 1,
                      adcRaw2.begin());
          }
        } else {
          adc1.assign(4, 0.1);
          adcRaw1 = adc1;
          adc2.assign(4, 0.1);
          adcRaw2 = adc2;
        }

        // Scale shared strip B ('adc') by ratio of peak of ADC counts from central strip A ('adc2')
        // to sum of A and neighbouring maxima C ('adc1')

        for (size_t k = 0; k < 4; ++k) {
          if (adc1[k] > 0 && adc2[k] > 0)
            adc[k] = adc[k] * adc2[k] / (adc1[k] + adc2[k]);
          if (adcRaw1[k] > 0 && adcRaw2[k] > 0)
            adcRaw[k] = adcRaw[k] * adcRaw2[k] / (adcRaw1[k] + adcRaw2[k]);
        }
        prelimData = CSCStripHitData(thisStrip, tmax_cluster, adcRaw, adc);
      }
    }
  }
  return prelimData;
}

runStrip()

std::vector< CSCStripHit > CSCHitFromStripOnly::runStrip	(	const CSCDetId &	id,
		const CSCLayer *	layer,
		const CSCStripDigiCollection::Range &	rstripd
	)

fact (20.10.09);

L1A (Begin looping) Attempt to redefine theStrips, to encode L1A phase bits


L1A (end Looping)

L1A

Print statement to check StripHit content w/ LA1

Definition at line 54 of file CSCHitFromStripOnly.cc.

References funct::abs(), CSCLayer::chamber(), clusterSize, fillPulseHeights(), findMaxima(), gainWeight, ganged(), CSCChamberSpecs::gangedStrips(), triggerObjects_cff::id, id_, createfilelist::int, isDeadStrip(), phase1PixelTopology::layer, layer_, LogTrace, makeCluster(), nstrips_, recoConditions_, CSCDetId::ring(), setGanged(), CSCChamber::specs(), strips_adc, strips_adcRaw, CSCRecoConditions::stripWeights(), theClosestMaximum, theClusterSize, theConsecutiveStrips, theMaxima, theStrips, tmax_cluster, and useCalib.

Referenced by CSCRecHitDBuilder::build().

                                                                                                   {
  std::vector<CSCStripHit> hitsInLayer;

  // cache layer info for ease of access
  id_ = id;
  layer_ = layer;
  nstrips_ = layer->chamber()->specs()->nStrips();

  setGanged(false);
  if (id_.ring() == 4 && layer_->chamber()->specs()->gangedStrips())
    setGanged(true);  //@@ ONLY ME1/1A CAN BE GANGED

  LogTrace("CSCHitFromStripOnly") << "[CSCHitFromStripOnly::runStrip] id= " << id_ << " nstrips= " << nstrips_
                                  << " ganged strips? " << ganged();

  tmax_cluster = 5;

  // Get gain correction weights for all strips in layer, and cache in gainWeight.
  // They're used in fillPulseHeights below.
  // When ME11a is ganged we only need the first 16 values of the 48 filled,
  // but 17-48 are just duplicates of 1-16 anyway

  if (useCalib) {
    recoConditions_->stripWeights(id, nstrips_, gainWeight);

    // *** START DUMP gainWeight
    //    std::cout << "gainWeight for id= " << id_ << " nstrips= " << nstrips_ << std::endl;
    //    for ( size_t i = 0; i!=10; ++i ) {
    //      for ( size_t j = 0; j!=8; ++j ) {
    //        std::cout << gainWeight[i*8 + j] << "   ";
    //      }
    //      std::cout << std::endl;
    //    }
    // *** END DUMP gainWeight
  }

  // Store pulseheights from SCA and find maxima (potential hits)
  fillPulseHeights(rstripd);
  findMaxima(id);

  // Make a Strip Hit out of each strip local maximum
  for (size_t imax = 0; imax != theMaxima.size(); ++imax) {
    // Initialize parameters entering the CSCStripHit
    clusterSize = theClusterSize;
    theStrips.clear();
    strips_adc.clear();
    strips_adcRaw.clear();

    // makeCluster calls findHitOnStripPosition to determine the centroid position

    // Remember, the array starts at 0, but the stripId starts at 1...
    float strippos = makeCluster(theMaxima[imax] + 1);

    //if ( strippos < 0 || tmax_cluster < 3 ){
    // the strippos (as calculated here) is not used later on in
    // with the negative charges allowed it can become negative
    if (tmax_cluster < 3) {
      theClosestMaximum.push_back(99);  // to keep proper vector size
      continue;
    }
    //---- If two maxima are too close the error assigned will be width/sqrt(12) - see CSCXonStrip_MatchGatti.cc
    int maximum_to_left = 99;  //---- If there is one maximum - the distance is set to 99 (strips)
    int maximum_to_right = 99;
    if (imax < theMaxima.size() - 1) {
      maximum_to_right = theMaxima.at(imax + 1) - theMaxima.at(imax);
    }
    if (imax > 0) {
      maximum_to_left = theMaxima.at(imax - 1) - theMaxima.at(imax);
    }
    if (std::abs(maximum_to_right) < std::abs(maximum_to_left)) {
      theClosestMaximum.push_back(maximum_to_right);
    } else {
      theClosestMaximum.push_back(maximum_to_left);
    }

    //---- Check if a neighbouring strip is a dead strip
    //bool deadStrip = isNearDeadStrip(id, theMaxima.at(imax));
    bool deadStripL = isDeadStrip(id, theMaxima.at(imax) - 1, nstrips_);
    bool deadStripR = isDeadStrip(id, theMaxima.at(imax) + 1, nstrips_);
    short int aDeadStrip = 0;
    if (!deadStripL && !deadStripR) {
      aDeadStrip = 0;
    } else if (deadStripL && deadStripR) {
      aDeadStrip = 255;
    } else {
      if (deadStripL) {
        aDeadStrip = theMaxima.at(imax) - 1;
      } else {
        aDeadStrip = theMaxima.at(imax) + 1;
      }
    }

    std::vector<int> theL1AStrips;
    for (int ila = 0; ila < (int)theStrips.size(); ila++) {
      bool stripMatchCounter = false;
      for (auto itl1 = rstripd.first; itl1 != rstripd.second; ++itl1) {
        int stripNproto = (*itl1).getStrip();
        if (!ganged()) {
          if (theStrips[ila] == stripNproto) {
            stripMatchCounter = true;
            auto sz = (*itl1).getOverlappedSample().size();
            int L1AbitOnPlace = 0;
            for (auto iBit = 0UL; iBit < sz; iBit++) {
              L1AbitOnPlace |= ((*itl1).getL1APhase(iBit) << (15 - iBit));
            }
            theL1AStrips.push_back(theStrips[ila] | L1AbitOnPlace);
          }
        } else {
          for (int tripl = 0; tripl < 3; ++tripl) {
            if (theStrips[ila] == (stripNproto + tripl * 16)) {
              stripMatchCounter = true;
              auto sz = (*itl1).getOverlappedSample().size();
              int L1AbitOnPlace = 0;
              for (auto iBit = 0UL; iBit < sz; iBit++) {
                L1AbitOnPlace |= ((*itl1).getL1APhase(iBit) << (15 - iBit));
              }
              theL1AStrips.push_back(theStrips[ila] | L1AbitOnPlace);
            }
          }
        }
      }
      if (!stripMatchCounter) {
        theL1AStrips.push_back(theStrips[ila]);
      }
    }

    CSCStripHit striphit(id,
                         strippos,
                         tmax_cluster,
                         theL1AStrips,
                         strips_adc,
                         strips_adcRaw,  
                         theConsecutiveStrips.at(imax),
                         theClosestMaximum.at(imax),
                         aDeadStrip);
    hitsInLayer.push_back(striphit);
  }

  /*   
      for(std::vector<CSCStripHit>::const_iterator itSHit=hitsInLayer.begin(); itSHit!=hitsInLayer.end(); ++itSHit){
         (*itSHit).print(); 
         }  
  */

  return hitsInLayer;
}

setConditions()

void CSCHitFromStripOnly::setConditions ( const CSCRecoConditions *  reco )

inline

Definition at line 46 of file CSCHitFromStripOnly.h.

References dt_dqm_sourceclient_common_cff::reco, and recoConditions_.

Referenced by CSCRecHitDBuilder::setConditions().

{ recoConditions_ = reco; }

setGanged()

void CSCHitFromStripOnly::setGanged ( bool  ig )

inline

Definition at line 49 of file CSCHitFromStripOnly.h.

References ganged_.

Referenced by runStrip().

{ ganged_ = ig; }

Member Data Documentation

calcped_

CSCPedestalChoice* CSCHitFromStripOnly::calcped_

private

Definition at line 87 of file CSCHitFromStripOnly.h.

Referenced by CSCHitFromStripOnly(), fillPulseHeights(), and ~CSCHitFromStripOnly().

clusterSize

int CSCHitFromStripOnly::clusterSize

private

Definition at line 104 of file CSCHitFromStripOnly.h.

Referenced by makeCluster(), makeStripData(), and runStrip().

gainWeight

float CSCHitFromStripOnly::gainWeight[80]

private

Definition at line 84 of file CSCHitFromStripOnly.h.

Referenced by fillPulseHeights(), and runStrip().

ganged_

bool CSCHitFromStripOnly::ganged_

private

Definition at line 109 of file CSCHitFromStripOnly.h.

Referenced by ganged(), and setGanged().

id_

CSCDetId CSCHitFromStripOnly::id_

private

Definition at line 78 of file CSCHitFromStripOnly.h.

Referenced by fillPulseHeights(), and runStrip().

layer_

const CSCLayer* CSCHitFromStripOnly::layer_

private

Definition at line 79 of file CSCHitFromStripOnly.h.

Referenced by runStrip().

nstrips_

unsigned CSCHitFromStripOnly::nstrips_

private

Definition at line 82 of file CSCHitFromStripOnly.h.

Referenced by findMaxima(), makeCluster(), and runStrip().

recoConditions_

const CSCRecoConditions* CSCHitFromStripOnly::recoConditions_

private

Definition at line 80 of file CSCHitFromStripOnly.h.

Referenced by fillPulseHeights(), isDeadStrip(), isNearDeadStrip(), runStrip(), and setConditions().

strips_adc

std::vector<float> CSCHitFromStripOnly::strips_adc

private

Definition at line 105 of file CSCHitFromStripOnly.h.

Referenced by findHitOnStripPosition(), and runStrip().

strips_adcRaw

std::vector<float> CSCHitFromStripOnly::strips_adcRaw

private

Definition at line 106 of file CSCHitFromStripOnly.h.

Referenced by findHitOnStripPosition(), and runStrip().

theClosestMaximum

std::vector<int> CSCHitFromStripOnly::theClosestMaximum

private

Definition at line 100 of file CSCHitFromStripOnly.h.

Referenced by findMaxima(), and runStrip().

theClusterSize

const int CSCHitFromStripOnly::theClusterSize = 3

staticprivate

Definition at line 91 of file CSCHitFromStripOnly.h.

Referenced by makeCluster(), and runStrip().

theConsecutiveStrips

std::vector<int> CSCHitFromStripOnly::theConsecutiveStrips

private

Definition at line 99 of file CSCHitFromStripOnly.h.

Referenced by findMaxima(), and runStrip().

theMaxima

std::vector<int> CSCHitFromStripOnly::theMaxima

private

Definition at line 98 of file CSCHitFromStripOnly.h.

Referenced by findMaxima(), makeStripData(), and runStrip().

thePulseHeightMap

PulseHeightMap CSCHitFromStripOnly::thePulseHeightMap

private

Definition at line 96 of file CSCHitFromStripOnly.h.

Referenced by fillPulseHeights(), findMaxima(), isPeakOK(), and makeStripData().

theStrips

std::vector<int> CSCHitFromStripOnly::theStrips

private

Definition at line 107 of file CSCHitFromStripOnly.h.

Referenced by makeCluster(), and runStrip().

theThresholdForAPeak

float CSCHitFromStripOnly::theThresholdForAPeak

private

Definition at line 92 of file CSCHitFromStripOnly.h.

Referenced by CSCHitFromStripOnly(), and isPeakOK().

theThresholdForCluster

float CSCHitFromStripOnly::theThresholdForCluster

private

Definition at line 93 of file CSCHitFromStripOnly.h.

Referenced by CSCHitFromStripOnly(), and isPeakOK().

tmax_cluster

int CSCHitFromStripOnly::tmax_cluster

private

Definition at line 103 of file CSCHitFromStripOnly.h.

Referenced by makeStripData(), and runStrip().

useCalib

bool CSCHitFromStripOnly::useCalib

private

Definition at line 90 of file CSCHitFromStripOnly.h.

Referenced by CSCHitFromStripOnly(), fillPulseHeights(), and runStrip().