July 1, 1988
* /NTUPLE/CUTS ICUT [ CHOPT FNAME ]
ICUT I 'Cut number' R=0:100
CHOPT C 'Options' D='P' R=' ,G,P,-,R,W,D'
FNAME C 'File name' D=' '
Define cut number ICUT for an ntuple. This cut can then be used
in subsequent commands NTUPLE/PLOT, PROJECT.
CHOPT='G' define a new cut ICUT using graphics input on the latest
1-DIM or 2-DIM projection of the Ntuple.
In case of a 1-DIM projection, give 2 points cutmin,cutmax.
In case of a 2-DIM projection, give up to 20 points to
delimit the selected area. The polygon will automatically
be closed by PAW.
CHOPT='P' Print definition of cut number ICUT.
'-' Reset cut ICUT
'R' read definition of cut ICUT from file FNAME.
'W' write definition of cut ICUT on file FNAME (text file).
'D' Draw cut contour.
Note that ICUT=0 means all cuts except for 'G' option.
If NUPD is negative then the histogram being filled has IDF=-NUPD. IDF may have been created with H/CREATE. Before filling IDF, the contents of IDF are reset if IDF already exists. Use NTUPLE/PROJECT to cumulate several passes into IDF. Note that NUPD negative is a convenient way to force user binning. This option must be used when options '+', 'U', 'S' are specified in CHOPT.
* NTUPLE/UWFUNC IDN FNAME [CHOPT]
IDN 'Ntuple identifier' I
FNAME 'File name' C
+
CHOPT 'Options' C D=' ' R=' ,E,P,T'
Generates automatically the Fortran skeleton of a selection function.
Example
If Ntuple ID=30 has variable names [X,Y,Z,ETOT,EMISS,etc] then
NTUPLE/UWFUNC 30 SELECT.FOR will generate the file SELECT.FOR with:
FUNCTION SELECT(XDUMMY)
COMMON/PAWIDN/IDNEVT,X,Y,Z,ETOT,EMISS,etc
*
SELECT=1.
END
Then using the command EDIT one can modify this file which could then
look something like: IDNEVT is the event number.
FUNCTION SELECT(XDUMMY)
COMMON/PAWIDN/IDNEVT,X,Y,Z,ETOT,EMISS,etc
*
IF(X**2+Y**2.GT.Z**2.OR.ETOT.GT.20.)THEN
SELECT=1.
ELSE
SELECT=0.
ENDIF
END
If CHOPT='E' then the local editor is invoked on FNAME.
='P' code to print events is generated.
='T' Tags for variables are generated in DATA.
Using the file VECTOR.DAT described in the documentation try:
VECTOR/CREATE X(3,10000).
VECTOR/READ X VECTOR.DAT.*
VECTOR/PLOT X(1)
VECTOR/PLOT X(1)%X(3)
![]() | Release Notes | Known bugs | FAQs | Contributions | Tutorial | Reference manual | Down load | Miscellaneous |