Skip to content

Commit

Permalink
Merge pull request #38841 from bsunanda/Run3-sim123
Browse files Browse the repository at this point in the history 
Run3-sim123 Update the validation scripts for HFNose in Validation/HGCalValidation and macros in SimG4CMS/Calo
  • Loading branch information
@cmsbuild
cmsbuild authored Jul 28, 2022
2 parents fa0034a + dec9747 commit 7178892
Show file tree
Hide file tree
Showing 10 changed files with 1,054 additions and 787 deletions.
403 changes: 232 additions & 171 deletions SimG4CMS/Calo/macros/MakeECPlot.C
View file

Large diffs are not rendered by default.

1,013 changes: 551 additions & 462 deletions SimG4CMS/Calo/macros/MakeGVPlots.C
View file

Large diffs are not rendered by default.

92 changes: 52 additions & 40 deletions SimG4CMS/Calo/macros/MakeHFNPlots.C
View file
Original file line number Diff line number Diff line change
Expand Up @@ -3,19 +3,23 @@
// Usage:
// .L MakeHFNPlots.C+g
//
// To make plot of layers from a file created using .SimHit/Digi/RecHitStudy
// To make plot of layers from a file created using SimHit/Digi/RecHit
// using the scripts runHFNoseSimHitStudy_cfg.py, runHFNoseDigiStudy_cfg.py,
// runHFNoseRecHitStudy_cfg.py in Validation/HGCalValidation
//
// makeLayerPlots(fname, type, todomin, todomax, tag, text, save)
//
// where
// fname std::string Name of the ROOT file (hfnSimHitD44tt.root)
// type int File type (0:SimHit; 2:Digi; 3:RecHit)
// todomin int Range of plots to make (minimum: -1
// todomax int maximum: 8)
// tag std::string To be added to the name of the canvas ("")
// text std::string To be added to the title of the histogram ("")
// save bool If the canvas is to be saved as jpg file (false)
// fname std::string Name of the ROOT file [hfnSimHitD94tt.root]
// type int File type (0:SimHit; 2:Digi; 3:RecHit) [0]
// todomin int Range of plots to make [0]
// todomax int (minimum: -1; maximum: 8) [7]
// tag std::string To be added to the name of the canvas [""]
// text std::string To be added to the title of the histogram [""]
// save bool If the canvas is to be saved as jpg file [false]
//
//////////////////////////////////////////////////////////////////////////////

#include <TCanvas.h>
#include <TChain.h>
#include <TFile.h>
Expand All @@ -37,51 +41,59 @@
#include <string>
#include <vector>

void makeLayerPlots(std::string fname="hfnSimHitD44tt.root", int type=0,
int todomin=0, int todomax=7, std::string tag="",
std::string text="", bool save=false) {

void makeLayerPlots(std::string fname = "hfnSimHitD94tt.root",
int type = 0,
int todomin = 0,
int todomax = 7,
std::string tag = "",
std::string text = "",
bool save = false) {
std::string dirnm[3] = {"hgcalSimHitStudy", "hfnoseDigiStudy", "hfnoseRecHitStudy"};
std::string units[3] = {" (mm)", " (cm)", " (cm)"};
std::string titlx[3] = {"SimHit", "Digi", "RecHit"};

gStyle->SetCanvasBorderMode(0); gStyle->SetCanvasColor(kWhite);
gStyle->SetPadColor(kWhite);
gStyle->SetCanvasBorderMode(0);
gStyle->SetCanvasColor(kWhite);
gStyle->SetPadColor(kWhite);
gStyle->SetOptStat(0);
if (type < 0 || type > 2) type = 0;
if (type < 0 || type > 2)
type = 0;
TFile *file = new TFile(fname.c_str());
if (file) {
TDirectory *dir = (TDirectory*)file->FindObjectAny(dirnm[type].c_str());
TDirectory *dir = (TDirectory *)file->FindObjectAny(dirnm[type].c_str());
char cname[100], name[100], title[100], xtitl[100], ytitl[100];
for (int i=todomin; i<=todomax; ++i) {
for (int i = todomin; i <= todomax; ++i) {
if (i < 0) {
sprintf (name, "RZ_HGCalHFNoseSensitive");
sprintf (title, "%s (%s)", text.c_str(), titlx[type].c_str());
sprintf (xtitl, "z %s", units[type].c_str());
sprintf (ytitl, "R %s", units[type].c_str());
sprintf(name, "RZ_HGCalHFNoseSensitive");
sprintf(title, "%s (%s)", text.c_str(), titlx[type].c_str());
sprintf(xtitl, "z %s", units[type].c_str());
sprintf(ytitl, "R %s", units[type].c_str());
} else {
sprintf (name, "XY_L%d", i+1);
sprintf (title, "%s (Layer %d %s)", text.c_str(), i+1, titlx[type].c_str());
sprintf (xtitl, "x %s", units[type].c_str());
sprintf (ytitl, "y %s", units[type].c_str());
sprintf(name, "XY_L%d", i + 1);
sprintf(title, "%s (Layer %d %s)", text.c_str(), i + 1, titlx[type].c_str());
sprintf(xtitl, "x %s", units[type].c_str());
sprintf(ytitl, "y %s", units[type].c_str());
}
TH2D* hist = (TH2D*)dir->FindObjectAny(name);
TH2D *hist = (TH2D *)dir->FindObjectAny(name);
std::cout << name << " read out at " << hist << std::endl;
if (hist != nullptr) {
sprintf (cname, "%s%s", name, tag.c_str());
TCanvas *pad = new TCanvas(cname,cname,500,500);
pad->SetRightMargin(0.10); pad->SetTopMargin(0.10);
hist->GetYaxis()->SetTitleOffset(1.2);
hist->GetYaxis()->SetTitle(ytitl);
hist->GetXaxis()->SetTitle(xtitl);
hist->SetTitle(title);
if (i < 0 && type == 0) hist->GetXaxis()->SetNdivisions(5);
hist->Draw("colz");
pad->Modified(); pad->Update();
if (save) {
sprintf (name, "c_%s.jpg", pad->GetName());
pad->Print(name);
}
sprintf(cname, "%s%s", name, tag.c_str());
TCanvas *pad = new TCanvas(cname, cname, 500, 500);
pad->SetRightMargin(0.10);
pad->SetTopMargin(0.10);
hist->GetYaxis()->SetTitleOffset(1.2);
hist->GetYaxis()->SetTitle(ytitl);
hist->GetXaxis()->SetTitle(xtitl);
hist->SetTitle(title);
if (i < 0 && type == 0)
hist->GetXaxis()->SetNdivisions(5);
hist->Draw("colz");
pad->Modified();
pad->Update();
if (save) {
sprintf(name, "c_%s.jpg", pad->GetName());
pad->Print(name);
}
}
}
}
Expand Down
222 changes: 166 additions & 56 deletions SimG4CMS/Calo/macros/MakeHitStudyPlots.C
View file
Original file line number Diff line number Diff line change
@@ -1,3 +1,28 @@
//////////////////////////////////////////////////////////////////////////////
//
// Usage:
// .L MakeHitStudyPlots.C+g
//
// To make plot of various quantities which are created by CaloSimHitStudy
// from one or two different settings
//
// makeHitStudyPlots(file1, tag1, file2, tag2, toDo, ratio, save, dirnm)
//
// where
// file1 std::string Name of the first ROOT file [analRun3Old.root]
// tag1 std::string Tag for the first file [Old]
// file2 std::string Name of the second ROOT file [analRun3New.root]
// tag2 std::string Tag for the second file [New]
// todo int The plot type to be made [0]
// if -1, 6 different types are plotted
// (3, 5, 8, 9, 10, 11)
// ratio bool if the ratio to be plotted [true]
// (works when both files are active)
// save bool If the canvas is to be saved as jpg file [false]
// dirnm std::string Name of the directory [CaloSimHitStudy]
//
//////////////////////////////////////////////////////////////////////////////

#include <TCanvas.h>
#include <TChain.h>
#include <TFile.h>
Expand All @@ -21,9 +46,10 @@ void makeHitStudyPlots(std::string file1 = "analRun3Old.root",
std::string tag1 = "Old",
std::string file2 = "analRun3New.root",
std::string tag2 = "New",
int todo = 0,
std::string dirnm = "CaloSimHitStudy",
bool save = false) {
int toDo = 0,
bool ratio = true,
bool save = false,
std::string dirnm = "CaloSimHitStudy") {
std::string names[18] = {"Edep",
"EdepEM",
"EdepHad",
Expand All @@ -44,12 +70,17 @@ void makeHitStudyPlots(std::string file1 = "analRun3Old.root",
"PtInc"};
int numb[18] = {9, 9, 9, 16, 9, 9, 9, 1, 1, 1, 16, 9, 9, 16, 1, 1, 1, 1};
int rebin[18] = {10, 10, 10, 1, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 1, 1, 1, 1};
int todos[6] = {3, 5, 8, 9, 10, 11};
bool debug(false);

gStyle->SetCanvasBorderMode(0);
gStyle->SetCanvasColor(kWhite);
gStyle->SetPadColor(kWhite);
gStyle->SetFillColor(kWhite);
gStyle->SetOptStat(1110);
if (ratio)
gStyle->SetOptStat(0);
else
gStyle->SetOptStat(1110);
TFile* file[2];
int nfile(0);
std::string tag(""), tags[2];
Expand All @@ -69,62 +100,141 @@ void makeHitStudyPlots(std::string file1 = "analRun3Old.root",
tag += tag2;
}
}
std::cout << "Use " << nfile << " files from " << file1 << " and " << file2 << std::endl;
for (int i1 = 0; i1 < numb[todo]; ++i1) {
int first(0);
double y1(0.90), dy(0.12);
double y2 = y1 - dy * nfile - 0.01;
TLegend* leg = new TLegend(0.74, y2 - nfile * 0.04, 0.89, y2);
leg->SetBorderSize(1);
leg->SetFillColor(10);
TCanvas* pad;
for (int ifile = 0; ifile < nfile; ++ifile) {
TDirectory* dir = (TDirectory*)file[ifile]->FindObjectAny(dirnm.c_str());
int todoMin = (toDo >= 0) ? 0 : 0;
int todoMax = (toDo >= 0) ? 0 : 5;
std::cout << "Use " << nfile << " files from " << file1 << " and " << file2 << " and look for " << todoMin << ":"
<< todoMax << std::endl;
for (int i0 = todoMin; i0 <= todoMax; ++i0) {
int todo = (todoMax == 0) ? toDo : todos[i0];
for (int i1 = 0; i1 < numb[todo]; ++i1) {
double y1(0.90), dy(0.12);
double y2 = y1 - dy * nfile - 0.01;
TLegend* leg = (ratio) ? (new TLegend(0.10, 0.86, 0.90, 0.90)) : (new TLegend(0.65, y2 - nfile * 0.04, 0.90, y2));
leg->SetBorderSize(1);
leg->SetFillColor(10);
TCanvas* pad;
TH1D* hist0[nfile];
char name[100], namec[100];
if (numb[todo] == 1) {
sprintf(name, "%s", names[todo].c_str());
sprintf(namec, "%s%s", names[todo].c_str(), tag.c_str());
} else {
sprintf(name, "%s%d", names[todo].c_str(), i1);
sprintf(namec, "%s%d%s", names[todo].c_str(), i1, tag.c_str());
}
TH1D* hist = static_cast<TH1D*>(dir->FindObjectAny(name));
std::cout << name << " read out at " << hist << " for " << tags[ifile] << std::endl;
if (hist != nullptr) {
hist->SetLineColor(first + 1);
hist->SetLineStyle(first + 1);
hist->GetYaxis()->SetTitleOffset(1.4);
if (rebin[todo] > 1)
hist->Rebin(rebin[todo]);
if (first == 0) {
pad = new TCanvas(namec, namec, 500, 500);
pad->SetRightMargin(0.10);
pad->SetTopMargin(0.10);
pad->SetLogy();
hist->Draw();
for (int ifile = 0; ifile < nfile; ++ifile) {
TDirectory* dir = (TDirectory*)file[ifile]->FindObjectAny(dirnm.c_str());
if (numb[todo] == 1) {
sprintf(name, "%s", names[todo].c_str());
sprintf(namec, "%s%s", names[todo].c_str(), tag.c_str());
} else {
hist->Draw("sames");
sprintf(name, "%s%d", names[todo].c_str(), i1);
sprintf(namec, "%s%d%s", names[todo].c_str(), i1, tag.c_str());
}
leg->AddEntry(hist, tags[ifile].c_str(), "lp");
pad->Update();
++first;
TPaveStats* st = ((TPaveStats*)hist->GetListOfFunctions()->FindObject("stats"));
if (st != NULL) {
st->SetLineColor(first);
st->SetTextColor(first);
st->SetY1NDC(y1 - dy);
st->SetY2NDC(y1);
st->SetX1NDC(0.65);
st->SetX2NDC(0.90);
y1 -= dy;
hist0[ifile] = static_cast<TH1D*>(dir->FindObjectAny(name));
if (debug)
std::cout << name << " read out at " << hist0[ifile] << " for " << tags[ifile] << std::endl;
}
if (!ratio) {
int first(0);
for (int ifile = 0; ifile < nfile; ++ifile) {
TH1D* hist(hist0[ifile]);
if (hist != nullptr) {
hist->SetLineColor(first + 1);
hist->SetLineStyle(first + 1);
hist->GetYaxis()->SetTitleOffset(1.4);
if (rebin[todo] > 1)
hist->Rebin(rebin[todo]);
if (first == 0) {
pad = new TCanvas(namec, namec, 500, 500);
pad->SetRightMargin(0.10);
pad->SetTopMargin(0.10);
pad->SetLogy();
hist->Draw();
} else {
hist->Draw("sames");
}
leg->AddEntry(hist, tags[ifile].c_str(), "lp");
pad->Update();
++first;
TPaveStats* st = ((TPaveStats*)hist->GetListOfFunctions()->FindObject("stats"));
if (st != NULL) {
st->SetLineColor(first);
st->SetTextColor(first);
st->SetY1NDC(y1 - dy);
st->SetY2NDC(y1);
st->SetX1NDC(0.65);
st->SetX2NDC(0.90);
y1 -= dy;
}
pad->Modified();
pad->Update();
leg->Draw("same");
pad->Update();
if (save) {
sprintf(name, "c_%s.pdf", pad->GetName());
pad->Print(name);
}
}
}
pad->Modified();
pad->Update();
leg->Draw("same");
pad->Update();
if (save) {
sprintf(name, "c_%s.pdf", pad->GetName());
pad->Print(name);
} else {
if (nfile == 2) {
int nbin = hist0[0]->GetNbinsX();
int nbinR = nbin / rebin[todo];
double xlow = hist0[0]->GetXaxis()->GetBinLowEdge(1);
double xhigh = hist0[0]->GetXaxis()->GetBinLowEdge(nbin) + hist0[0]->GetXaxis()->GetBinWidth(nbin);
;
if (numb[todo] == 1) {
sprintf(name, "%sRatio", names[todo].c_str());
sprintf(namec, "%sRatio%s", names[todo].c_str(), tag.c_str());
} else {
sprintf(name, "%s%dRatio", names[todo].c_str(), i1);
sprintf(namec, "%s%dRatio%s", names[todo].c_str(), i1, tag.c_str());
}
pad = new TCanvas(namec, namec, 500, 500);
TH1D* histr = new TH1D(name, hist0[0]->GetTitle(), nbinR, xlow, xhigh);
sprintf(name, "Ratio (%s/%s)", tags[0].c_str(), tags[1].c_str());
histr->GetXaxis()->SetTitle(hist0[0]->GetXaxis()->GetTitle());
histr->GetYaxis()->SetTitle(name);
histr->GetXaxis()->SetLabelOffset(0.005);
histr->GetXaxis()->SetTitleOffset(1.00);
histr->GetYaxis()->SetLabelOffset(0.005);
histr->GetYaxis()->SetTitleOffset(1.20);
histr->GetYaxis()->SetRangeUser(0.0, 2.0);
double sumNum(0), sumDen(0), maxDev(0);
for (int j = 0; j < nbinR; ++j) {
double tnum(0), tden(0), rnum(0), rden(0);
for (int i = 0; i < rebin[todo]; ++i) {
int ib = j * rebin[todo] + i + 1;
tnum += hist0[0]->GetBinContent(ib);
tden += hist0[1]->GetBinContent(ib);
rnum += ((hist0[0]->GetBinError(ib)) * (hist0[0]->GetBinError(ib)));
rden += ((hist0[1]->GetBinError(ib)) * (hist0[1]->GetBinError(ib)));
}
if (tden > 0 && tnum > 0) {
double rat = tnum / tden;
double err = rat * sqrt((rnum / (tnum * tnum)) + (rden / (tden * tden)));
histr->SetBinContent(j + 1, rat);
histr->SetBinError(j + 1, err);
double temp1 = (rat > 1.0) ? 1.0 / rat : rat;
double temp2 = (rat > 1.0) ? err / (rat * rat) : err;
sumNum += (fabs(1 - temp1) / (temp2 * temp2));
sumDen += (1.0 / (temp2 * temp2));
if (fabs(1 - temp1) > maxDev)
maxDev = fabs(1 - temp1);
}
}
histr->Draw();
sprintf(name, "%s vs %s", tag1.c_str(), tag2.c_str());
leg->AddEntry(histr, name, "lp");
leg->Draw("same");
pad->Update();
TLine* line = new TLine(xlow, 1.0, xhigh, 1.0);
line->SetLineColor(2);
line->SetLineWidth(2);
line->SetLineStyle(2);
line->Draw("same");
pad->Modified();
pad->Update();
sumNum = (sumDen > 0) ? (sumNum / sumDen) : 0;
sumDen = (sumDen > 0) ? 1.0 / sqrt(sumDen) : 0;
if (sumNum == 0)
sumDen = 0;
std::cout << tag1 << " vs " << tag2 << " " << hist0[0]->GetXaxis()->GetTitle() << " Mean deviation " << sumNum
<< " +- " << sumDen << " maximum " << maxDev << std::endl;
}
}
}
Expand Down
Loading

0 comments on commit 7178892

Please sign in to comment.