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