SiStripGainsPCLHarvester.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    CalibTracker/SiStripChannelGain
// Class:      SiStripGainsPCLHarvester
// 
//
//  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 "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 "DQMServices/Core/interface/DQMEDHarvester.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

// user includes
#include "CalibTracker/SiStripChannelGain/interface/APVGainStruct.h"

// ROOT includes
#include "TH1F.h"
#include "TH2S.h"
#include "TProfile.h"
#include "TF1.h"

// System includes
#include <unordered_map>

class SiStripGainsPCLHarvester : public  DQMEDHarvester {
   public:
      explicit SiStripGainsPCLHarvester(const edm::ParameterSet& ps);
      virtual void beginRun(edm::Run const& run, edm::EventSetup const & isetup);
      virtual void endRun(edm::Run const& run, edm::EventSetup const & isetup);
      static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

   private:

      virtual void checkBookAPVColls(const edm::EventSetup& setup);
      virtual void dqmEndJob(DQMStore::IBooker& ibooker_, DQMStore::IGetter& igetter_);

      int statCollectionFromMode(const char* tag) const;

      void algoComputeMPVandGain(const MonitorElement* Charge_Vs_Index);
      void getPeakOfLandau(TH1* InputHisto, double* FitResults, double LowRange=50, double HighRange=5400);
      bool IsGoodLandauFit(double* FitResults); 

      bool produceTagFilter(const MonitorElement* Charge_Vs_Index);
      std::unique_ptr<SiStripApvGain> getNewObject(const MonitorElement* Charge_Vs_Index);

      bool doStoreOnDB;
      unsigned int GOOD;
      unsigned int BAD;
      unsigned int MASKED;

      double tagCondition_NClusters;
      double tagCondition_GoodFrac;

      int NStripAPVs;
      int NPixelDets;
      double MinNrEntries;

      std::string m_Record;
      std::string m_DQMdir;
      std::string m_calibrationMode;
      std::vector<std::string> dqm_tag_;  

      int CalibrationLevel;

      edm::ESHandle<TrackerGeometry> tkGeom_;
      const TrackerGeometry *bareTkGeomPtr_;  // ugly hack to fill APV colls only once, but checks

      std::vector<std::shared_ptr<stAPVGain> > APVsCollOrdered;
      std::unordered_map<unsigned int, std::shared_ptr<stAPVGain> > APVsColl; 

};