Skip to content

B
Be_libraries
Switch branch/tag
Find file
Normal viewHistoryPermalink
figPresCorrGS.cxx 3.21 KB
Newer Older
Vratislav Chudoba's avatar
New scripts for presentation use.
3d649d46
 
Vratislav Chudoba committed
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 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 
//#include "TFile.h"
//#include "TCanvas.h"
//#include "TH1F.h"

void figPresCorrGS(const Bool_t saveFigs = 0)
{

	TFile *fr1 = new TFile("figures/PRCepsilonT.root", "READ");
//	TFile *fr2 = new TFile("figures/PRCepsilonT.root", "READ");

	TH1F *hExp[4];
	TH1F *hSim[4];
	TH1F *hInput[4];

	TH1F *hCurr = 0;

	hCurr = (TH1F*)fr1->Get("heeT0_0_4560");
	hExp[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("heeT0_5_6075");
	hExp[1] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("heeT0_5_7590");
	hExp[2] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("heeT0_0_90120");
	hExp[3] = new TH1F(*hCurr);

	hCurr = (TH1F*)fr1->Get("hseT0_0_4560");
	hSim[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("hseT0_5_6075");
	hSim[1] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("hseT0_5_7590");
	hSim[2] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("hseT0_0_90120");
	hSim[3] = new TH1F(*hCurr);

	hCurr = (TH1F*)fr1->Get("hieT0_0_4560");
	hInput[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("hieT0_5_6075");
	hInput[1] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("hieT0_5_7590");
	hInput[2] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr1->Get("hieT0_0_90120");
	hInput[3] = new TH1F(*hCurr);



	TCanvas *c1 = new TCanvas();
	c1->ToggleEditor();
//	c1->SetWindowSize(800, 1200);
	c1->SetWindowSize(1000, 1000);
	c1->Divide(2, 2, .001, .001);


	for (Int_t i = 0; i < 4; i++) {
		c1->cd(i+1);

//		cout << i << endl;

		hSim[i]->SetTitle("");
		hSim[i]->SetLineColor(kGray+1);
		hSim[i]->SetFillColor(kGray+1);
		hSim[i]->GetXaxis()->SetNdivisions(405);
		hSim[i]->GetYaxis()->SetTitle("");
		if (i == 0 || i == 1) {
			hSim[i]->GetYaxis()->SetNdivisions(1005);
		}
//		hInput[i]->SetLineWidth(2);
		hInput[i]->SetLineColor(kRed);

		hSim[i]->Draw("");
		hExp[i]->Draw("E same");
		hInput[i]->Draw("same");

		hExp[i]->Chi2Test(hSim[i], "P");
	}
	c1->Update();


	TCanvas *c2 = new TCanvas();
	c2->ToggleEditor();
	c2->SetWindowSize(1000, 500);
	c2->Divide(2, 1, .001, .001);

	TFile *fr3 = new TFile("figures/PRCepsilonY.root");
	TFile *fr4 = new TFile("figures/PRCcosThetaT.root");

	TH1F *hCorrExp[2];
	TH1F *hCorrSim[2];
	TH1F *hCorrInput[2];

	hCurr = (TH1F*)fr3->Get("heeY0_0");
	hCorrExp[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr4->Get("hecoskT0_0_7590");
	hCorrExp[1] = new TH1F(*hCurr);

	hCurr = (TH1F*)fr3->Get("hseY0_0");
	hCorrSim[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr4->Get("hscoskT0_0_7590");
	hCorrSim[1] = new TH1F(*hCurr);

	hCurr = (TH1F*)fr3->Get("hieY0_0");
	hCorrInput[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr4->Get("hicoskT0_0_7590");
	hCorrInput[1] = new TH1F(*hCurr);

	for (Int_t i = 0; i < 2; i++) {
		c2->cd(i+1);

		hCorrSim[i]->SetTitle("");
		hCorrSim[i]->SetLineColor(kGray+1);
		hCorrSim[i]->SetFillColor(kGray+1);
		hCorrSim[i]->GetXaxis()->SetNdivisions(505);
		hCorrSim[i]->GetYaxis()->SetTitle("");
		if (i == 1) {
			hCorrSim[i]->GetXaxis()->SetTitleSize(0.08);
			hCorrSim[i]->GetXaxis()->SetTitleOffset(0.73);
		}
//		hCorrInput[i]->SetLineWidth(2);
		hCorrInput[i]->SetLineColor(kRed);

		hCorrSim[i]->Draw("");
		hCorrExp[i]->Draw("E same");
		hCorrInput[i]->Draw("same");

		hCorrExp[i]->Chi2Test(hCorrSim[i], "P");
	}
	c2->Update();
	c1->Update();

	if (saveFigs) {
		c1->SaveAs("figures/PresGSepsilonT.png");
		c2->SaveAs("figures/PresGSotherCorrs.png");
	}

}