Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-10-25 09:36:33

0001 #include "CondCore/Utilities/interface/PayloadInspectorModule.h"
0002 #include "CondCore/Utilities/interface/PayloadInspector.h"
0003 #include "CondCore/CondDB/interface/Time.h"
0004 #include "DataFormats/EcalDetId/interface/EBDetId.h"
0005 #include "DataFormats/EcalDetId/interface/EEDetId.h"
0006 #include "CondCore/EcalPlugins/plugins/EcalDrawUtils.h"
0007 
0008 // the data format of the condition to be inspected
0009 #include "CondFormats/EcalObjects/interface/EcalTBWeights.h"
0010 
0011 #include "TH2F.h" // a 2-D histogram with four bytes per cell (float)
0012 #include "TCanvas.h"
0013 #include "TLine.h"
0014 #include "TStyle.h"
0015 #include "TLatex.h"//write mathematical equations.
0016 #include "TPave.h"
0017 #include "TPaveStats.h"
0018 #include <string>
0019 #include <fstream>
0020 
0021 namespace {
0022 
0023 /*******************************************************
0024  2d plot of Ecal TBWeights of 1 IOV
0025  *******************************************************/
0026 class EcalTBWeightsPlot: public cond::payloadInspector::PlotImage<EcalTBWeights>{
0027     public:
0028         EcalTBWeightsPlot():
0029             cond::payloadInspector::PlotImage<EcalTBWeights>("Ecal TBWeights - map "){
0030                 setSingleIov(true);
0031         }
0032 
0033     bool fill(const std::vector<std::tuple<cond::Time_t, cond::Hash> >& iovs)override {
0034         auto iov=iovs.front();
0035         std::shared_ptr <EcalTBWeights> payload = fetchPayload(std::get<1> (iov));
0036         unsigned int run=std::get<0> (iov);
0037         TH2F* align;
0038         int NbRows;
0039 
0040         if(payload.get()){
0041             EcalTBWeights::EcalTBWeightMap map=(*payload).getMap();
0042             NbRows=map.size();
0043             align=new TH2F("Ecal PTM Temperatures","EcalXtalGroupId          EcalTDCId          EcalWeightSet",
0044                 3,0,3,NbRows,0,NbRows);
0045         
0046             double row = NbRows-0.5;
0047 
0048             for (EcalTBWeights::EcalPTMTemperatureMap::const_iterator it =map.begin(); 
0049                 it != map.end();it++) {
0050                 
0051                 //uint32_t mapKey=it->first;
0052                 //float val=it->second;
0053 
0054                 //align->Fill(0.5,row,mapKey);
0055                 //align->Fill(1.5,row,val);
0056 
0057                 row--;
0058             }
0059             
0060         }else
0061             return false;
0062 
0063         gStyle->SetPalette(1);
0064         gStyle->SetOptStat(0);
0065         TCanvas canvas("CC map", "CC map", 1000, 1000);
0066         TLatex t1;
0067         t1.SetNDC();
0068         t1.SetTextAlign(26);
0069         t1.SetTextSize(0.05);
0070         t1.SetTextColor(2);
0071         t1.DrawLatex(0.5, 0.96,Form("Ecal TBWeights, IOV %i", run));
0072 
0073 
0074         TPad* pad = new TPad("pad", "pad", 0.0, 0.0, 1.0, 0.94);
0075         pad->Draw();
0076         pad->cd();
0077         align->Draw("TEXT");
0078         TLine* l = new TLine;
0079         l->SetLineWidth(1);
0080 
0081         for (int i = 1; i < NbRows; i++) {
0082           double y = (double) i;
0083           l = new TLine(0., y, 3., y);
0084           l->Draw();
0085         }
0086 
0087         for (int i = 1; i < 3; i++) {
0088           double x = (double) i;
0089           double y = (double) NbRows;
0090           l = new TLine(x, 0., x, y);
0091           l->Draw();
0092         }
0093 
0094         align->GetXaxis()->SetTickLength(0.);
0095         align->GetXaxis()->SetLabelSize(0.);
0096         align->GetYaxis()->SetTickLength(0.);
0097         align->GetYaxis()->SetLabelSize(0.);
0098 
0099         std::string ImageName(m_imageFileName);
0100         canvas.SaveAs(ImageName.c_str());
0101 
0102         return true;
0103     }
0104 };
0105 
0106 } //close namespace
0107 
0108 // Register the classes as boost python plugin
0109 PAYLOAD_INSPECTOR_MODULE(EcalTBWeights) {
0110   PAYLOAD_INSPECTOR_CLASS(EcalTBWeightsPlot);
0111 }