UInt_t width,height;
void AutoRMS()
{
  Int_t j,i,n;
  for(n=0;n<=20;n++)
  {
  	 //TImage *signal=TImage::Open("temp.jpg");  

 	 if(n==0){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_0Amp.jpg");}
 	 if(n==1){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_1Amp.jpg");}
 	 if(n==2){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_2Amp.jpg");}
 	 if(n==3){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_3Amp.jpg");}
 	 if(n==4){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_4Amp.jpg");}
 	 if(n==5){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_5Amp.jpg");}
 	 if(n==6){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_6Amp.jpg");}
 	 if(n==7){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_7Amp.jpg");}
 	 if(n==8){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_8Amp.jpg");}
 	 if(n==9){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_9Amp.jpg");}
 	 if(n==10){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_10Amp.jpg");}
 	 if(n==11){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_11Amp.jpg");}
 	 if(n==12){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_12Amp.jpg");}
 	 if(n==13){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_13Amp.jpg");}
	 if(n==14){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_14Amp.jpg");}
 	 if(n==15){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_15Amp.jpg");}
 	 if(n==16){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_16Amp.jpg");}
 	 if(n==17){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_17Amp.jpg");}
 	 if(n==18){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_18Amp.jpg");}
 	 if(n==19){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_19Amp.jpg");}
 	 if(n==20){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_First_20Amp.jpg");}
 	 UInt_t *pixaray=signal->GetRgbaArray();

 	 width=signal->GetWidth();
 	 height=signal->GetHeight();
 	 //TH2F *PictPixl=new TH2F("PictPixl","PictPixl",width,0,width,height,0,height); //whole picture
	 //TH2F *PictPixl=new TH2F("PictPixl","PictPixl",height,0,height,width,0,width); //whole picture
 	 TH2F *PictPixl=new TH2F("PictPixl","PictPixl",200,260,460,250,150,400); //picture around spot

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

 	 for( j=0;j<=height;j++)
    		{
     		 for( i=0;i<=width;i++)
			{
	 		  index=720*480-j*width+i;

          		 if(290<i && i<440 && 255<j && j<350) //6

          		 //if(((i-355)*(i-355))+((j-280)*(j-280))<(6400)) //picture around spot

          		 if(((i-336)*(i-336))+((j-272)*(j-272))>(81)) //remove filament spot

	 	   	 PictPixl->Fill(i,j,pixaray[index]&0xffffff);
	}
    }

 	 TH1D *px = PictPixl->ProjectionX("x-projection");
 	 TH1D *py = PictPixl->ProjectionY("y-projection");

         Int_t Entriesx,Entriesy; Double_t RMSx,RMSy;
	 Entriesx=px->GetEntries(); RMSx=px->GetRMS(1);
	 Entriesy=py->GetEntries(); RMSy=py->GetRMS(1);

 	 printf("%g\t%g\t%g\t\t%g\t%g\n",n,Entriesx,RMSx,Entriesy,RMSy);
	 PictPixl->Delete();
/*
 	 TH1D *px = PictPixl->ProjectionX("x-projection");
 	 printf("X projection Entries=%g\n",px->GetEntries());
 	 printf("X projection RMS=%g\n",  px->GetRMS(1));

 	 TH1D *py = PictPixl->ProjectionY("y-projection");
 	 printf("Y projection Entries=%g\n",py->GetEntries());
 	 printf("Y projection RMS=%g\n",  py->GetRMS(1));
*/

/*
 	 TCanvas *sp= new TCanvas("sp","sp", 10,10, 770,750);//picture around spot
 	 sp->Divide(2,2);
 	 sp->cd(1);
 	 signal->Draw();

 	 sp->cd(2);
 	 PictPixl->Draw("colz");

 	 sp->cd(3);
 	 TH1D *px = PictPixl->ProjectionX("x-projection");
 	 px->Draw();
 	 px->Fit("gaus");
 	 printf("X projection Entries=%g\n",px->GetEntries());
 	 printf("X projection RMS=%g\n",  px->GetRMS(1));

 	 sp->cd(4);
 	 TH1D *py = PictPixl->ProjectionY("y-projection");
 	 py->Draw();
 	 py->Fit("gaus");

 	 printf("Y projection Entries=%g\n",py->GetEntries());
 	 printf("Y projection RMS=%g\n",  py->GetRMS(1));

*/
 }
}