Re: [ROOT] SetKFIN in TPythia6

From: Rene Brun (Rene.Brun@cern.ch)
Date: Mon Sep 13 2004 - 07:39:20 MEST


Hi Suyong,

Thanks for this suggestion. Now implemented in the CVS version.

Rene Brun

On Sun, 12 
Sep 2004, Suyong Choi wrote:

> Hi,
> 
> I just found that SetKFIN doesn't work as I had
> expected. This is because in fortran KFIN(2,-40,40)
> whereas in TPythia6Calls.h, it's
> 
> struct Pysubs_t {
>   int    MSEL;
>   int    MSELPD;
>   int    MSUB[500];
>   int    KFIN[81][2];//
>   double CKIN[200];
> };
> 
>  void        SetKFIN(int i, int j, Int_t kfin=1) {
> fPysubs->KFIN[j-1][i-1] = kfin; }
> 
> so, if one follows the pythia manual and calls SetKFIN
> with negative j, then one ends up doing something
> unintended.
> 
> So, I'd like to suggest that GetKFIN and SetKFIN be
> written
> 
> Int_t GetKFIN(int i, int j) { return
> fPysubs->KFIN[j+40][i-1]; }
> 
> void SetKFIN(int i, int j, Int_t kfin=1) {
> fPysubs->KFIN[j+40][i-1] = kfin; }
> 
> Regards,
> Suyong
> 
> 
> 
> 		
> __________________________________
> Do you Yahoo!?
> Yahoo! Mail Address AutoComplete - You start. We finish.
> http://promotions.yahoo.com/new_mail 
> 



This archive was generated by hypermail 2b29 : Sun Jan 02 2005 - 05:50:09 MET