void parms2TeX_tavg(void)
{
	
	int Nangles = 4;
	int angles[] = {12.0, 16.0, 20.0, 90.0};
	
	int Nschemes = 6;
	const char* schemes[] = {
		"LAYERS_FINE",
		"LAYERS_1234_SECTORS1111",
		"LAYERS_22222_SECTORS11222",
		"LAYERS_244_SECTORS111",
		"LAYERS_3322_SECTORS1122",
		"LAYERS_334_SECTORS111"
	};
	
	const char* schemes_short[] = {
		"FINE",
		"1234",
		"22222",
		"244",
		"3322",
		"334"
	};
	
	// Loop over angles
	for(int iangle=0; iangle<Nangles; iangle++)
	//for(int iangle=0; iangle<1; iangle++)
	{
		//--------------- tavg ----------------
		double p0[6], p0err[6];
		double p1[6], p1err[6];
		double cov[6]; // off-diagonal covariance term
	
		// Loop over schemes
		for(int ischeme=0; ischeme<Nschemes; ischeme++){
			// Open file with tavg parameters for this angle/scheme
			char fname[256];
			sprintf(fname, "%ddegrees/%s/tavg_parms.dat", angles[iangle], schemes[ischeme]);
			ifstream ifs(fname);
			
			if(!ifs.is_open()){
				cout<<"Can't open \""<<fname<<"\"!"<<endl;
				continue;
			}
			
			string seg_scheme;
			int angle;
			
			ifs >> seg_scheme;
			ifs >> angle;
			ifs >> p0[ischeme] >> p0err[ischeme];
			ifs >> p1[ischeme] >> p1err[ischeme];
			ifs >> cov[ischeme];
			ifs.close();

		}
		
		cout<<endl;
		cout<<"% "<<angles[iangle]<<" degrees  tavg"<<endl;
		
		// Calculate 334 resolutions so we can find percentage change from that
		double E = 0.5;
		double norm500 = sqrt(pow(p0[Nschemes-1]/sqrt(E), 2.0) + pow(p1[Nschemes-1], 2.0));
		double E = 1.0;
		double norm1000 = sqrt(pow(p0[Nschemes-1]/sqrt(E), 2.0) + pow(p1[Nschemes-1], 2.0));

		// Loop over all segementation schemes
		for(int ischeme=0; ischeme<Nschemes; ischeme++){
			double p_0 = fabs(p0[ischeme]);
			double p_1 = fabs(p1[ischeme]);
			double p_0_err = fabs(p0err[ischeme]);
			double p_1_err = fabs(p1err[ischeme]);

			double E = 0.5;
			double r500= sqrt(p_0*p_0/E + p_1*p_1);
			double r500_err = p_0/E*p_0_err/r500 + p_1*p_1_err/r500;

			double E = 1.0;
			double r1000= sqrt(p_0*p_0/E + p_1*p_1);
			double r1000_err = p_0/E*p_0_err/r500 + p_1*p_1_err/r500;
			
			char str[256];
			sprintf(str, "%s\t& ($%d\\pm %d$)ps\t& ($%d\\pm %d$)ps\t& %3.1f\\%%\t& %3.1f\\%% \\\\"
				,schemes_short[ischeme]
				,(int)floor(0.5+r500)
				,(int)floor(0.5+r500_err)
				,(int)floor(0.5+r1000)
				,(int)floor(0.5+r1000_err)
				,100.0*(norm500 - r500)/norm500
				,100.0*(norm1000 - r1000)/norm1000);
			cout<<str<<endl;
			cout<<"\\hline"<<endl;
		}
	}

}