171,7 → 171,7 |
Init(); |
double izkoristek = Single(detector, parameters, TVector3(xh,yh,zh),theta, phi); |
PrintGuideHead(); |
PrintGuideStat(detector, 1); |
PrintGuideStat(izkoristek); |
DrawData(detector, parameters, theta, izkoristek); |
} |
//------------------------------------------------------------------------------------------ |
183,7 → 183,7 |
double izkoristek = Grid(detector, parameters, NN, theta); |
//printf("izkoristek = %.3lf\n", izkoristek); |
PrintGuideHead(); |
PrintGuideStat(detector, (NN+1)*(NN+1)); |
PrintGuideStat(izkoristek); |
DrawData(detector, parameters, theta, izkoristek); |
} |
//------------------------------------------------------------------------------------------ |
194,7 → 194,7 |
Init(); |
double izkoristek = RandYZ(detector, parameters, NN, theta, phi); |
//printf("izkoristek = %.3lf\n", izkoristek); |
PrintGuideHead(); PrintGuideStat(detector, NN); |
PrintGuideHead(); PrintGuideStat(izkoristek); |
DrawData(detector, parameters, theta, izkoristek); |
} |
//------------------------------------------------------------------------------------------ |
206,12 → 206,11 |
//CDetector detector = new CDetector(); |
double izkoristek = RandIso(detector, parameters, NN, theta, nnrays, showr); |
//printf("izkoristek = %.3lf\n", izkoristek); |
PrintGuideHead(); PrintGuideStat(detector, NN); |
PrintGuideHead(); PrintGuideStat(izkoristek); |
DrawData(detector, parameters, theta, izkoristek); |
//TCanvas *canvasDetector = new TCanvas("canvasDetector","canvasDetector",500,500); |
//canvasDetector->cd(); |
//detector->Draw(); |
//delete detector; |
} |
|
|
233,7 → 232,7 |
Init(); |
acc[i] = RandIso(detector, parameters, NN, theta); |
//printf("%6.2lf |%6.1lf\n", step[i], acc[i]*100.0); |
PrintGuideStat(detector, NN); |
PrintGuideStat(acc[i]); |
} |
|
//char sbuff[256]; |
263,8 → 262,8 |
CDetector *detector = new CDetector(TVector3(-2,0,0), parameters); |
acc[i] = RandYZ(detector, parameters, NN, step[i], phi); |
//printf("%6.1lf |%6.1lf\n", step[i], acc[i]*100.0); |
PrintGuideStat(detector, NN); |
delete detector; |
PrintGuideStat(acc[i]); |
//delete detector; |
} |
|
//char sbuff[256]; |
313,8 → 312,8 |
//double acceptance = -1.0; |
hAcceptance->Fill(d, a, acceptance); |
//printf("%.2lf | %.2lf | ", x,y); |
//PrintGuideStat(detector, NN); |
//delete detector; //works fine, 50x50 grid takes ~4MB of RAM |
//PrintGuideStat(izkoristek); |
delete detector; //works fine, 50x50 grid takes ~4MB of RAM |
} |
} |
|
330,7 → 329,7 |
TFile *file = new TFile("acceptance.root","RECREATE"); |
hAcceptance->Write(); |
file->Close(); |
delete file; |
//delete file; |
|
//hAcceptance->SetContour(100); |
//gStyle->SetPalette(1,0); |
370,7 → 369,7 |
printf("%.1lf | %.2lf | ", step[i], sipm_d); |
Init(); |
acc[i] = RandIso(detector, parameters, NN, theta); |
PrintGuideStat(detector, NN); |
PrintGuideStat(acc[i]); |
sprintf(sbuff, "d%2d.gif", i); |
//c3dview->SaveAs(sbuff); |
} |
429,7 → 428,7 |
printf("%.1lf | ", step[i]); |
Init(); |
acc[i] = RandIso(detector, parameters, NN, theta); |
PrintGuideStat(detector, NN); |
PrintGuideStat(acc[i]); |
delete detector; |
} |
|
461,7 → 460,7 |
CDetector *detector = new CDetector(CENTER, parameters); |
Init(); |
acc[i] = RandIso(detector, parameters, NN, theta); |
PrintGuideStat(detector, NN); |
PrintGuideStat(acc[i]); |
} |
|
//char sbuff[256]; |