the current t2d is chi's 9th order poly iterative fit

From: Chris Crawford (chris2@lns.mit.edu)
Date: Thu Nov 10 2005 - 00:49:29 EST


   A while back it was reported that the recrunch was using the 6th
order Garfield T2D, but actually it is Chi's. See the following lines
(don't know how the logic got so twisted!):

~blast/pro/Blast_Params/blastrc
# ---- Hits -----
*.Hit.T0: 6400
*.Hit.Conv: 0.5
*.Hit.Vel: 0.00182 # 2.02/1100.
*.Hit.MaxHits: 1000 # max hits / event
*.Hit.Cal: on # use t0 calibrations from file
*.Hit.Dist: on # use nonlinear tdc->distance function
*.Hit.SixthOrder: on # on = 6th order Garfield, off = 3rd
*.Hit.Hyperbolic: off # on = hyperbolic t2d
*.Hit.TdcMin: 2000 # 2000=8us (normal range=2us)
*.Hit.TdcMax: 6600 # full range (close to the wire)
*.Hit.DistMin: 0 # closest to wire
*.Hit.DistMax: 9999 # farthest from wire
*.Hit.Offset: 0.05 # stagger multiple hits
*.Hit.Color: 2
*.Hit.Width: 2

*.Hit.9th_cal: true

TBLRecon.cc:
   // Linking
   fWire = new TBLWc1WireCal();

   if (gOpt->GetValue("Hit.Dist",true)) { // nonlinear t2d
     if(gOpt->GetValue("Hit.9th_cal", true))
       ft2d = new t2d_cal_9th_Order(*fWire);
     else if(gOpt->GetValue("Hit.Hyperbolic",true))
       ft2d = new t2d_hyperbolic(*fWire);
     else if(gOpt->GetValue("Hit.SixthOrder",true))
       ft2d = new t2d_GarField_6th_Order(*fWire);
     else
       ft2d = new t2d_GarField_3rd_Order(*fWire);
   }
   else if(gOpt->GetValue("Hit.Cal", true)) { // improved linear t2d
     ft2d = new t2d_cubic(*fWire);
   }
   else { // simple linear t2d
     ft2d = new t2d_linear();
   }



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