Print

Print


Hi all,
After a very quick look at the code in fitNtp.cc I see that the b0 mixing
correction is still done at mes level (maybe I'm wrong but I see that
piece of code in :
theFit() :

     for(int k=1;k<40;k++){
       tempbinchb = ((TH1D*)gDirectory->Get(namebch))->GetBinContent(k);
       tempbinb0os = ((TH1D*)gDirectory->Get(nameb0os))->GetBinContent(k);
       tempbinb0ss = ((TH1D*)gDirectory->Get(nameb0ss))->GetBinContent(k);
       temperrchb = ((TH1D*)gDirectory->Get(namebch))->GetBinError(k);
       temperrb0os = ((TH1D*)gDirectory->Get(nameb0os))->GetBinError(k);
       temperrb0ss = ((TH1D*)gDirectory->Get(nameb0ss))->GetBinError(k);
       tempbin = tempbinchb + ((1-chid)/(1-2*chid)) * tempbinb0os -
(chid/(1-2*chid)) * tempbinb0ss;
       if(MIXCORR==0){
	 tempbin = tempbinchb + tempbinb0os;
	 temperr = sqrt(temperrchb * temperrchb + temperrb0os *
temperrb0os);
       }else{
	 if(tempbin<0) tempbin=0;
	 temperr = sqrt(temperrchb*temperrchb + ((1-chid)/(1-2*chid)) *
((1-chid)/(1-2*chid)) * temperrb0os * temperrb0os + (chid/(1-2*chid)) *
(chid/(1-2*chid)) * temperrb0ss * temperrb0ss);
       }
       ((TH1D*)gDirectory->Get(name))->SetBinContent(k, tempbin);
       ((TH1D*)gDirectory->Get(name))->SetBinError(k, temperr);
    }

This is the same code as that i'm using right now to make comparisons...

Maybe I'm lloking at the wrong piece of code..
Can please Daniele point me to the right one?
Alessio

______________________________________________________
Alessio Sarti     Universita' & I.N.F.N. Ferrara
 tel  +39-0532-781928  Ferrara
roma  +39-06-49914338
SLAC +001-650-926-2972

"... e a un Dio 'fatti il culo' non credere mai..."
(F. De Andre')