Skip to content

Commit

Permalink
update BeamSpotOnline payload inspector to include parameters added in …
Browse files Browse the repository at this point in the history 
…cms-sw#35338
  • Loading branch information
@mmusich
mmusich committed Sep 28, 2021
1 parent c1cb2f1 commit 2ddb4d7
Show file tree
Hide file tree
Showing 2 changed files with 145 additions and 33 deletions.
47 changes: 29 additions & 18 deletions CondCore/BeamSpotPlugins/interface/BeamSpotPayloadInspectorHelper.h
View file
Expand Up @@ -28,21 +28,32 @@ namespace BeamSpotPI {
}

enum parameters {
X,
Y,
Z,
sigmaX,
sigmaY,
sigmaZ,
dxdz,
dydz,
lastLumi,
lastRun,
lastFill,
nTracks,
nPVs,
creationTime,
END_OF_TYPES
X = 0, // 0 regular BS methods
Y = 1, // 1
Z = 2, // 2
sigmaX = 3, // 3
sigmaY = 4, // 4
sigmaZ = 5, // 5
dxdz = 6, // 6
dydz = 7, // 7
lastLumi = 8, // 8 additional int's
lastRun = 9, // 9
lastFill = 10, // 10
nTracks = 11, // 11
nPVs = 12, // 12
nUsedEvents = 13, // 13
maxPVs = 14, // 14
meanPV = 15, // 15 additional float's
meanErrorPV = 16, // 16
rmsPV = 17, // 17
rmsErrorPV = 18, // 18
creationTime = 19, // 19 additional cond::Time_t
startTimeStamp = 20, // 20
endTimeStamp = 21, // 21
startTime = 22, // 22 additional std::string
endTime = 23, // 23
lumiRange = 24, // 24
END_OF_TYPES = 25,
};

/************************************************/
Expand Down Expand Up @@ -233,7 +244,7 @@ namespace BeamSpotPI {
canvas.cd(1)->SetTopMargin(0.05);
canvas.cd(1)->SetBottomMargin(0.06);
canvas.cd(1)->SetLeftMargin(0.15);
canvas.cd(1)->SetRightMargin(0.03);
canvas.cd(1)->SetRightMargin(0.01);
canvas.cd(1)->Modified();
canvas.cd(1)->SetGrid();

Expand Down Expand Up @@ -331,8 +342,8 @@ namespace BeamSpotPI {
canvas.cd(2);
canvas.cd(2)->SetTopMargin(0.05);
canvas.cd(2)->SetBottomMargin(0.06);
canvas.cd(2)->SetLeftMargin(0.15);
canvas.cd(2)->SetRightMargin(0.03);
canvas.cd(2)->SetLeftMargin(0.25);
canvas.cd(2)->SetRightMargin(0.01);
canvas.cd(2)->Modified();
canvas.cd(2)->SetGrid();

Expand Down
131 changes: 116 additions & 15 deletions CondCore/BeamSpotPlugins/plugins/BeamSpotOnline_PayloadInspector.cc
View file
Expand Up @@ -62,11 +62,17 @@ namespace {
public:
std::shared_ptr<TH2F> fillTheExtraHistogram() const override {
gStyle->SetHistMinimumZero();
auto h2_ExtraBSParameters = std::make_shared<TH2F>("ExtraParameters", "", 1, 0.0, 1.0, 6, 0, 6.);

const auto span = parameters::startTime - parameters::lastLumi;
edm::LogInfo("BeamSpotOnlineParameters") << "the span is" << span;

auto h2_ExtraBSParameters =
std::make_shared<TH2F>("ExtraParameters", "", 1, 0.0, 1.0, span, 0, static_cast<float>(span));
h2_ExtraBSParameters->SetStats(false);

std::function<int(parameters)> mycutFunctor = [this](parameters my_param) {
int ret(-999.);
//_____________________________________________________________________________
std::function<int(parameters)> myIntFunctor = [this](parameters my_param) {
int ret(-999);
switch (my_param) {
case lastLumi:
return this->m_payload->GetLastAnalyzedLumi();
Expand All @@ -78,28 +84,85 @@ namespace {
return this->m_payload->GetNumTracks();
case nPVs:
return this->m_payload->GetNumPVs();
case END_OF_TYPES:
case nUsedEvents:
return this->m_payload->GetUsedEvents();
case maxPVs:
return this->m_payload->GetMaxPVs();
default:
return ret;
}
};

//_____________________________________________________________________________
std::function<float(parameters)> myFloatFunctor = [this](parameters my_param) {
float ret(-999.);
switch (my_param) {
case meanPV:
return this->m_payload->GetMeanPV();
case meanErrorPV:
return this->m_payload->GetMeanErrorPV();
case rmsPV:
return this->m_payload->GetRmsPV();
case rmsErrorPV:
return this->m_payload->GetRmsErrorPV();
default:
return ret;
}
};

//_____________________________________________________________________________
std::function<std::string(parameters)> myStringFunctor = [this](parameters my_param) {
std::string ret("");
switch (my_param) {
case startTime:
return this->m_payload->GetStartTime();
case endTime:
return this->m_payload->GetEndTime();
case lumiRange:
return this->m_payload->GetLumiRange();
default:
return ret;
}
};

unsigned int yBin = 6;
for (int foo = parameters::lastLumi; foo != parameters::END_OF_TYPES; foo++) {
//_____________________________________________________________________________
std::function<cond::Time_t(parameters)> myTimeFunctor = [this](parameters my_param) {
cond::Time_t ret(1);
switch (my_param) {
case creationTime:
return this->m_payload->GetCreationTime();
case startTimeStamp:
return this->m_payload->GetStartTimeStamp();
case endTimeStamp:
return this->m_payload->GetEndTimeStamp();
default:
return ret;
}
};

unsigned int yBin = span;
for (int foo = parameters::lastLumi; foo != parameters::startTime; foo++) {
parameters param = static_cast<parameters>(foo);
std::string theLabel = this->getStringFromTypeEnum(param);
h2_ExtraBSParameters->GetYaxis()->SetBinLabel(yBin, theLabel.c_str());

edm::LogInfo("BeamSpotOnline_PayloadInspector")
<< theLabel.c_str() << " : " << mycutFunctor(param) << std::endl;

if (param == BeamSpotPI::creationTime) {
h2_ExtraBSParameters->SetBinContent(1, yBin, m_payload->GetCreationTime());
} else {
h2_ExtraBSParameters->SetBinContent(1, yBin, mycutFunctor(param));
if (foo <= parameters::maxPVs) {
const auto output = try_<int, std::out_of_range>(std::bind(myIntFunctor, param), print_error);
edm::LogInfo("BeamSpotOnline_PayloadInspector") << theLabel.c_str() << " : " << output << std::endl;
h2_ExtraBSParameters->SetBinContent(1, yBin, output);
} else if (foo <= parameters::rmsErrorPV) {
const auto output = try_<float, std::out_of_range>(std::bind(myFloatFunctor, param), print_error);
edm::LogInfo("BeamSpotOnline_PayloadInspector") << theLabel.c_str() << " : " << output << std::endl;
h2_ExtraBSParameters->SetBinContent(1, yBin, output);
} else if (foo <= parameters::endTimeStamp) {
const auto output = try_<cond::Time_t, std::out_of_range>(std::bind(myTimeFunctor, param), print_error);
edm::LogInfo("BeamSpotOnline_PayloadInspector") << theLabel.c_str() << " : " << output << std::endl;
h2_ExtraBSParameters->SetBinContent(1, yBin, output);
//} else if( foo <=parameters::lumiRange){
// const auto output = try_<std::string,std::out_of_range>(std::bind(myStringFunctor, param), print_error);
//edm::LogInfo("BeamSpotOnline_PayloadInspector") << theLabel.c_str() << " : " << output << std::endl;
//h2_ExtraBSParameters->SetBinContent(1, yBin, output);
}

yBin--;
}

Expand Down Expand Up @@ -141,12 +204,50 @@ namespace {
return "# tracks";
case nPVs:
return "# PVs";
case nUsedEvents:
return "# events";
case maxPVs:
return "max PVs";
case meanPV:
return "#LT # PV #GT";
case meanErrorPV:
return "#LT # PV #GT error";
case rmsPV:
return "RMS(# PV)";
case rmsErrorPV:
return "RMS(# PV) error";
case creationTime:
return "time";
return "creation time";
case startTimeStamp:
return "start timestamp";
case endTimeStamp:
return "end timestamp";
case startTime:
return "startTime";
case endTime:
return "endTime";
case lumiRange:
return "lumiRange";
default:
return "should never be here";
}
}

//slightly better error handler
static void print_error(const std::exception& e) { edm::LogError("BeamSpotOnlineParameters") << e.what() << '\n'; }

// method to catch exceptions
template <typename T, class Except, class Func, class Response>
T try_(Func f, Response r) const {
try {
LogDebug("BeamSpotOnlineParameters") << "I have tried" << std::endl;
return f();
} catch (Except& e) {
LogDebug("BeamSpotOnlineParameters") << "I have caught!" << std::endl;
r(e);
return static_cast<T>(-999);
}
}
};

} // namespace
Expand Down

0 comments on commit 2ddb4d7

Please sign in to comment.