21-Aug-01 M. Yosoi rev. 31-Aug-04 M. Yosoi rev. 23-Aug-10 M. Yosoi rev. 08-Feb-15 M. Yosoi last rev. 07-Sep-20 M. Yosoi ************************************************************************ Brief manual for analyzing data in the UNIX workstation which are obtained by means of the TAMIDAQ system using the Grand Raiden FPP,LAS, etc. (Most of programs were originally written for the FACOM computer with LANA graphical packages to analyze data taken by the Q- or IQ-system. When the main machine was switched to the UNIX system, some parts were revised and extended to have wider applicability.) ************************************************************************ This analyzer can be used for both the online and offline almost in the same manner. The PAW in the CERN Library is employed for the histograming. ** new (Apr-2003)** Ntuple is newly supported, so that you can do the two-dimensional cut. ** new (Aug-2010)** Slight modification was made for the new computer system (HP, Redhat Enterprise Linux 5.4). ** new (Feb-2015)** Modification was made for the new DAQ system with V1190. I put modified programs in saho:/home/yosoi/Rcnpdaq/new-ana/src/ . (I changed "j11main4.f", "j11anna_mm.f", "vdc_m.f", and "inidef.f". You need to make the builder file as "LeCroy3377 compatible".) ** new (Sep-2020)** Slight modification was made for the new computer system (saho: CentOS8.2). See makefile and tcplib.c. ********* << Necessaries >> (the last general version is 7.0, in which some changes are made for the e200 experiment. All files are in the disk of [miho]->[saho].) source programs ; e.g., in the directory of ~yosoi/Rcnpdaq/e200/src. (The contents may be changed without notice.) You should copy a set of files (***.f, ***.c, ***.h and makefile). Compile and link is done by "make target" or "make j11vme#mm_ex". (# is one of 4, 3, 2, 1 and none. which depends on the blocksize of data. (4:128kB,3:80kB,2:64kB,1:48kB,nothing:16kB)) If your machine is a PC/Linux, some files must be exchanged by files in ~yosoi/Rcnpdaq/e200/src/linux. ** new (Aug-2010)** When running on the new Linux system (saho), copy tcplib.c.saho to tcplib.c. j11main4.f is also updated. See makefile. job control file ; ~yosoi/Rcnpdaq/e200/gr.in is an example. You need to change the content for every run. (As for the BLP data, see ~yosoi/Rcnpdaq/e156/Blp.) spectrum definition file ; There are some examples in the directory of ~yosoi/Rcnpdaq/e200/config or e156/config. Please write your own definition file referring to "grsingl" etc. other files ; ~yosoi/Rcnpdaq/e200/off_run (off_run0, tcp_run). (You need "chmod u+x" for above shell-script files.) macro files for PAW, i.e., hcopy.kumac, sum.kumac, etc. (They are in ~yosoi/macro or ~yosoi/Rcnpdaq/e200/macro. "hcopy.ray" is also convenient to make a postscript file of particle rays using gnuplot.) << Example of Job Execution >> % off_run0 ( : return ) Definition-File ? config/grsingl Scaler-File ? gr.sca log-File ? gr.log (if you want to stop,) ctrl-c (if you want to change it to background job,) ctrl-z % bg (N.B., in the [miho], batch jobs are recommended, then I make "off_run.bat" with "offrun.sh" and "run.cmd".) Output files after job is successfully finished ; gr.log log file of its job. gr.sca final scaler data. ***.hb HBOOK file for PAW. (*** is the name defined in the "gr.in") ***.kumac macro file for PAW. (exe ***#f: read HBOOK file, exe ***#1:plot spectrum ID=1 with x,y-axis titles. A macro file includes only 90 macros. ***a.kumac is for the ID more than 90, and ***b.kumac is for more than 180.) ***n.kumac macro file for PAW.(macro to plot a spectrum with its name. exe ***n#XPOS1 etc. The name list can be shown by exe ***#d. ***na.kumac, ***nb.kumac ... are similar to above.) ***-RAYXZ, ***.RAYXZ etc. files to plot particle rays by gnuplot. ( gnuplot> load "***-RAYXZ") fort.31 for the correction of zero timing of VDC. (Since the LRS3377 TDC became to be used, this file is of less importance.) << Example of Job Control File (gr.in) >> ---------------------------------------------------------------------------- gr r1005 ; TITLE GRS0 ; interactive(4capital)/no(blank) #repeat r05 ; OUTPUT FILE (HBOOK) 1 1 0 0 0 0 ; 1ST & LAST FILE / MIN. & MAX. BLOCK (>0)/dump/nobuild /home/yosoi/tmp1/e200/run1005.bld ; INPUT FILE NAME(s) ---------------------------------------------------------------------------- line 1: title (less than 50 characters) line 2: global file name (four capital letters) for PAW in the case of interactive run. Shared memory does not released after job-end. You should delete it. (check ID by 'ipcs -m' then ipcrm -m ID (or ipcrm shm ID for Linux)) If you write the number of Block (I4) successively, output of block scaler is the sampling one. (default is 1. See PRINTSCA in the spectrum definition file.) Blank is the case without using global file line 3: name attached prior to the some output files line 4: range of files and range of blocks to be analyzed (0 0 means from first to last) and dump option (off(0),on(1): when off, error and warning messages more than 100 are cut). In the online use, this line is -1 8 0 0 1 0. (** 2011Oct ** Note of range of block: MIN. =< n < MAX.) line 5: data file to be analyzed. When you want to analyze more than two files successively, other data files should be written from the next line. A blank line is needed in the end of files. << Content of Each Source File >> @ j11anna_mm.f ; main part of the analyzer. It consists of EVENT1 (dealing with each event), DEFINE (definition of spectra), SPTBL (table of spectrum name), etc. This file should be changed corresponding to each different experiment. (e.g., different module configuration, addition of non-standard spectrum ...) A simple correction for multi-cluster of wire chambers can be made by parameters of multig(l), etc. @ j11main4.f ; group of programs needed to decompose each event. (It usually doesn't require any modifications.) @ hout.f, gate1.f ; group of programs for histograming and gating. (It usually doesn't require any modifications.) @ inidef.f ; group of programs of man-machine interface utility between the spectrum definition file and analyzer. (It usually doesn't require any modifications.) @ champrm_ex.f ; initial data of chamber configurations etc. A few parts of this file should be changed for different experiment. Especially the time zero of each wire of VDC must be checked first. (see WIRE-TDC spectra.) You can choose 10 points for the z-position of focal plane (ZFV1, ZFV2). ZFV1(10) is important because the correction of focal plane (curvature etc.) can be done for the 10th plane. @ vdc_m.f ; group of programs to make spectra of hit pattern, position, cluster distribution etc. and to calculate efficiency of each VDC plane. You can choose how to calculate the position by the parameter LCAL1, LCAL2 in subroutine POSIT. @ rayvdc_ex.f ; group of programs for ray-tracing from all VDC's. @ mwpc.f ; group of programs to make position spectrum etc. of each MWPC of FPP. @ raypc_m.f ; group of programs for ray-tracing from more than four MWPC's. @ option_m.f ; group of programs to make non-standard spectra. For example, FMIST is used to make missing mass spectrum from the focal plane position spectra of both Grand Raiden and LAS in the (p,2p) experiments. @ mtread4.c, mtdump.c ; c-programs to read and dump raw data. (mtread3.c, etc. are substitutions for different block sizes.) @ tcplib.c ; libraries used in the online. << Example of Spectrum Definition File >> --------------------------------------------------------------------------- ! ! G-RAIDEN ! SAMPLE DEFINITION FILE ! ! ***** AVAILABLE OPERAND NAMES ***** ! ** STANDARD G-RAIDEN ** ! - 'NIM-IN ', 'NIM-IN-P ', 'HIT-NUMBER ', 'CLUST-SIZE ', ! - 'CLUST-NUM ', 'POSITION-ID ', 'POSITION ', 'SINGLE-ANGLE', ! - 'DIFF-MEAN ', 'DELTA-POS ', 'V-DELTA ', 'X-POSITION ', ! - 'Y-POSITION ', 'Z-POSITION ', 'GRAD-X ', 'GRAD-Y ', ! - 'GRAD-U ', 'GRAD-V ', 'GRAD-X-CORR ', 'GRAD-Y-CORR ', ! - 'ADC ', 'TDC ', 'TIME-DIFF-LR', 'ADC-LR-MEAN ', ! - 'RAY-ID ', 'VDC-MWPC ', 'WIRE-TDC-ALL', 'DRIFT-ALL ', ! - 'WIRE-TDC ', 'COIN-REG ', 'ADCP ', 'TDCP ', ! ** WITH FPP ** ! - 'X-POS-VDC ', 'Y-POS-VDC ', 'Z-POS-VDC ', 'RAY-ID-PC ', ! - 'HIT-NUM-PC ', 'CLUST-SIZ-PC', 'CLUST-NUM-PC', 'POS-ID-PC ', ! - 'POS-PC-1CL ', 'POS-PC-2CL ', 'POS-PC-3CL ', 'POS-PC-4CL ', ! - 'X-POS-PC ', 'Y-POS-PC ', 'Z-POS-PC ', 'GRAD-X-PC ', ! - 'GRAD-Y-PC ', 'GRAD-U-PC ', 'GRAD-V-PC ', 'GRAD-X-FPP ', ! - 'GRAD-Y-FPP ', 'PC-MWPC ', 'X-POS-DIFF ', 'Y-POS-DIFF ', ! - 'POS-PC-SNGL ', 'NIM-IN-ALL ', 'INT-REG ', 'RAY-ID-LAS ', ! ** WITH LAS & OTHERS ** ! - 'X-POS-LAS ', 'Y-POS-LAS ', 'Z-POS-LAS ', 'GRAD-X-LAS ', ! - 'GRAD-Y-LAS ', 'GRAD-U-LAS ', 'GRAD-V-LAS ', 'GRAD-XL-CORR', ! - 'GRAD-YL-CORR', 'NIM-BIT ', 'COIN-BIT ', 'COIN2-BIT ' ! - 'MIN-WTDC ', 'MIN+1-WTDC ', 'MIN-1-WTDC ', 'COIN2-REG ', ! FILTER1 RAY-ID ;00 0.0 50.0 GWIRE FILTER1 ADC-LR-MEAN ;01 100.0 250.0 GPI1P FILTER1 X-POSITION ;10 170.0 225.0 GELAS FILTER1 GRAD-X-CORR ;00 -0.6 0.6 GTHETA FILTER1 NIM-IN ;00 0.0 8.0 HATAOUT FILTER1 NIM-IN ;00 8.0 16.0 HATAIN ANDFILT GPI1P HATAIN $ HATA1 ANDFILT GPI1P HATAOUT $ HATA0 OR-FILT HATA1 HATA0 $ HATAOR OPTION1 PRINTCNF OPTION1 PRINTVDC OPTION1 PRINTSCA ! ! ------------ FOCAL PLANE SCINTILLATOR --------- ADDONE1 ADC ;01 0. 1024. 512 ADC01 'DE1-L ' 'CHANNNEL' 'COUNTS/CHANNEL' GATE-NAME GWIRE $ ADDONE1 TDC ;01 300. 500. 200 TDC01 'TDC1-L' 'CHANNNEL' 'COUNTS/CHANNEL' GATE-NAME GWIRE $ ADDONE1 TDC ;05 0. 1024. 512 TOF 'TOF ' 'CHANNEL' 'COUNTS/CHANNEL' GATE-NAME GWIRE GELAS $ ADDONE1 TDC ;93 0. 1024. 512 TOFC 'TOF-C ' 'CHANNEL' 'COUNTS/CHANNEL' GATE-NAME GWIRE GELAS $ ADDONE1 ADC-LR-MEAN ;01 0. 800. 400 PI01 'PARTICLE-ID-1' 'MEAN-LR' 'COUNTS' GATE-NAME GWIRE $ ! --------------- ID CHECK -------------- ADDONE1 NIM-IN ;00 -1. 31. 32 NIMIN 'NIM-IN ' 'CHANNEL' 'COUNTS/CHANNEL' ADDONE1 RAY-ID ;00 -10. 50. 60 RAYID 'RAY-ID ' 'CHANNEL' 'COUNTS/CHANNEL' ! ------------ FOCAL PLANE VDC ------------- ADDONE1 WIRE-TDC-ALL;01 0. 512. 512 WTDC1 'VDC1-TDC' 'CHANNEL' 'COUNTS/CHANNEL' ADDONE1 VDC-MWPC ;01 -1. 192. 192 MWPC1 'MWPC1 ' 'CHANNEL' 'COUNTS/CHANNEL' ADDONE1 POSITION ;01 -1. 1200. 1200 POS01 'POS-1 ' 'CHANNEL' 'COUNTS/CHANNEL' ADDONE1 X-POSITION ;01 -600. 600. 2400 XPOS1 'X-POSITION-1' 'X (MM)' 'COUNTS/0.5MM' GATE-NAME HATAOR $ ADDONE1 Y-POSITION ;01 -50. 50. 200 YPOS1 'Y-POSITION-1' 'Y (MM)' 'COUNTS/0.5MM' GATE-NAME HATAOR $ ADDONE1 GRAD-X-CORR ;00 -4. 4. 800 GRADXC 'GRAD-X-CORR' 'THETA-C (DEG)' 'COUNTS/0.01DEG' ADDONE1 X-POSITION ;10 -600. 600. 1200 XX10U 'XX-1OU' ' X-1O ' 'COUNTS/CHANNEL' SPIN-MODE UP GATE-NAME HATA0 GTHETA $ ADDONE1 X-POSITION ;10 -600. 600. 1200 XX10D 'XX-1OD' ' X-10 ' 'COUNTS/CHANNEL' SPIN-MODE DOWN GATE-NAME HATA0 GTHETA $ ADDONE2 X-POSITION ;10 GRAD-X-CORR ;00 0. 600. 120 -2.5 2.5 100 X10&AC 'X10 VS. THETA-C' ' XX-10 ' ' THETA-C ' ADDONE2 X-POSITION ;10 GRAD-X-CORR ;00 0. 600. 120 -2.5 2.5 100 X10&ACG 'X10 VS. THETA-C GATED' ' XX-10 ' ' THETA-C ' GATE-NAME GTHETA $ ADDONE2 TDC ;93 GRAD-X-CORR ;00 0. 2000. 200 -10. 10. 200 TOFC&AC 'TOF-C VS. THETA-C' ' TOF-C ' ' THETA-C ' GATE-NAME GELAS $ ADDRAY4 X-POSITION ;01 Z-POSITION ;01 X-POSITION ;04 Z-POSITION ;04 0. 600. -400. 600. 400 RAY&XZ 'RAY XZ-PLANE' ' X ' ' Z ' ADDRAY4 Y-POSITION ;01 Z-POSITION ;01 Y-POSITION ;04 Z-POSITION ;04 -25. 25. -400. 600. 200 RAY&YZ 'RAY YZ-PLANE' ' Y ' ' Z ' ! --------------- For Ntupl -------------- NTUPLE1 X-POSITION ;01 -600. 600. NXPOS10 'X-POSITION-1' 'X (MM)' 'COUNTS/0.5MM' GATE-NAME GPI1P GTHETA $ NTUPLE1 Y-POSITION ;01 -50. 50. NYPOS10 'Y-POSITION-1' 'Y (MM)' 'COUNTS/0.5MM' GATE-NAME GPI1P GTHETA $ NTUPLE1 GRAD-X-CORR ;00 -4. 4. NGRADXC 'GRAD-X-CORR' 'THETA-C (DEG)' 'COUNTS/0.01DEG' GATE-NAME GPI1P $ ENDFILE --------------------------------------------------------------------------- Commands which can be used in the definition file are only nine. FILTER1, FILTER2, OR-FILT, ANDFILT, ADDONE1, ADDONE2, ADDRAY4, OPTION1, NTUPLE1, ENDFILE Each command requires some operands except for ENDFILE. All commands must start from the first column of the line, while a series of operands may be written with free format and only the turn is meaningful. More than one blank or carriage return separates them. The index of the array operand is written after ";". "$" is the delimiter of a row of gates. The line whose first column is "!" is a comment. ENDFILE must be put at the last of significant part of the file. (N.B., FILTER2 command isn't now supported. Don't use TAB in the definition file.) FILTER1 OperandName;n Xmin Xmax GateName (1-dim gate) ( if Xmin < Xmax, Xmin=< X < Xmax is the gate region. if Xmin > Xmax, Xmax=< X < Xmin is the veto region.) OR-FILT GateName-1, GateName-2 .... $ GateName (OR of gates) ANDFILT GateName-1, GateName-2 .... $ GateName (AND of gates) ADDONE1 OperandName ;n Xmin Xmax ndiv MemberName (1-dim hist.) ### Flags ADDONE2 OperandName-x;n1 OperandName-y;n2 Xmin Xmax nXdiv Ymin Ymax nYdiv MemberName (2-dim hist.) ### Flags ADDRAY4 OperandName-x1;n1 OperandName-y1;n2 OperandName-x2;n1 OperandName-y2;n2 Xmin Xmax Ymin Ymax npoint MemberName (ray-plot) ### Flags OPTION1 PRINTCNF, PRINTVDC, or PRINTSCA (output control) (print configuration file, vdc information or scaler of each block. You'd better not use PRINTSCA when output unit is a log file.) NTUPLE1 OperandName;n Xmin Xmax $B!!!!(B MemberName (for Ntuple) ### Flags (ID of the Ntuple is 1000. How to plot a spectrum: for example, nt/plot 1000.nxpos10 (when the MemberName is NXPOS10) Set x=Xmin-1 for xXmax. There is no restriction for Xmin=Xmax. The event-loop of Ntuple is done for the event satisfying the common gate(s) in all NTUPLE1-definitions.) N.B.) The following eight operands can not be used in NTUPLE1, because more than one data are filled in the histgram for the same event. 'VDC-MWPC ', 'WIRE-TDC-ALL', 'DRIFT-ALL ', 'PC-MWPC ', 'CLUST-SIZ-PC', 'NIM-BIT ', 'COIN-BIT ', 'COIN2-BIT ' ENDFILE (This is indispensable.) Xmin/Xmax/Ymin/Ymax ; 4-byte real constatt/variable n/n1-n4/nXdiv/nYdiv/Npoint ; 4-byte integer constant/variable ****Name*** ; character (MemberName is less than 8 characters. Others are less than 12 characters.) ### = 'SubTitle' 'Xcomment' 'Ycomment' Flags ; SPIN-MODE UP (or DOWN) GATE-NAME GateName-1 GateName-2 ...... $ (### and Flags are optional.) Now 76 operands for the spectrum definition are given, which are shown in the first comment part of the above sample file. If the operand corresponds to array variable, the index is written after ';'. If not, the index is 0. @@@ Some limits in the spectrum definition file @@@ Maximum number of histograms defined in a file; 200 for 1-dimensional hist. 50 for 2-dimensional hist. 200 for Ntuple Maximum number of GateNames in a series ; 20 (before "$") Maximum number of gates defined by FILTER1 ; 200 In a series of gate names in the OR-FILT or ANDFILT, those defined by other OR-FILT or ANDFILT may be included. In that case, the parent OR-FILT (or ANDFILT) must only consist of gate names defined by FILTER1. @@@ Explanation of each operand @@@ Operand name Index Content ## NIMIN, ADC, TDC etc. ## NIM-IN-ALL 00 NIMIN (Input Register) NIM-IN 00 lower 4 bits of NIMIN except for BLP event NIM-IN-P 00 lower 4 bits of NIMIN for BLP event NIM-BIT 00 bit pattern of NIMIN INT-REG 00 Interupt Register (now using for 2nd-level) COIN-REG 00 Coincidence Register COIN-BIT 00 bit pattern of Coin. Reg. COIN2-REG 00 Coincidence Register-2 (range of current integrator) COIN2-BIT 00 bit pattern of Coin. Reg.-2 ADC 01-200 ADC {default 1:DE1L,2:DE1R,3:DE2L,4:DE2R} TDC 01-200 TDC {default 1:DE1L,2:DE1R,3:DE2L,4:DE2R,5:RF} ADCP 01-20 ADC for BLP {default 1:LF,2:RB,3:RF,4:LB,5:UF,6:DB,7:DF,8:UB} TDCP 01-20 TDC for BLP {default 1:LF,2:RB,3:RF,4:LB,5:UF,6:DB,7:DF,8:UB} TIME-DIFF-LR 01-50 time difference between left and right photomul signals ADC-LR-MEAN 01-50 geometrical mean of left and right photomul signals (PI) ## VDC ## HIT-NUMBER 01-20 number of hit wires in each chamber plane of VDC CLUST-SIZE 01-20 size of each cluster in each chamber plane of VDC CLUST-NUM 01-20 number of clusters per event in each chamber plane of VDC POSITION-ID 01-20 error ID of each chamber plane of VDC (see below) POSITION 01-20 position spectrum of each chamber plane of VDC SINGLE-ANGLE 01-20 angle spectrum obtained only by each plane of VDC DIFF-MEAN 01-20 (D(i+1)-D(i-1))/2 [D: drift length] in each plane of VDC DELTA-POS 01-20 (D(i+1)-D(i-1))/2-D(i) : a kind of resolution V-DELTA 01-20 (difference between V-plane position and ray tracing one) Now, difference between GRAD-X(U) and SINGLE-ANGLE except for 3,6,9,and 12. GR(1:X1,2:U1,4:x2,5:U2) LAS(7:X1,8:U1,10:X2,11:U2) VDC-MWPC 01-20 hit pattern spectrum in each chamber plane of VDC WIRE-TDC-ALL 01-20 TDC spectrum of all wires in each chamber plane of VDC DRIFT-ALL 01-20 drift length spectrum in each chamber plane of VDC WIRE-TDC 01-2000 TDC spectrum of each wire in each chamber plane of VDC RAY-ID{-LAS} 00 error ID of ray-tarcing (see below) ({ } is for LAS VDC) X-POSITION 01-10 | position spectra obtained by ray-tracing (focal plane Y-POSITION 01-10 | coordinate). 10 z-positions are given in ZFV1(ZFV2) Z-POSITION 01-10 | in champrm_ex.f. 10th is special for correction. {for LAS VDC, these are X -POS-LAS, Y-POS-LAS, Z-POS-LAS} GRAD-X{-LAS} 00 | angle spectra obtained by ray-tracing. GRAD-Y{-LAS} 00 | GRAD-U{-LAS} 00 | theta, phi, chi+(U), chi-(V) GRAD-V{-LAS} 00 | GRAD-X{L}-CORR 00 | incident angles to the spectrometer (theta-in, phi-in) GRAD-Y{L}-CORR 00 | calculated from theta, phi, and x-position. X-POS-VDC 01-10 | position spectra in the central-ray coordinate. Y-POS-VDC 01-10 | (These are used in the FPP experiment.) Z-POS-VDC 01-10 | 10 z-positions are given in ZFPC MIN-WTDC 01-20 TDC spectrum of minimum drift-time wire in each plane MIN+1-WTDC 01-20 TDC spectrum of minimum drift-time wire+1 in each plane MIN-1-WTDC 01-20 TDC spectrum of minimum drift-time wire-1 in each plane ## MWPC ## HIT-NUM-PC 01-20 number of hit wires in each chamber plane of MWPC CLUST-SIZ-PC 01-20 size of each cluster in each chamber plane of MWPC CLUST-NUM-PC 01-20 number of clusters per event in each plane of MWPC POS-ID-PC 01-20 error ID of each chamber plane of MWPC (see below) POS-PC-SNGL 01-20 position spectrum of each plane of MWPC (good event) POS-PC-1CL 01-20 position spectrum of each plane of MWPC (1st cluster) POS-PC-2CL 01-20 position spectrum of each plane of MWPC (2nd cluster) POS-PC-3CL 01-20 position spectrum of each plane of MWPC (3rd cluster) POS-PC-4CL 01-20 position spectrum of each plane of MWPC (4th cluster) PC-MWPC 01-20 hit pattern spectrum in each chamber plane of MWPC RAY-ID-PC 00 error ID of ray-tracing (see below) X-POS-PC 01-10 | position spectra obtained by ray-tracing (central- Y-POS-PC 01-10 | ray coordinate). 10 z-positions are given in ZFPC Z-POS-PC 01-10 | in champrm_ex.f. GRAD-X-PC 00 | angle spectra obtained by ray-tracing. GRAD-Y-PC 00 | GRAD-U-PC 00 | theta, phi, chi+(+45deg), chi-(-45deg) GRAD-V-PC 00 | GRAD-X-FPP 00 theta(PC) - theta(VDC) GRAD-Y-FPP 00 phi(PC) - phi(VDC) X-POS-DIFF 01-10 (X-POS-PC(i)) - (X-POS-VDC(i)) Y-POS-DIFF 01-10 (Y-POS-PC(i)) - (Y-POS-VDC(i)) N.B.-1) POSITION-ID 0; good event (1 cluster(+single hits)) 1; number of hits =< 1 2; edge of cluster has minimum drift time (If correction is done with LCAL1=1, then POSITION-ID=0 and change the sign of CLUST(I0(BSIZE) 4; multi single hits (--> CLUST-SIZE=1) 5; multi ( >3) clusters (1 cluster: >= 2 hits) 6; two clusters; both are bad (1 cluster: >= 2 hits) 7; two clusters; both are good (1 cluster: >= 2 hits) 8; Two Clusters; only one is good (1 cluster: >= 2 hits) 9; Three clusters (1 cluster: >= 2 hits) (POSITION-ID=5-9 --> CLUST-SIZE=-1) 10; number of hits >= MAXHIT (POSITION-ID=10 --> CLUST-SIZE=10) POS-ID-PC 0; good event (1 cluster (include corrected event)) 2; no hit 5; two-cluster event (can not be corrected)) (now event in which the space is one between two clusters is corrected. (POS-PC-2CL-->-1)) 7; multi cluster event ( > 2) 10; number of hits >= MAXHIT 15; bad chamber event RAY-ID 0; ray-tarcing is succeeded (RAY-ID-PC) i=1-10; POSITION-ID(POS-ID-PC) is not 0 in i-plane. (in the case of multi cluster correction (multi*=1), POSITION-ID=RAYID=0 and change the sign of CLUST-NUM) i=11-100; POSITION-ID(POS-ID-PC)'s of two planes are not 0. (e.g., 24 indicate i=2 and 4 are not good) i=-6 to -1; POSITION-ID(POS-ID-PC) of only one plane is 0. i=-7; 4 planes have good two clusters. (multi-ray) i=-8; POSITION-ID's of all 4 planes are not 0. i=-9; bad vdc(mwpc) event i=-10; no region ID N.B.-2) Spectra other than simple ADC, TDC and hit pattern of wire chamber have some parameters to be optimized for each different experiment. Please pay attention in using those spectra. Mon-standard spectra defined by user may be assigned to array of ADC or TDC which is not used as ADC and TDC. << How to Display Histograms using PAW >> in the same directory as analyzer has run, % paw (interactive case) PAW> glob GRS0 PAW> h/li ; show the list of spectra PAW> h/pl 1 ; display histogram ID=1 : : (in using HBOOK file after analyzer job finished) (when the 3rd line of gr.in is r02) PAW> exe r02#f : read HBOOK file (= h/file 1 r02.hb) PAW> exe r02#d : list of spectrum ID and MemberName PAW> h/pl 10 (or exe r01#10 etc.) (If you want to make postscript file of histogram(s)) PAW> exe hcopy#on PAW> h/pl 20 : : PAW> exe hcopy#off (A file named paw.ps is made. If you want to print out) PAW> sh lpr -Prcphaser1a paw.ps For a detailed description, please refer to the manual of PAW. << How to Display particle Rays >> % gnuplot gnuplot> load "r02-RAYXZ" $B!J(Br02-RAYYZ, etc. $B!K(B (if change the range) gnuplot> set yrange [-200:200] gnuplot> replot @ how to make postscript file gnuplot> load "hcopy.ray" gnuplot> replot (postscript file name is ray.ps) gnuplot> quit (exit from gnuplot) (END)