Fix for FMD DA
[u/mrichter/AliRoot.git] / ITS / AliITSLegoPlot.C
1 void AliITSLegoPlot(Float_t nchtheta, Float_t nchphi, const char *inFile = "galice.root") {
2 // macro to visualize the lego plots generated by gAlive->RunLego
3    
4    gROOT->Reset();
5    TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject(inFile);
6    if (!file->IsOpen()) {
7        cerr<<"Can't open "<<inFile<<" !" << endl;
8        return 1;
9    } // end if !file
10     
11
12    Float_t theta = 10;
13    Float_t phi   = 170;
14    Int_t ncont   = 50;
15
16    TCanvas *cradl = new TCanvas("cradl","radl",10,10,800,600);
17    cradl->SetTheta(theta);
18    cradl->SetPhi(phi);
19    TH2F *hradl = (TH2F*)file->Get("hradl");
20    hradl->SetStats(kFALSE);
21    hradl->GetXaxis()->SetTitle("Phi (degrees)");
22    hradl->GetYaxis()->SetTitle("Theta (degrees)");
23    hradl->SetFillColor(2);
24    hradl->SetContour(ncont);
25    hradl->Draw("colz");
26
27
28    TCanvas *cradlx = new TCanvas("cradl1","radl",50,50,800,600);
29    hradl->ProjectionX();
30    hradl_px->SetStats(kFALSE);
31 //   hradl_px->SetOptLogY();
32 //   hradl_px->SetMinimum(0.001);   
33 //   hradl_px->SetMaximum(1);
34    hradl_px->Scale(1./nchtheta);   
35    hradl_px->GetXaxis()->SetTitle("Phi (degrees)");
36    hradl_px->GetYaxis()->SetTitle("X/X0");
37    hradl_px->Draw();
38 //   cout << "Average over Phi: " << hradl_px->GetSumOfWeights()/nchphi << " X/X0" << endl;
39    
40    
41    TCanvas *cradly = new TCanvas("cradl2","radl",100,100,800,600);
42    hradl->ProjectionY();
43    hradl_py->SetStats(kFALSE);
44 //   hradl_py->SetOptLogY();
45 //   hradl_py->SetMinimum(0.001);   
46 //   hradl_py->SetMaximum(1);   
47    hradl_py->Scale(1./nchphi);
48    hradl_py->GetXaxis()->SetTitle("Theta (degrees)");
49    hradl_py->GetYaxis()->SetTitle("X/X0");
50    hradl_py->Draw();   
51 //   cout << "Average over Theta: " << hradl_py->GetSumOfWeights()/nchtheta << " X/X0" << endl;  
52    cout << "Average: " << hradl_py->GetSumOfWeights()/nchtheta << " X/X0" << endl;     
53 }