SCORE and energy conservation

From: maestro (paolo.maestro@pi.infn.it)
Date: Mon Jul 16 2007 - 18:46:19 CEST

  • Next message: Joseph Comfort: "Particle branching ratios"

    Hi all,
    I am using the SCORE and EVENDAT cards to
    score the total and electromagnetic energy deposited in a calorimeter.
    My geometry consists of 70 regions and describes
    a sampling calorimeter made of scintillating plates and tungsten
    and preceded by a carbon pre-shower target.
    I attached the xxxEVT.SCOR.txt file, I obtained by running the Eventdat.f
    routine on the output binary file xxxEVT.SCOR, produced by simulating
    100 GeV protons impinging the calorimeter.
    I noticed that in some events (for istance Ncase=2) the sum of the
    energy deposits
    (quantity scored 208) in each region exceeds the beam particle energy,
    that is the total
    scored energy is > 100 GeV. Instead, the energy balance, I mean the sum
    of all the Endist variables,
    is equal to 100 GeV.
    Maybe there is something I have not understood properly, since I
    thought that summing all the scores
    should give the beam particle energy. Can you explain this ?
    Thanks
             Paolo


             Cream Monte Carlo
     DATE: 7/16/ 7, TIME: 13:57:44
    Number of regions: 70 Number of scored quantities: 2
    The scored quantities are: 208 211

    NCASE, WEIPRU, ENETOT: 1 1.000 100.0
    ENDIST:
        1 0.3187 2 0.3505
        3 5.5645E-03 4 0.000
        5 99.25 6 0.000
        7 0.000 8 3.2074E-03
        9 0.000 10 0.000
       11 0.000 12 7.2369E-02
    Quantity n. 1 (208):
     Scoring per region:
          1 99.25
          2 0.000
          3 3.0407E-04
          4 1.0989E-04
          5 3.2804E-04
          6 2.5839E-04
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 4.2160E-02
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 3.7514E-02
         21 2.0740E-04
         22 3.9773E-04
         23 7.4964E-04
         24 5.5082E-04
         25 8.3320E-03
         26 1.1679E-04
         27 7.4096E-03
         28 9.4905E-05
         29 9.5583E-03
         30 6.8891E-05
         31 1.0835E-02
         32 2.2770E-03
         33 1.5522E-02
         34 6.8151E-05
         35 7.7159E-03
         36 8.5441E-05
         37 8.0855E-03
         38 7.6108E-05
         39 3.0912E-02
         40 8.1245E-04
         41 1.5576E-02
         42 7.7392E-05
         43 1.3169E-02
         44 1.3763E-04
         45 1.1894E-02
         46 5.4802E-04
         47 9.0640E-03
         48 7.2751E-05
         49 1.1391E-02
         50 6.6458E-05
         51 1.2986E-02
         52 7.7892E-05
         53 2.4095E-02
         54 1.3664E-04
         55 1.3461E-02
         56 1.1325E-03
         57 5.6466E-02
         58 6.1753E-04
         59 2.4365E-02
         60 4.7813E-04
         61 5.8515E-02
         62 2.6660E-04
         63 8.5070E-02
         64 3.4489E-04
         65 9.4300E-03
         66 0.1655
         67 0.000
         68 0.000
         69 0.000
         70 0.000
    Quantity n. 2 (211):
     Scoring per region:
          1 1.475
          2 0.000
          3 0.000
          4 0.000
          5 0.000
          6 1.2082E-04
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 1.1439E-02
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 7.1020E-03
         21 0.000
         22 0.000
         23 3.9038E-04
         24 3.3353E-04
         25 1.7411E-03
         26 0.000
         27 1.0166E-03
         28 0.000
         29 2.7401E-03
         30 0.000
         31 4.4681E-03
         32 7.3909E-04
         33 8.8102E-03
         34 0.000
         35 1.0219E-03
         36 0.000
         37 1.5451E-03
         38 0.000
         39 2.4438E-02
         40 7.3035E-04
         41 8.6394E-03
         42 0.000
         43 6.1278E-03
         44 0.000
         45 9.1399E-04
         46 4.7165E-04
         47 2.3534E-03
         48 0.000
         49 4.9805E-03
         50 0.000
         51 6.2883E-03
         52 0.000
         53 1.7207E-02
         54 0.000
         55 6.6730E-03
         56 1.0394E-03
         57 4.9677E-02
         58 5.2665E-04
         59 1.7617E-02
         60 3.9978E-04
         61 3.4695E-02
         62 2.0104E-04
         63 7.8000E-02
         64 2.6874E-04
         65 5.9498E-03
         66 4.1860E-02
         67 0.000
         68 0.000
         69 0.000
         70 0.000
     838774 0 0. 0. 0. 0.

    NCASE, WEIPRU, ENETOT: 2 1.000 100.0
    ENDIST:
        1 10.19 2 47.75
        3 0.1152 4 0.000
        5 36.65 6 1.145
        7 0.000 8 0.3071
        9 0.000 10 0.000
       11 0.000 12 3.839
    Quantity n. 1 (208):
     Scoring per region:
          1 36.65
          2 0.000
          3 5.7585E-04
          4 1.4132E-04
          5 6.2736E-04
          6 1.0852E-04
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 0.1859
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 0.6604
         21 1.3986E-02
         22 3.6599E-02
         23 2.6425E-02
         24 1.7098E-02
         25 2.265
         26 4.6511E-02
         27 1.691
         28 2.8882E-02
         29 7.008
         30 3.1685E-02
         31 5.307
         32 4.5538E-02
         33 7.090
         34 4.6543E-02
         35 7.704
         36 3.8386E-02
         37 3.612
         38 2.111
         39 4.786
         40 4.1907E-02
         41 3.171
         42 4.2814E-02
         43 4.135
         44 2.6847E-02
         45 5.439
         46 3.1408E-02
         47 5.811
         48 2.5929E-02
         49 5.956
         50 2.8457E-02
         51 6.169
         52 4.7328E-02
         53 2.868
         54 3.5282E-02
         55 12.46
         56 0.9592
         57 14.00
         58 3.7665E-02
         59 8.749
         60 3.6232E-02
         61 5.347
         62 2.5944E-02
         63 3.918
         64 2.0768E-02
         65 0.2543
         66 4.285
         67 4.3594E-04
         68 7.5659E-04
         69 2.2790E-03
         70 4.7840E-05
    Quantity n. 2 (211):
     Scoring per region:
          1 4.064
          2 0.000
          3 1.7592E-04
          4 3.2545E-05
          5 2.4096E-04
          6 6.8245E-07
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 8.7181E-02
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 0.2316
         21 1.8903E-03
         22 6.7434E-03
         23 8.9887E-03
         24 7.3170E-03
         25 0.5296
         26 1.1512E-02
         27 0.9188
         28 1.8408E-02
         29 1.742
         30 2.1967E-02
         31 2.316
         32 3.0163E-02
         33 2.752
         34 3.3572E-02
         35 3.045
         36 3.2390E-02
         37 3.282
         38 3.5963E-02
         39 3.147
         40 3.7377E-02
         41 2.893
         42 3.1214E-02
         43 2.701
         44 2.1443E-02
         45 2.242
         46 2.5679E-02
         47 2.209
         48 2.1236E-02
         49 2.234
         50 2.2955E-02
         51 2.216
         52 2.0499E-02
         53 2.115
         54 2.3037E-02
         55 2.104
         56 2.1768E-02
         57 2.267
         58 1.9505E-02
         59 2.312
         60 2.1239E-02
         61 2.121
         62 1.6467E-02
         63 1.613
         64 1.2243E-02
         65 0.1712
         66 1.999
         67 4.3534E-04
         68 7.5360E-04
         69 4.4302E-04
         70 4.6766E-05
     110984797 0 0. 0. 0. 0.

    NCASE, WEIPRU, ENETOT: 3 1.000 100.0
    ENDIST:
        1 11.65 2 50.36
        3 0.1612 4 0.000
        5 32.83 6 0.6628
        7 0.000 8 0.3726
        9 0.000 10 0.000
       11 0.000 12 3.964
    Quantity n. 1 (208):
     Scoring per region:
          1 32.83
          2 0.000
          3 5.2931E-03
          4 2.1252E-03
          5 7.5216E-03
          6 1.799
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 1.121
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 2.688
         21 6.1996E-02
         22 6.1664E-02
         23 8.6381E-02
         24 3.2272E-02
         25 3.135
         26 2.5389E-02
         27 3.358
         28 3.7889E-02
         29 7.445
         30 3.7769E-02
         31 4.362
         32 4.9459E-02
         33 5.468
         34 7.5438E-02
         35 11.96
         36 6.2592E-02
         37 8.117
         38 5.7400E-02
         39 6.183
         40 7.0193E-02
         41 11.15
         42 6.1975E-02
         43 8.662
         44 4.9206E-02
         45 11.64
         46 1.121
         47 3.260
         48 3.2448E-02
         49 3.937
         50 2.3011E-02
         51 2.385
         52 2.3289E-02
         53 4.813
         54 2.9355E-02
         55 1.629
         56 2.1607E-02
         57 5.343
         58 1.5323E-02
         59 3.947
         60 1.3569E-02
         61 2.045
         62 2.1816E-02
         63 1.749
         64 1.2195E-02
         65 0.1288
         66 1.846
         67 2.9360E-04
         68 3.1983E-05
         69 8.9359E-02
         70 9.6389E-07
    Quantity n. 2 (211):
     Scoring per region:
          1 1.683
          2 0.000
          3 3.0789E-04
          4 4.1809E-04
          5 2.7095E-03
          6 3.8835E-04
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 0.4020
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 0.8680
         21 8.3273E-03
         22 1.5042E-02
         23 2.2933E-02
         24 1.1461E-02
         25 0.7326
         26 8.7970E-03
         27 1.483
         28 2.0950E-02
         29 2.113
         30 2.5813E-02
         31 2.923
         32 3.7046E-02
         33 3.799
         34 4.8123E-02
         35 4.371
         36 5.1205E-02
         37 4.775
         38 4.6054E-02
         39 4.656
         40 5.0476E-02
         41 4.719
         42 4.5670E-02
         43 3.667
         44 3.8388E-02
         45 3.102
         46 2.8683E-02
         47 2.470
         48 2.1351E-02
         49 1.986
         50 1.6283E-02
         51 1.591
         52 1.6889E-02
         53 1.287
         54 9.2059E-03
         55 0.7751
         56 1.0071E-02
         57 0.8665
         58 6.1078E-03
         59 0.7276
         60 6.7784E-03
         61 0.8296
         62 9.9595E-03
         63 0.7843
         64 6.0556E-03
         65 7.3719E-02
         66 0.7908
         67 2.9259E-04
         68 0.000
         69 3.5670E-03
         70 0.000
     227896448 0 0. 0. 0. 0.

    NCASE, WEIPRU, ENETOT: 4 1.000 100.0
    ENDIST:
        1 0.2207 2 9.7711E-02
        3 5.3415E-06 4 0.000
        5 99.68 6 0.000
        7 0.000 8 0.000
        9 0.000 10 0.000
       11 0.000 12 -3.4674E-12
    Quantity n. 1 (208):
     Scoring per region:
          1 99.68
          2 0.000
          3 3.1749E-04
          4 1.5069E-04
          5 3.1246E-04
          6 1.4335E-04
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 4.2972E-02
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 6.7489E-02
         21 8.6736E-04
         22 1.0493E-03
         23 1.3409E-03
         24 8.7531E-04
         25 1.6688E-02
         26 7.7991E-05
         27 9.2980E-03
         28 7.7888E-05
         29 7.3109E-03
         30 7.6047E-05
         31 1.0001E-02
         32 1.9739E-03
         33 8.0407E-03
         34 7.7023E-05
         35 7.9952E-03
         36 1.0386E-04
         37 7.3361E-03
         38 1.4341E-04
         39 7.7219E-03
         40 7.7055E-05
         41 6.8106E-03
         42 7.3536E-05
         43 7.9812E-03
         44 6.7479E-05
         45 9.0389E-03
         46 1.4041E-04
         47 1.3386E-02
         48 7.2736E-05
         49 8.4895E-03
         50 8.3796E-05
         51 7.8622E-03
         52 7.6832E-05
         53 7.9792E-03
         54 1.2075E-04
         55 9.3883E-03
         56 1.7500E-04
         57 7.2401E-03
         58 1.1968E-04
         59 1.0525E-02
         60 6.6645E-05
         61 8.0093E-03
         62 1.6385E-04
         63 7.8860E-03
         64 6.4816E-05
         65 1.2720E-03
         66 2.4178E-02
         67 0.000
         68 0.000
         69 0.000
         70 0.000
    Quantity n. 2 (211):
     Scoring per region:
          1 0.000
          2 0.000
          3 0.000
          4 0.000
          5 0.000
          6 0.000
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 1.1775E-02
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 3.6779E-02
         21 5.8920E-04
         22 7.3005E-04
         23 9.6846E-04
         24 6.1086E-04
         25 9.8802E-03
         26 0.000
         27 2.4225E-03
         28 0.000
         29 4.9819E-04
         30 0.000
         31 3.0972E-03
         32 1.8983E-03
         33 1.4068E-03
         34 0.000
         35 1.2809E-03
         36 0.000
         37 4.1165E-04
         38 0.000
         39 1.0847E-03
         40 0.000
         41 2.8996E-04
         42 0.000
         43 1.1652E-03
         44 0.000
         45 2.3024E-03
         46 0.000
         47 1.3129E-03
         48 0.000
         49 2.0214E-03
         50 0.000
         51 1.1965E-03
         52 0.000
         53 1.4207E-03
         54 0.000
         55 2.4020E-03
         56 0.000
         57 7.8652E-04
         58 0.000
         59 3.6473E-03
         60 0.000
         61 1.1287E-03
         62 0.000
         63 1.2174E-03
         64 0.000
         65 0.000
         66 5.3863E-03
         67 0.000
         68 0.000
         69 0.000
         70 0.000
     228122554 0 0. 0. 0. 0.

    NCASE, WEIPRU, ENETOT: 5 1.000 100.0
    ENDIST:
        1 0.2204 2 7.1961E-02
        3 1.5010E-03 4 0.000
        5 99.71 6 0.000
        7 0.000 8 0.000
        9 0.000 10 0.000
       11 0.000 12 -2.3874E-12
    Quantity n. 1 (208):
     Scoring per region:
          1 99.71
          2 0.000
          3 2.7347E-04
          4 1.2552E-04
          5 3.2813E-04
          6 1.0892E-04
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 4.5230E-02
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 5.0251E-02
         21 6.7499E-04
         22 3.1855E-04
         23 4.6554E-04
         24 2.3980E-04
         25 8.2915E-03
         26 7.5524E-05
         27 8.6313E-03
         28 7.1596E-05
         29 7.0957E-03
         30 7.4343E-05
         31 7.2953E-03
         32 7.9026E-05
         33 7.7854E-03
         34 6.9678E-05
         35 8.6358E-03
         36 1.7555E-04
         37 7.9896E-03
         38 6.5576E-05
         39 7.2686E-03
         40 1.0456E-04
         41 7.5676E-03
         42 6.9101E-05
         43 7.1405E-03
         44 1.1465E-04
         45 8.6676E-03
         46 7.4218E-05
         47 1.0499E-02
         48 9.4636E-05
         49 7.6924E-03
         50 6.7763E-05
         51 6.8346E-03
         52 1.0927E-04
         53 9.5359E-03
         54 8.3071E-05
         55 1.2886E-02
         56 6.9841E-05
         57 8.2160E-03
         58 6.4256E-05
         59 7.7121E-03
         60 8.8284E-05
         61 9.1505E-03
         62 8.0693E-05
         63 9.7925E-03
         64 1.7773E-04
         65 1.2749E-03
         66 2.4046E-02
         67 0.000
         68 0.000
         69 0.000
         70 0.000
    Quantity n. 2 (211):
     Scoring per region:
          1 7.2889E-04
          2 0.000
          3 0.000
          4 0.000
          5 0.000
          6 0.000
          7 0.000
          8 0.000
          9 0.000
         10 0.000
         11 0.000
         12 0.000
         13 0.000
         14 1.4004E-02
         15 0.000
         16 0.000
         17 0.000
         18 0.000
         19 0.000
         20 1.7602E-02
         21 4.6024E-04
         22 0.000
         23 1.0301E-04
         24 0.000
         25 1.4840E-03
         26 0.000
         27 1.7210E-03
         28 0.000
         29 5.5017E-04
         30 0.000
         31 9.2747E-04
         32 0.000
         33 9.0390E-04
         34 0.000
         35 2.0470E-03
         36 1.0978E-04
         37 1.1098E-03
         38 0.000
         39 8.9333E-04
         40 0.000
         41 6.9699E-04
         42 0.000
         43 6.6667E-04
         44 0.000
         45 1.9205E-03
         46 0.000
         47 3.3061E-03
         48 0.000
         49 8.6386E-04
         50 0.000
         51 1.5019E-04
         52 0.000
         53 2.8008E-03
         54 0.000
         55 6.3359E-03
         56 0.000
         57 1.2196E-03
         58 0.000
         59 9.1465E-04
         60 0.000
         61 2.7407E-03
         62 0.000
         63 2.8699E-03
         64 1.0918E-04
         65 1.8946E-04
         66 5.2604E-03
         67 0.000
         68 0.000
         69 0.000
         70 0.000
     228305941 0 0. 0. 0. 0.

     End of a run of 5 particles


  • Next message: Joseph Comfort: "Particle branching ratios"

    This archive was generated by hypermail 2.1.6 : Mon Jul 16 2007 - 18:11:02 CEST