UInt_t width,height,NP,N;
Int_t Entriesx,Entriesy; 
Double_t RMSx[100],RMSy[100], B[100];
NP=20;

void AutoRMS()
{
  Int_t j,i,n,m;
  for(n=0;n<=2*NP;n++)
  {
	 if(0<=n && n<=NP){m=n;}
	 if(NP<n){m=-n+NP;}
	 //B[n]=0.00355145 + 0.019446*n;
	 B[n]=1.0*m;
 	 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");}

 	 if(n==NP+1){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_1Amp.jpg");}
 	 if(n==NP+2){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_2Amp.jpg");}
 	 if(n==NP+3){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_3Amp.jpg");}
 	 if(n==NP+4){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_4Amp.jpg");}
 	 if(n==NP+5){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_5Amp.jpg");}
 	 if(n==NP+6){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_6Amp.jpg");}
 	 if(n==NP+7){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_7Amp.jpg");}
 	 if(n==NP+8){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_8Amp.jpg");}
 	 if(n==NP+9){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_9Amp.jpg");}
 	 if(n==NP+10){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_10Amp.jpg");}
 	 if(n==NP+11){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_11Amp.jpg");}
 	 if(n==NP+12){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_12Amp.jpg");}
 	 if(n==NP+13){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_13Amp.jpg");}
	 if(n==NP+14){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_14Amp.jpg");}
 	 if(n==NP+15){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_15Amp.jpg");}
 	 if(n==NP+16){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_16Amp.jpg");}
 	 if(n==NP+17){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_17Amp.jpg");}
 	 if(n==NP+18){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_18Amp.jpg");}
 	 if(n==NP+19){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_19Amp.jpg");}
 	 if(n==NP+20){TImage *signal=TImage::Open("Rotated_HRRL_Emit_test_Quad_Scan_Second_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");

	 Entriesx=px->GetEntries(); RMSx[n]=px->GetRMS(1);
	 Entriesy=py->GetEntries(); RMSy[n]=py->GetRMS(1);

 	 printf("%g\t%g\t%g\t%g\t\t%g\t%g\n",n,m,Entriesx,RMSx[n],Entriesy,RMSy[n]);
	 PictPixl->Delete();
 }

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

	TGraph *Hpx = new TGraph ( n, B, RMSx);
	Hpx ->SetTitle("X Projection RMS vs. Magnet Current");
	Hpx ->GetXaxis()->SetTitle("Current From Shunt Resistance (Amps)");
	Hpx ->GetYaxis()->SetTitle("RMS (Pixels)");
	Hpx ->GetYaxis()->CenterTitle();
	Hpx ->SetMarkerColor(4);
	Hpx ->SetMarkerStyle(21);
	Hpx ->SetMarkerSize(0.5);
	Hpx ->Draw("AP");

 	sp->cd(2);
	TGraph *Hpy = new TGraph ( n, B, RMSy);
	Hpy ->SetTitle("Y Projection RMS vs. Magnet Current");
	Hpy ->GetXaxis()->SetTitle("Current From Shunt Resistance (Amps)");
	Hpy ->GetYaxis()->SetTitle("RMS (Pixels)");
	Hpy ->GetYaxis()->CenterTitle();
	Hpy ->SetMarkerColor(4);
	Hpy ->SetMarkerStyle(21);
	Hpy ->SetMarkerSize(0.5);
	Hpy ->Draw("AP");

}