/[PAMELA software]/PamVMC/trk/src/f77/gpdspe.F
ViewVC logotype

Annotation of /PamVMC/trk/src/f77/gpdspe.F

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.5 - (hide annotations) (download)
Fri Jun 12 18:39:58 2009 UTC (15 years, 7 months ago) by pam-rm2
Branch: MAIN
CVS Tags: v1r0, HEAD
Changes since 1.1: +0 -0 lines
- Introduced user-defined names of output files and random seeds number.
Users can do it use options of PamVMCApplication constructor:
PamVMCApplication(const char* name,  const char *title, const char*
filename="pamtest", Int_t seed=0).
The Random object that I use is TRandom3 object which has astronomical
large period (in case of default initialization 0). All random generators
in the code use this object by calling of gRandom singleton which keeps
it.

- Corrected TOF digitization routine. No problems with TDC hits due to
hadronic interactions anymore.

- Some small changes was done to compile code under Root 5.23. +
geant4_vmc v. 2.6 without any warnings

- Some classes of PamG4RunConfiguartion was changed for geant4_vmc v.
2.6.Some obsolete classes was deleted as soon as developers implemented
regions.

- Navigation was changed from "geomRootToGeant4" to "geomRoot", because on
VMC web page written that as soon as Geant4 has no option ONLY/MANY
translation of overlapped geometry to Geant4 through VGM could be wrong.
I'd like to stay with Root navigation:
http://root.cern.ch/root/vmc/Geant4VMC.html. This should be default
option.

- New Tracker digitization routine written by Sergio was implemented

- PamVMC again became compatible with geant4_vmc v.2.5 and ROOT 5.20.
 The problem was that ROOT developers introduced in TVirtualMC class a new
method SetMagField and new base class:TVirtualMagField from which
user-defined classes shoukd be derived

1 nikolas 1.1 *
2    
3     SUBROUTINE GPDSPE
4     *************************************************************************
5     * *
6     * To perform SPE digitization *
7     * *
8     * Called by: GUDIGI *
9     * Author: Francesco Cafagna, 23/02/96 12.50.38 *
10     * *
11     *************************************************************************
12     PARAMETER (RADDEG=57.295779)
13    
14     *#include "gppdigi.inc"
15     c +SEQ,GPDSPE.
16     *#include "gpsed.inc"
17     *et. 04/11/2005
18     #include "gpstripspe.inc"
19     *end et.
20     *
21     * INTEGER IS,ID,NV,NVM,GPNHIT,NTHIT,NHFRST
22     * CHARACTER*4 CSET
23     * LOGICAL CONT,HF
24     *
25     * Data block for digitalization
26     *
27     * INTEGER NTRA,IDIG,LTRA,NPLANE,LKDN,NTDIGI,JJ
28     * REAL XI,YI,XO,YO,PATH,DE,GPGAUS
29     * REAL VPOS(3),BMAGNET(3),ANGLEX,ANGLEY,XMULT,YMULT,DXPOS,DYPOS
30     *
31     * DATA NTRA,LTRA/1,1/
32    
33     *et. 04/11/2005
34     *
35     *
36     * Vector zeroed
37     *
38     CALL VZERO(ISTRIPX,MAXSTR)
39     CALL VZERO(NPSTRIPX,MAXSTR)
40     CALL VZERO(NTSTRIPX,MAXSTR)
41     CALL VZERO(QSTRIPX,MAXSTR)
42     CALL VZERO(XSTRIPX,MAXSTR)
43     CALL VZERO(ISTRIPY,MAXSTR)
44     CALL VZERO(NPSTRIPY,MAXSTR)
45     CALL VZERO(NTSTRIPY,MAXSTR)
46     CALL VZERO(QSTRIPY,MAXSTR)
47     CALL VZERO(YSTRIPY,MAXSTR)
48     NSTRPX=0
49     NSTRPY=0
50    
51     *
52     do ip=1,nplanes
53     do iip=1,nsensors
54     isignyesx=0
55     isignyesy=0
56     do iiip=1,nstripx
57     if(proxtanti(ip,iip,iiip).gt.0) then
58     isignyesx=1
59     GO TO 456
60     endif
61     enddo
62     do iiip=1,nstripy
63     if(proytanti(ip,iip,iiip).gt.0) then
64     isignyesy=1
65     GO TO 456
66     endif
67     enddo
68     456 if(isignyesx.eq.1.or.isignyesy.eq.1) then
69     call gpucollection(ip,iip)
70     endif
71     enddo
72     enddo
73    
74    
75     * PRINT*,'NSTRPX', NSTRPX, 'NSTRPY', NSTRPY
76     c do i=1, nstrpx
77     c PRINT*,i,' NPX ',NPSTRIPX(i), ' NTX ', NTSTRIPX(i)
78     c PRINT*,i,' ISX ', ISTRIPX(i)
79     c PRINT*,i,' QX ',QSTRIPX(i),' XSX ', XSTRIPX(i)
80     c enddo
81     c do i=1, nstrpy
82     c PRINT*,i,' NPY ',NPSTRIPY(i), ' NTY ', NTSTRIPY(i)
83     c PRINT*,i,' ISY ', ISTRIPY(i)
84     c PRINT*,i,' QY ',QSTRIPY(i),' YSY ', YSTRIPY(i)
85     c enddo
86    
87    
88     20 RETURN
89     END

  ViewVC Help
Powered by ViewVC 1.1.23