| 90,7 → 90,7 |
|---|
| char sbuff[256]; |
| sprintf(sbuff, "SiPM = %.1lf, L.y. = %.2lf, d = %.1lf | gap = (%.1lf, %.1lf, %.1lf), #theta = %.1lf | acceptance = %.3lf", |
| parameters.getA(), parameters.getLightYield()*acc, parameters.getD(), |
| parameters.getGap().x(), parameters.getGap().y(), parameters.getGap().z(), theta, acc); |
| parameters.getGap().x(), parameters.getGap().y(), parameters.getGap().z(), theta, acc*100); |
| |
| if(!only2d) { |
| RTCanvas *cdata = new RTCanvas((char*)"Data", sbuff, 350, 0, 950, 950); |
| 98,12 → 98,11 |
|---|
| int cpc = 1; |
| cdata->cd(cpc++); |
| TH2F* generated = detector->GetGenerated(); |
| int nGenerated = generated->GetEntries(); |
| int minimum = generated->GetBinContent(generated->GetMinimumBin()); |
| double nGenerated = generated->GetEntries(); |
| //int minimum = generated->GetBinContent(generated->GetMinimumBin()); |
| int maximum = generated->GetBinContent(generated->GetMaximumBin()); |
| generated->GetZaxis()->SetRangeUser(0, 1.05*maximum); |
| double variation = (maximum-minimum)/(double)nGenerated; |
| printf("Statistical variation (max-min)/all = %f perc. \n", variation*100); |
| //double variation = (maximum-minimum)/(double)nGenerated; |
| generated->SetTitle("Generated"); |
| generated->Draw("colz"); |
| cdata->cd(cpc++); //((detector->GetLG())->GetHIn())->Draw("COLZ"); |
| 115,6 → 114,10 |
|---|
| |
| cdata->cd(cpc++); |
| TH2F* histoActive = (detector->GetHActive()); |
| //double nAccepted = histoActive->GetEntries(); |
| //double variation = 1 / sqrt(nAccepted); |
| double variation = sqrt(acc*(1-acc)/nGenerated); |
| printf("Statistical error = %f perc.\n", variation*100); |
| histoActive->Draw("COLZ"); |
| cdata->cd(cpc++); |
| TH2F* histoLaser = (detector->GetHLaser()); |