| 1 |
cafagna |
3.1 |
* |
| 2 |
|
|
* $Id$ |
| 3 |
|
|
* |
| 4 |
|
|
* $Log$ |
| 5 |
|
|
* |
| 6 |
|
|
*CMZ : 2.01/00 06/03/2000 13.07.03 by Francesco Cafagna |
| 7 |
|
|
*CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna |
| 8 |
|
|
*CMZ : 1.02/00 22/10/99 11.10.56 by Francesco Cafagna |
| 9 |
|
|
*CMZ : 1.01/00 21/05/96 08.53.57 by Francesco Cafagna |
| 10 |
|
|
*-- Author : Paolo (The Magic) Papini 23/04/96 |
| 11 |
|
|
#if defined(GPAMELA_PAOLO) |
| 12 |
|
|
SUBROUTINE GPSPEHD(NHIT,NPLANE,XI,YI,ZI,XO,YO,ZO,PATH,DE,LKDN) |
| 13 |
|
|
#endif |
| 14 |
|
|
#if !defined(GPAMELA_PAOLO) |
| 15 |
|
|
SUBROUTINE GPSPEHD(NHIT,NUMV,NLEV,XI,YI,ZI,XO,YO,ZO,PATH,DE, |
| 16 |
|
|
+ LKDN) |
| 17 |
|
|
#endif |
| 18 |
|
|
************************************************************************ |
| 19 |
|
|
* * |
| 20 |
|
|
* To perform hit/digit conversion for spectrometer. * |
| 21 |
|
|
* This routine fills the LKD1,LKD2,LKD3,LKD4 array * |
| 22 |
|
|
* and it up-date LKND as total number of digit * |
| 23 |
|
|
* * |
| 24 |
|
|
* Variables definition: * |
| 25 |
|
|
* IN: * |
| 26 |
|
|
* NHIT - hit number * |
| 27 |
|
|
* NPLANE - plane number * |
| 28 |
|
|
* XI,YI,ZI,XO,YO,ZO - input and output hit coordinates * |
| 29 |
|
|
* PATH,DE - hit path lenght and dE/dx * |
| 30 |
|
|
* * |
| 31 |
|
|
* OUT: * |
| 32 |
|
|
* LKDN - total number of digit * |
| 33 |
|
|
* * |
| 34 |
|
|
* Called by: GPDSPE * |
| 35 |
|
|
* Author: Paolo (The Magic) Papini, 23/04/96 15.45.26 * |
| 36 |
|
|
* * |
| 37 |
|
|
************************************************************************ |
| 38 |
|
|
#include "gcunit.inc" |
| 39 |
|
|
#include "gpphit.inc" |
| 40 |
|
|
#include "gppdigi.inc" |
| 41 |
|
|
PARAMETER (NH=NHSPE,NHM=NHMSPE,NVOL=4) |
| 42 |
|
|
#include "gpspe.inc" |
| 43 |
|
|
#include "gpsed.inc" |
| 44 |
|
|
#include "gpgeo.inc" |
| 45 |
|
|
#include "gpdspe.inc" |
| 46 |
|
|
* |
| 47 |
|
|
* Data block for digitalization |
| 48 |
|
|
* |
| 49 |
|
|
INTEGER NHIT,NPLANE,LKDN,NSIL,NSTRIP,NADC,NEWD |
| 50 |
|
|
REAL XSIL,YSIL,XBORD,YBORD |
| 51 |
|
|
REAL XI,YI,XO,YO,PATH,DE,CH_DIG(2048) |
| 52 |
|
|
REAL RR(3),BB(3) |
| 53 |
|
|
#if !defined(GPAMELA_PAOLO) |
| 54 |
|
|
INTEGER LNAM,LNUM |
| 55 |
|
|
DIMENSION NUMV(NLEV),LNAM(NVOL),LNUM(NVOL) |
| 56 |
|
|
REAL XYZISI(3),XYZOSI(3),XYZIMA(3),XYZOMA(3),XYZMO(3),XYZSO(3) |
| 57 |
|
|
#endif |
| 58 |
|
|
* |
| 59 |
|
|
BB(1)=0. |
| 60 |
|
|
BB(2)=0. |
| 61 |
|
|
BB(3)=0. |
| 62 |
|
|
#if defined(GPAMELA_MAGN) |
| 63 |
|
|
RR(1)=(XI+XO)/2. |
| 64 |
|
|
RR(2)=(YI+YO)/2. |
| 65 |
|
|
RR(3)=(ZI+ZO)/2. |
| 66 |
|
|
CALL GUFLD(RR,BB) |
| 67 |
|
|
#endif |
| 68 |
|
|
NSEPX=1+INT(ABS(XO-XI)/STEPX) |
| 69 |
|
|
NSEPY=1+INT(ABS(YO-YI)/STEPY) |
| 70 |
|
|
NSEPZ=1+INT(ABS(ZO-ZI)/STEPZ) |
| 71 |
|
|
#if defined(GPAMELA_PAOLO) |
| 72 |
|
|
C* |
| 73 |
|
|
C* X-view digitalization |
| 74 |
|
|
C* |
| 75 |
|
|
C NSIL=INT(((XI+XO)/2.+3.*TSPA(1))/(2.*TSPA(1)))+1 |
| 76 |
|
|
C* go way if no in the silicon |
| 77 |
|
|
C IF(NSIL.LT.1.OR.NSIL.GT.3) GOTO 10 |
| 78 |
|
|
C* n.division of hit path |
| 79 |
|
|
C NSEP=MAX0(NSEPX,NSEPZ) |
| 80 |
|
|
C* energy loss in each path division |
| 81 |
|
|
C DDE=DE/FLOAT(NSEP) |
| 82 |
|
|
C* XI-ccor. in the silicon |
| 83 |
|
|
C XSIL_I=MOD(XI+3.*TSPA(1),2.*TSPA(1)) |
| 84 |
|
|
C* XO-coor. in the silicon |
| 85 |
|
|
C XSIL_O=MOD(XO+3.*TSPA(1),2.*TSPA(1)) |
| 86 |
|
|
C* x-coor. of sensitive edge |
| 87 |
|
|
C XBORD=(2.*TSPA(1)-NXSTR*WXSTR)/2. |
| 88 |
|
|
C* XI of the strip |
| 89 |
|
|
C NSTR_I=INT((XSIL_I-XBORD)/WXSTR)+1 |
| 90 |
|
|
C* XO of the strip |
| 91 |
|
|
C NSTR_O=INT((XSIL_O-XBORD)/WXSTR)+1 |
| 92 |
|
|
C NSTR_MIN=MIN0(NSTR_I,NSTR_O)-10 |
| 93 |
|
|
C* minimum strip to be calculated |
| 94 |
|
|
C IF(NSTR_MIN.LT.1) NSTR_MIN=1 |
| 95 |
|
|
C NSTR_MAX=MAX0(NSTR_I,NSTR_O)+10 |
| 96 |
|
|
C* maximum strip to be calculated |
| 97 |
|
|
C IF(NSTR_MAX.GT.NXSTR) NSTR_MAX=NXSTR |
| 98 |
|
|
C+SELF |
| 99 |
|
|
#endif |
| 100 |
|
|
#if !defined(GPAMELA_PAOLO) |
| 101 |
|
|
* |
| 102 |
|
|
* Fill up the GCVOLU common block |
| 103 |
|
|
* |
| 104 |
|
|
* CALL GFPATH (ISSPE,IDTSPA,NUMV,NLEV,LNAM,LNUM) |
| 105 |
|
|
c LNAM(1) =4HPAME |
| 106 |
|
|
c LNAM(2) =4HSPEB |
| 107 |
|
|
LNAM(1) =4HTRPB |
| 108 |
|
|
LNAM(2) =4HTRPL |
| 109 |
|
|
LNAM(3) =4HTRSL |
| 110 |
|
|
LNAM(4) =4HTSPA |
| 111 |
|
|
CALL GLVOLU (NLEV,LNAM,NUMV,IER) |
| 112 |
|
|
* |
| 113 |
|
|
* X-view digitalization |
| 114 |
|
|
* |
| 115 |
|
|
NSIL=NUMV(NLEV-1) |
| 116 |
|
|
* go way if no in the silicon |
| 117 |
|
|
IF(NSIL.LT.1.OR.NSIL.GT.3) GOTO 10 |
| 118 |
|
|
* |
| 119 |
|
|
* Calculate the stripes |
| 120 |
|
|
* |
| 121 |
|
|
NPLANE = NUMV(NLEV-3) |
| 122 |
|
|
* n.division of hit path |
| 123 |
|
|
NSEP=MAX0(NSEPX,NSEPZ) |
| 124 |
|
|
* energy loss in each path division |
| 125 |
|
|
DDE=DE/FLOAT(NSEP) |
| 126 |
|
|
* trasform coordinates from MARS to silicon center |
| 127 |
|
|
XYZIMA(1) = XI |
| 128 |
|
|
XYZIMA(2) = YI |
| 129 |
|
|
XYZIMA(3) = ZI |
| 130 |
|
|
CALL GMTOD(XYZIMA,XYZISI,1) |
| 131 |
|
|
XYZOMA(1) = XO |
| 132 |
|
|
XYZOMA(2) = YO |
| 133 |
|
|
XYZOMA(3) = ZO |
| 134 |
|
|
CALL GMTOD(XYZOMA,XYZOSI,1) |
| 135 |
|
|
* xy-coor. of sensitive edge |
| 136 |
|
|
XBORD=(2.*TSPA(1)-NXSTR*WXSTR)/2. |
| 137 |
|
|
YBORD=(2.*TSPA(2)-NYSTR*WYSTR)/2. |
| 138 |
|
|
* xyz-coord. of the paddle origin in MARS |
| 139 |
|
|
CALL VZERO(XYZSO,3) |
| 140 |
|
|
CALL GDTOM(XYZSO,XYZMO,1) |
| 141 |
|
|
* Move origin of coordinates to the paddle edge, a guard zone added |
| 142 |
|
|
XYZISI(1) = XYZISI(1) + (TSPA(1) - XBORD) |
| 143 |
|
|
XYZISI(2) = XYZISI(2) + (TSPA(2) - YBORD) |
| 144 |
|
|
* |
| 145 |
|
|
XYZOSI(1) = XYZOSI(1) + (TSPA(1) - XBORD) |
| 146 |
|
|
XYZOSI(2) = XYZOSI(2) + (TSPA(2) - YBORD) |
| 147 |
|
|
* Calculated the strip number now |
| 148 |
|
|
NSTR_I=INT(XYZISI(1)/WXSTR)+1 |
| 149 |
|
|
NSTR_O=INT(XYZOSI(1)/WXSTR)+1 |
| 150 |
|
|
NSTR_MIN=MIN0(NSTR_I,NSTR_O)-10 |
| 151 |
|
|
* minimum strip to be calculated |
| 152 |
|
|
IF(NSTR_MIN.LT.1) NSTR_MIN=1 |
| 153 |
|
|
NSTR_MAX=MAX0(NSTR_I,NSTR_O)+10 |
| 154 |
|
|
* maximum strip to be calculated |
| 155 |
|
|
IF(NSTR_MAX.GT.NXSTR) NSTR_MAX=NXSTR |
| 156 |
|
|
#endif |
| 157 |
|
|
C PRINT *,'X DIGIT, N.HIT:',NHIT |
| 158 |
|
|
C PRINT *,'NSTR_MIN,NSTR_MAX,NSEP',NSTR_MIN,NSTR_MAX,NSEP |
| 159 |
|
|
|
| 160 |
|
|
* loop on the strip (for the diffusion of Q) |
| 161 |
|
|
DO IS=NSTR_MIN,NSTR_MAX |
| 162 |
|
|
* x-center of the strip IS |
| 163 |
|
|
XSTR=XBORD+(IS-0.5)*WXSTR+(2*NSIL-5)*TSPA(1) |
| 164 |
|
|
C+SELF,IF=PAOLO |
| 165 |
|
|
C* z-center of the strip |
| 166 |
|
|
C ZSTR=ZSPEC+(3.5-NPLANE)*2.*MGPL(3)+(3-NPLANE)*2.*TRPB(3)+ |
| 167 |
|
|
C + 2.*TRCP(3)+2.*TRRP(3)+2.*TSKP(3) |
| 168 |
|
|
C+SELF |
| 169 |
|
|
#if !defined(GPAMELA_PAOLO) |
| 170 |
|
|
* XSTR=XYZMO(1) |
| 171 |
|
|
* Z-coord. of the implanted strip. The upper is the Y view |
| 172 |
|
|
ZSTR=XYZMO(3)-TSPA(3) |
| 173 |
|
|
#endif |
| 174 |
|
|
* magnetic field effect |
| 175 |
|
|
XI=XI-ABS(ZI-ZSTR)*0.2475*BB(2)/15. |
| 176 |
|
|
XO=XO-ABS(ZO-ZSTR)*0.2475*BB(2)/15. |
| 177 |
|
|
* |
| 178 |
|
|
C PRINT *,ZSPEC,NPLANE,MGPL(3),TRPB(3),TRCP(3),TRRP(3),TSKP(3) |
| 179 |
|
|
* charge definition for the strip IS |
| 180 |
|
|
CH_DIG(IS)=0. |
| 181 |
|
|
* loop on the hit path division |
| 182 |
|
|
DO IH=1,NSEP |
| 183 |
|
|
* x-coor. hit division |
| 184 |
|
|
XH=XI+(IH-0.5)*(XO-XI)/FLOAT(NSEP) |
| 185 |
|
|
* Z-coor. hit division |
| 186 |
|
|
ZH=ZI+(IH-0.5)*(ZO-ZI)/FLOAT(NSEP) |
| 187 |
|
|
IX=1+INT(ABS(XH-XSTR)/STEPX) |
| 188 |
|
|
IZ=1+INT(ABS(ZH-ZSTR)/STEPZ) |
| 189 |
|
|
C PRINT *,'IX,IZ',IX,IZ,XH,XSTR,ZH,ZSTR |
| 190 |
|
|
C PRINT *,IS,IH,IX,IZ,CHARX(IX,IZ) |
| 191 |
|
|
IF(IX.LE.NPCHXY.AND.IZ.LE.NPCHZ) THEN |
| 192 |
|
|
IF(MOD(IS,2).EQ.1) THEN |
| 193 |
|
|
CH_DIG(IS)=CH_DIG(IS)+DDE*CHARX(IX,IZ) |
| 194 |
|
|
ELSE |
| 195 |
|
|
CH_DIG(IS)=CH_DIG(IS)+DDE*CHARX_I(IX,IZ) |
| 196 |
|
|
ENDIF |
| 197 |
|
|
ENDIF |
| 198 |
|
|
ENDDO |
| 199 |
|
|
END DO |
| 200 |
|
|
* loop for C-coupling |
| 201 |
|
|
DO IS=NSTR_MIN,NSTR_MAX |
| 202 |
|
|
* X-view has double strip (ISR is the read strip) |
| 203 |
|
|
ISR=INT((IS+1)/2) |
| 204 |
|
|
CHARGE=0. |
| 205 |
|
|
DO ISSR=ISR-2,ISR+2,1 |
| 206 |
|
|
IF(2*ISSR-2.GE.NSTR_MIN.AND.2*ISSR.LE.NSTR_MAX) CHARGE= |
| 207 |
|
|
+ CHARGE+CCX(1+IABS(ISSR-ISR))* (CH_DIG(2*ISSR-1)+CH_DIG(2* |
| 208 |
|
|
+ ISSR)/2+CH_DIG(2*ISSR-2)/2) |
| 209 |
|
|
ENDDO |
| 210 |
|
|
C CHARGE=CH_DIG(IS) |
| 211 |
|
|
* |
| 212 |
|
|
* digitaz record |
| 213 |
|
|
NEWD=1 |
| 214 |
|
|
* X-view has double strip (ISR is the read strip) |
| 215 |
|
|
IF(LKDN.NE.0) THEN |
| 216 |
|
|
DO I=1,LKDN |
| 217 |
|
|
IF(NPLANE.EQ.LKD1(I).AND.NSIL.EQ.LKD2(I) |
| 218 |
|
|
+ .AND.ISR.EQ.LKD3(I)) THEN |
| 219 |
|
|
NEWD=0 |
| 220 |
|
|
RLKD4(I)=RLKD4(I)+CHARGE |
| 221 |
|
|
ENDIF |
| 222 |
|
|
ENDDO |
| 223 |
|
|
ENDIF |
| 224 |
|
|
IF(NEWD.EQ.1) THEN |
| 225 |
|
|
LKDN=LKDN+1 |
| 226 |
|
|
IF (LKDN.GT.NDMSPE) THEN |
| 227 |
|
|
WRITE(10000,CHMAIL) LKDN |
| 228 |
|
|
CALL GMAIL(1,0) |
| 229 |
|
|
GO TO 30 |
| 230 |
|
|
ENDIF |
| 231 |
|
|
LKD1(LKDN)=NPLANE |
| 232 |
|
|
LKD2(LKDN)=NSIL |
| 233 |
|
|
LKD3(LKDN)=ISR |
| 234 |
|
|
RLKD4(LKDN)=CHARGE |
| 235 |
|
|
LKDH(LKDN)=NHIT |
| 236 |
|
|
ENDIF |
| 237 |
|
|
ENDDO |
| 238 |
|
|
10 CONTINUE |
| 239 |
|
|
* |
| 240 |
|
|
* Y-view digitalization |
| 241 |
|
|
* |
| 242 |
|
|
NSEP=MAX0(NSEPY,NSEPZ) |
| 243 |
|
|
DDE=DE/NSEP |
| 244 |
|
|
C+SELF,IF=PAOLO |
| 245 |
|
|
C NSIL=INT(((YI+YO)/2.+2.*TSPA(2))/(2.*TSPA(2)))+4 |
| 246 |
|
|
C IF(NSIL.LT.4.OR.NSIL.GT.5) GOTO 20 |
| 247 |
|
|
C YSIL_I=MOD(YI+2.*TSPA(2),2.*TSPA(2)) |
| 248 |
|
|
C YSIL_O=MOD(YO+2.*TSPA(2),2.*TSPA(2)) |
| 249 |
|
|
C YBORD=(2.*TSPA(2)-NYSTR*WYSTR)/2. |
| 250 |
|
|
C NSTR_I=INT((YSIL_I-YBORD)/WYSTR)+1 |
| 251 |
|
|
C NSTR_O=INT((YSIL_O-YBORD)/WYSTR)+1 |
| 252 |
|
|
C NSTR_MIN=MIN0(NSTR_I,NSTR_O)-5 |
| 253 |
|
|
C IF(NSTR_MIN.LT.1) NSTR_MIN=1 |
| 254 |
|
|
C NSTR_MAX=MAX0(NSTR_I,NSTR_O)+5 |
| 255 |
|
|
C IF(NSTR_MAX.GT.NYSTR) NSTR_MAX=NYSTR |
| 256 |
|
|
C+SELF. |
| 257 |
|
|
#if !defined(GPAMELA_PAOLO) |
| 258 |
|
|
NSIL = 5 - (NUMV(NLEV)-1) |
| 259 |
|
|
IF(NSIL.LT.4.OR.NSIL.GT.5) GOTO 20 |
| 260 |
|
|
* Calculate the strip number now |
| 261 |
|
|
NSTR_I=INT((XYZISI(2))/WYSTR)+1 |
| 262 |
|
|
NSTR_O=INT((XYZOSI(2))/WYSTR)+1 |
| 263 |
|
|
NSTR_MIN=MIN0(NSTR_I,NSTR_O)-5 |
| 264 |
|
|
* minimum strip to be calculated |
| 265 |
|
|
IF(NSTR_MIN.LT.1) NSTR_MIN=1 |
| 266 |
|
|
NSTR_MAX=MAX0(NSTR_I,NSTR_O)+5 |
| 267 |
|
|
* maximum strip to be calculated |
| 268 |
|
|
IF(NSTR_MAX.GT.NYSTR) NSTR_MAX=NYSTR |
| 269 |
|
|
#endif |
| 270 |
|
|
DO IS=NSTR_MIN,NSTR_MAX |
| 271 |
|
|
YSTR=YBORD+(IS-0.5)*WYSTR+(NSIL-5)*2.*TSPA(2) |
| 272 |
|
|
C+SELF,IF=PAOLO |
| 273 |
|
|
C ZSTR=ZSPEC+(3.5-NPLANE)*2.*MGPL(3)+(3-NPLANE)*2.*TRPB(3)+ |
| 274 |
|
|
C + 2.*TRCP(3)+2.*TRRP(3)+2.*TSKP(3)+2.*TSPA(3) |
| 275 |
|
|
C+SELF |
| 276 |
|
|
#if !defined(GPAMELA_PAOLO) |
| 277 |
|
|
* YSTR=XYZMO(2) |
| 278 |
|
|
* Z-coord. of the implanted strip. The lower is the X view |
| 279 |
|
|
ZSTR=XYZMO(3) + TSPA(3) |
| 280 |
|
|
#endif |
| 281 |
|
|
* magnetic field effect |
| 282 |
|
|
XI=XI+ABS(ZI-ZSTR)*0.04716*BB(2)/15. |
| 283 |
|
|
XO=XO+ABS(ZO-ZSTR)*0.04716*BB(2)/15. |
| 284 |
|
|
* |
| 285 |
|
|
CH_DIG(IS)=0. |
| 286 |
|
|
DO IH=1,NSEP |
| 287 |
|
|
YH=YI+(IH-0.5)*(YO-YI)/NSEP |
| 288 |
|
|
ZH=ZI+(IH-0.5)*(ZO-ZI)/NSEP |
| 289 |
|
|
IY=1+INT(ABS(YH-YSTR)/STEPY) |
| 290 |
|
|
IZ=1+INT(ABS(ZH-ZSTR)/STEPZ) |
| 291 |
|
|
C PRINT *,IS,IH,IY,IZ,CHARY(IY,IZ) |
| 292 |
|
|
IF(IY.LE.NPCHXY.AND.IZ.LE.NPCHZ) CH_DIG(IS)=CH_DIG(IS)+DDE* |
| 293 |
|
|
+ CHARY(IY,IZ) |
| 294 |
|
|
ENDDO |
| 295 |
|
|
ENDDO |
| 296 |
|
|
DO IS=NSTR_MIN,NSTR_MAX |
| 297 |
|
|
CHARGE=0. |
| 298 |
|
|
DO ISS=IS-2,IS+2,1 |
| 299 |
|
|
IF(ISS.GE.NSTR_MIN.AND.ISS.LE.NSTR_MAX) CHARGE= |
| 300 |
|
|
+ CHARGE+CCY(1+IABS(ISS-IS))*CH_DIG(ISS) |
| 301 |
|
|
ENDDO |
| 302 |
|
|
* |
| 303 |
|
|
* digitaz record |
| 304 |
|
|
NEWD=1 |
| 305 |
|
|
* the Y view paddles are connected |
| 306 |
|
|
IF(NSIL.EQ.5) NSIL=4 |
| 307 |
|
|
IF(LKDN.NE.0) THEN |
| 308 |
|
|
DO I=1,LKDN |
| 309 |
|
|
IF(NPLANE.EQ.LKD1(I).AND.NSIL.EQ.LKD2(I) |
| 310 |
|
|
+ .AND.IS.EQ.LKD3(I)) THEN |
| 311 |
|
|
NEWD=0 |
| 312 |
|
|
RLKD4(I)=RLKD4(I)+CHARGE |
| 313 |
|
|
ENDIF |
| 314 |
|
|
ENDDO |
| 315 |
|
|
ENDIF |
| 316 |
|
|
IF(NEWD.EQ.1) THEN |
| 317 |
|
|
LKDN=LKDN+1 |
| 318 |
|
|
IF (LKDN.GT.NDMSPE) THEN |
| 319 |
|
|
WRITE(CHMAIL,10000) LKDN |
| 320 |
|
|
CALL GMAIL(1,0) |
| 321 |
|
|
GO TO 30 |
| 322 |
|
|
ENDIF |
| 323 |
|
|
LKD1(LKDN)=NPLANE |
| 324 |
|
|
LKD2(LKDN)=NSIL |
| 325 |
|
|
LKD3(LKDN)=IS |
| 326 |
|
|
RLKD4(LKDN)=CHARGE |
| 327 |
|
|
LKDH(LKDN)=NHIT |
| 328 |
|
|
ENDIF |
| 329 |
|
|
ENDDO |
| 330 |
|
|
20 CONTINUE |
| 331 |
|
|
* |
| 332 |
|
|
10000 FORMAT(' GPSPEHD: too much digitalization= ',I5) |
| 333 |
|
|
30 RETURN |
| 334 |
|
|
END |
Parent Directory
| 
Revision Log