Yes, this is the wrong piece of code. You should look at the method Fitmessub() and not Fitmes() in fitUtil.icc Daniele > 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') > >