[Clas_offline] problem with particle decay in gsim

[Maurik Holtrop]

Zhiwen,

I am not sure I understand correctly. Are you stating that if you do not let particles decay, i.e. you send gsim banks with just pions, electrons and protons, you DO get the expected results?

Are you setting up the cooking to cook gsim output?

Best,
	Maurik


On Mar 29, 2010, at 12:08 PM, Zhiwen Zhao wrote:

> 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
> _______________________________________________
> Clas_offline mailing list
> Clas_offline at jlab.org
> https://mailman.jlab.org/mailman/listinfo/clas_offline