#include "TpcPadPlane.h"
#include "TpcGem.h"
#include "TpcPadShapePool.h"
#include <iostream>
#include <vector>
#include "TCanvas.h"
#include "TVector2.h"
#include "TWbox.h"
#include "TPad.h"



void plotPadPlane(){

  TpcGem* _gem=new TpcGem(5000,           // Gain
                  0.02);          // Spread
  
  TpcPadShapePool* _padShapes = new TpcPadShapePool("2mmPads.dat",
				   *_gem,
				   0.5, // lookup range
				   0.02, // Lookup Step
				   0.01); // LookupIntegrationStep
  
  
  TpcPadPlane* _padPlane= new TpcPadPlane(100,100,
			     1.,1.,
			     -50.,-50.);
  
  // Put Pads into PadPlane
  _padPlane->ReadFromFile("padplane.dat", _padShapes);
  
  
  std::cout<<*_padPlane<<std::endl;

  TCanvas* c=new TCanvas("c","regions",10,10,600,600);
  gPad->Range(-50,-50,50,50);

  // plot regions
  int npads=_padPlane->GetNPads();
  std::vector<TVector2> regmin(1000);
  std::vector<TVector2> regmax(1000);

  for(int k=0;k<1000;k++){
    regmin[k].Set(9999.,9999.);
    regmax[k].Set(-9999.,-9999.);
  }


  for(int i=0;i<npads;++i){
    TpcPad* apad=_padPlane->GetPad(i);
    int sectorid=apad->sectorId();
    regmin[sectorid].Set(min(apad->x(),regmin[sectorid].X()),
			 min(apad->y(),regmin[sectorid].Y()));
    regmax[sectorid].Set(max(apad->x(),regmax[sectorid].X()),
			 max(apad->y(),regmax[sectorid].Y()));
  }

  int counter=0;

  for(int j=0;j<1000;++j){
    
    if(regmin[j].Mod()>100 || regmax[j].Mod()>100) continue;
    std::cout<<"Sector["<<j<<"] x="<<regmin[j].X()<<".."<<regmax[j].X()
	     <<" y="<<regmin[j].Y()<<".."<<regmax[j].Y()<<std::endl;
    TWbox* b= new TWbox(regmin[j].X(),regmin[j].Y(),
		      regmax[j].X(),regmax[j].Y());
    b->SetBorderSize(0);
    b->SetFillColor(j%5);
    b->Draw();
    ++counter;
  }

  std::cout<<counter<<" Sectors found"<<std::endl;

  c->Update();
  c->Draw();

  c->SaveAs("regions.eps");
}