Subversion Repositories f9daq

Compare Revisions

Ignore whitespace Rev 46 → Rev 47

/sipmScan/examples/sipm.C
6,6 → 6,8
#include "TCanvas.h"
#include "TStyle.h"
#include "TPad.h"
#include "TPaveText.h"
//#include "TLabel.h"
#include "TF1.h"
#include "TGraph.h"
#include "TSpectrum.h"
15,11 → 17,9
 
double getNoise(TH2F*, int, int);
 
int sipm(char filename[256] = "test", char plopt[256]="all", int chXstart=0, int chXend=7, int chYstart=0, int chYend=7, bool debug = false)
int sipm(char filename[256] = "test", char plopt[256]="all", int parameter1=0, int parameter2=7, int chYstart=0, int chYend=7, bool debug = false)
{
const int c_nChannels = 64;
const double c_xOffset = 0; // mm
const double c_yOffset = 0;
int map[8][8]={{32,34,53,55,40,42,61,63},
{48,50,37,39,56,58,45,47},
43,6 → 43,16
 
// set draw style
RTSetStyle(gStyle);
// Print results to .pdf (option "p")
if(strstr(plopt, "p") != NULL) {
TCanvas* canvasFirst = new TCanvas("canvasFirst","",500,500);
TPaveText* text = new TPaveText(.05,.1,.95,.8);
text->AddText("SiPM Scan results");
canvasFirst->cd();
text->Draw();
canvasFirst->Print("result.pdf(");
}
if( strstr(plopt, "all") != NULL ) {
TCanvas *canvas2 = new TCanvas("canvas2","Hits x;;",2000,2000);
68,10 → 78,16
TCanvas *canvas10 = new TCanvas("canvas10","Ch x;;",500,500);
TH1F* h_hitsx;
canvas10->cd();
int chPosition = map[parameter1][parameter2];
char hname[128];
sprintf(hname, "hnhitsx%d", chXstart);
sprintf(hname, "hnhitsx%d", chPosition);
h_hitsx = (TH1F*)rootfile->Get(hname);
h_hitsx->Draw();
h_hitsx->Draw();
if(strstr(plopt, "p") != NULL) {
canvas10->Print("result.pdf");
}
}
if( strstr(plopt, "y") != NULL ) {
78,10 → 94,15
TCanvas *canvas11 = new TCanvas("canvas11","Ch x;;",500,500);
TH1F* h_hitsy;
canvas11->cd();
int chPosition = map[parameter1][parameter2];
char hname[128];
sprintf(hname, "hnhitsy%d", chXstart);
sprintf(hname, "hnhitsy%d", chPosition);
h_hitsy = (TH1F*)rootfile->Get(hname);
h_hitsy->Draw();
if(strstr(plopt, "p") != NULL) {
canvas11->Print("result.pdf");
}
}
if( strstr(plopt, "share") != NULL ) {
102,7 → 123,7
TCanvas *canvas4 = new TCanvas("canvas4","canvas4",500,500);
canvas4->cd();
for(int i=chXstart; i<=chXend; i++) {
for(int i=parameter1; i<=parameter2; i++) {
TH1F* h_hitsx;
char hname[128];
int chPosition = map[i][chYstart];
111,7 → 132,7
h_hitsx->SetTitle("Scan X;x [mm]; Entries");
h_hitsx->GetYaxis()->SetTitleOffset(1.3);
h_hitsx->SetStats(0);
if (i == chXstart)
if (i == parameter1)
h_hitsx->Draw();
else {
h_hitsx->SetLineColor(i+1);
126,7 → 147,7
for(int i=chYstart; i<=chYend; i++) {
TH1F* h_hitsy;
char hname[128];
int chPosition = map[chXstart][i];
int chPosition = map[parameter1][i];
sprintf(hname, "hnhitsy%d", chPosition);
h_hitsy = (TH1F*)rootfile->Get(hname);
h_hitsy->SetTitle("Scan Y;y [mm]; Entries");
153,12 → 174,12
// Get the filled histogram from scan
TH1F* h_x;
char name[32];
sprintf(name, "hnhitsx%d", chXstart);
sprintf(name, "hnhitsx%d", parameter1);
h_x = (TH1F*)rootfile->Get(name);
h_x->DrawCopy();
TH1F* h_y;
sprintf(name, "hnhitsy%d", chXstart);
sprintf(name, "hnhitsy%d", parameter1);
h_y = (TH1F*)rootfile->Get(name);
// Create and fill corrected histogram
229,7 → 250,7
h_line[j] = new TH1F("h_line", "h_line", binsX, xLowUser, xUpUser);
}
for(int j=chXstart; j<=chXend; j++) {
for(int j=parameter1; j<=parameter2; j++) {
int chPosition = map[j][chYstart];
char hname[128];
sprintf(hname, "h2d%d", chPosition);
245,10 → 266,10
double eta = -log(1 - signal);
double x = xLowUser + k*(xUpUser-xLowUser)/double(binsX);
//double y = l*(yUpUser-yLowUser)/double(binsY);
h_line[j]->Fill(x-c_xOffset, signal);
h_line[j]->Fill(x, signal);
//}
}
if (j == chXstart) {
if (j == parameter1) {
h_line[j]->SetTitle("SiPM#2 w/o noise subtraction;x[mm];Hits");
//h_line[j]->GetYaxis()->SetRangeUser(-0.05, 0.3);
//h_line[j]->GetYaxis()->SetRangeUser(-50, 2500);
265,9 → 286,9
TCanvas* canvas61 = new TCanvas("canvas61","canvas61",8*200,300);
canvas61->cd();
TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX,xLowUser,xUpUser, binsY,yLowUser,yUpUser);
for(int i=chXstart; i<=chXend; i++) {
//int canvasPosition = nX*(i-chYstart)+(j-chXstart)+1;
//int canvasPosition = nX*(chYend-i)+chXstart+1;
for(int i=parameter1; i<=parameter2; i++) {
//int canvasPosition = nX*(i-chYstart)+(j-parameter1)+1;
//int canvasPosition = nX*(chYend-i)+parameter1+1;
//if (debug) printf("canvas %d\n",canvasPosition);
int chPosition = map[i][chYstart];
char hname[128];
292,7 → 313,7
//double y = l*(yUpUser-yLowUser)/double(binsY);
double x = h[histogram]->GetXaxis()->GetBinCenter(k);
double y = h[histogram]->GetYaxis()->GetBinCenter(l);
h_corrected->Fill(x-c_xOffset, y-c_yOffset, eta);
h_corrected->Fill(x, y, eta);
}
}
}
315,10 → 336,10
TCanvas* canvas7 = new TCanvas("canvas7","canvas7", nX*400,nY*400);
printf("nx %d ny %d\n",nX,nY);
canvas7->Divide(nX,nY);
for(int i=chXend-1; i<=chXend+1; i++) {
for(int j=chXstart-1; j<=chXstart+1; j++) {
//int canvasPosition = nX*(i-chYstart)+(j-chXstart)+1;
int canvasPosition = nX*(chXend+1-i) + (j-chXstart+1) +1;
for(int i=parameter2-1; i<=parameter2+1; i++) {
for(int j=parameter1-1; j<=parameter1+1; j++) {
//int canvasPosition = nX*(i-chYstart)+(j-parameter1)+1;
int canvasPosition = nX*(parameter2+1-i) + (j-parameter1+1) +1;
if (debug) printf("canvas %d\n",canvasPosition);
canvas7->cd(canvasPosition);
char hname[128];
334,7 → 355,7
canvas8->cd();
gStyle->SetOptStat(0);
char hname[128];
int chPosition = map[chXstart][chXend];
int chPosition = map[parameter1][parameter2];
sprintf(hname, "h2d%d", chPosition);
TH2F* h_2d = (TH2F*)rootfile->Get(hname);
 
350,6 → 371,9
Double_t yUpUser = h_2d->GetYaxis()->GetXmax();
if (debug) printf("xLow %f xUp %f\n",xLowUser,xUpUser);
TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX, xLowUser, xUpUser, binsY, yLowUser, yUpUser);
double diffX = xUpUser - xLowUser;
double diffY = yUpUser - yLowUser;
//TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX, 3.00, 3.00+diffX, binsY, 8.00, 8.00+diffY);
//TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX, minX, maxX, binsY, minY, maxY);
double noise = getNoise(h_2d, 1, 70);
359,24 → 383,28
double signal = h_2d->GetBinContent(k,j); // detected
//double p = ((signal - noise) > 1) ? (signal-noise) : 1;
double p = signal - noise;
p /= 10000.;
/*
p /= 1000.;
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);
double y = h_2d->GetYaxis()->GetBinCenter(j);
if (debug) printf("x=%f y=%f\n", x, y);
h_corrected->Fill(x + c_xOffset, y + c_yOffset, eta);
h_corrected->Fill(x, y, p);
}
}
h_corrected->SetTitle("n_pe = - ln P(0);x[mm];y[mm]");
//h_corrected->SetTitle(";x[mm];y[mm]");
gStyle->SetPalette(52,0);
h_corrected->GetZaxis()->SetRangeUser(-0.05,0.7);
h_corrected->SetContour(50);
//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();
h_corrected->Draw("colz");
// collection efficiency
391,6 → 419,11
}
}
printf("Signal sum = %f\n # of points = %d\n",efficiency,nPoints);
if(strstr(plopt, "p") != NULL) {
canvas8->Print("result.pdf");
}
}
/** Draws the sum of channel signals
398,7 → 431,7
* Suitable for 8x8 chs scan
*/
if( strstr(plopt, "sum") != NULL ) {
int nX = chXend - chXstart + 1;
int nX = parameter2 - parameter1 + 1;
int nY = chYend - chYstart + 1;
TCanvas* canvas12 = new TCanvas("canvas12","c2",8*200, 8*200);
canvas12->cd();
422,11 → 455,11
// 2d histogram noise subtraction and poisson scaling
for(int i=chYstart; i<=chYend; i++) {
for(int j=chXstart; j<=chXend; j++) {
for(int j=parameter1; j<=parameter2; j++) {
int chPosition = map[j][i];
char hname[128];
sprintf(hname, "h2d%d", chPosition);
int histogram = nX*(i-chYstart)+(j-chXstart);
int histogram = nX*(i-chYstart)+(j-parameter1);
h[histogram] = (TH2F *) rootfile->Get(hname);
int noise = getNoise(h[histogram], 1, 170);
if (debug) printf("noise: %d\n",noise);
442,7 → 475,7
//printf("p=%f p0=%f log(p0)=%f eta=%f\n",p,p0,log(p0),eta);
double x = k*(xUpUser-xLowUser)/double(binsX);
double y = l*(yUpUser-yLowUser)/double(binsY);
h_corrected->Fill(x-c_xOffset, y-c_yOffset, eta);
h_corrected->Fill(x, y, eta);
}
}
469,12 → 502,12
TCanvas* canvas9 = new TCanvas("canvas9","canvas9", 500,500);
canvas9->cd();
char hname[128];
sprintf(hname, "hnhitsx%d", 36);
sprintf(hname, "hnhitsx%d", 49);
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])",17.02,17.30);
TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",parameter1,parameter2);
err->SetParameter(0,2500);
err->SetParameter(1, h_laser->GetMaximum());
err->SetParameter(2, h_laser->GetBinCenter(h_laser->GetMaximumBin()));
500,6 → 533,12
gStyle->SetOptStat(0);
h_map->Draw("text");
}
if(strstr(plopt, "p") != NULL) {
TCanvas* canvasLast = new TCanvas("canvasLast","",500,500);
canvasLast->Print("result.pdf)");
}
return(0);
}