[Brahms-dev-l] absorption parameters

I am puzzled about the parameters in the fstree2datamap for absorption.
The code reads (see below) that absorption is smaller of FS than FFS.
Has this been used for a long time. I hope I read oit wrong or it has not been used.

In any case I am not able to reproduce these numbers (for p-bar) and get more absorption.
will post some plot for tomorrow.

flemming



//_____________________________________________________________________________
TF1* GetAbsorpFunc(Int_t pid, Char_t pol, Int_t angle)
{
  // set correction parameters
  TF1* absorpFunc = new TF1("absorpFunc", "[0]-[1]*exp(-[2]*x)", 0.0, 30.0);
  absorpFunc->SetParameters(1.0, 0.0, 0.0);

  if (pid != Selector::kProton) return absorpFunc;

  if (gPidMethod == "H1") // means only FFS
    {
      if (angle <= 4 && pol == 'B') // proton
        absorpFunc->SetParameters(0.9544, 0.1038, 0.193);
      if (angle <= 4 && pol == 'A') // anti-proton
        absorpFunc->SetParameters(0.9386, 0.1029, 0.1784);
      
      if (angle > 4 && pol == 'B') // proton
        absorpFunc->SetParameters(0.9824, 0.0217, 0.06053);
      if (angle > 4 && pol == 'A') // anti-proton
        absorpFunc->SetParameters(0.9704, 0.02623, 0.1277);
     absorpFunc->SetParameters(0.9704, 0.02623, 0.1277);
    }
  else
    {
      if (angle <= 4 && pol == 'B') // proton
        absorpFunc->SetParameters( 0.995, 0.2289, 0.1105);
      if (angle <= 4 && pol == 'A') // anti-proton
        absorpFunc->SetParameters(0.9637, 0.2710, 0.148 );
          
      if (angle > 4 && pol == 'B') // proton
        absorpFunc->SetParameters(0.9627, 0.05202, 0.1181);
      if (angle > 4 && pol == 'A') // anti-proton
        absorpFunc->SetParameters(0.9516, 0.06252, 0.1741);
    }

--------------------------------------------
Flemming Videbaek
Physics Department 
Bldg 510-D
Brookhaven National Laboratory
Upton, NY11973

phone: 631-344-4106
fax:        631-344-1334
e-mail: videbaek @ bnl.gov


_______________________________________________
Brahms-dev-l mailing list
Brahms-dev-l@lists.bnl.gov
http://lists.bnl.gov/mailman/listinfo/brahms-dev-l