[BLAST_ANAWARE] 6th-order time-to-distance relations added to reconstruction

From: Aaron Joseph Maschinot (ajmasch@MIT.EDU)
Date: Tue Jan 20 2004 - 14:33:48 EST


hi all:

the 6th-order time-to-distance relations have just been added to the
reconstruction. try them out, and let me know if you see any errors.

some points:

1) the 6th-order relations are now the default relations in the
   reconstruction. thus, if you CVS update your /BlastLib2 and
   /Blast_Params directories and do absolutely nothing else, you will be
   using the 6th-order relations.

2) in general, in order to use the 6th-order relations in the
   reconstruction, open the blastrc file and set the *.Hit.Dist entry to
   'on' and the *.Hit.SixthOrder entry to 'on':

    *.Hit.Dist: on # use nonlinear tdc->distance function
    *.Hit.SixthOrder: on # on = 6th order Garfield, off = 3rd

   if both flags are 'on', then the 6th-order time-to-distance relations
   (as stored within the blastrc file entry *.Time2Dis6thOrd.File) will
   be used. once again, this is currently the default.

3) in order to use the 3rd-order (i.e. cubic) relations (which worked
   just fine in the reconstruction but were not completely invertible and
   thus didn't work in the monte carlo), open the blastrc file and set the
   *.Hit.Dist entry to 'on' and the *.Hit.SixthOrder entry to 'off':

    *.Hit.Dist: on # use nonlinear tdc->distance function
    *.Hit.SixthOrder: off # on = 6th order Garfield, off = 3rd

   then the 3rd-order time-to-distance relations (as stored within the
   blastrc file entry *.Time2Dis3rdOrd.File) will be used. i kept this
   option in the reconstruction for backward compatibility purposes.

4) in order to use the linear relations, set the *.Hit.Dist entry to
   'off'; the *.Hit.SixthOrder entry will then not be checked:

    *.Hit.Dist: off # use nonlinear tdc->distance function

5) for the moment, the monte carlo can only invert the linear relations.
   within the next day or so, i will add the 6th-order relations, so the
   monte carlo will then be able invert either the 6th-order or the linear
   relations. THE CUBIC RELATIONS ARE NOT INVERTIBLE AND THUS WILL NEVER
   BE ENTERED INTO THE MONTE CARLO.

aaron



This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:30 EST