1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
#include "catch.hpp"
#include <iostream>
#include <numeric> // std::accumulate
#include "TCanvas.h"
#include "DQM/TrackerRemapper/interface/SiStripTkMaps.h"
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
TEST_CASE("SiStripTkMaps testing", "[SiStripTkMaps]") {
//_____________________________________________________________
SECTION("Check SiStrip Tk Maps plotting") {
gStyle->SetOptStat(0);
SiStripTkMaps theMap("COLZA L");
theMap.bookMap("testing SiStripTkMaps", "counts");
const auto detInfo =
SiStripDetInfoFileReader::read(edm::FileInPath(SiStripDetInfoFileReader::kDefaultFile).fullPath());
unsigned int count = 0;
for (const auto& it : detInfo.getAllData()) {
count++;
theMap.fill(it.first, count);
}
const auto filledIds = theMap.getTheFilledIds();
TCanvas c = TCanvas("c", "c");
theMap.drawMap(c, "");
c.SaveAs("SiStripsTkMaps.png");
std::cout << "SiStripTkMaps filled " << filledIds.size() << " DetIds" << std::endl;
REQUIRE(filledIds.size() == count);
}
//_____________________________________________________________
SECTION("Check empty SiStripTkMaps") {
SiStripTkMaps theMap("COLZA L");
theMap.bookMap("testing SiStripTkMaps", "counts");
TCanvas c = TCanvas("c", "c");
theMap.drawMap(c, "");
c.SaveAs("SiStripsEmptyTkMaps.png");
REQUIRE(true);
}
}
|