{

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


TCanvas* c1=new TCanvas("Riemann","Riemann",10,10,600,400);
TCanvas* c2=new TCanvas("Riemann2","Riemann",20,20,600,400);
TFile* infile=new TFile("RiemannFinder.root");

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

int i=0;
while(true){
  char name[200];
  sprintf(name,"plane%i",i);
  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==0)plane.DrawClone();
  else plane.DrawClone("same");
  ++i; 
}



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);

i=0;
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");

}