macros/HistoBaseFiles/ShSizeHist.C

0001
0002 void ShSizeHist(const char *part = "e")
0003 {
0004   int rad[11] = {0,1,2,3,4,5,6,7,8,9,100};
0005   int pz[8] = {1,2,4,8,16,32,40,50};
0006   char *calo[] = {"E","HO","HI","BH","MAG"};
0007   double x[11];
0008   double y[11];
0009   double xer[11] = {0};
0010   double yer[11];
0011   char hname[100];
0012   char fname[100];
0013   sprintf(fname,"%spercent.root",part);
0014   char area[100];
0015   sprintf(area,"%s",part);
0016   TFile *fin = TFile::Open(fname);
0017   gROOT->cd();
0018   for (int i=0; i<8; i++)
0019     {
0020       for (int j=0; j<5; j++)
0021     {
0022       for (int k=0; k<11;k++)
0023         {
0024           // TFile *fin = TFile::Open(fname);
0025           x[k] = rad[k];
0026           sprintf(hname,"%sSH%dRad%d%s",part,pz[i],rad[k],calo[j]);
0027           TH1 *h1 = (TH1F *) fin->Get(hname);
0028           if (! fin)
0029         {
0030           cout << "DEATH BY STUPIDITY" << endl;
0031         }
0032           cout << hname << endl;
0033           y[k] = h1->GetMean(1);
0034           yer[k] = h1->GetRMS(1);
0035         }
0036       gr = new TGraphErrors(10,x,y,xer,yer);
0037       char gname[100];
0038       sprintf(gname, "%s%d%sRADERR",part,pz[i],calo[j]);
0039       gr->SetTitle(gname);
0040       gr->SetName(gname);
0041       gr->Draw("ALP");
0042       gr->SetMarkerStyle(20);
0043       gr->SetMarkerColor(2);
0044       char fname[100];
0045       sprintf(fname,"%sShEdepErr.root",part);
0046       TFile *fout = TFile::Open(fname, "UPDATE");
0047       gr->Write();
0048       fout->Write();
0049       fout->Close();
0050       for (int m=0; m<11;m++)
0051         {
0052           y[m] = 0;
0053           yer[m] = 0;
0054           x[m] = 0;
0055         }
0056     }
0057     }
0058  fin->Close();
0059 }