Gverify test: geometry tag Y2006

Gverify : package designed to check a GEANT geometry for two particular types of errors: daughters not entirely contained in their mother volumes, and daughters which intersected. Is was written by Saul Youssef (1990). See also the following presentation given by James Bellinger (he has ported Gverify to Linux)

• GVerify talk for CDFII Tracking meeting (16 Mar. 2005) by James N. Bellinger (CDF)

So, what problems could STAR user expect from broken geometry? Many, actually, because detector material is not located where user expect it to be. Wrong efficiency estimations, noise checks, even subsystem acceptance could be very different...

Whoops... 1180 daugthers outside of mother volume + 529 daughters found to be intersecting...

Results/details could be found here (11 chunks, click y2006_0XX for details) :

Chunk	Daughter volumes protruded from Mother volume	Daughter volumes intersecting each other	Daughters close to int.
y2006_000	371	324	1
y2006_001	159	45	0
y2006_002	228	0	0
y2006_003	76	28	4
y2006_004	10	6	2
y2006_005	103	0	0
y2006_006	122	16	0
y2006_007	7	0	0
y2006_008	29	49	9
y2006_009	46	40	0
y2006_010	29	21	0
TOTAL: 1180 protrusions, 529 intersections, 18 close to intersection

Well, STAR geometry is very symmetric, so total number of unique cases should be at least one order of magnitude lower than table shows (my guess).

1. Run the scripts written by James. Those with f77 inside (not ifc), usually end with .com extension - to get dfunc, gvols, types and ming compiled
2. GVERIFY.com needs to be modified slightly to compile at RCAS:
   1. CERNSTUF='-L/cern/2004/lib/ -lgeant321 -lmathlib /cern/2004/lib/libpacklib_noshift.a /cern/2004/lib/libkernlib.a'
   2. FORTSTUF='/usr/lib/gcc-lib/i386-redhat-linux/3.2.3/libg2c.a -lF77 -lnsl'

Enjoy!

1. starsim -w 0, start starsim framework
2. gdrop all
3. detp geom Y200Xy, e.g.: Y2004c, Y2006 etc.
4. make geometry
5. gclose all
6. grfile y200Xy.rz
7. quit, quit starsim

1. Convert standard STAR .rz file to .fz file
2. Run Extractor: ./gvextrtd.exe(gvextrtx hacked to read .rz files), input file is not in exchange mode. Extractor creates "y2005.txt" file
3. Run Reducer: ./gvreducx.exe y2005.ext, y2005.red is produced
4. Run Splitter: ./gvsplitx.exe y2005.red, I've chosen chunks of 500
5. Run Gverify: ./gverifyx.exe y2005_008.red, y2005.rec is produced
6. Run Summarizer: ./gvsumryx.exe y2005_008.rec > y2005_008.sum, voila, all Errors instantly reported
7. OPTIONALLY run volume name converter: ./conv.php y2005_008.sum > y2005_008.conv.txt, to revert volume names back to normal

WARNING: it is surely NOT the best way to convert .rz to .fz formats ;)

What we do: read .rz file by standard GRFILE routine, dump it back by GFOUT as required by Gverify:

1. starsim -w 1
2. CALL GRFILE(50,'y2005.rz','I')
3. CALL GOPEN(60,'O',1024,0)
4. CALL GFOUT(60,'INIT',1,0,' ',0)
5. quit
6. mv fort.60 y2005.fz

Done!

Note: one can also use gvextrtd.exe which is hacked to use STAR standard .rz files...