Subversion Repositories f9daq

Compare Revisions

Ignore whitespace Rev 84 → Rev 85

/lightguide/trunk/src/raySimulator.cpp
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());