void fdc_report(const char *fname)
{
	gROOT->Reset();
	
	// Coordinates of center of active area
	double x_center=56.0;
	double y_center=51.5;
	double r_active=48.5;
	
	// Get pointers to histos and tree from file
	TFile *f = new TFile(fname);
	TH2D *depth_mean = (TH2D*)gROOT->FindObject("depth_mean");
	TH2D *depth_rms = (TH2D*)gROOT->FindObject("depth_rms");
	TTree *scan = (TTree*)gROOT->FindObject("scan");

	// Define variables
	double mean_min=10000.0, mean_max=0.0, mean_mean=0.0;
	double rms_min=10000.0, rms_max=0.0, rms_mean=0.0;
	double rms_min_x, rms_min_y, rms_max_x, rms_max_y;
	double mean_min_x, mean_min_y, mean_max_x, mean_max_y;
	int N_inside_active = 0;
	

	// Loop over bins in mean and rms histos that are in active area
	// of cathode and calculate mean, min, and max
	int i,j;
	for(i=1; i<=depth_mean->GetNbinsX(); i++){
		double x = depth_mean->GetXaxis()->GetBinCenter(i);
		for(j=1; j<=depth_mean->GetNbinsY(); j++){
			double y = depth_mean->GetYaxis()->GetBinCenter(j);
			
			double r2 = pow(x-x_center, 2.0) + pow(y-y_center, 2.0);
			double r = sqrt(r2);
			
			if(r < (r_active-2.54/2.0)){
				double mean = depth_mean->GetBinContent(i,j);
				double rms = depth_rms->GetBinContent(i,j);
				
				if(mean==0.0)continue; // skip zeros
				
				mean_mean += mean;
				rms_mean += rms;
				N_inside_active++;

				if(rms<rms_min){
					rms_min=rms;
					rms_min_x = x;
					rms_min_y = y;
				}
				if(rms>rms_max){
					rms_max = rms;
					rms_max_x = x;
					rms_max_y = y;
				}
				if(mean<mean_min){
					mean_min = mean;
					mean_min_x = x;
					mean_min_y = y;
				}
				if(mean>mean_max){
					mean_max = mean;
					mean_max_x = x;
					mean_max_y = y;
				}
			}
		}
	}
	
	mean_mean/=(double)N_inside_active;
	rms_mean/=(double)N_inside_active;
	
	cout<<"======================================================="<<endl;
	cout<<" Report for \""<<fname<<"\""<<endl;
	cout<<"--------------------------------------"<<endl;
	cout<<"Active area radius used: "<<r_active<<" cm"<<endl;
	cout<<"Cathode center used: x="<<x_center<<" cm   y="<<y_center<<" cm"<<endl;
	cout<<endl;
	cout<<" Values below reported for active area only:"<<endl;
	cout<<endl;
	cout<<"mean depth: "<<mean_mean<<" mm"<<endl;
	cout<<"min. depth: "<<mean_min<<" mm     "<<"   location: x="<<mean_min_x-x_center<<" y="<<mean_min_y-y_center<<endl;
	cout<<"max. depth: "<<mean_max<<" mm     "<<"   location: x="<<mean_max_x-x_center<<" y="<<mean_max_y-y_center<<endl;
	cout<<"     range: "<<(mean_max-mean_min)*1000.0<<" microns"<<endl;
	cout<<"  mean RMS: "<<rms_mean<<" microns"<<endl;
	cout<<"  min. RMS: "<<rms_min <<" microns"<<"   location: x="<<rms_min_x-x_center<<" y="<<rms_min_y-y_center<<endl;
	cout<<"  max. RMS: "<<rms_max <<" microns"<<"   location: x="<<rms_max_x-x_center<<" y="<<rms_max_y-y_center<<endl;
	cout<<endl;
	cout<<"Report produced on: ";
	cout.flush();
	system("date");

	
}