R__LOAD_LIBRARY ($PLUTOLIBDIR/libPluto.so)

int dielectron_sources() {
  PFireball fireOmega ("w"  , 2.0, 0.100, 0., 1., 0., 0., 0., 0., 0. );
  PFireball firePhi   ("phi", 2.0, 0.100, 0., 1., 0., 0., 0., 0., 0. );
  fireOmega.SetW (0.7);
  firePhi  .SetW (0.3);
  fireOmega.Print ();
  firePhi  .Print ();

  PParticle omega ("w") , ep2 ("e+") , em2 ("e-") ;
  PParticle* arr_fireOmega_omega[] = { &fireOmega , &omega  };
  PParticle* arr_omega_ee       [] = { &omega , &ep1 , &em1 };
  PChannel * ch_fireOmega_omega    = new PChannel ( arr_fireOmega_omega , 1 );
  PChannel * ch_omega_ee           = new PChannel ( arr_omega_ee        , 2 );

  PParticle phi   ("phi") , ep1 ("e+") , em1 ("e-") ;
  PParticle* arr_firePhi_phi    [] = { &firePhi , &phi   };
  PParticle* arr_phi_ee         [] = { &phi, &ep2 , &em2 };
  PChannel * ch_firePhi_phi        = new PChannel ( arr_firePhi_phi     , 1 );
  PChannel * ch_phi_ee             = new PChannel ( arr_phi_ee          , 2 );

  PChannel *  tab_chans[] = { ch_fireOmega_omega, ch_omega_ee ,
                              ch_firePhi_phi    , ch_phi_ee    };

  PReaction* r_TwoFBs = new PReaction ( tab_chans, "omegaphi_ee", 4, 1, 0, 0, 0);
  r_TwoFBs->Print();
  r_TwoFBs->loop (1000, 1);
  return 0;
}