Changeset 1337 for trunk/examples/extended/electromagnetic/TestEm2/geant3/src

Timestamp:

Sep 30, 2010, 2:47:17 PM (14 years ago)

Author:

garnier

Message:

tag geant4.9.4 beta 1 + modifs locales

Location:

trunk/examples/extended/electromagnetic/TestEm2/geant3/src

Files:

• guout.F (modified) (2 diffs)
• uglast.F (modified) (2 diffs)
• uhinit.F (modified) (3 diffs)

trunk/examples/extended/electromagnetic/TestEm2/geant3/src/guout.F

@@ -37 +37 @@
          CALL HFILL(4,BIN,100*DEDL(I)/(DLX0*PKINE(3)),1.)
          BIN = FLOAT(I)*DLX0
-         CALL HFILL(5,BIN,100*DLC    /PKINE(3),1.)                                                                                                                 
+         CALL HFILL(6,BIN,100*DLC    /PKINE(3),1.)                                                                                                                 
    2  CONTINUE                                                                  
 
@@ -50 +50 @@
          
          BIN = (FLOAT(I)-0.5)*DRX0
-         CALL HFILL(7,BIN,100*DEDR(I)/(DRX0*PKINE(3)),1.)
+         CALL HFILL(8,BIN,100*DEDR(I)/(DRX0*PKINE(3)),1.)
          BIN = FLOAT(I)*DRX0
-         CALL HFILL(8,BIN,100*DRC    /PKINE(3),1.)                                                  
+         CALL HFILL(10,BIN,100*DRC    /PKINE(3),1.)                                                  
    3  CONTINUE                                                                  
 *

trunk/examples/extended/electromagnetic/TestEm2/geant3/src/uglast.F

@@ -30 +30 @@
          XSEL2C(I) = CNORM*SQRT(ABS(XEVENT*SEL2C(I) - SEL1C(I)**2))
     2 CONTINUE
-         CALL HPAK (6,XSEL2C)                                                                  
+         CALL HPAK (5,XSEL2)    
+         CALL HPAK (7,XSEL2C)
 *
 * *** radial profile                                                                                
@@ -39 +40 @@
          XSER2C(I) = CNORM*SQRT(ABS(XEVENT*SER2C(I) - SER1C(I)**2))             
     3 CONTINUE
-         CALL HPAK (9,XSER2C)
+         CALL HPAK ( 9,XSER2)    
+         CALL HPAK (11,XSER2C)
 *
 * *** total track length

trunk/examples/extended/electromagnetic/TestEm2/geant3/src/uhinit.F

@@ -25 +25 @@
       ZMAX=NLTOT*DLX0
                                                                                            
-      CALL HBPROF(4,'longit energy profile    (in percent of E inc)'          
+      CALL HBPROF(4,'longit energy profile    (in percent of E inc)'
+     *, NLTOT, 0.,ZMAX, 0., 1000.,' ')
+      CALL HBPROF(5,'rms longit energy profile (in percent of E inc)'
      *, NLTOT, 0.,ZMAX, 0., 1000.,' ')
 *      
@@ -31 +33 @@
       ZMAX = ZMIN + NLTOT*DLX0
                                                                                                                           
-      CALL HBPROF(5,'cumul longit energy dep. (in percent of E inc)'          
+      CALL HBPROF(6,'cumul longit energy dep. (in percent of E inc)'          
      *, NLTOT, ZMIN,ZMAX, 0., 100.,' ')
-      CALL HBOOK1(6,'resolution: cumul L energy dep. (% of E inc)'          
+      CALL HBOOK1(7,'resolution: cumul L energy dep. (% of E inc)'          
      *, NLTOT, ZMIN,ZMAX, 0.0)                                  
 *
@@ -39 +41 @@
       RMAX=NRTOT*DRX0
       
-      CALL HBPROF(7,'radial energy profile    (in percent of E inc)'          
+      CALL HBPROF(8,'radial energy profile    (in percent of E inc)'          
      *, NRTOT, 0.,RMAX, 0., 1000.,' ')
+      CALL HBPROF(9,'rms radial energy profile (in percent of E inc)'          
+     *, NRTOT, 0.,RMAX, 0., 1000.,' ')     
 *
       RMIN = 0.5*DRX0
       RMAX = RMIN + NRTOT*DRX0
       
-      CALL HBPROF(8,'cumul radial energy dep. (in percent of E inc)'          
+      CALL HBPROF(10,'cumul radial energy dep. (in percent of E inc)'          
      *, NRTOT, RMIN,RMAX, 0., 100.,' ')
-      CALL HBOOK1(9,'resolution: cumul R energy dep. (% of E inc)'          
+      CALL HBOOK1(11,'resolution: cumul R energy dep. (% of E inc)'          
      *, NRTOT, RMIN,RMAX, 0.0)