macros/HistoBaseFiles/Surf3DHistoMaker.C

0001 void Surf3DHistoMaker(int ncol, int gcol)
0002 {
0003   TFile *fin = TFile::Open("neutron_gamma_edep.root");
0004   gROOT->cd();
0005   TCanvas *c1 = new TCanvas();
0006   c1->SetLogz();
0007   TH2 *h2 = (TH2F *) fin->Get("neutron")->Clone();
0008   h2->SetXTitle("dtheta");
0009   h2->SetYTitle("dphi");
0010   h2->SetZTitle("Energy");
0011   h2->GetXaxis()->SetTitleOffset(2);
0012   h2->GetYaxis()->SetTitleOffset(2);
0013   h2->GetZaxis()->SetTitleOffset(1.2);
0014   TH2 *h2a = (TH2F *) fin->Get("gamma")->Clone();
0015   h2a->SetXTitle("dtheta");
0016   h2a->SetYTitle("dphi");
0017   h2a->SetZTitle("Energy");
0018   h2a->GetXaxis()->SetTitleOffset(2);
0019   h2a->GetYaxis()->SetTitleOffset(2);
0020   h2a->GetZaxis()->SetTitleOffset(1.2);
0021   h2->SetLineColor(ncol);
0022   h2a->SetLineColor(gcol);
0023   h2->Draw("surf");
0024   h2a->Draw("surf same");
0025   // TFile *fout = TFile::Open("finalhistos.root","UPDATE");
0026   c1->Print("Neutron_vs_Gamma.png");
0027   fin->Close();
0028   //fout->Close();
0029 }