---

CSCHitFromStripOnly.h

Go to the documentation of this file.

#ifndef CSCRecHitD_CSCHitFromStripOnly_h
#define CSCRecHitD_CSCHitFromStripOnly_h

//---- Possible changes from Stoyan Stoynev - NU
#include <RecoLocalMuon/CSCRecHitD/src/CSCStripData.h>
#include <RecoLocalMuon/CSCRecHitD/src/CSCStripHitData.h>
#include <RecoLocalMuon/CSCRecHitD/src/CSCStripHit.h>
#include <RecoLocalMuon/CSCRecHitD/src/CSCRecoConditions.h>

#include <DataFormats/CSCDigi/interface/CSCStripDigiCollection.h>

#include <FWCore/Framework/interface/Frameworkfwd.h>
#include <FWCore/ParameterSet/interface/ParameterSet.h>

#include <vector>

class CSCLayer;
class CSCChamberSpecs;
class CSCLayerGeometry;
class CSCStripDigi;
class CSCPeakBinOfStripPulse;


class CSCHitFromStripOnly 
{
  
 public:

  typedef std::vector<CSCStripData> PulseHeightMap;
  
  explicit CSCHitFromStripOnly( const edm::ParameterSet& ps );
  
  ~CSCHitFromStripOnly();
  
  std::vector<CSCStripHit> runStrip( const CSCDetId& id, const CSCLayer* layer, const CSCStripDigiCollection::Range& rstripd );

  void setConditions( const CSCRecoConditions* reco ) {
    recoConditions_ = reco;
  }

 protected:
  
  void fillPulseHeights( const CSCStripDigiCollection::Range& rstripd );  

  void findMaxima();    

  float makeCluster( int centerStrip );

  std::vector<int> theMaxima;
  std::vector<int>  theConsecutiveStrips;//... with charge for a given maximum
  std::vector<int> theClosestMaximum; // this is number of strips to the closest other maximum

  PulseHeightMap thePulseHeightMap;
  
  float findHitOnStripPosition( const std::vector<CSCStripHitData>& data, const int& centerStrip );
  
  CSCDetId id_;    
  const CSCLayer * layer_;
  const CSCLayerGeometry * layergeom_;
  const CSCChamberSpecs * specs_;
  
 private:
  
  CSCStripHitData makeStripData( int centerStrip, int offset );
  bool isNearDeadStrip(const CSCDetId& id, int centralStrip); 

  // Variables entering the CSCStripHit construction:
  int tmax_cluster;
  int clusterSize;
  std::vector<float> strips_adc;
  std::vector<float> strips_adcRaw;
  std::vector<int> theStrips;
  
  // The cuts for forming the strip hits are described in the data/.cfi file
  bool useCalib;
  static const int theClusterSize = 3;
  float theThresholdForAPeak;
  float theThresholdForCluster;
  //bool isData;

  float gainWeight[80];

  // Peaking time for strip hit
  int TmaxOfCluster;            // in time bins;
  // Number of strips in layer
  unsigned Nstrips;

  const CSCRecoConditions* recoConditions_;

  CSCPeakBinOfStripPulse* pulseheightOnStripFinder_;
};

#endif

---