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