macro/plot/CalcDecile_RhoFluct.cc

0001 root -l hadd_RhoFluct.root <<EOF
0002
0003   TH1D* hparamt = new TH1D("hparamt", ";impactparam;N_{events}", 10000, 0.,22);
0004   TH1D* gross   = new TH1D("gross", ";impact param [fm];N_{events}", 100, 0.,22);
0005   double x[11];
0006   double* q = new double[11];
0007   for (int i=0;i<11;++i) {
0008     x[i] = .1*i;
0009   }
0010
0011   T->Draw("impactparam>>hparamt");
0012   hparamt->GetQuantiles(11,q,x);
0013   for (int i=0;i<11;++i) {
0014     cout << " i("<<i<<") x("<<x[i]<<") q("<<q[i]<<")"<<endl;
0015   }
0016
0017   T->Draw("impactparam>>gross");
0018 TCanvas*c = new TCanvas("c","c");
0019   hparamt->Draw("PE");
0020   /* c->WaitPrimitive(); */
0021   TFile *fout = new TFile("CalcDecile_RhoFluct.root","recreate");
0022   TH1D* keep = (TH1D*) gross->Clone("h_impactparam");
0023   TH1D* bounds = new TH1D("bound",";[fm];", 10, q);
0024   keep->Write();
0025   bounds->Write();
0026   fout->Close();
0027
0028
0029 EOF