{

//gROOT->ProcessLine(".x tpc/TestChamber/macro/christian_style.C");


TCanvas* c1=new TCanvas("Riemann","Riemann",10,10,600,400);
TCanvas* c2=new TCanvas("Riemann2","Riemann2",20,20,600,400);
TFile* infile=new TFile("RiemannInit.root");
if(infile==NULL)return;
infile->Print();

c1->cd();
TF3* plane=new TF3("plane","x*[0]+y*[1]+z*[2]+[3]",-1,1,-1,1,0,1.2);
plane->SetNpz(20);
plane->SetTitle("Planes");

int i=1;
while(true){
  char name[200];
  sprintf(name,"plane%i",i);
  std::cout<<"Plane "<<i<<std::endl;
  TLorentzVector* par=(TLorentzVector*)infile->Get(name);
  if(par==NULL)break;
  plane->SetParameters(par->X(),par->Y(),par->Z(),par->T());
  plane->SetLineColor(kMagenta+i+100);
  if(i==1)plane->DrawClone();
  else plane->DrawClone("same");
  ++i; 
}

std::cout<<"Drawn planes"<<std::endl;

return;
}

//TF3 sphere("sphere","sqrt(x*x+y*y+(z-0.5)*(z-0.5))-0.499",-1,1,-1,1,0,1.2);
//sphere.SetFillStyle(0);
//sphere.SetLineColor();
//sphere.SetNpz(30);
//sphere.SetTitle("Riemann Sphere");
//sphere.SetLineColor(17);
//sphere.Draw("same");

//c2->cd();
//sphere.DrawClone();

//TCanvas* c=new TCanvas("Tracks","Tracks",10,10,600,600);

//return;

i=1;
while(true){
  char name[200];
  sprintf(name,"track%i",i);
  TPolyMarker3D* marker=(TPolyMarker3D*)infile->Get(name);
  if(marker==NULL){std::cout<<"Marker not found"<<std::endl;break;}
  sprintf(name,"clusters%i",i);
  TPolyMarker3D* clust=(TPolyMarker3D*)infile->Get(name);
  if(clust==NULL){std::cout<<"Clust not found"<<std::endl;break;}
  
  marker->SetMarkerColor(kMagenta+i);marker->SetMarkerStyle(20);
  marker->SetMarkerSize(0.75);
  clust->SetMarkerColor(1);clust->SetMarkerStyle(5);clust->SetMarkerSize(0.5);
  
  //if(i==0)marker->Draw();
  c1->cd();
  marker->Draw("same");
  clust->Draw("same");

  //c2->cd();
  //marker->DrawClone("same");
  //clust->DrawClone("same");

  //c->cd();
  //clust->DrawClone("AP");

  ++i;
}
//delete c;
//c1->SaveAs("riemann2.eps");
//c2->SaveAs("riemann1.eps");

}