From Graeme.Watt@cern.ch Sat Jan 14 00:44:09 2012 Date: Sat, 14 Jan 2012 06:43:58 +0100 From: Graeme Watt To: Xiaochao Zheng Subject: Re: request for MSTW structure function code Dear Xiaochao, I attach a file ("mstw2008structurefunctions.f") containing the Fortran code for F_2 and F_L used in the MSTW2008 fit. To use it, link to "mstwpdf.f" and "alphaS.f", and call the MSTWNC() subroutine in your main program. There is an example of usage given in the attached "example_sf.f". Unfortunately, only the photon contribution is currently implemented and a version extended to include the gamma-Z interference is not a high priority for me at the moment. Best regards, Graeme. On Jan 13, 2012, at 10:31 PM, Xiaochao Zheng wrote: > > Dear Dr. Watt: > > I am working on the data analysis of an experiment recently completed at Jefferson Lab. During the experiment we measured the parity violation asymmetry in the electron-deuteron DIS. In order to interpret our asymmetry, I need to calculate the expected asymmetry using world's best fit of structure functions such as R(gamma)=sigma_L/sigma_T, F1(gamma), F1(gamma-Z) and F3(gamma-Z), where the latter two are the interference terms. > > According to the MSTW2008 website, structure function code (neutral exchange, photon) is available upon request. May I have the code at the LO and the NLO? NNLO would be helpful too if available. > > Also, is it possible to extend these structure function codes to include the gamma-Z interference structure functions F1,3(gamma-Z)? These would be crucial in interpreting our data. > > (Fortran version would be the best). > > Best regards, > > Xiaochao > --------------------------------------------------------------- > Assistant Professor, Department of Physics, Univ. of Virginia > UVa office (during academic year): > 134 Physics Building, 382 McCormick Rd. > PO Box 400714, Charlottesville, VA 22904-4714 > Phone: 434-243-4032 > > JLab office (during summer months, May-August): > Rm 12/F357, MS 12H5, 12000 Jefferson Ave, Newport News, VA 23606 > Phone: 757-269-5433 > --------------------------------------------------------------- [ Part 2, Application/OCTET-STREAM (Name: ] [ "mstw2008structurefunctions.f") 139KB. ] [ Unable to print this part. ] [ Part 3, Application/OCTET-STREAM (Name: "example_sf.f") 8KB. ] [ Unable to print this part. ]