Back to home page

Project CMSSW displayed by LXR

 
 

    


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 }