EcalTPGFineGrainTowerEE_PayloadInspector.cc

Go to the documentation of this file.

#include "CondCore/Utilities/interface/PayloadInspectorModule.h"
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
#include "CondCore/EcalPlugins/plugins/EcalDrawUtils.h"

// the data format of the condition to be inspected
#include "CondFormats/EcalObjects/interface/EcalTPGFineGrainTowerEE.h"

#include "TH2F.h"
#include "TCanvas.h"
#include "TStyle.h"
#include "TLine.h"
#include "TLatex.h"

#include <string>

namespace {

/***********************************************
   2d plot of EcalTPGFineGrainTowerEE of 1 IOV
************************************************/
class EcalTPGFineGrainTowerEEPlot : public cond::payloadInspector::PlotImage<EcalTPGFineGrainTowerEE> {

  public:

    EcalTPGFineGrainTowerEEPlot() : cond::payloadInspector::PlotImage<EcalTPGFineGrainTowerEE>("EcalTPGFineGrainTowerEE - map ") {
      setSingleIov(true);
    }

    bool fill( const std::vector<std::tuple<cond::Time_t,cond::Hash> >& iovs ) override{
      
     TH2F* endc_p = new TH2F("EE+","EE+ Tower TPG FineGrain",22, 0, 22, 22, 0, 22);
     TH2F* endc_m = new TH2F("EE-","EE- Tower TPG FineGrain",22, 0, 22, 22, 0, 22);

     auto iov = iovs.front();
     std::shared_ptr<EcalTPGFineGrainTowerEE> payload = fetchPayload( std::get<1>(iov) );
     unsigned int run = std::get<0>(iov);
     double minEE = 0, maxEE = 1;


    if( payload.get() ){
        const EcalTPGFineGrainTowerEEMap &towerMap = (*payload).getMap();

        EcalTPGFineGrainTowerEEMapIterator it;
        for(it = towerMap.begin(); it != towerMap.end(); ++it) {
          if(EcalScDetId::validHashIndex((*it).first)) {
            EcalScDetId ttId = EcalScDetId::unhashIndex((*it).first); 

            int ix = ttId.ix();
            int iy = ttId.iy();
            int zside = ttId.zside();
            
            uint32_t weight = (uint32_t)((*it).second);

            if(zside == -1)
              endc_m->Fill(ix, iy, weight);
            else
              endc_p->Fill(ix, iy, weight);

            if(maxEE < weight)
              maxEE = weight;

            if(minEE>weight)
              minEE=weight;
          }
        }//tower map
    }//payload


    TCanvas canvas("CC map","CC map",800,800);
    TLatex t1;
    t1.SetNDC();
    t1.SetTextAlign(26);
    t1.SetTextSize(0.05);
    t1.DrawLatex(0.5, 0.96, Form("Ecal TPGFineGrain Tower EE, IOV %i", run));

    
    TPad* padem = new TPad("padem","padem", 0., 0.3, 0.45, 0.75);
    padem->Draw();
    TPad* padep = new TPad("padep","padep", 0.55, 0.3, 1., 0.75);
    padep->Draw();

    TLine* l = new TLine(0., 0., 72., 0.);
    l->Draw();

    padem->cd();
    DrawEE_Tower(endc_m,l,minEE,maxEE);

 
    padep->cd();
    DrawEE_Tower(endc_p,l,minEE,maxEE);

    std::string ImageName(m_imageFileName);
    canvas.SaveAs(ImageName.c_str());

    return true;
  }// fill method

};


}


// Register the classes as boost python plugin
PAYLOAD_INSPECTOR_MODULE(EcalTPGFineGrainTowerEE){
  PAYLOAD_INSPECTOR_CLASS(EcalTPGFineGrainTowerEEPlot);
}