UInt_t width,height;
void Sign()
{

  TImage *signal=TImage::Open("HRRLBeamSpot_Color.png");  
  UInt_t *pixaray=signal->GetArgbArray();
  width=signal->GetWidth();
  height=signal->GetHeight();

  Float_t Calibration;

  Calibration = 500/4; // 5 mm / 40 pixels

  TH2F *PictPixl=new TH2F("PictPixl","PictPixl",height,0,height,width,-0.5,width-0.5);
  TH1F *PictPixlX=new TH1F("PictPixlX","PictPixlX",height,0-0.5,height-0.5);
  TH1F *PictPixlY=new TH1F("PictPixlY","PictPixlY",width,-0.5,width-0.5);
  TH1F *PictPixlXcal=new TH1F("PictPixlXcal","PictPixlXcal",height,-0.5*Calibration,(height-0.5)*Calibration);
  TH1F *PictPixlYcal=new TH1F("PictPixlYcal","PictPixlYcal",width,-0.5*Calibration,(width-0.5)*Calibration);
  Int_t index;

  printf("Height=%d\tWidth=%d\n",height,width);

  //  for(Int_t j=0;j<height*2;j++)
  for(Int_t j=0;j<height;j++)
    {
      for(Int_t i=0;i<width;i++)
	{
	  index=i*width+j;
	  //	  Radius=sqrt((i-352)*(i-352) + (j-230)*(j-230));
	  //	  Radius=sqrt((j-352)*(j-352) + (i-230)*(i-230));
	  /*
	  if(pixaray[index]&0xffffff>0 && i > 600)
	    {
	      printf("%d\t%d\t%d\n",i,j,pixaray[index]&0xffffff);
	      printf("%d\t%d\t%d\n",i,j,pixaray[index]&0xffffff);
	      printf("%d\t%d\t%d\n",i,j,pixaraybckg[index]);
	      printf("Radius=%g\n",Radius);
	      printf("test1\n");
	    }
	  if((i>200 && i < 280)&& (j>420 && j<580))
	  if((i>200 && i < 340)&& (j>280 && j<420))
	  */
	    {
	      PictPixl->Fill(i,j,pixaray[index]&0xffffff);
	      PictPixlX->Fill(i,abs(pixaray[index]&0xffffff-24000000));
	      PictPixlY->Fill(j,abs((pixaray[index]&0xffffff)-21800000));
	      PictPixlXcal->Fill(i*Calibration,abs(pixaray[index]&0xffffff-24000000));
	      PictPixlYcal->Fill(j*Calibration,abs((pixaray[index]&0xffffff)-21800000));
	      /*
	      PictPixlY->Fill(j,abs((pixaray[index]&0xffffff)));

	      if(j==250)
		{
		  SigSumY250->Fill(i,pixaray[index]&0xffffff);
		  BackSumY250->Fill(i,pixaraybckg[index]&0xffffff);
		}

	      
	      PictPixlSub->Fill(i,j,(pixaray[index]&0xffffff)-1.0*(pixaraybckg[index]&0xffffff));
	      PictPixlSubX->Fill(i,(pixaray[index]&0xffffff)ffff));
	      PictPixlSubY->Fill(j,(pixaray[index]&0xffffff)-1.0*(pixaraybckg[index]&0xffffff));
	      */
	    }
	}
    }

  

  TCanvas *c1= new TCanvas("c1","c1", 10,10, 770,750);
  c1->Divide(2,2);
  c1->cd(1);
  signal->Draw();
  c1->cd(2);
  //  PictPixlX->Draw();;
  PictPixlXcal->Draw();;
  c1->cd(3);
  //  PictPixlY->Draw();
  PictPixlYcal->Draw();
}