CSCHitFromStripOnly.h

Go to the documentation of this file.

#ifndef CSCRecHitD_CSCHitFromStripOnly_h
#define CSCRecHitD_CSCHitFromStripOnly_h

#include "DataFormats/MuonDetId/interface/CSCDetId.h"

#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>
#include <array>

class CSCLayer;
class CSCStripDigi;
class CSCPedestalChoice;

class CSCHitFromStripOnly 
{
  
 public:

  typedef std::array<CSCStripData,100> 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;
  } 
 
  bool ganged() { return ganged_;}
  void setGanged( bool ig ) { ganged_ = ig;}

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

  void findMaxima(const CSCDetId& id);  
  // What we call a peak
  bool isPeakOK(int iStrip, float heightCluster);  

  float makeCluster( int centerStrip );

  CSCStripHitData makeStripData( int centerStrip, int offset );

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

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

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

// MEMBER DATA

  // Hold pointers to current layer, conditions data
  CSCDetId id_;    
  const CSCLayer * layer_;
  const CSCRecoConditions* recoConditions_;
  // Number of strips in layer
  unsigned nstrips_;
  // gain correction weights and crosstalks read in from conditions database.
  float gainWeight[80];

  // The specific pedestal calculator
  CSCPedestalChoice* calcped_;

  // The cuts for forming the strip hits are described in the config file
  bool useCalib;
  static const int theClusterSize = 3;
  float theThresholdForAPeak;
  float theThresholdForCluster;


  // working buffer for sca pulseheights
  PulseHeightMap thePulseHeightMap;

  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

  // Variables entering the CSCStripHit construction:
  int tmax_cluster; // Peaking time for strip hit, in time bin units
  int clusterSize;
  std::vector<float> strips_adc;
  std::vector<float> strips_adcRaw;
  std::vector<int> theStrips;

  bool ganged_; // only True if ME1/1A AND it is ganged
  
};

#endif