File indexing completed on 2024-04-06 11:58:57
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[41]={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.};
0012 float y[41]={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.,0.};
0013
0014 y[0] = h1etacoefmin1->GetRMS();
0015 y[1] = h1etacoefmin2->GetRMS();
0016 y[2] = h1etacoefmin3->GetRMS();
0017 y[3] = h1etacoefmin4->GetRMS();
0018 y[4] = h1etacoefmin5->GetRMS();
0019 y[5] = h1etacoefmin6->GetRMS();
0020 y[6] = h1etacoefmin7->GetRMS();
0021 y[7] = h1etacoefmin8->GetRMS();
0022 y[8] = h1etacoefmin9->GetRMS();
0023 y[9] = h1etacoefmin10->GetRMS();
0024 y[10] = h1etacoefmin11->GetRMS();
0025 y[11] = h1etacoefmin12->GetRMS();
0026 y[12] = h1etacoefmin13->GetRMS();
0027 y[13] = h1etacoefmin14->GetRMS();
0028 y[14] = h1etacoefmin15->GetRMS();
0029 y[15] = h1etacoefmin16_1->GetRMS();
0030
0031 h1etacoefmin1->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin1.gif");
0032 h1etacoefmin2->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin2.gif");
0033 h1etacoefmin3->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin3.gif");
0034 h1etacoefmin4->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin4.gif");
0035 h1etacoefmin5->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin5.gif");
0036 h1etacoefmin6->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin6.gif");
0037 h1etacoefmin7->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin7.gif");
0038 h1etacoefmin8->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin8.gif");
0039 h1etacoefmin9->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin9.gif");
0040 h1etacoefmin10->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin10.gif");
0041 h1etacoefmin11->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin11.gif");
0042 h1etacoefmin12->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin12.gif");
0043 h1etacoefmin13->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin13.gif");
0044 h1etacoefmin14->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin14.gif");
0045 h1etacoefmin15->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin15.gif");
0046 h1etacoefmin16_1->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin16_1.gif");
0047
0048 h2etacoefmin1->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin1.gif");
0049 h2etacoefmin2->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin2.gif");
0050 h2etacoefmin3->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin3.gif");
0051 h2etacoefmin4->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin4.gif");
0052 h2etacoefmin5->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin5.gif");
0053 h2etacoefmin6->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin6.gif");
0054 h2etacoefmin7->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin7.gif");
0055 h2etacoefmin8->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin8.gif");
0056 h2etacoefmin9->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin9.gif");
0057 h2etacoefmin10->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin10.gif");
0058 h2etacoefmin11->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin11.gif");
0059 h2etacoefmin12->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin12.gif");
0060 h2etacoefmin13->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin13.gif");
0061 h2etacoefmin14->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin14.gif");
0062 h2etacoefmin15->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin15.gif");
0063 h2etacoefmin16_1->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin16_1.gif");
0064
0065 TFile f1("coefficients_219_val_endcap_minus_8.9mln.root");
0066
0067 y[16] = h1etacoefmin17->GetRMS();
0068 y[17] = h1etacoefmin18->GetRMS();
0069 y[18] = h1etacoefmin19->GetRMS();
0070 y[19] = h1etacoefmin20->GetRMS();
0071 y[20] = h1etacoefmin21->GetRMS();
0072 y[21] = h1etacoefmin22->GetRMS();
0073 y[22] = h1etacoefmin23->GetRMS();
0074 y[23] = h1etacoefmin24->GetRMS();
0075 y[24] = h1etacoefmin25->GetRMS();
0076 y[25] = h1etacoefmin26->GetRMS();
0077 y[26] = h1etacoefmin27->GetRMS();
0078 y[27] = h1etacoefmin28->GetRMS();
0079 y[28] = h1etacoefmin29->GetRMS();
0080
0081
0082 h1etacoefmin17->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin17.gif");
0083 h1etacoefmin18->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin18.gif");
0084 h1etacoefmin19->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin19.gif");
0085 h1etacoefmin20->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin20.gif");
0086 h1etacoefmin21->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin21.gif");
0087 h1etacoefmin22->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin22.gif");
0088 h1etacoefmin23->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin23.gif");
0089 h1etacoefmin24->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin24.gif");
0090 h1etacoefmin25->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin25.gif");
0091 h1etacoefmin26->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin26.gif");
0092 h1etacoefmin27->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin27.gif");
0093 h1etacoefmin28->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin28.gif");
0094 h1etacoefmin29->Draw();c2h->Print("GIF/ALL/8.9mln/h1etacoefmin29.gif");
0095
0096 h2etacoefmin17->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin17.gif");
0097 h2etacoefmin18->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin18.gif");
0098 h2etacoefmin19->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin19.gif");
0099 h2etacoefmin20->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin20.gif");
0100 h2etacoefmin21->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin21.gif");
0101 h2etacoefmin22->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin22.gif");
0102 h2etacoefmin23->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin23.gif");
0103 h2etacoefmin24->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin24.gif");
0104 h2etacoefmin25->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin25.gif");
0105 h2etacoefmin26->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin26.gif");
0106 h2etacoefmin27->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin27.gif");
0107 h2etacoefmin28->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin28.gif");
0108 h2etacoefmin29->Draw("box");c2h->Print("GIF/ALL/8.9mln/h2etacoefmin29.gif");
0109
0110
0111 TFile f2("coefficients_219_val_hf_minus_8.9mln.root");
0112
0113 y[29] = h1etacoefmin30->GetRMS();
0114 y[30] = h1etacoefmin31->GetRMS();
0115 y[31] = h1etacoefmin32->GetRMS();
0116 y[32] = h1etacoefmin33->GetRMS();
0117 y[33] = h1etacoefmin34->GetRMS();
0118 y[34] = h1etacoefmin35->GetRMS();
0119 y[35] = h1etacoefmin36->GetRMS();
0120 y[36] = h1etacoefmin37->GetRMS();
0121 y[37] = h1etacoefmin38->GetRMS();
0122 y[38] = h1etacoefmin39->GetRMS();
0123 y[39] = h1etacoefmin40->GetRMS();
0124 y[40] = h1etacoefmin41->GetRMS();
0125
0126
0127 h1etacoefmin30->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin30.gif");
0128 h1etacoefmin31->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin31.gif");
0129 h1etacoefmin32->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin32.gif");
0130 h1etacoefmin33->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin33.gif");
0131 h1etacoefmin34->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin34.gif");
0132 h1etacoefmin35->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin35.gif");
0133 h1etacoefmin36->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin36.gif");
0134 h1etacoefmin37->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin37.gif");
0135 h1etacoefmin38->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin38.gif");
0136 h1etacoefmin39->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin39.gif");
0137 h1etacoefmin40->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin40.gif");
0138 h1etacoefmin41->Draw(); c2h->Print("GIF/ALL/8.9mln/h1etacoefmin41.gif");
0139
0140 h2etacoefmin30->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin30.gif");
0141 h2etacoefmin31->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin31.gif");
0142 h2etacoefmin32->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin32.gif");
0143 h2etacoefmin33->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin33.gif");
0144 h2etacoefmin34->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin34.gif");
0145 h2etacoefmin35->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin35.gif");
0146 h2etacoefmin36->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin36.gif");
0147 h2etacoefmin37->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin37.gif");
0148 h2etacoefmin38->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin38.gif");
0149 h2etacoefmin39->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin39.gif");
0150 h2etacoefmin40->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin40.gif");
0151 h2etacoefmin41->Draw("box"); c2h->Print("GIF/ALL/8.9mln/h2etacoefmin41.gif");
0152
0153
0154 cout<<y[27]<<endl;
0155
0156 TH2F* h3= new TH2F(" "," ",2,0.,42.,100,0.,0.15);
0157 h3->GetXaxis()->SetTitle("ieta");
0158 h3->GetYaxis()->SetTitle("RMS");
0159 h3->Draw();
0160 Int_t np=41;
0161
0162 TGraph* gr1 = new TGraph(np,x,y);
0163 gr1->Draw("P");
0164 c2h->Print("GIF/accuracy_minus_8.9mln.gif");
0165 c2h->Print("GIF/accuracy_minus_8.9mln.C");
0166
0167 }