c--------------------------------------------------------------------------- c Copyright (c) 1996 Southeastern Universities Research Association, c Thomas Jefferson National Accelerator Facility c c This software was developed under a United States Government c license described in the NOTICE file included as part of this c distribution. c c Eddy A.J.M. Offermann, 12000 Jefferson Ave., Newport News, VA 23606 c Email: offerman@cebaf.gov Tel: (757) 249-7544 Fax: (757) 249-5800 c--------------------------------------------------------------------------- c include 'coolhands/parameter.h' integer MAXRUNS ! maximum number of runs that can be analyzed ! at once. parameter(MAXRUNS=100) c---spectra parameter c NSGOOD number of spectra defined in ESPACE integer NSGOOD parameter(NSGOOD=300) c---event parameter integer MAX_FLAG ! Maximum size of event identification array parameter (MAX_FLAG=4) c---transport parameter integer NMATR_FOC_TAR1, ! # matrix elements for the focal plane > NMATR_FOC_TAR2, ! # matrix elements converting focal-plane ! coordinates to target > NMATR_TAR_FOC, ! # matrix elements converting target ! coordinates to focal plane > NMATR_PATH, ! # matrix elements for the path length diff > PORDER, ! order of polynomial describing the focal- ! plane dependence of the matrix elements > POWSM_FOC_TAR, ! maximum value of j+k+l in ! M_jkl*thfp**j*yfp**k*phfp**l > POWSM_TAR_FOC, ! maximum value of i+j+k+l+m in ! M_ijklm*xtg**i*thtg**j*ytg**k*phtg**l*dp/p**m > POWSM_PATH ! maximum value of i+j+k+l in ! M_ijkl*xfp**i*thfp**j*yfp**k*phfp**l parameter(NMATR_FOC_TAR1=3,NMATR_FOC_TAR2=200,NMATR_TAR_FOC=100, > NMATR_PATH=100,PORDER=7,POWSM_FOC_TAR=10, > POWSM_TAR_FOC=1,POWSM_PATH=4) c c R. Suleiman 3/25/1998 c Changed NMATR_PATH=33 to NMATR_PATH=100 to calculate the path c correction to higher orders. c c---detector parameter c NPLANE_VDC number of vdc planes per spectrometer c NMAX_WC maximum number of wires in a plane c NHITS_WC maximum number of wires excited by one particle in c a plane c NHITS_WC_MIN c minimum number of wires excited by one particle in c a plane c NGAP_CLUS The minimum amount of unexcited wires between two c clusters c NCLUS_WC maximum number of clusters in one plane c NHITS_CLUS maximum number of wires excited per cluster c NHITS_CLUS_MIN c minimum number of wires excited by one particle in c a cluster to be valid c NMAX_S0 max. number of paddles in scintillator, S0 c NMAX_S1 max. number of paddles in scintillator, S1 c NMAX_S2 max. number of paddles in scintillator, S2 c NMAX_S3 max. number of paddles in scintillator, S3 c NHITS_PD maximum # hits per paddle c NMAX_AERO number of mirrors in AEROGEL cherenkov c NMAX_GAS number of mirrors in GAS cherenkov c NHITS_MIR maximum # hits per Cherenkov mirror c NMAX_PRESHOW c number of blocks in PRESHOWER absorber c NMAX_SHOW number of blocks in SHOWER absorber c NHITS_BLOCK maximum # hits per shower block integer NPLANE_VDC,NMAX_WC,NHITS_WC_MIN,NHITS_WC,NGAP_CLUS, > NCLUS_WC,NHITS_CLUS,NHITS_CLUS_MIN,NMAX_S0,NMAX_S1, > NMAX_S2,NMAX_S3,NHITS_PD,NMAX_AERO,NMAX_GAS, > NHITS_MIR,NMAX_PRESHOW,NMAX_SHOW,NHITS_BLOCK parameter(NPLANE_VDC=4,NMAX_WC=368,NHITS_WC_MIN=2,NGAP_CLUS=2, > NCLUS_WC=4,NHITS_CLUS=16,NHITS_CLUS_MIN=2, > NHITS_WC=NCLUS_WC*NHITS_CLUS,NMAX_S0=1,NMAX_S1=6, > NMAX_S2=6,NMAX_S3=9, > NHITS_PD=2,NMAX_AERO=26,NMAX_GAS=10, > NHITS_MIR=2,NMAX_PRESHOW=48,NMAX_SHOW=96, > NHITS_BLOCK=2) c NPLANE_SCINT number of scintillator planes c NFPPWIRES maximum number of wire groups per plane c NFPPMAX maximum number of hits per plane per event c NFPP_PLANES number of FPP planes c IMAXGEOM number of geometry parameters per detector integer NPLANE_SCINT,NFPPWIRES,NFPPMAX,NFPP_PLANES,IMAXGEOM parameter (NPLANE_SCINT=4,NFPPWIRES=35,NFPPMAX=35, > NFPP_PLANES=24,IMAXGEOM=12) c c MWK980218: spares added. c integer*4 NMAX_SPARE1 ! not used? integer*4 NMAX_SPARE2 ! not used? integer*4 NMAX_SPARE3 ! not used? integer*4 NMAX_SPARE4 ! not used? integer*4 NMAX_SPARE5 ! not used? integer*4 NMAX_SPARE6 ! not used? integer*4 NMAX_SPARE7 ! Spare7 is used to read the latched ! trigger pattern. The trigger ! supervisor had 12 output corresponding ! to 12 inputs AFTER prescaling. integer*4 NMAX_SPARE8 ! not used? integer*4 NMAX_SPARE9 ! not used? integer*4 NMAX_SPARE10 ! not used? parameter ( NMAX_SPARE1 = 0 ) parameter ( NMAX_SPARE2 = 0 ) parameter ( NMAX_SPARE3 = 0 ) parameter ( NMAX_SPARE4 = 0 ) parameter ( NMAX_SPARE5 = 0 ) parameter ( NMAX_SPARE6 = 0 ) parameter ( NMAX_SPARE7 = 12 ) parameter ( NMAX_SPARE8 = 0 ) parameter ( NMAX_SPARE9 = 0 ) parameter ( NMAX_SPARE10 = 0 ) c c TDC_AVER average scintillator TDC value after correction real*8 TDC_AVER parameter (TDC_AVER=2000.d0) c laser hitpattern definition integer PULSER_HIT parameter (PULSER_HIT=50) c detector hitpattern definition c++++spectrometer E integer GAS_HIT_E,S0_HIT_E,S1_HIT_E,S2_HIT_E, > VDC1U_HIT_E,VDC1V_HIT_E,VDC2U_HIT_E,VDC2V_HIT_E, > AERO_HIT_E,PRESHOW_HIT_E,SHOW_HIT_E parameter(VDC1U_HIT_E = 1, > VDC1V_HIT_E = 2, > VDC2U_HIT_E = 3, > VDC2V_HIT_E = 4, > S1_HIT_E = 5, > S2_HIT_E = 6, > AERO_HIT_E = 7, > GAS_HIT_E = 8, > PRESHOW_HIT_E = 9, > SHOW_HIT_E = 10, > S0_HIT_E = 11 ) c++++spectrometer H integer GAS_HIT_H,S0_HIT_H,S1_HIT_H,S2_HIT_H, > VDC1U_HIT_H,VDC1V_HIT_H,VDC2U_HIT_H,VDC2V_HIT_H, > AERO_HIT_H,FPPFU1_HIT_H, > FPPFU2_HIT_H,FPPFU3_HIT_H,FPPFU4_HIT_H,FPPFU5_HIT_H, > FPPFU6_HIT_H,FPPRU1_HIT_H,FPPRU2_HIT_H,FPPRU3_HIT_H, > FPPRU4_HIT_H,FPPRU5_HIT_H,FPPFV1_HIT_H,FPPFV2_HIT_H, > FPPFV3_HIT_H,FPPFV4_HIT_H,FPPFV5_HIT_H,FPPFV6_HIT_H, > FPPRV1_HIT_H,FPPRV2_HIT_H,FPPRV3_HIT_H,FPPRV4_HIT_H, > FPPRV5_HIT_H,FPPRX1_HIT_H,FPPRX2_HIT_H,S3_HIT_H parameter(VDC1U_HIT_H = 17, > VDC1V_HIT_H = 18, > VDC2U_HIT_H = 19, > VDC2V_HIT_H = 20, > S1_HIT_H = 21, > S2_HIT_H = 22, > AERO_HIT_H = 23, > GAS_HIT_H = 24, > FPPFV1_HIT_H = 25, > FPPFV2_HIT_H = 26, > FPPFV3_HIT_H = 27, > FPPFU1_HIT_H = 28, > FPPFU2_HIT_H = 29, > FPPFU3_HIT_H = 30, > FPPFV4_HIT_H = 31, > FPPFV5_HIT_H = 32, > FPPFV6_HIT_H = 33, > FPPFU4_HIT_H = 34, > FPPFU5_HIT_H = 35, > FPPFU6_HIT_H = 36, > FPPRU1_HIT_H = 37, > FPPRU2_HIT_H = 38, > FPPRV1_HIT_H = 39, > FPPRV2_HIT_H = 40, > FPPRX1_HIT_H = 41, > FPPRX2_HIT_H = 42, > FPPRU3_HIT_H = 43, > FPPRU4_HIT_H = 44, > FPPRU5_HIT_H = 45, > FPPRV3_HIT_H = 46, > FPPRV4_HIT_H = 47, > FPPRV5_HIT_H = 48, > S3_HIT_H = 49, > S0_HIT_H = 50 ) c error modes for the VDC's c ERR1: a nohit or all clusters are of ERR1 type. c ERR2: wp.hits > NHITS_WC c ERR3: some drift times resulting in drift distances larger c than 1.25 times distance between wire plane and cathode c ERR4: some drift times are correcting for the stop offset c negative c error modes for the cluster's c ERR1: NHITS_CLUS_MIN > wp.clus(j).nhits c ERR2: wp.clus(j).nhits > NHITS_CLUS c ERR3: The 1st or wp.clus(j).nhits wire had the shortest c drift time. (smallest time is external) c ERR4: TDC value out-of-range c ERR5: fit of cluster failed c error modes for the track's c ERR1: Tracks share the same link c ERR2: Track projection to the scintillator plane corresponds c to a paddle that did not fire. integer VDC_OK,VDC_ERR1,VDC_ERR2,VDC_ERR3,VDC_ERR4 parameter(VDC_OK = 1) parameter(VDC_ERR1 = 2, > VDC_ERR2 = 3, > VDC_ERR3 = 4, > VDC_ERR4 = 5) integer CLUS_OK,CLUS_ERR1,CLUS_ERR2,CLUS_ERR3,CLUS_ERR4,CLUS_ERR5 parameter(CLUS_OK = 1) parameter(CLUS_ERR1 = 2, > CLUS_ERR2 = 3, > CLUS_ERR3 = 4, > CLUS_ERR4 = 5, > CLUS_ERR5 = 6) integer TRACK_OK,TRACK_ERR1,TRACK_ERR2 parameter(TRACK_OK = 1) parameter(TRACK_ERR1 = 2, > TRACK_ERR2 = 3) integer ID_TAB_OFF ! HBOOK ID offset used for tx lookup ! tables for each VDC wire plane parameter(ID_TAB_OFF=4900) c---optimize parameter integer MAXPEAKS, ! maximum number of peaks that can used in ! optimizing the matrix elements for thtg,ytg, ! phtg,delta > NCAL, ! maximum number of reference points used in ! the calibration, f.i. sieve slit holes or ! y positions on the target > MAXBND, ! maximum number of boundary conditions ! in optimization > NWALK, ! maximum number of coeff used to correct for ! timewalk in the TOF paddle signal > MAXPNTS, ! maximum number of points used in the ! minimization of aberation functions > ID_FIT_OFF,! MBOOK ID offset used for fit NTUPLES > FOC_EVENT_STRUCT_LEN, ! number of REAL*8 elements in the structure ! FOC_EVENT > TAR_EVENT_STRUCT_LEN, ! number of REAL*8 elements in the structure ! TAR_EVENT > EVENT_STRUCT_LEN, ! number of REAL*8 elements in the structure ! EVENT > VDC_PLANE_STRUCT_LEN ! number of INTEGER elements in the structure ! VDC_PLANE that is changed for each event parameter(MAXPEAKS=50,NCAL=50,NWALK=1, > MAXBND=2, > MAXPNTS=100000, > ID_FIT_OFF=ID_TAB_OFF+2*NPLANE_VDC, > FOC_EVENT_STRUCT_LEN=32,TAR_EVENT_STRUCT_LEN=17, > EVENT_STRUCT_LEN=3+4*FOC_EVENT_STRUCT_LEN+ > TAR_EVENT_STRUCT_LEN, > VDC_PLANE_STRUCT_LEN=3+6*NHITS_WC+NCLUS_WC+ > NCLUS_WC*(11*NHITS_CLUS+19)) c c R. Suleiman 3/25/1998 c Changed MAXPNTS=10000 to MAXPNTS=100000. This parameter is also c used in coincidence timing optimization which limits the number of c events you can use in the optimization. c c---target parameter integer MAXISO, ! maximum number of isotopes in target > NFOIL ! maximum number of target foils parameter(MAXISO=3,NFOIL=3) c added 04 Mar 1998 by L.Todor for sieve coordinate reconstruction real DSIEVE_E, DSIEVE_H parameter(DSIEVE_E=1.1834,DSIEVE_H=1.1745) c added 13 OCt 1999 by M. Jones for ext_tar_cor.f c factors detemined from 1st order optics matrix c for correcting delta and theta_tg for x-target dependenc real DEL_XTG_COR,THT_XTG_COR parameter(DEL_XTG_COR=5.18,THT_XTG_COR=.61)