trainbow.py : Rainbow deviation angle comparison between Opticks and Geant4ΒΆ

To simulate the rainbow events:

trainbow-

trainbow-- --spol
trainbow-- --ppol

trainbow-- --spol --tcfg4
trainbow-- --ppol --tcfg4

Expected output, is a scattering angle plot and history comparison table:

In [28]: run rainbow_cfg4.py
WARNING:opticks.ana.evt:init_index S-Pol G4 finds too few (ps)phosel uniques : 1
WARNING:opticks.ana.evt:init_index S-Pol G4 finds too few (rs)recsel uniques : 1
WARNING:opticks.ana.evt:init_index S-Pol G4 finds too few (rsr)reshaped-recsel uniques : 1
                       5:rainbow   -5:rainbow           c2
                8ccd        819160       819654             0.15  [4 ] TO BT BT SA
                 8bd        102089       101615             1.10  [3 ] TO BR SA
               8cbcd         61869        61890             0.00  [5 ] TO BT BR BT SA
              8cbbcd          9617         9577             0.08  [6 ] TO BT BR BR BT SA
             8cbbbcd          2604         2687             1.30  [7 ] TO BT BR BR BR BT SA
            8cbbbbcd          1056         1030             0.32  [8 ] TO BT BR BR BR BR BT SA
               86ccd          1014         1000             0.10  [5 ] TO BT BT SC SA
           8cbbbbbcd           472          516             1.96  [9 ] TO BT BR BR BR BR BR BT SA
                 86d           498          473             0.64  [3 ] TO SC SA
          bbbbbbbbcd           304          294             0.17  [10] TO BT BR BR BR BR BR BR BR BR
          8cbbbbbbcd           272          247             1.20  [10] TO BT BR BR BR BR BR BR BT SA
          cbbbbbbbcd           183          161             1.41  [10] TO BT BR BR BR BR BR BR BR BT
                 4cd           161          139             1.61  [3 ] TO BT AB
                86bd           138          142             0.06  [4 ] TO BR SC SA
               8c6cd           126          106             1.72  [5 ] TO BT SC BT SA
                4ccd           100          117             1.33  [4 ] TO BT BT AB
              86cbcd            88          110             2.44  [6 ] TO BT BR BT SC SA
                  4d            51           54             0.09  [2 ] TO AB
               8cc6d            38           40             0.05  [5 ] TO SC BT BT SA
             8cc6ccd            19           33             3.77  [7 ] TO BT BT SC BT BT SA
                         1000000      1000000         1.09