SiStripThreshold_PayloadInspector.cc

Go to the documentation of this file.

#include "CondCore/Utilities/interface/PayloadInspectorModule.h"
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "CondCore/CondDB/interface/Time.h"

// the data format of the condition to be inspected
#include "CondFormats/SiStripObjects/interface/SiStripThreshold.h"
#include "DataFormats/DetId/interface/DetId.h"                     
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"   
#include "CondFormats/SiStripObjects/interface/SiStripDetSummary.h"  
#include "FWCore/MessageLogger/interface/MessageLogger.h"      

// needed for the tracker map
#include "CommonTools/TrackerMap/interface/TrackerMap.h"   

// auxilliary functions
#include "CondCore/SiStripPlugins/interface/SiStripPayloadInspectorHelper.h"
#include "CalibTracker/StandaloneTrackerTopology/interface/StandaloneTrackerTopology.h"

#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"

#include <memory>
#include <sstream>
#include <iostream>

// include ROOT 
#include "TH2F.h"
#include "TLegend.h"
#include "TCanvas.h"
#include "TLine.h"
#include "TStyle.h"
#include "TLatex.h"
#include "TPave.h"
#include "TPaveStats.h"

using namespace std;

namespace {

  /************************************************
    test class
  *************************************************/

  class SiStripThresholdTest : public cond::payloadInspector::Histogram1D<SiStripThreshold> {

  public:
    SiStripThresholdTest() : cond::payloadInspector::Histogram1D<SiStripThreshold>("SiStrip Threshold test",
                                           "SiStrip Threshold test", 10,0.0,10.0),
                 m_trackerTopo{StandaloneTrackerTopology::fromTrackerParametersXMLFile(edm::FileInPath("Geometry/TrackerCommonData/data/trackerParameters.xml").fullPath())} 
    {
      Base::setSingleIov( true );
    }
    
    bool fill( const std::vector<std::tuple<cond::Time_t,cond::Hash> >& iovs ) override{ 
      for ( auto const & iov: iovs) {
    std::shared_ptr<SiStripThreshold> payload = Base::fetchPayload( std::get<1>(iov) ); 
    if( payload.get() ){
     
      fillWithValue(1.);
     
      std::stringstream ss;
      ss << "Summary of strips threshold:" << std::endl;

      payload->printSummary(ss,&m_trackerTopo);

      std::vector<uint32_t> detid;
      payload->getDetIds(detid);
          
      std::cout<<ss.str()<<std::endl;
 
    }
      }
      return true;
    }// fill
  private:
    TrackerTopology m_trackerTopo;
  };

  /************************************************************
    1d histogram of SiStripThresholds of 1 IOV - High Threshold 
  *************************************************************/

  class SiStripThresholdValueHigh : public cond::payloadInspector::Histogram1D<SiStripThreshold> {
    
  public:
    SiStripThresholdValueHigh() : cond::payloadInspector::Histogram1D<SiStripThreshold>("SiStrip High threshold values (checked per APV)",
                                            "SiStrip High threshold values (cheched per APV)", 10,0.0,10){
      Base::setSingleIov( true );
    }
    
    bool fill( const std::vector<std::tuple<cond::Time_t,cond::Hash> >& iovs ) override{

      edm::FileInPath fp_ = edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat");
      SiStripDetInfoFileReader* reader = new SiStripDetInfoFileReader(fp_.fullPath());

      for ( auto const & iov: iovs) {
    std::shared_ptr<SiStripThreshold> payload = Base::fetchPayload( std::get<1>(iov) );
    if( payload.get() ){
      std::vector<uint32_t> detid;
      payload->getDetIds(detid);

      for (const auto & d : detid) {
        //std::cout<<d<<std::endl;
        SiStripThreshold::Range range=payload->getRange(d);

        int nAPVs = reader->getNumberOfApvsAndStripLength(d).first;
        
        for( int it=0; it<nAPVs; ++it ){

          auto hth = payload->getData(it*128,range).getHth();
          //std::cout<<hth<<std::endl;
          fillWithValue(hth);
        }
      }
    }
      }

      delete reader;
      return true;
    }
  };

  /************************************************************
    1d histogram of SiStripThresholds of 1 IOV - Low Threshold 
  *************************************************************/

  class SiStripThresholdValueLow : public cond::payloadInspector::Histogram1D<SiStripThreshold> {
    
  public:
    SiStripThresholdValueLow() : cond::payloadInspector::Histogram1D<SiStripThreshold>("SiStrip Low threshold values (checked per APV)",
                                            "SiStrip Low threshold values (cheched per APV)", 10,0.0,10){
      Base::setSingleIov( true );
    }
    
    bool fill( const std::vector<std::tuple<cond::Time_t,cond::Hash> >& iovs ) override{

      edm::FileInPath fp_ = edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat");
      SiStripDetInfoFileReader* reader = new SiStripDetInfoFileReader(fp_.fullPath());

      for ( auto const & iov: iovs) {
    std::shared_ptr<SiStripThreshold> payload = Base::fetchPayload( std::get<1>(iov) );
    if( payload.get() ){
      std::vector<uint32_t> detid;
      payload->getDetIds(detid);

      for (const auto & d : detid) {
        //std::cout<<d<<std::endl;
        SiStripThreshold::Range range=payload->getRange(d);

        int nAPVs = reader->getNumberOfApvsAndStripLength(d).first;
            
        for( int it=0; it<nAPVs; ++it ){

          auto lth = payload->getData(it*128,range).getLth();
          //std::cout<<hth<<std::endl;
          fillWithValue(lth);
        }
      }
    }
      }

      delete reader;
      return true;
    }
  };





} 

PAYLOAD_INSPECTOR_MODULE(SiStripThreshold){
  PAYLOAD_INSPECTOR_CLASS(SiStripThresholdTest);
  PAYLOAD_INSPECTOR_CLASS(SiStripThresholdValueHigh);
  PAYLOAD_INSPECTOR_CLASS(SiStripThresholdValueLow);

}