#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

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

Definition at line 36 of file CSCHitFromStripOnly.h.

Constructor & Destructor Documentation

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 ( )

Definition at line 46 of file CSCHitFromStripOnly.cc.

References calcped_.

46 { delete calcped_; }

CSCHitFromStripOnly::calcped_

CSCPedestalChoice * calcped_

Definition: CSCHitFromStripOnly.h:87

Member Function Documentation

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;
       }
     }
   }
 }

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, 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;
 }

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, isotrackApplyRegressor::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

bool CSCHitFromStripOnly::ganged ( )

inline

Definition at line 48 of file CSCHitFromStripOnly.h.

References ganged_.

Referenced by fillPulseHeights(), and runStrip().

48 { return ganged_; }

CSCHitFromStripOnly::ganged_

bool ganged_

Definition: CSCHitFromStripOnly.h:109

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);
 }

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);
 }

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;
 }

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;
 }

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

private

makeStripData

Definition at line 241 of file CSCHitFromStripOnly.cc.

References gpuClustering::adc, SplitLinear::begin, clusterSize, filterCSVwithJSON::copy, spr::find(), mps_fire::i, isotrackApplyRegressor::k, LogTrace, min(), hltrates_dqm_sourceclient-live_cfg::offset, jetcorrextractor::sign(), theMaxima, thePulseHeightMap, tmax, and tmax_cluster.

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;
 }

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(), gpuClustering::id, id_, isDeadStrip(), phase1PixelTopology::layer, layer_, LogTrace, makeCluster(), CSCChamberSpecs::nStrips(), 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;
 }