[G8b_run] Photon Corrections (taggerSag)

Charles Hanretty hanretty at hadron.physics.fsu.edu
Thu Oct 29 14:14:36 EDT 2009


Hello all,
      Below is the C++ translation from the Fortran code that Mike wrote to 
correct the photon energies.

   float E0 = 4.55;
   double egdivE0;
   egdivE0 = pho_energy/E0 + 0.00134;

   if (egdivE0 < 0.401) {
     egdivE0 = egdivE0 - 0.00599597 + 0.0390912*egdivE0 -
0.0573218*egdivE0*egdivE0;
   }
   if ((egdivE0 >= 0.401) && (egdivE0 <= 0.674)) {
          egdivE0 = egdivE0 - 0.0181293 + 0.0675599*egdivE0 - 
0.0541268*egdivE0*egdivE0;
   }
   if (egdivE0 > 0.674) {
          egdivE0 = egdivE0 - 0.0511162 + 0.138208*egdivE0 - 
0.0867757*egdivE0*egdivE0;
   }

   pho_energy = egdivE0*E0 + pho_energy*(-0.0035);

Where "pho_energy" is the energy of the photon as taken from the data 
file and "E0" is the (e-)beam energy.

    A CLAS-note for the g8b Tagger Sag study is almost finished.

-Chuck


More information about the G8b_run mailing list