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')
|