SiStripGainsPCLWorker.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    CalibTracker/SiStripChannelGain
// Class:      SiStripGainsPCLWorker
//
//
//  Original Author: L. Quertermont (calibration algorithm)
//  Contributors:    M. Verzetti    (data access)
//                   A. Di Mattia   (PCL multi stream processing and monitoring)
//                   M. Delcourt    (monitoring)
//                   M. Musich      (migration to thread-safe DQMStore access)
//
//  Created:  Wed, 12 Apr 2017 14:46:48 GMT
//
 
// CMSSW includes
#include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
#include "CalibFormats/SiStripObjects/interface/SiStripGain.h"
#include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
#include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
#include "CalibTracker/Records/interface/SiStripGainRcd.h"
#include "CalibTracker/Records/interface/SiStripQualityRcd.h"
#include "CondFormats/SiStripObjects/interface/SiStripApvGain.h"
#include "DQMServices/Core/interface/DQMGlobalEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "DataFormats/GeometrySurface/interface/RectangularPlaneBounds.h"
#include "DataFormats/GeometrySurface/interface/TrapezoidalPlaneBounds.h"
#include "DataFormats/SiStripCluster/interface/SiStripClusterCollection.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2D.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit1D.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/Framework/interface/ESWatcher.h"
#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "Geometry/CommonDetUnit/interface/TrackingGeometry.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/PixelGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
#include "TrackingTools/PatternTools/interface/TrajTrackAssociation.h"
#include "TrackingTools/PatternTools/interface/Trajectory.h"
 
#include "CalibTracker/SiStripChannelGain/interface/APVGainStruct.h"
#include "CalibTracker/SiStripChannelGain/interface/APVGainHelpers.h"
 
// System includes
#include <unordered_map>
 
//
// class declaration
//
 
class SiStripGainsPCLWorker : public DQMGlobalEDAnalyzer<APVGain::APVGainHistograms> {
public:
  explicit SiStripGainsPCLWorker(const edm::ParameterSet &);
 
  void bookHistograms(DQMStore::IBooker &,
                      edm::Run const &,
                      edm::EventSetup const &,
                      APVGain::APVGainHistograms &) const override;
  void dqmAnalyze(edm::Event const &, edm::EventSetup const &, APVGain::APVGainHistograms const &) const override;
 
  static void fillDescriptions(edm::ConfigurationDescriptions &descriptions);
 
private:
  void beginJob() override;
  void dqmBeginRun(edm::Run const &, edm::EventSetup const &, APVGain::APVGainHistograms &) const override;
  void endJob() override;
  void checkBookAPVColls(const TrackerGeometry *bareTkGeomPtr, APVGain::APVGainHistograms &histograms) const;
 
  std::vector<std::string> dqm_tag_;
 
  int statCollectionFromMode(const char *tag) const;
 
  double MinTrackMomentum;
  double MaxTrackMomentum;
  double MinTrackEta;
  double MaxTrackEta;
  unsigned int MaxNrStrips;
  unsigned int MinTrackHits;
  double MaxTrackChiOverNdf;
  int MaxTrackingIteration;
  bool AllowSaturation;
  bool FirstSetOfConstants;
  bool Validation;
  bool OldGainRemoving;
  bool useCalibration;
  bool doChargeMonitorPerPlane; 
  std::string m_DQMdir;                  
  std::string m_calibrationMode;         
  std::vector<std::string> VChargeHisto; 
  //Data members for processing
 
  edm::EDGetTokenT<std::vector<bool> > TrigTech_token_;
  edm::EDGetTokenT<std::vector<double> > trackchi2ndof_token_;
  edm::EDGetTokenT<std::vector<float> > trackp_token_;
  edm::EDGetTokenT<std::vector<float> > trackpt_token_;
  edm::EDGetTokenT<std::vector<double> > tracketa_token_;
  edm::EDGetTokenT<std::vector<double> > trackphi_token_;
  edm::EDGetTokenT<std::vector<unsigned int> > trackhitsvalid_token_;
  edm::EDGetTokenT<std::vector<int> > trackalgo_token_;
  edm::EDGetTokenT<std::vector<int> > trackindex_token_;
  edm::EDGetTokenT<std::vector<unsigned int> > rawid_token_;
  edm::EDGetTokenT<std::vector<double> > localdirx_token_;
  edm::EDGetTokenT<std::vector<double> > localdiry_token_;
  edm::EDGetTokenT<std::vector<double> > localdirz_token_;
  edm::EDGetTokenT<std::vector<unsigned short> > firststrip_token_;
  edm::EDGetTokenT<std::vector<unsigned short> > nstrips_token_;
  edm::EDGetTokenT<std::vector<bool> > saturation_token_;
  edm::EDGetTokenT<std::vector<bool> > overlapping_token_;
  edm::EDGetTokenT<std::vector<bool> > farfromedge_token_;
  edm::EDGetTokenT<std::vector<unsigned int> > charge_token_;
  edm::EDGetTokenT<std::vector<double> > path_token_;
  edm::EDGetTokenT<std::vector<double> > chargeoverpath_token_;
  edm::EDGetTokenT<std::vector<unsigned char> > amplitude_token_;
  edm::EDGetTokenT<std::vector<double> > gainused_token_;
  edm::EDGetTokenT<std::vector<double> > gainusedTick_token_;
 
  std::string EventPrefix_;  //("");
  std::string EventSuffix_;  //("");
  std::string TrackPrefix_;  //("track");
  std::string TrackSuffix_;  //("");
  std::string CalibPrefix_;  //("GainCalibration");
  std::string CalibSuffix_;  //("");
 
  // maps histograms index to topology
  std::map<unsigned int, APVloc> theTopologyMap;
};
 
inline int SiStripGainsPCLWorker::statCollectionFromMode(const char *tag) const {
  std::vector<std::string>::const_iterator it = dqm_tag_.begin();
  while (it != dqm_tag_.end()) {
    if (*it == std::string(tag))
      return it - dqm_tag_.begin();
    it++;
  }
 
  if (std::string(tag).empty())
    return 0;  // return StdBunch calibration mode for backward compatibility
 
  return None;
}