---

SiStripClusterInfo.h

Go to the documentation of this file.

#ifndef ANALYSISDATAFORMATS_SISTRIPCLUSTERINFO_H
#define ANALYSISDATAFORMATS_SISTRIPCLUSTERINFO_H
// -*- C++ -*-
//
// Package:     AnalysisDataFormats
// Class  :     SiStripClusterInfo
// 
//
// Original Author:  Evelyne Delmeire
//         Created:  Sun Jan 27 16:07:51 CET 2007
//
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "DataFormats/Common/interface/DetSet.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/DetSetNew.h"
#include "DataFormats/Common/interface/DetSetVectorNew.h"
#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"

#include "CondFormats/SiStripObjects/interface/SiStripPedestals.h"
#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"

#include "CalibFormats/SiStripObjects/interface/SiStripGain.h"

#include "CommonTools/SiStripZeroSuppression/interface/SiStripPedestalsSubtractor.h"
#include "CommonTools/SiStripZeroSuppression/interface/SiStripCommonModeNoiseSubtractor.h"

#include <vector>
#include <sstream>
#include <string>
#include "boost/cstdint.hpp"

class SiStripRawDigi;

class SiStripDigi;

class SiStripClusterInfo {
  
 public:
  
  // constructor with the detid (in case it will be eliminated from the SiStripCluster)
  SiStripClusterInfo(const uint32_t cluster_detId,
                     const SiStripCluster& cluster, 
                     const edm::EventSetup& es, 
                     std::string CMNSubtractionMode="Median");
  
  // constructor without the detid (taken from SiStripCluster)
  SiStripClusterInfo(const SiStripCluster& cluster, 
                     const edm::EventSetup& es, 
                     std::string CMNSubtractionMode="Median");
  
  ~SiStripClusterInfo();
    
  uint32_t getDetId() const {return cluster_detId_;};
  
  uint16_t                           getFirstStrip() const {return cluster_->firstStrip();}
  float                                getPosition() const {return cluster_->barycenter();}
  const std::vector<uint8_t>&  getStripAmplitudes() const {return cluster_->amplitudes();} 
  
  
  float                                   getWidth() const {return cluster_->amplitudes().size();}
  float                                  getCharge() const;
  uint16_t                          getMaxPosition() const;
  uint16_t                          getMaxPosition(const SiStripCluster*) const;
  float                               getMaxCharge() const;  
  std::pair<float,float>               getChargeLR() const;
  std::pair<float,float> getChargeLRFirstNeighbour() const;
  //
  std::vector<float> getRawChargeCLR( const edm::DetSet<SiStripRawDigi>&,
                                      edmNew::DetSet<SiStripCluster>&,
                                      std::string );
  
  
  std::vector<float>               getStripNoises() const; 
  std::vector<float> getStripNoisesRescaledByGain() const; 
  float                                  getNoise()  ;
  float                    getNoiseRescaledByGain()  ;
  float                        getNoiseForStripNb(uint16_t istrip ) const;
  
  std::vector<float> getApvGains() const; 
  float        getGainForStripNb(uint16_t istrip ) const;
  
  float               getSignalOverNoise() ;
  float getSignalOverNoiseRescaledByGain() ;

  
  std::pair< std::vector<float>,std::vector<float> > getRawDigiAmplitudesLR(      uint32_t&                          neighbourStripNr, 
                                                                            const edm::DetSetVector<SiStripRawDigi>&    rawDigis_dsv_, 
                                                                                  edmNew::DetSetVector<SiStripCluster>& clusters_dsv_,
                                                                                  std::string                           rawDigiLabel);

  std::pair< std::vector<float>,std::vector<float> > getRawDigiAmplitudesLR(      uint32_t&                   neighbourStripNr, 
                                                                            const edm::DetSet<SiStripRawDigi>&    rawDigis_ds_, 
                                                                                  edmNew::DetSet<SiStripCluster>& clusters_ds_,
                                                                                  std::string                     rawDigiLabel); 


  
  std::pair< std::vector<float>,std::vector<float> >  getDigiAmplitudesLR(      uint32_t&                          neighbourStripNr,
                                                                          const edm::DetSetVector<SiStripDigi>&          digis_dsv_,
                                                                                edmNew::DetSetVector<SiStripCluster>& clusters_dsv_);

  std::pair< std::vector<float>,std::vector<float> > getDigiAmplitudesLR(      uint32_t&                   neighbourStripNr, 
                                                                         const edm::DetSet<SiStripDigi>&          digis_ds_, 
                                                                               edmNew::DetSet<SiStripCluster>& clusters_ds_);
 

 private:
  
  
  void rawdigi_algorithm(const edm::DetSet<SiStripRawDigi>    rawDigis_ds_,
                               edmNew::DetSet<SiStripCluster> clusters_ds_,
                               std::string                    rawDigiLabel);
  
  void digi_algorithm(const edm::DetSet<SiStripDigi>         digis_ds_,
                            edmNew::DetSet<SiStripCluster> cluster_ds_);
  
  void findNeigh(char*                                  mode,
                 edmNew::DetSet<SiStripCluster> clusters_ds_,
                 std::vector<int16_t>&                  vadc,
                 std::vector<int16_t>&                vstrip);
 
  const edm::EventSetup&      es_; 
  edm::ESHandle<SiStripNoises>        noiseHandle_;
  edm::ESHandle<SiStripGain>           gainHandle_;
  edm::ESHandle<SiStripPedestals> pedestalsHandle_;
    
  const SiStripCluster *cluster_;
  uint32_t cluster_detId_;
   
  uint32_t neighbourStripNr_;
  SiStripCommonModeNoiseSubtractor* SiStripCommonModeNoiseSubtractor_;
  std::string CMNSubtractionMode_;
  bool validCMNSubtraction_;  
  SiStripPedestalsSubtractor* SiStripPedestalsSubtractor_;
  float                amplitudeC_;
  std::vector<float>   amplitudesL_;
  std::vector<float>   amplitudesR_;
  
};

#endif // ANALYSISDATAFORMATS_SISTRIPCLUSTER_H

---