File indexing completed on 2024-04-06 12:09:09
0001
0002 #include "DQM/TrackingMonitorClient/plugins/TrackingDQMClientHeavyIons.h"
0003
0004 #include "DQMServices/Core/interface/DQMStore.h"
0005 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0006 #include "FWCore/ServiceRegistry/interface/Service.h"
0007 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0008
0009 #include <TH1F.h>
0010 #include <TClass.h>
0011 #include <TString.h>
0012 #include <string>
0013 #include <cmath>
0014 #include <climits>
0015 #include <boost/tokenizer.hpp>
0016
0017 using namespace std;
0018 using namespace edm;
0019
0020 typedef dqm::harvesting::MonitorElement ME;
0021
0022 TrackingDQMClientHeavyIons::TrackingDQMClientHeavyIons(const edm::ParameterSet& pset) {
0023 TopFolder_ = pset.getParameter<std::string>("FolderName");
0024 }
0025
0026 void TrackingDQMClientHeavyIons::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
0027 typedef vector<string> vstring;
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044 ibooker.cd();
0045 ibooker.setCurrentFolder(TopFolder_);
0046
0047 histName = "DCAStats_";
0048 DCAStats = ibooker.book2D(histName, histName, 2, 0, 2, 4, 0, 4);
0049 DCAStats->getTH2F()->GetYaxis()->SetBinLabel(1, "Mean");
0050 DCAStats->getTH2F()->GetYaxis()->SetBinLabel(2, "RMS, #sigma");
0051 DCAStats->getTH2F()->GetYaxis()->SetBinLabel(3, "Skewness ,#gamma_{1}");
0052 DCAStats->getTH2F()->GetYaxis()->SetBinLabel(4, "Kurtosis, #gamma_{2}");
0053 DCAStats->setBinLabel(1, "Longitudinal");
0054 DCAStats->setBinLabel(2, "Transverse");
0055 DCAStats->setOption("text");
0056
0057 histName = "LongDCASig_HeavyIonTk";
0058 ME* element = igetter.get(TopFolder_ + "/" + histName);
0059
0060 DCAStats->setBinContent(1, 1, element->getTH1F()->GetMean());
0061 DCAStats->setBinContent(1, 2, element->getTH1F()->GetRMS());
0062 DCAStats->setBinContent(1, 3, element->getTH1F()->GetSkewness());
0063 DCAStats->setBinContent(1, 4, element->getTH1F()->GetKurtosis());
0064
0065 histName = "TransDCASig_HeavyIonTk";
0066 ME* element1 = igetter.get(TopFolder_ + "/" + histName);
0067
0068 DCAStats->setBinContent(2, 1, element1->getTH1F()->GetMean());
0069 DCAStats->setBinContent(2, 2, element1->getTH1F()->GetRMS());
0070 DCAStats->setBinContent(2, 3, element1->getTH1F()->GetSkewness());
0071 DCAStats->setBinContent(2, 4, element1->getTH1F()->GetKurtosis());
0072 }