StRoot: StEmcPool/StPhotonAnalysis/macros/binwidth/doBin.C Source File

00001 Double_t doBin(Double_t x,Double_t d)
00002 {
00003   Double_t alpha=9.125;
00004   //right
00005   x+=d/2.;
00006   cout<<"r: "<<x<<endl;
00007   Double_t retR=x*(1./(1.-alpha)) * pow(1.+x,1.-alpha);
00008   cout<<retR<<endl;
00009   retR=retR - (1./((1.-alpha)*(2.-alpha)))*pow(1.+x,1.-alpha); 
00010   cout<<retR<<endl;
00011   //left
00012   x-=d;
00013   cout<<"l: "<<x<<endl;
00014   Double_t retL=x*(1./(1.-alpha)) * pow(1.+x,2.-alpha);
00015   cout<<retL<<endl;
00016   retL=retL - (1./((1.-alpha)*(2.-alpha)))*pow(1.+x,2.-alpha); 
00017   cout<<retL<<endl;
00018   //ratio
00019   x+=d/2.;
00020   cout<<"c: "<<x<<endl;
00021   Double_t pqcd=pow(1.+x,-1.0*alpha);
00022   pqcd/=(retR-retL);
00023   return pqcd*d*x;
00024 }