Hi, I have observed the same thing... When I saw the numbers the first time I thought the same ting... I searched through Djamels thesis and other places but were unable to find the number written elsewhere. It surely looks odd!!! Unless people using their own versions of method1 has done anything to these numbers they could well end up in the correction... Cheers, Truls On Wed, 2005-12-21 at 19:46 -0500, Flemming Videbaek wrote: > 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-lReceived on Thu Dec 22 05:41:44 2005
This archive was generated by hypermail 2.1.8 : Thu Dec 22 2005 - 05:41:57 EST