50,16 → 50,17 |
canvas9->cd(); |
char hname[128]; |
int channel = map[chX][chY]; |
sprintf(hname, "hnhitsx%d", channel); |
sprintf(hname, "hnhitsy%d", channel); |
TH1F* h_laser = (TH1F*)rootfile->Get(hname); |
h_laser->Draw(); |
h_laser->SetStats(1); |
|
TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",parameter1,parameter2); |
err->SetParameter(0,2500); |
err->SetParameter(1, h_laser->GetMaximum()); |
//TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",parameter1,parameter2); |
TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",h_laser->GetBinCenter(3), h_laser->GetBinCenter(h_laser->GetMaximumBin())); |
err->SetParameter(0, h_laser->GetMinimum()); |
err->SetParameter(1, h_laser->GetMaximum() - h_laser->GetMinimum()); |
err->SetParameter(2, h_laser->GetBinCenter(h_laser->GetMaximumBin())); |
err->SetParameter(3, 0.001); |
err->SetParameter(3, 0.025); |
h_laser->Fit(err,"qr"); |
h_laser->Fit(err,"lr"); |
double sigma = err->GetParameter(3); |