[Frost] What do the intercepts mean in ELoss Correction?

[Sungkyun Park]

Hi Eugene,

Now I check the eloss code event by event and I think I find something.
Can you check it instead of me because FORTRAN is not a familiar job for me.

I think in the FORTRAN code, momcor.F (which is located in /u/home/supark/clas/packages/eloss)
the initialization of the variables is not worked. 

In the following, the Initial part is the place to check the initial values of the variables.
the targ (0,0,0.2666166) do not have the zero values. In the following, the same event part is between # and #.
The first event in the following has Vertex (in Butanol) and a1 (in Butanol). That is, targ(1) only has some value and the others should be zero.
But in the final process, targ has the value (0.3915622  , 0.000000  , 0.2666166 )
0.3915622 is from targets. 0.2666166 is the initial value.
In this event, eloss calculate energy loss in the butanol target and CH2 target.

And in the second event, the initial values of targ is the last values of targ in the first event.

I think this is the reason why we get the reduction of the number of events in CH2 target after applying the eloss correction.

Eugene, let me know this is reasonable or not.
I attach the output of 100 runs together.

Sung
Florida State University

  ###################################################
 Initial    : Vertex ( -0.4297619     -0.2481232       1.174964     )
 Initial    : le (   0.000000     ) pout ( -6.0856879E-02  0.2732520      8.9600444E-02 )
 Initial    : targ (   0.000000       0.000000      0.2666166     )
 Initial    : L (   0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       2.959804      3.0873390E-02  7.4647874E-02  0.8575045       0.000000      7.5778753E-02   1.448219       0.000000     )
 Initial    : st_offset (   2.000000     )
 Initial    : nwrite (           1 )

 Initial FROST : pout ( -6.0856879E-02  0.2732520      8.9600444E-02 )
 CALL FROST function
 CALL targets function
               : targ (   0.000000       0.000000      0.2666166     )
  CALL CyInt function
   vertex ( -0.4297619     -0.2481232       1.174964     )
   cdir   ( -0.2070411      0.9296305      0.3048296     )
   center (   0.000000       0.000000       0.000000     )
   dir    (   0.000000       0.000000       1.000000     )
   RC     (  0.7620000     )

  After CyInt  : a1( -0.5919008      0.4798931       1.413683     ) a2(   0.000000       0.000000       0.000000     ) ierr(           1 )
               : targ (   0.000000       0.000000      0.2666166     )
  Vertex -> A2 : in Butanol
  a1 -> A2 : in Butanol
 Variable L and targ

 For the targets function
  Vertex ( -0.4297619     -0.2481232       1.174964     )
  pout   ( -6.1271627E-02  0.2751143      9.0211079E-02 )
  targ   (  0.3915622       0.000000      0.2666166     )
 Butanol    : le  (  0.3915622     ) pout ( -6.1326720E-02  0.2753617      9.0292193E-02 )
 Butanol    : pin ( -6.1374489E-02  0.2755761      9.0362519E-02 )
 Carbon     : le  (   0.000000     ) pout ( -6.1374489E-02  0.2755761      9.0362519E-02 )
 Carbon     : pin ( -6.1374489E-02  0.2755761      9.0362519E-02 )
 CH2        : le  (  0.2666166     ) pout ( -6.1374489E-02  0.2755761      9.0362519E-02 )
 CH2        : pin ( -6.1561316E-02  0.2764150      9.0637587E-02 )
 ###################################################


 ###################################################
 Initial    : Vertex ( -0.4793870      0.2767742       6.049476     )
 Initial    : le (   0.000000     ) pout (  0.2033056      0.3828902      0.4422177     )
 Initial    : targ (  0.3915622       0.000000      0.2666166     )
 Initial    : L (  0.7831244       0.000000      8.3122686E-02  0.3085811       0.000000      7.3703863E-02   1.330781       0.000000      9.5213510E-02  1.0752519E-07   0.000000      5.9778359E-02  0.6805524       0.000000      6.0655493E-02   1.448219       0.000000     )
 Initial    : st_offset (   2.000000     )
 Initial    : nwrite (           1 )

 Initial FROST : pout (  0.2033056      0.3828902      0.4422177     )
 CALL FROST function
 CALL targets function
               : targ (  0.3915622       0.000000      0.2666166     )
  CALL CyInt function
   vertex ( -0.4793870      0.2767742       6.049476     )
   cdir   (  0.3282992      0.6182934      0.7140959     )
   center (   0.000000       0.000000       0.000000     )
   dir    (   0.000000       0.000000       1.000000     )
   RC     (  0.7620000     )

  After CyInt  : a1( -0.2428396      0.7222693       6.563999     ) a2(   0.000000       0.000000       0.000000     ) ierr(           1 )
               : targ (  0.3915622       0.000000      0.2666166     )
  Vertex -> A4 in Carbon
  a1 -> A4 in Carbon
   d   0.5191545
   d   0.7205238
   x   0.7140958
   targ   0.3915622      0.1049299      0.2666166
 Variable L and targ

 For the targets function
  Vertex ( -0.4793870      0.2767742       6.049476     )
  pout   (  0.2037915      0.3838053      0.4432746     )
  targ   (  0.3915622      0.1049299      0.2666166     )
 Butanol    : le  (  0.3915622     ) pout (  0.2037915      0.3838053      0.4432746     )
 Butanol    : pin (  0.2038594      0.3839331      0.4434223     )
 Carbon     : le  (  0.1049299     ) pout (  0.2038594      0.3839331      0.4434223     )
 Carbon     : pin (  0.2039871      0.3841736      0.4437000     )
 CH2        : le  (  0.2666166     ) pout (  0.2039871      0.3841736      0.4437000     )
 CH2        : pin (  0.2042536      0.3846756      0.4442798     )
 ###################################################

-------------- next part --------------
A non-text attachment was scrubbed...
Name: 100_Ruuns_Output
Type: application/octet-stream
Size: 169573 bytes
Desc: not available
Url : https://mailman.jlab.org/pipermail/frost/attachments/20100508/0aad5061/attachment-0001.obj