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

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

Definition at line 39 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_(0), calcped_(0), ganged_(0) {
  
  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 50 of file CSCHitFromStripOnly.cc.

References calcped_.

                                          {
  delete calcped_;
}

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 360 of file CSCHitFromStripOnly.cc.

References calcped_, popcon2dropbox::copy(), gainWeight, ganged(), id_, CSCPedestalChoice::pedestal(), recoConditions_, thePulseHeightMap, tmax, and useCalib.

Referenced by runStrip(), and setGanged().

                                                                                       {
  
  // 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 510 of file CSCHitFromStripOnly.cc.

References popcon2dropbox::copy(), mps_fire::i, LogTrace, digi_MixPreMix_cfi::strip, strips_adc, strips_adcRaw, and w.

Referenced by makeCluster(), and setGanged().

                                                                                                                  {
  
  float strippos = -1.;
  
  if ( data.size() < 1 ) 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 424 of file CSCHitFromStripOnly.cc.

References mps_fire::i, isDeadStrip(), isPeakOK(), gen::k, checklumidiff::l, nstrips_, theClosestMaximum, theConsecutiveStrips, theMaxima, and thePulseHeightMap.

Referenced by runStrip(), and setGanged().

                                                       {
  
  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)){
      theMaxima.push_back(i);
      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)){
      theMaxima.push_back(i);
      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)){
      theMaxima.push_back(i);
      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;
        }
      }
      theConsecutiveStrips.push_back(numberOfConsecutiveStrips);
    }
  }
}

bool CSCHitFromStripOnly::ganged ( )

inline

Definition at line 51 of file CSCHitFromStripOnly.h.

References ganged_.

Referenced by fillPulseHeights(), and runStrip().

{ return ganged_;}

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

private

Is the strip 'bad'?

Definition at line 568 of file CSCHitFromStripOnly.cc.

References CSCRecoConditions::badStrip(), recoConditions_, and theClusterSize.

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

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

} 

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

private

Is either neighbour 'bad'?

Definition at line 560 of file CSCHitFromStripOnly.cc.

References CSCRecoConditions::nearBadStrip(), and recoConditions_.

Referenced by setGanged().

                                                                                           {

  //@@ 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 496 of file CSCHitFromStripOnly.cc.

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

Referenced by findMaxima(), and setGanged().

                                                                 {
  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 218 of file CSCHitFromStripOnly.cc.

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

Referenced by runStrip(), and setGanged().

                                                        {
  
  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 252 of file CSCHitFromStripOnly.cc.

References ecalMGPA::adc(), begin, clusterSize, popcon2dropbox::copy(), spr::find(), mps_fire::i, gen::k, LogTrace, min(), PFRecoTauDiscriminationByIsolation_cfi::offset, Validation_hcalonly_cfi::sign, theMaxima, thePulseHeightMap, tmax, and tmax_cluster.

Referenced by makeCluster(), and setGanged().

                                                                              {
  
  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 61 of file CSCHitFromStripOnly.cc.

References funct::abs(), CSCLayer::chamber(), clusterSize, fillPulseHeights(), findMaxima(), gainWeight, ganged(), CSCChamberSpecs::gangedStrips(), hcalTTPDigis_cfi::id, id_, createfilelist::int, isDeadStrip(), 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;
}

void CSCHitFromStripOnly::setConditions ( const CSCRecoConditions *  reco )

inline

Definition at line 47 of file CSCHitFromStripOnly.h.

References trackingPlots::reco, and recoConditions_.

Referenced by CSCRecHitDBuilder::setConditions().

                                                      {
    recoConditions_ = reco;
  } 

void CSCHitFromStripOnly::setGanged ( bool  ig )

inline

Definition at line 52 of file CSCHitFromStripOnly.h.

References data, fillPulseHeights(), findHitOnStripPosition(), findMaxima(), ganged_, isDeadStrip(), isNearDeadStrip(), isPeakOK(), makeCluster(), makeStripData(), and PFRecoTauDiscriminationByIsolation_cfi::offset.

Referenced by runStrip().

{ ganged_ = ig;}

Member Data Documentation

CSCPedestalChoice* CSCHitFromStripOnly::calcped_

private

Definition at line 92 of file CSCHitFromStripOnly.h.

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

int CSCHitFromStripOnly::clusterSize

private

Definition at line 110 of file CSCHitFromStripOnly.h.

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

float CSCHitFromStripOnly::gainWeight[80]

private

Definition at line 89 of file CSCHitFromStripOnly.h.

Referenced by fillPulseHeights(), and runStrip().

bool CSCHitFromStripOnly::ganged_

private

Definition at line 115 of file CSCHitFromStripOnly.h.

Referenced by ganged(), and setGanged().

CSCDetId CSCHitFromStripOnly::id_

private

Definition at line 83 of file CSCHitFromStripOnly.h.

Referenced by fillPulseHeights(), and runStrip().

const CSCLayer* CSCHitFromStripOnly::layer_

private

Definition at line 84 of file CSCHitFromStripOnly.h.

Referenced by runStrip().

unsigned CSCHitFromStripOnly::nstrips_

private

Definition at line 87 of file CSCHitFromStripOnly.h.

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

const CSCRecoConditions* CSCHitFromStripOnly::recoConditions_

private

Definition at line 85 of file CSCHitFromStripOnly.h.

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

std::vector<float> CSCHitFromStripOnly::strips_adc

private

Definition at line 111 of file CSCHitFromStripOnly.h.

Referenced by findHitOnStripPosition(), and runStrip().

std::vector<float> CSCHitFromStripOnly::strips_adcRaw

private

Definition at line 112 of file CSCHitFromStripOnly.h.

Referenced by findHitOnStripPosition(), and runStrip().

std::vector<int> CSCHitFromStripOnly::theClosestMaximum

private

Definition at line 106 of file CSCHitFromStripOnly.h.

Referenced by findMaxima(), and runStrip().

const int CSCHitFromStripOnly::theClusterSize = 3

staticprivate

Definition at line 96 of file CSCHitFromStripOnly.h.

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

std::vector<int> CSCHitFromStripOnly::theConsecutiveStrips

private

Definition at line 105 of file CSCHitFromStripOnly.h.

Referenced by findMaxima(), and runStrip().

std::vector<int> CSCHitFromStripOnly::theMaxima

private

Definition at line 104 of file CSCHitFromStripOnly.h.

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

PulseHeightMap CSCHitFromStripOnly::thePulseHeightMap

private

Definition at line 102 of file CSCHitFromStripOnly.h.

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

std::vector<int> CSCHitFromStripOnly::theStrips

private

Definition at line 113 of file CSCHitFromStripOnly.h.

Referenced by makeCluster(), and runStrip().

float CSCHitFromStripOnly::theThresholdForAPeak

private

Definition at line 97 of file CSCHitFromStripOnly.h.

Referenced by CSCHitFromStripOnly(), and isPeakOK().

float CSCHitFromStripOnly::theThresholdForCluster

private

Definition at line 98 of file CSCHitFromStripOnly.h.

Referenced by CSCHitFromStripOnly(), and isPeakOK().

int CSCHitFromStripOnly::tmax_cluster

private

Definition at line 109 of file CSCHitFromStripOnly.h.

Referenced by makeStripData(), and runStrip().

bool CSCHitFromStripOnly::useCalib

private

Definition at line 95 of file CSCHitFromStripOnly.h.

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