Ignore:
Timestamp:
Nov 5, 2010, 3:45:55 PM (14 years ago)
Author:
garnier
Message:

update ti head

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/hadronic/models/cascade/cascade/src/G4WatcherGun.cc

    r1337 r1340  
    2323// * acceptance of all terms of the Geant4 Software license.          *
    2424// ********************************************************************
     25// $Id: G4WatcherGun.cc,v 1.15 2010/10/22 20:41:05 mkelsey Exp $
     26// GEANT4 tag: $Name: hadr-casc-V09-03-85 $
    2527//
    26 // $Id: G4WatcherGun.cc,v 1.13.2.1 2010/06/25 09:45:10 gunter Exp $
    27 // GEANT4 tag: $Name: geant4-09-04-beta-01 $
     28// To include numerous high-Z watchers (e.g., for proton-lead collisions)
     29// set the preprocessor flag G4CASCADE_WATCHER_HIGHZ
    2830//
    2931// 20100407  M. Kelsey -- Replace std::vector<>::resize(0) with ::clear(),
    3032//              and create vectors pre-sized to maximum needed.
    31 
    32 //#define PB
     33//
    3334
    3435#include "G4WatcherGun.hh"
     
    6263  cs.push_back(217.4);
    6364  errs.push_back(1.22);
    64   watchers.push_back(G4NuclWatcher(0.0, as, cs, errs, false, false));
     65  watchers.push_back(G4NuclWatcher(0, as, cs, errs, false, false));
    6566
    6667  // Z = 1
     
    8081  cs.push_back(338.6);
    8182  errs.push_back(3.1);
    82   watchers.push_back(G4NuclWatcher(1.0, as, cs, errs, false, false));
     83  watchers.push_back(G4NuclWatcher(1, as, cs, errs, false, false));
    8384
    8485  // Z = 1
     
    9293  cs.push_back(338.6);
    9394  errs.push_back(3.1);
    94   watchers.push_back(G4NuclWatcher(1.0, as, cs, errs, false, true));
     95  watchers.push_back(G4NuclWatcher(1, as, cs, errs, false, true));
    9596 
    9697  // Z = -1
     
    101102  cs.push_back(198.3);
    102103  errs.push_back(1.0);
    103   watchers.push_back(G4NuclWatcher(-1.0, as, cs, errs, false, false));
     104  watchers.push_back(G4NuclWatcher(-1, as, cs, errs, false, false));
    104105
    105106  // Z = 2
     
    113114  cs.push_back(781.0);
    114115  errs.push_back(5.0);
    115   watchers.push_back(G4NuclWatcher(2.0, as, cs, errs, false, true));
    116 
    117 #ifdef PB
     116  watchers.push_back(G4NuclWatcher(2, as, cs, errs, false, true));
     117
     118#ifdef G4CASCADE_WATCHER_HIGHZ
    118119  // Z = 22
    119120  // watchers for pb208 + 1 GeV p
     
    139140  cs.push_back(0.06);
    140141  errs.push_back(0.023);
    141   watchers.push_back(G4NuclWatcher(22.0, as, cs, errs, false, true));
     142  watchers.push_back(G4NuclWatcher(22, as, cs, errs, false, true));
    142143
    143144  // Z = 23
     
    166167  cs.push_back(0.05);
    167168  errs.push_back(0.012);
    168   watchers.push_back(G4NuclWatcher(23.0, as, cs, errs, false, true));
     169  watchers.push_back(G4NuclWatcher(23, as, cs, errs, false, true));
    169170
    170171  // Z = 24
     
    193194  cs.push_back(0.092);
    194195  errs.push_back(0.011);
    195   watchers.push_back(G4NuclWatcher(24.0, as, cs, errs, true, true));
     196  watchers.push_back(G4NuclWatcher(24, as, cs, errs, true, true));
    196197
    197198  // Z = 25
     
    223224  cs.push_back(0.051);
    224225  errs.push_back(0.015);
    225   watchers.push_back(G4NuclWatcher(25.0, as, cs, errs, false, true));
     226  watchers.push_back(G4NuclWatcher(25, as, cs, errs, false, true));
    226227
    227228  // Z = 26
     
    253254  cs.push_back(0.09);
    254255  errs.push_back(0.014);
    255   watchers.push_back(G4NuclWatcher(26.0, as, cs, errs, false, true));
     256  watchers.push_back(G4NuclWatcher(26, as, cs, errs, false, true));
    256257
    257258  // Z = 27
     
    286287  cs.push_back(0.108);
    287288  errs.push_back(0.046);
    288   watchers.push_back(G4NuclWatcher(27.0, as, cs, errs, false, true));
     289  watchers.push_back(G4NuclWatcher(27, as, cs, errs, false, true));
    289290
    290291  // Z = 28
     
    322323  cs.push_back(0.054);
    323324  errs.push_back(0.036);
    324   watchers.push_back(G4NuclWatcher(28.0, as, cs, errs, true, true));
     325  watchers.push_back(G4NuclWatcher(28, as, cs, errs, true, true));
    325326
    326327  // Z = 29
     
    361362  cs.push_back(0.041);
    362363  errs.push_back(0.024);
    363   watchers.push_back(G4NuclWatcher(29.0, as, cs, errs, false, true));
     364  watchers.push_back(G4NuclWatcher(29, as, cs, errs, false, true));
    364365
    365366  // Z = 30
     
    397398  cs.push_back(0.082);
    398399  errs.push_back(0.034);
    399   watchers.push_back(G4NuclWatcher(30.0, as, cs, errs, false, true));
     400  watchers.push_back(G4NuclWatcher(30, as, cs, errs, false, true));
    400401
    401402  // Z = 31
     
    436437  cs.push_back(0.061);
    437438  errs.push_back(0.014);
    438   watchers.push_back(G4NuclWatcher(31.0, as, cs, errs, true, true));
     439  watchers.push_back(G4NuclWatcher(31, as, cs, errs, true, true));
    439440
    440441  // Z = 32
     
    478479  cs.push_back(0.04);
    479480  errs.push_back(0.018);
    480   watchers.push_back(G4NuclWatcher(32.0, as, cs, errs, false, true));
     481  watchers.push_back(G4NuclWatcher(32, as, cs, errs, false, true));
    481482
    482483  // Z = 33
     
    520521  cs.push_back(0.023);
    521522  errs.push_back(0.01);
    522   watchers.push_back(G4NuclWatcher(33.0, as, cs, errs, false, true));
     523  watchers.push_back(G4NuclWatcher(33, as, cs, errs, false, true));
    523524
    524525  // Z = 34
     
    562563  cs.push_back(0.045);
    563564  errs.push_back(0.012);
    564   watchers.push_back(G4NuclWatcher(34.0, as, cs, errs, false, true));
     565  watchers.push_back(G4NuclWatcher(34, as, cs, errs, false, true));
    565566
    566567  // Z = 35
     
    607608  cs.push_back(0.033);
    608609  errs.push_back(0.02);
    609   watchers.push_back(G4NuclWatcher(35.0, as, cs, errs, false, true));
     610  watchers.push_back(G4NuclWatcher(35, as, cs, errs, false, true));
    610611
    611612  // spallation part
     
    656657  cs.push_back(0.023);
    657658  errs.push_back(0.016);
    658   watchers.push_back(G4NuclWatcher(61.0, as, cs, errs, true, true));
     659  watchers.push_back(G4NuclWatcher(61, as, cs, errs, true, true));
    659660
    660661  // Z = 69
     
    716717  cs.push_back(0.014);
    717718  errs.push_back(0.01);
    718   watchers.push_back(G4NuclWatcher(69.0, as, cs, errs, true, true));
     719  watchers.push_back(G4NuclWatcher(69, as, cs, errs, true, true));
    719720
    720721  // Z = 73
     
    791792  cs.push_back(0.016);
    792793  errs.push_back(0.009);
    793   watchers.push_back(G4NuclWatcher(73.0, as, cs, errs, true, true));
     794  watchers.push_back(G4NuclWatcher(73, as, cs, errs, true, true));
    794795
    795796  // Z = 77
     
    878879  cs.push_back(0.037);
    879880  errs.push_back(0.019);
    880   watchers.push_back(G4NuclWatcher(77.0, as, cs, errs, true, true));
     881  watchers.push_back(G4NuclWatcher(77, as, cs, errs, true, true));
    881882
    882883  // Z = 81
     
    956957  cs.push_back(17.525);
    957958  errs.push_back(0.526);
    958   watchers.push_back(G4NuclWatcher(81.0, as, cs, errs, true, true));
     959  watchers.push_back(G4NuclWatcher(81, as, cs, errs, true, true));
    959960
    960961  // Z = 82
     
    10251026  cs.push_back(63.653);
    10261027  errs.push_back(9.573);
    1027   watchers.push_back(G4NuclWatcher(82.0, as, cs, errs, true, true));
     1028  watchers.push_back(G4NuclWatcher(82, as, cs, errs, true, true));
    10281029
    10291030#else
     
    10881089  cs.push_back(0.003);
    10891090  errs.push_back(0.001);
    1090   watchers.push_back(G4NuclWatcher(80.0, as, cs, errs, true, true));
     1091  watchers.push_back(G4NuclWatcher(80, as, cs, errs, true, true));
    10911092
    10921093  // Z = 77
     
    11631164  cs.push_back(3.54);
    11641165  errs.push_back(0.45);
    1165   watchers.push_back(G4NuclWatcher(77.0, as, cs, errs, true, true));
     1166  watchers.push_back(G4NuclWatcher(77, as, cs, errs, true, true));
    11661167
    11671168  // Z = 73
     
    12321233  cs.push_back(0.01);
    12331234  errs.push_back(0.008);
    1234   watchers.push_back(G4NuclWatcher(73.0, as, cs, errs, true, true));
     1235  watchers.push_back(G4NuclWatcher(73, as, cs, errs, true, true));
    12351236
    12361237  // Z = 61
     
    12741275  cs.push_back(0.02);
    12751276  errs.push_back(0.003);
    1276   watchers.push_back(G4NuclWatcher(61.0, as, cs, errs, true, true));
     1277  watchers.push_back(G4NuclWatcher(61, as, cs, errs, true, true));
    12771278
    12781279  // Z = 79
     
    13431344  cs.push_back(0.36);
    13441345  errs.push_back(0.05);
    1345   watchers.push_back(G4NuclWatcher(79.0, as, cs, errs, true, true));
     1346  watchers.push_back(G4NuclWatcher(79, as, cs, errs, true, true));
    13461347
    13471348  // Z = 78
     
    14151416  cs.push_back(5.49);
    14161417  errs.push_back(0.7);
    1417   watchers.push_back(G4NuclWatcher(78.0, as, cs, errs, true, true));
     1418  watchers.push_back(G4NuclWatcher(78, as, cs, errs, true, true));
    14181419
    14191420  // Z = 72
     
    14751476  cs.push_back(0.14);
    14761477  errs.push_back(0.02);
    1477   watchers.push_back(G4NuclWatcher(72.0, as, cs, errs, true, true));
     1478  watchers.push_back(G4NuclWatcher(72, as, cs, errs, true, true));
    14781479
    14791480  // Z = 66
     
    15141515  cs.push_back(0.55);
    15151516  errs.push_back(0.1);
    1516   watchers.push_back(G4NuclWatcher(66.0, as, cs, errs, true, true));
     1517  watchers.push_back(G4NuclWatcher(66, as, cs, errs, true, true));
    15171518
    15181519  // Z = 65
     
    15511552  errs.push_back(0.1);
    15521553
    1553   watchers.push_back(G4NuclWatcher(65.0, as, cs, errs, true, true));
     1554  watchers.push_back(G4NuclWatcher(65, as, cs, errs, true, true));
    15541555#endif 
    15551556
Note: See TracChangeset for help on using the changeset viewer.