[Clas_offline] problem with particle decay in gsim

[Zhiwen Zhao]

Dear All,

I am having some problems with particle decay in gsim
For example, if I let K0 decay in gsim, I will have very few Pip and Pim 
reconstructed after cooking, even charged particles reconstructed are very few.
The branching ratio is not even close what it should be.
Some other channels have been tested also, like Lambda1520, Sigma1385, Phi1020.
None of them has correct result after decay in gsim and cooking.
The same cooking code is use to cook production runs for eg3.
So the problem seems related to gsim only.

gsim.ffread
======================================
MLIST
CUTS   5.e-3 5.e-3 5.e-3 5.e-3 5.e-3
CCCUTS 1.e-3 1.e-3 1.e-3 1.e-3 1.e-3
DCCUTS 1.e-4 1.e-4 1.e-4 1.e-4 1.e-4
ECCUTS 1.e-4 1.e-4 1.e-4 1.e-4 1.e-4
SCCUTS 1.e-4 1.e-4 1.e-4 1.e-4 1.e-4
MAGTYPE 2
MAGSCALE -0.5000 0
GEOM 'ALL' 'ST'
NOGEOM 'PTG ' 'MINI'
NOMCDATA 'ALL'
NOSEC 'TORU' 'MINI'
TARGET 'g11a'
TGMATE 'DEUT'
STTYPE 1
STZOFF -41.
TGPOS 0. 0. -9.
SCHDLEN +10.81
SIGBEAM 0.5
ATLEN=376.
BEAM 5.776
RUNG 10
STOP
END
======================================

One MCTK bank generated by fsgen for reachtion gamma + d -> K0 + Lambda1520
and Lambda1520 is allowed to decay to P K- in fsgen and K0 should decay in gsim
==========================================================================
Group:  MCTK    Sector: 0       Nhits:  6       Next ind:       0
cx: 0.000  cy: 0.000  cz: 1.000  pmom: 2.608  mass: 0.000  charge: 0.000  id: 
22  flag:   1  beg_vtx:   0
   end_vtx:   1  parent:   0
cx: 0.626  cy: -0.570  cz: -0.533  pmom: 0.121  mass: 0.938  charge: 1.000  id: 
2212  flag:   1  beg_vtx:
   2  end_vtx:   0  parent:   0
cx: -0.061  cy: 0.301  cz: 0.952  pmom: 1.524  mass: 1.533  charge: 0.000  id: 
53  flag:  11  beg_vtx:
3  end_vtx:   1  parent:   0
cx: 0.014  cy: -0.303  cz: 0.953  pmom: 1.283  mass: 0.498  charge: 0.000  id: 
311  flag:   1  beg_vtx:
4  end_vtx:   0  parent:   0
cx: 0.076  cy: 0.526  cz: 0.847  pmom: 0.849  mass: 0.938  charge: 1.000  id: 
2212  flag:   1  beg_vtx:
5  end_vtx:   0  parent:   3
cx: -0.211  cy: 0.015  cz: 0.977  pmom: 0.748  mass: 0.494  charge: -1.000  id: 
-321  flag:   1  beg_vtx:
   6  end_vtx:   0  parent:   3
==========================================================================

Any suggestion is greatly appreciated.


Zhiwen Zhao