1 | README for SurfaceChecker |
---|
2 | |
---|
3 | |
---|
4 | |
---|
5 | History: |
---|
6 | ---------------------------------- |
---|
7 | |
---|
8 | 19.Dec 05 |
---|
9 | - Code is now using GetPointOnSurface() |
---|
10 | - Added new elliptic solid EllipticalCone |
---|
11 | - Added Tet |
---|
12 | |
---|
13 | 26.July 05 |
---|
14 | - Added support for segmented torus (fPhi,fPhiSegment) |
---|
15 | |
---|
16 | 25.July 05 |
---|
17 | - Sphere, Shell and other G4Sphere objects are now using the same |
---|
18 | code of SurfacePoint generation. |
---|
19 | - Added code for |
---|
20 | * small shells: Q2Shell Q3Shell Q4Shell Q5Shell Q6Shell Q7Shell Q8Shell |
---|
21 | * special regular spheres: HalfSphere HollowSphere HalfHollowSphere |
---|
22 | |
---|
23 | 4.July 05 |
---|
24 | - Added Ellipsoid to the check |
---|
25 | |
---|
26 | 1.July 05 |
---|
27 | - SurfaceChecker supports now G4Sphere("aSphere",rmin,rmax,phi,dphi,theta,dtheta) without |
---|
28 | restricting the parameters. |
---|
29 | - corrected error in the surface equation of the twisted tubs solid. (The edge of the hyperbolic |
---|
30 | surface was not compatible with the implementation in G4TwistedTubs). |
---|
31 | |
---|
32 | 31.May 05 |
---|
33 | - change in singlesolid.kumac: |
---|
34 | added an argument in order to specifiy the name of |
---|
35 | the hbook file directly, e.g |
---|
36 | PAW>> exec singlesolid hbk/yourfile.hbk |
---|
37 | - corrected a misleading information in extract.pl |
---|
38 | |
---|
39 | 27.May 05 |
---|
40 | - change in testall.pl: additionally executes at the end the paw macro single.kumac |
---|
41 | the default number of events is set to 100000 |
---|
42 | - added a paw macro single.kumac: |
---|
43 | analyses systematically all hbook files found in folder hbk ( hbk/run_cone_10000.hbk ) |
---|
44 | and produces two plots |
---|
45 | - solid-full.ps : comparing all solids in a full range (1e-12 to 1e5 mm) |
---|
46 | - solid-zoom.ps : comparing all solids in a reduced range (1e-12 to 1e-8 mm) |
---|
47 | Attention: only files respecting the naming rule run_solidname_nevents.hbk will be considered! |
---|
48 | - added a paw macro singlesolid.kumac: |
---|
49 | Can be used to to the same analysis for one single file only. No naming rule for the bhook file. |
---|
50 | |
---|
51 | 23.May 05 |
---|
52 | - changed naming of general surface parameters (hbk only): now u,v insted of u,phi |
---|
53 | - Detector is directly selectable from macro file for the following solids: |
---|
54 | Torus (360 deg only, test of outer surface) |
---|
55 | Box (one side only) |
---|
56 | Sphere (solid sphere, 360 deg) |
---|
57 | Tube (solid tub, 360 deg) |
---|
58 | Orb |
---|
59 | Cone |
---|
60 | TwistedTubs (test of hyperbolic surface ; all parameters free) |
---|
61 | TwistedBox (test of 0deg side) |
---|
62 | TwistedTrd (test of 0deg side) |
---|
63 | TwistedTrap (regular case: equal endcaps, alpha=0, phi=0, theta=0 ; test of 0deg side) |
---|
64 | TwistedTrap2 (equal endcaps ; test of 0deg side ) |
---|
65 | TwistedTrap3 (general case ; test of 0deg side ) |
---|
66 | - Added testall.pl: a script to systematically test a list of solids with SurfaceChecker. |
---|
67 | Produces hbook files in folder hbk, eg. for 500 events and G4Box file will be named hbk/run_box_500.hbk |
---|
68 | |
---|
69 | Author: |
---|
70 | Oliver Link, Oliver.Link@cern.ch |
---|
71 | Thu May 19 14:50:39 CEST 2005 |
---|
72 | |
---|
73 | |
---|
74 | SurfaceChecker |
---|
75 | ==================================== |
---|
76 | |
---|
77 | The following two steps have to be executed in order: |
---|
78 | |
---|
79 | 1) SurfaceChecker run1.mac > data/yourfile.data |
---|
80 | |
---|
81 | Comments: |
---|
82 | - in run1.mac: select your detector (Box, Torus, Cons, Sphere, Orb, TwistedTrap..) |
---|
83 | - run1.mac can be replaced by some longer tests,e.g. runLong.mac |
---|
84 | - runspecial.mac is used to save the random number at the end of |
---|
85 | a given event. Restart with runSingle.mac at the previous random |
---|
86 | number position. |
---|
87 | - You have to redirect the output to a file in the directory |
---|
88 | "data", with the filename extension ".data". |
---|
89 | |
---|
90 | 2) extract.pl |
---|
91 | |
---|
92 | This perl script parses the output file (located in the "data" |
---|
93 | directory) and converts the information into a ntuple. The file will |
---|
94 | have the same name as your data file (e.g. yourfile.hbk) and will |
---|
95 | be stored in the directory "hbk". |
---|
96 | The format of the ntuple (with Id 1) is the following: |
---|
97 | |
---|
98 | u,v : general surface parameters |
---|
99 | x,y,z: generated (true) surface point in cartesian coordinates |
---|
100 | px,py,pz: vertex position in cartesian coordinates |
---|
101 | dist: (true) distance between vertex and intersection |
---|
102 | vx,vy,vz: momentum of particle at vertex |
---|
103 | rx,ry,rz: reconstructed intersection |
---|
104 | delta : distance between reconstructed and true intersection |
---|
105 | theta: angle between momentum and surface normal at intersection |
---|
106 | |
---|