This code enables you to choose one of 100 different binned jet-correction functions, obtained from the "standard" binned function by fluctuating each bin according to a gaussian as wide as the corresponding statistical error. The jet correction file and the corresponding errors can now be chosen (they should be the same as used to correct the MC in your default analysis)

Files you need:

• fbin_r.f function
• random.file: set of 100 arrays of random numbers
• iseq.inc, an 'include' file containing a DATA statement, allowing the user to specify the selected sequence number ISEQ (1 to 100)

  ISEQ = 0 corresponds to the central values.

  ISEQ = 10 and ISEQ = 81 correspond to the recommended values for a reasonable estimate of the statistical error from jet corrections (choose one of the two or take the maximum observed variation)

  (I know the include file containing just one DATA statement is not elegant, it was done for practical purposes. You can modify the code yourself if you wish, giving the ISEQ number to the function the way you like best, either as a second argument, or as a parameter, etc..)