Message: How to add some files in my code | Not Logged In (login) |
Dear all GEANT4 Users, Hi ! I need some help to built some files in my code; I have built the code in GEANT4, I am getting the hits values in two sensitive gas gaps, As Sensitive gas gap 1 and Sensitive Gas gap 2
I want to make my calculations in Gas gap 1
Gas gap 2
Gas gap 1 and 2
and Gas gap 1 Or gas gap 2. I found in my G321 code I have these defined in GUSTEP.f and GUOUT.f Which are givem as; ================================================================ SUBROUTINE GUSTEP common / iprimary / iflag_prim common / iflag_1 / iflag1,iflag2 common /iflag_2 /iflag_and,iflag_or call hf2(21,vect(1),vect(2),1.) call hf2(22,vect(1),vect(3),1.) call hf2(23,vect(2),vect(3),1.) if(ngkine.GT.0) call gsking(0) ! traccia tutti i secondari if(ntmult.eq.1.and.ipart.eq.ikine) then ! primary particle gamma iflag_prim = 1 endif c particlle secondarie c if(ipart.eq.ikine) then if(inwvol.eq.1.and.numed.eq.14) then call hf1(30,GEKIN,1.) endif endif if(abs(charge).gt.0) then ! solo particelle secondarie cariche c I gap c if(inwvol.eq.1.and.numed.eq.1) then iflag1=1 endif ! fine I gap c II gap c if(inwvol.eq.1.and.numed.eq.2) then iflag2=1 endif ! fine II gap c endif ! fine solo particelle secondarie cariche c C I AND II C if (iflag1.eq.1.and.iflag2.eq.1) then iflag_and=1 endif C C C I OR II if(iflag1.eq.1) then iflag_or = 1 endif c call gdebug END =================================================== c---------------------------------------------------------------------------- subroutine guout COMMON /GCNUMX/ NALIVE,NTMSTO common / iprimary / iflag_prim common / iflag_1 / iflag1,iflag2 common /iflag_2 /iflag_and,iflag_or common / itracce / nprim,isign1,isign2,isign_and,isign_or if (iflag_prim.eq.1) then nprim = nprim + 1 iflag_prim = 0 endif if (iflag1.eq.1) then isign1=isign1+1 iflag1=0 call hf2(24,vect(1),vect(2),1.) call hf2(25,vect(1),vect(3),1.) call hf2(26,vect(2),vect(3),1.) endif C if (iflag2.eq.1) then isign2=isign2+1 iflag2=0 call hf2(27,vect(1),vect(2),1.) call hf2(28,vect(1),vect(3),1.) call hf2(29,vect(2),vect(3),1.) endif C if(iflag_and.eq.1)then isign_and=isign_and+1 iflag_and=0 endif C if(iflag_or.eq.1)then isign_or=isign_or+1 iflag_or=0 endif return end I hope you will help me, I am waiting for your reply, jamil |
Inline Depth: | Outline Depth: | Add message: |
to: |