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/TrackerGeometryBuilder/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::ConcurrentBooker &, 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;
}