Functions
makeContourPlots.C File Reference
#include "Utilities.h"

Go to the source code of this file.

Functions

TCanvas * InitialiseCanvas ()
 
void DrawLegend (LegendMap &legobj, int textCol=kBlack)
 
void DrawDelChisqAndContours (TH2D *delchisq, ContourCollection &contours)
 
void DrawUniverses (TH2D *background, std::vector< ContourCollection > &uniContours, ContourCollection &asimov, ContourCollection &median, CLNumber cln)
 
void DrawHeatMaps (std::vector< TH2F * > heatMaps, ContourCollection asimovContours, ContourCollection medianContours, CLNumber cln, int max)
 
void DrawHiddenParameter (TH2D *h, std::pair< float, float > const &range)
 
void makeContourPlots (std::string file, std::vector< int > contours, std::string xlabel, std::string ylabel)
 

Function Documentation

void DrawDelChisqAndContours ( TH2D *  delchisq,
ContourCollection &  contours 
)

Definition at line 55 of file makeContourPlots.C.

References canv, DrawContour(), DrawLegend(), GetBestFitPoint(), InitialiseCanvas(), kAsimov, kEmpty, kPTRed, string, xaxis, and yaxis.

Referenced by makeContourPlots().

57 {
58  TCanvas* canv = InitialiseCanvas();
59  delchisq->Draw("colz");
60  DrawContour(contours, CLNumber::kEmpty, kWhite);
61 
62  LegendMap legmap = {
63  {contours[0].clnLatex, contours[0].contours[0]},
64  {contours[1].clnLatex, contours[1].contours[0]},
65  {contours[2].clnLatex, contours[2].contours[0]},
66  };
67 
68  TGraph* bfp = nullptr;
69  if (contours[0].clt == CLType::kAsimov){
70  bfp = GetBestFitPoint(CLType::kAsimov);
71  bfp->SetMarkerStyle(29);
72  //bfp->SetMarkerSize(0.5);
73  bfp->SetMarkerColor(kPTRed);
74  bfp->Draw("same p");
75  legmap["Best Fit"] = bfp;
76  }
77 
78  DrawLegend(legmap, kWhite);
79 
80  std::string name =
81  "DelChiSqAndContours" +
82  contours[0].cltName +
83  xaxis + yaxis;
84 
85 
86  canv->SaveAs((std::string(name+".png")).c_str());
87  canv->SaveAs((std::string(name+".pdf")).c_str());
88  delete canv;
89 }
name
const XML_Char * name
Definition: expat.h:151
yaxis
std::string yaxis
Definition: makeContourPlotsOld.C:19
canv
TCanvas * canv
Definition: checkMeasured_Target.C:142
LegendMap
std::map< std::string, TObject * > LegendMap
Definition: Utilities.h:135
xaxis
std::string xaxis
Definition: makeContourPlotsOld.C:18
GetBestFitPoint
TGraph * GetBestFitPoint(CLType clt, int universeNum=-1)
Definition: Utilities.h:436
DrawLegend
void DrawLegend(LegendMap &legobj, int textCol=kBlack)
Definition: makeContourPlots.C:38
kPTRed
Int_t kPTRed
Definition: Utilities.h:50
DrawContour
void DrawContour(ContourCollection cc, CLNumber cln=CLNumber::kEmpty, int colorOverride=-1, int styleOverride=-1)
Definition: Utilities.h:180
InitialiseCanvas
TCanvas * InitialiseCanvas()
Definition: makeContourPlots.C:15
string
enum BeamMode string
void DrawHeatMaps ( std::vector< TH2F * >  heatMaps,
ContourCollection  asimovContours,
ContourCollection  medianContours,
CLNumber  cln,
int  max 
)

Definition at line 132 of file makeContourPlots.C.

References canv, CLName, DrawContour(), DrawLegend(), InitialiseCanvas(), kPTOrange, kPTRed, string, xaxis, and yaxis.

Referenced by makeContourPlots().

137 {
138  TCanvas* canv = InitialiseCanvas();
139  heatMaps[cln]->SetMaximum(max);
140  heatMaps[cln]->Draw("colz");
141  DrawContour(asimovContours, cln, kPTOrange);
142  DrawContour(medianContours, cln, kPTRed);
143 
144  LegendMap legmap = {
145  {"Asimov", asimovContours[0].contours[0]},
146  {"Median", medianContours[0].contours[0]}
147  };
148  DrawLegend(legmap);
149 
150  std::string name = xaxis + yaxis + CLName[cln] + "HeatMaps";
151 
152  canv->SaveAs((std::string(name + ".png")).c_str());
153  canv->SaveAs((std::string(name + ".pdf")).c_str());
154  delete canv;
155 
156 }
name
const XML_Char * name
Definition: expat.h:151
yaxis
std::string yaxis
Definition: makeContourPlotsOld.C:19
canv
TCanvas * canv
Definition: checkMeasured_Target.C:142
LegendMap
std::map< std::string, TObject * > LegendMap
Definition: Utilities.h:135
xaxis
std::string xaxis
Definition: makeContourPlotsOld.C:18
DrawLegend
void DrawLegend(LegendMap &legobj, int textCol=kBlack)
Definition: makeContourPlots.C:38
CLName
std::vector< std::string > CLName
Definition: Utilities.h:15
kPTRed
Int_t kPTRed
Definition: Utilities.h:50
DrawContour
void DrawContour(ContourCollection cc, CLNumber cln=CLNumber::kEmpty, int colorOverride=-1, int styleOverride=-1)
Definition: Utilities.h:180
InitialiseCanvas
TCanvas * InitialiseCanvas()
Definition: makeContourPlots.C:15
cet::sqlite::max
T max(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:68
kPTOrange
Int_t kPTOrange
Definition: Utilities.h:50
string
enum BeamMode string
void DrawHiddenParameter ( TH2D *  h,
std::pair< float, float > const &  range 
)

Definition at line 159 of file makeContourPlots.C.

References canv, InitialiseCanvas(), string, xaxis, and yaxis.

Referenced by makeContourPlots().

160  {
161  gStyle->SetPalette(kCandy);
162  TCanvas* canv = InitialiseCanvas();
163  h->SetMaximum(range.second);
164  h->SetMinimum(range.first);
165  h->SetContour(1000);
166  h->Draw("colz");
167 
168  std::string name = (std::string)h->GetName();
169  canv->SaveAs((std::string(xaxis + yaxis + name + ".png")).c_str());
170  canv->SaveAs((std::string(xaxis + yaxis + name + ".pdf")).c_str());
171 }
name
const XML_Char * name
Definition: expat.h:151
yaxis
std::string yaxis
Definition: makeContourPlotsOld.C:19
canv
TCanvas * canv
Definition: checkMeasured_Target.C:142
xaxis
std::string xaxis
Definition: makeContourPlotsOld.C:18
make_syst_table_plots.h
h
Definition: make_syst_table_plots.py:240
InitialiseCanvas
TCanvas * InitialiseCanvas()
Definition: makeContourPlots.C:15
string
enum BeamMode string
void DrawLegend ( LegendMap &  legobj,
int  textCol = kBlack 
)

Definition at line 38 of file makeContourPlots.C.

References MECModelEnuComparisons::leg, and lo().

Referenced by ana::DataMCPair::CName(), DrawDelChisqAndContours(), DrawHeatMaps(), DrawUniverses(), DataMCPair::OverlayDataMCSyst(), and DataMCPair::OverlayDataMCSystNorm().

38  {
39  TLegend* leg = new TLegend(0.65, 0.15, 0.85, 0.15+(0.05*legobj.size()));
40  leg->SetLineWidth(0);
41  leg->SetFillStyle(0);
42  leg->SetTextColor(textCol);
43 
44  for (auto& lo : legobj){
45  if (lo.first == "Best Fit")
46  leg->AddEntry(lo.second, lo.first.c_str(), "p");
47  else
48  leg->AddEntry(lo.second, lo.first.c_str(), "l");
49  }
50 
51  leg->Draw("same");
52 }
lo
TSpline3 lo("lo", xlo, ylo, 12,"0")
void DrawUniverses ( TH2D *  background,
std::vector< ContourCollection > &  uniContours,
ContourCollection &  asimov,
ContourCollection &  median,
CLNumber  cln 
)

Definition at line 92 of file makeContourPlots.C.

References canv, CLName, CLNameLaTeX, DrawContour(), DrawLegend(), MECModelEnuComparisons::i, InitialiseCanvas(), string, xaxis, and yaxis.

Referenced by makeContourPlots().

97 {
98  // set the title for the background histogram to reflect the confidence level
99  background->SetTitle(("Universe Contours - " + CLNameLaTeX[cln]).c_str());
100 
101  TCanvas* canv = InitialiseCanvas();
102  background->Draw();
103  for (int i = 0; i < uniContours.size(); ++i){
104  DrawContour(uniContours[i], cln, -1, 1);
105  }
106  DrawContour(asimov, cln, -1, 1);
107  DrawContour(median, cln, -1, 1);
108 
109  //TGraph* bfp = nullptr;
110  //bfp = GetUniverseBestFitPoints();
111  //bfp->SetMarkerStyle(29);
112  //bfp->SetMarkerColor(kPTRed);
113  //bfp->Draw("same p");
114 
115  LegendMap legmap = {
116  {"Asimov", asimov[0].contours[0]},
117  {"Median", median[0].contours[0]},
118  {"Universes", uniContours[0][0].contours[0]},
119  //{"Best Fit", bfp}
120  };
121 
122  DrawLegend(legmap);
123 
124  std::string name = xaxis + yaxis + CLName[cln] + "UniverseContours";
125 
126  canv->SaveAs((std::string(name+".png")).c_str());
127  canv->SaveAs((std::string(name+".pdf")).c_str());
128  delete canv;
129 }
name
const XML_Char * name
Definition: expat.h:151
yaxis
std::string yaxis
Definition: makeContourPlotsOld.C:19
canv
TCanvas * canv
Definition: checkMeasured_Target.C:142
LegendMap
std::map< std::string, TObject * > LegendMap
Definition: Utilities.h:135
xaxis
std::string xaxis
Definition: makeContourPlotsOld.C:18
CLNameLaTeX
std::vector< std::string > CLNameLaTeX
Definition: Utilities.h:24
DrawLegend
void DrawLegend(LegendMap &legobj, int textCol=kBlack)
Definition: makeContourPlots.C:38
CLName
std::vector< std::string > CLName
Definition: Utilities.h:15
median
double median(TH1D *h)
Definition: absCal.cxx:524
DrawContour
void DrawContour(ContourCollection cc, CLNumber cln=CLNumber::kEmpty, int colorOverride=-1, int styleOverride=-1)
Definition: Utilities.h:180
InitialiseCanvas
TCanvas * InitialiseCanvas()
Definition: makeContourPlots.C:15
string
enum BeamMode string
TCanvas* InitialiseCanvas ( )

Definition at line 15 of file makeContourPlots.C.

References canv, xaxis, and yaxis.

Referenced by DrawDelChisqAndContours(), DrawHeatMaps(), DrawHiddenParameter(), and DrawUniverses().

15  {
16 
17  TCanvas* canv = new TCanvas("canv", "", 100, 100);
18  int canvx = 500;
19  int canvy = 500;
20  if (xaxis == "Th24" && yaxis == "Dmsq41"){
21  canvx = 400;
22  canvy = 600;
23  canv->SetLogx();
24  canv->SetLogy();
25  }
26  if (xaxis == "Th23" && yaxis == "Dmsq32"){
27  canvx = 600;
28  canvy = 400;
29  }
30 
31  canv->SetCanvasSize(canvx, canvy);
32  canv->SetRightMargin(0.18);
33 
34  return canv;
35 }
yaxis
std::string yaxis
Definition: makeContourPlotsOld.C:19
canv
TCanvas * canv
Definition: checkMeasured_Target.C:142
xaxis
std::string xaxis
Definition: makeContourPlotsOld.C:18
void makeContourPlots ( std::string  file,
std::vector< int >  contours,
std::string  xlabel,
std::string  ylabel 
)

Definition at line 175 of file makeContourPlots.C.

References DrawDelChisqAndContours(), DrawHeatMaps(), DrawHiddenParameter(), DrawUniverses(), e, GetBestFitPoint(), GetContours(), GetDeltaChiSqr(), GetHeatMaps(), GetHiddenParameter(), GetUniverseBestFitPoints(), GetUniverseContours(), kAsimov, kMedian, kNoType, kOneSigma, kThreeSigma, kTwoSigma, make_pair(), Save1DContour(), SetCMFColorPalette(), SetGlobalVariables(), xaxis, and yaxis.

178  {
179 
180  SetGlobalVariables(file, xlabel, ylabel, contours);
181  SetCMFColorPalette();
182 
183  // get the contours and other plots we're going to want to use
184  TH2D* background = GetDeltaChiSqr(CLType::kNoType);
185  TH2D* asimovDeltaChisq = GetDeltaChiSqr(CLType::kAsimov);
186  TH2D* medianDeltaChisq = GetDeltaChiSqr(CLType::kMedian);
187  ContourCollection asimovContours = GetContours(CLType::kAsimov);
188  ContourCollection medianContours = GetContours(CLType::kMedian);
189  std::vector<ContourCollection> uniContours = GetUniverseContours();
190  std::vector<TH2F*> heatMaps = GetHeatMaps();
191  TGraph* uniBestFitPoints = GetUniverseBestFitPoints();
192  TGraph* asimovBestFit = GetBestFitPoint(CLType::kAsimov);
193  std::vector<TH2D*> hiddenParVec;
194  std::vector<std::pair<float, float>> hiddenParRange;
195 
196  if(xlabel == "Th23" && ylabel == "Dmsq32"){
197  hiddenParVec.emplace_back(GetHiddenParameter("Th13"));
198  hiddenParRange.emplace_back(std::make_pair(0.14, 0.15));
199  hiddenParVec.emplace_back(GetHiddenParameter("dCP"));
200  hiddenParRange.emplace_back(std::make_pair(0, 6.283185));
201  }
202  else if(xlabel == "Th24" && ylabel == "Dmsq41"){
203  hiddenParVec.emplace_back(GetHiddenParameter("Th34"));
204  hiddenParRange.emplace_back(std::make_pair(0., 0.76));
205  hiddenParVec.emplace_back(GetHiddenParameter("Th23"));
206  hiddenParRange.emplace_back(std::make_pair(0.4, 0.65));
207  hiddenParVec.emplace_back(GetHiddenParameter("Dmsq32"));
208  hiddenParRange.emplace_back(std::make_pair(2.45, 2.57));
209  }
210  else if(xlabel == "Th24" && ylabel == "Th34"){
211  hiddenParVec.emplace_back(GetHiddenParameter("Dmsq41"));
212  hiddenParRange.emplace_back(std::make_pair(0., 100.));
213  hiddenParVec.emplace_back(GetHiddenParameter("Th23"));
214  hiddenParRange.emplace_back(std::make_pair(0.4, 0.65));
215  hiddenParVec.emplace_back(GetHiddenParameter("Dmsq32"));
216  hiddenParRange.emplace_back(std::make_pair(2.45e-3, 2.57e-3));
217  }
218  Save1DContour(kAsimov, xaxis);
219  Save1DContour(kAsimov, yaxis);
220 
221  DrawDelChisqAndContours(asimovDeltaChisq, asimovContours);
222  DrawDelChisqAndContours(medianDeltaChisq, medianContours);
223  DrawUniverses(background, uniContours, asimovContours, medianContours, kOneSigma);
224  DrawUniverses(background, uniContours, asimovContours, medianContours, kTwoSigma);
225  DrawUniverses(background, uniContours, asimovContours, medianContours, kThreeSigma);
226  DrawHeatMaps(heatMaps, asimovContours, medianContours, kOneSigma, uniContours.size());
227  DrawHeatMaps(heatMaps, asimovContours, medianContours, kTwoSigma, uniContours.size());
228  DrawHeatMaps(heatMaps, asimovContours, medianContours, kThreeSigma, uniContours.size());
229 
230  for(size_t hp = 0; hp < hiddenParVec.size(); ++hp)
231  DrawHiddenParameter(hiddenParVec[hp], hiddenParRange[hp]);
232 
233 }
DrawDelChisqAndContours
void DrawDelChisqAndContours(TH2D *delchisq, ContourCollection &contours)
Definition: makeContourPlots.C:55
ContourCollection
std::vector< Contour > ContourCollection
Definition: Utilities.h:134
GetDeltaChiSqr
TH2D * GetDeltaChiSqr(CLType clt, int universeNum=-1)
Definition: Utilities.h:381
SetCMFColorPalette
void SetCMFColorPalette()
Definition: Utilities.h:138
yaxis
std::string yaxis
Definition: makeContourPlotsOld.C:19
GetHiddenParameter
TH2D * GetHiddenParameter(std::string hiddenParName)
Definition: Utilities.h:539
make_pair
std::pair< Spectrum *, CheatDecomp * > make_pair(SpectrumLoaderBase &loader_data, SpectrumLoaderBase &loader_mc, HistAxis *axis, Cut *cut, const SystShifts &shift, const Var &wei)
Definition: DataMCLoad.C:336
xaxis
std::string xaxis
Definition: makeContourPlotsOld.C:18
GetUniverseContours
std::vector< ContourCollection > GetUniverseContours()
Definition: Utilities.h:509
DrawUniverses
void DrawUniverses(TH2D *background, std::vector< ContourCollection > &uniContours, ContourCollection &asimov, ContourCollection &median, CLNumber cln)
Definition: makeContourPlots.C:92
GetHeatMaps
std::vector< TH2F * > GetHeatMaps()
Definition: Utilities.h:523
GetBestFitPoint
TGraph * GetBestFitPoint(CLType clt, int universeNum=-1)
Definition: Utilities.h:436
GetUniverseBestFitPoints
TGraph * GetUniverseBestFitPoints()
Definition: Utilities.h:456
DrawHiddenParameter
void DrawHiddenParameter(TH2D *h, std::pair< float, float > const &range)
Definition: makeContourPlots.C:159
DrawHeatMaps
void DrawHeatMaps(std::vector< TH2F * > heatMaps, ContourCollection asimovContours, ContourCollection medianContours, CLNumber cln, int max)
Definition: makeContourPlots.C:132
file
TFile * file
Definition: cellShifts.C:17
SetGlobalVariables
void SetGlobalVariables(std::string f, std::string xax, std::string yax, std::vector< int > contours)
Definition: Utilities.h:52
e
Float_t e
Definition: plot.C:35
Save1DContour
void Save1DContour(CLType clt, std::string var, int universeNum=-1)
Definition: Utilities.h:217
GetContours
ContourCollection GetContours(CLType clt, int universeNum=-1)
Definition: Utilities.h:469
Generated on Thu Feb 25 2021 14:16:20 for NOvA by   doxygen 1.8.11