| 347,6 → 347,7 |
|---|
| sprintf(hname, "h2d%d", chPosition); |
| TH2F* h_2d = (TH2F*)rootfile->Get(hname); |
| h_2d->Draw("colz"); |
| if ( !((i == parameter2) and (j == parameter1)) ) h_2d->GetZaxis()->SetRangeUser(0,200); |
| } //x |
| } |
| |
| 377,6 → 378,7 |
|---|
| //TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX, minX, maxX, binsY, minY, maxY); |
| |
| double noise = getNoise(h_2d, 1, 70); |
| noise = 50; |
| if(debug) printf("Noise = %f\n", noise); |
| for(int k=minX; k<=maxX; k++) { |
| for(int j=minY; j<=maxY; j++) { |
| 383,13 → 385,11 |
|---|
| double signal = h_2d->GetBinContent(k,j); // detected |
| //double p = ((signal - noise) > 1) ? (signal-noise) : 1; |
| double p = signal - noise; |
| /* |
| p /= 1000.; |
| p /= 10000.; |
| double p0 = 1.0 - p; // events with zero photons |
| //double eta = (-log(p0) * p0 ) / (1-p0-0.00001); |
| double eta = -log(p0); // constant of the poissonian statistics |
| if (debug) printf("p=%f p0=%f log(p0)=%f eta=%f\n",p,p0,log(p0),eta); |
| */ |
| //double x = xLowUser + k*(xUpUser - xLowUser) / double(binsX); |
| double x = h_2d->GetXaxis()->GetBinCenter(k); |
| //double y = yLowUser + j*(yUpUser-yLowUser)/double(binsY); |
| 401,10 → 401,10 |
|---|
| //h_corrected->SetTitle("n_pe = - ln P(0);x[mm];y[mm]"); |
| h_corrected->SetTitle(";x[mm];y[mm]"); |
| //gStyle->SetPalette(52,0); // black and white for print |
| //h_corrected->GetZaxis()->SetRangeUser(-0.05,1); |
| //h_corrected->SetContour(100); |
| //gStyle->SetPalette(1); |
| //SetGS(); |
| //gStyle->SetPalette(1); //black and white 2nd option |
| //SetGS(); // inverse |
| h_corrected->GetZaxis()->SetRangeUser(-0.05,0.30); |
| h_corrected->SetContour(30); |
| h_corrected->Draw("colz"); |
| |
| // collection efficiency |