HcalCalibAlgos/macros/print_acc_minus_8.9mln_wcut.C

0001 {
0002   gROOT->Reset();
0003   gStyle->SetPalette(1);
0004   gStyle                  -> SetOptStat(1111);
0005
0006   TCanvas c2h("c2h","2-d options",10,10,800,600);
0007   c2h.SetFillColor(10);
0008   gPad->SetDrawOption("e1p");
0009   TFile f0("coefficients_219_val_barrel_minus_8.9mln.root");
0010
0011       float x[40]={1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,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.};
0012       float y[40]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.};
0013        y[0] = h1etacoefmin1l->GetRMS();
0014        y[1] = h1etacoefmin2l->GetRMS();
0015        y[2] = h1etacoefmin3l->GetRMS();
0016        y[3] = h1etacoefmin4l->GetRMS();
0017        y[4] = h1etacoefmin5l->GetRMS();
0018        y[5] = h1etacoefmin6l->GetRMS();
0019        y[6] = h1etacoefmin7l->GetRMS();
0020        y[7] = h1etacoefmin8l->GetRMS();
0021        y[8] = h1etacoefmin9l->GetRMS();
0022        y[9] = h1etacoefmin10l->GetRMS();
0023        y[10] = h1etacoefmin11l->GetRMS();
0024        y[11] = h1etacoefmin12l->GetRMS();
0025        y[12] = h1etacoefmin13l->GetRMS();
0026        y[13] = h1etacoefmin14l->GetRMS();
0027        y[14] = h1etacoefmin16_1l->GetRMS();
0028
0029        h1etacoefmin1l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin1l.gif");
0030        h1etacoefmin2l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin2l.gif");
0031        h1etacoefmin3l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin3l.gif");
0032        h1etacoefmin4l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin4l.gif");
0033        h1etacoefmin5l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin5l.gif");
0034        h1etacoefmin6l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin6l.gif");
0035        h1etacoefmin7l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin7l.gif");
0036        h1etacoefmin8l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin8l.gif");
0037        h1etacoefmin9l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin9l.gif");
0038        h1etacoefmin10l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin10l.gif");
0039        h1etacoefmin11l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin11l.gif");
0040        h1etacoefmin12l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin12l.gif");
0041        h1etacoefmin13l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin13l.gif");
0042        h1etacoefmin14l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin14l.gif");
0043        h1etacoefmin15l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin15l.gif");
0044        h1etacoefmin16_1l->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin16_1l.gif");
0045
0046        TFile f1("coefficients_219_val_endcap_minus_8.9mln.root");
0047        h1etacoefmin17a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin17a.gif");
0048        h1etacoefmin18a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin18a.gif");
0049        h1etacoefmin19a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin19a.gif");
0050        h1etacoefmin20a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin20a.gif");
0051        h1etacoefmin21a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin21a.gif");
0052        h1etacoefmin22a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin22a.gif");
0053        h1etacoefmin23a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin23a.gif");
0054        h1etacoefmin24a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin24a.gif");
0055        h1etacoefmin25a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin25a.gif");
0056        h1etacoefmin26a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin26a.gif");
0057        h1etacoefmin27a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin27a.gif");
0058        h1etacoefmin28a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin28a.gif");
0059        h1etacoefmin29a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin29a.gif");
0060
0061        h1etacoefmin17b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin17b.gif");
0062        h1etacoefmin18b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin18b.gif");
0063        h1etacoefmin19b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin19b.gif");
0064        h1etacoefmin20b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin20b.gif");
0065        h1etacoefmin21b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin21b.gif");
0066        h1etacoefmin22b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin22b.gif");
0067        h1etacoefmin23b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin23b.gif");
0068        h1etacoefmin24b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin24b.gif");
0069        h1etacoefmin25b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin25b.gif");
0070        h1etacoefmin26b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin26b.gif");
0071        h1etacoefmin27b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin27b.gif");
0072        h1etacoefmin28b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin28b.gif");
0073        h1etacoefmin29b->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin29b.gif");
0074
0075
0076        TFile f2("coefficients_219_val_hf_minus_8.9mln.root");
0077
0078        y[28] = h1etacoefmin30a->GetRMS();
0079        y[29] = h1etacoefmin31a->GetRMS();
0080        y[30] = h1etacoefmin32a->GetRMS();
0081        h1etacoefmin30a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin30a.gif");
0082        h1etacoefmin31a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin31a.gif");
0083        h1etacoefmin32a->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin32a.gif");
0084
0085
0086        /*
0087        y[15] = h1etacoefmin17_1->GetRMS();
0088        y[16] = h1etacoefmin18_1->GetRMS();
0089        y[17] = h1etacoefmin19_1->GetRMS();
0090        y[18] = h1etacoefmin20_1->GetRMS();
0091        y[19] = h1etacoefmin21_1->GetRMS();
0092        y[20] = h1etacoefmin22_1->GetRMS();
0093        y[21] = h1etacoefmin23_1->GetRMS();
0094        y[22] = h1etacoefmin24_1->GetRMS();
0095        y[23] = h1etacoefmin25_1->GetRMS();
0096        y[24] = h1etacoefmin26_1->GetRMS();
0097        y[25] = h1etacoefmin27_1->GetRMS();
0098        y[26] = h1etacoefmin28_1->GetRMS();
0099        y[27] = h1etacoefmin29_1->GetRMS();
0100        y[28] = h1etacoefmin30_1->GetRMS();
0101        y[29] = h1etacoefmin31_1->GetRMS();
0102        y[30] = h1etacoefmin32_1->GetRMS();
0103        y[31] = h1etacoefmin33_1->GetRMS();
0104        y[32] = h1etacoefmin34_1->GetRMS();
0105        y[33] = h1etacoefmin35_1->GetRMS();
0106        y[34] = h1etacoefmin36_1->GetRMS();
0107        y[35] = h1etacoefmin37_1->GetRMS();
0108        y[36] = h1etacoefmin38_1->GetRMS();
0109        y[37] = h1etacoefmin39_1->GetRMS();
0110        y[38] = h1etacoefmin40_1->GetRMS();
0111        y[39] = h1etacoefmin41_1->GetRMS();
0112     */
0113
0114   TH2F* h3= new TH2F(" "," ",2,0.,42.,100,0.,0.15);
0115   h3->GetXaxis()->SetTitle("ieta");
0116   h3->GetYaxis()->SetTitle("RMS");
0117   h3->Draw();
0118   Int_t np=40;
0119
0120   TGraph*  gr1 = new TGraph(np,x,y);
0121   gr1->Draw("P");
0122   c2h->Print("GIF/accuracy_minus_8.9mlnl.gif");
0123   c2h->Print("GIF/accuracy_minus_8.9mlnl.C");
0124
0125
0126 //           h1etacoefmin2D1->Draw("box");c2h->Print("GIF/h1etacoefminD1.gif");
0127 //           h1etacoefmin2D1_noise->Draw("box");c2h->Print("GIF/h1etacoefminD1_noise.gif");
0128 //           h1etacoefmin2D1->Draw("box");c2h->Print("GIF/h1etacoefminD1.gif");
0129 //           h1etacoefmin2D1_noise->Draw("box");c2h->Print("GIF/h1etacoefminD1_noise.gif");
0130
0131 //           h1etacoefmin2D16->Draw("box");c2h->Print("GIF/h1etacoefminD16.gif");
0132 //           h1etacoefmin2D16_noise->Draw("box");c2h->Print("GIF/h1etacoefminD16_noise.gif");
0133 //           h1etacoefmin2D16->Draw("box");c2h->Print("GIF/h1etacoefminD16.gif");
0134 //           h1etacoefmin2D16_noise->Draw("box");c2h->Print("GIF/h1etacoefminD16_noise.gif");
0135
0136 }