00001 aperELEMENT(char *gmadFile, int x=1, int y=1)
00002 {
00003 gStyle->SetOptTitle(0);
00004 BDSCanvas *b = new BDSCanvas(gmadFile);
00005
00006 TH1F *aper_x = new TH1F("Apertures X","Apertures X",b->GetNumberElement(),b->GetEndPositionElement());
00007 TH1F *aper_y = new TH1F("Apertures Y","Apertures Y",b->GetNumberElement(),b->GetEndPositionElement());
00008
00009 for(Int_t i=1;i<b->GetNumberElement();i++)
00010 {
00011 aper_x->SetBinContent(i,b->GetElementXAper()[i]);
00012 aper_y->SetBinContent(i,b->GetElementYAper()[i]);
00013 }
00014 aper_x->SetStats(0);
00015 aper_x->SetLineColor(2);
00016 aper_x->SetLineWidth(1);
00017 aper_y->SetLineColor(4);
00018 aper_y->SetLineWidth(1);
00019 if(x) aper_x->Draw();
00020 aper_x->GetXaxis()->SetTitle("Distance from IP [m]");
00021 aper_x->GetYaxis()->SetTitle("Apertures [m]");
00022 if(x && y )aper_y->Draw("same");
00023 if(!x && y) aper_y->Draw();
00024 if (x && y)
00025 {
00026 leg = new TLegend(0.9,0.87,0.99,0.99);
00027 char title[20];
00028 leg->AddEntry(aper_x,"Horiz.","l");
00029 leg->AddEntry(aper_y,"Vert.","l");
00030 leg->Draw();
00031 }
00032 b->cd();
00033 }
00034
