[Rivet] Final particles with negative squared mass

Paolo Francavilla paolo.francavilla at pi.infn.it
Tue Sep 8 14:36:48 BST 2009


Hi,
Thanks for the help.
I am using Rivet, but I found that in the projection for the final 
state, from time to time I have some negative mass2(). What does it means?
Usually the value is of the order of 10e-9.
I am using Pythia6.421

Thanks
with best regards,
Paolo Francavilla

Andy Buckley wrote:
> Hi Paolo,
>
> Sorry about the delay: I meant to look into this and got side-tracked. 
> We assume in the bootstrap script that the GSL and its headers are 
> installed in a system-wide location and so no special flag needs to be 
> passed to the Rivet build system.
>
> I just had a look at the configure script in Rivet, and I think there 
> is a bug in the GSL detection which I think means that the results 
> from running gsl-config don't get propagated correctly through the 
> build. In most cases this isn't a problem because the GSL is already 
> found in a standard place, but I'll have to fix it for you.
>
> Sorry for the inconvenience: in the meantime I suggest building Rivet 
> by hand from the point where the bootstrap fails: that way you can 
> pass the GSL location explicitly to the configure script.
>
> Best wishes,
> Andy
>
>
>
> Paolo Francavilla wrote:
>> Dear all,
>> I have a problem installing Rivet1.1.3 on a machine at CERN (with afs)
>> The problem seems to be related to the headers of the GSL Library.
>> Could you help me?
>>
>> This is what I am doing:
>> export  
>> LD_LIBRARY_PATH=/afs/cern.ch/sw/lcg/external/GSL/1.10/slc4_amd64_gcc34/lib:$LD_LIBRARY_PATH 
>>
>> export  
>> PATH=/afs/cern.ch/sw/lcg/external/GSL/1.10/slc4_amd64_gcc34/bin:$PATH
>> ./rivet-bootstrap --prefix=$PWD/istall
>>
>> the error is:
>>
>> In file included from ../../include/Rivet/Math/Matrices.hh:8,
>>
>>                  from ../../include/Rivet/Rivet.hh:78,
>>
>>                  from ChargedFinalState.cc:2:
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:7:28: gsl/gsl_vector.h: No 
>> such file or directory
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:8:28: gsl/gsl_matrix.h: No 
>> such file or directory
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:9:27: gsl/gsl_eigen.h: No such 
>> file or directory
>>
>> In file included from ../../include/Rivet/Math/Matrices.hh:8,
>>
>>                  from ../../include/Rivet/Rivet.hh:78,
>>
>>                  from ChargedFinalState.cc:2:
>>
>> ../../include/Rivet/Math/MatrixDiag.hh: In function 
>> `Rivet::EigenSystem<N> Rivet::diagonalize(const Rivet::Matrix<N>&)':
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: `gsl_matrix' was 
>> not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: `A' was not 
>> declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: there are no 
>> arguments to `gsl_matrix_alloc' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: (if you use 
>> `-fpermissive', G++ will accept your code, but allowing the use of an 
>> undeclared name is deprecated)
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:124: error: `vecs' was not 
>> declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:124: error: there are no 
>> arguments to `gsl_matrix_alloc' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:125: error: `gsl_vector' was 
>> not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:125: error: `vals' was not 
>> declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:125: error: there are no 
>> arguments to `gsl_vector_alloc' that depend on a template parameter, 
>> so a declaration of `gsl_vector_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:126: error: 
>> `gsl_eigen_symmv_workspace' was not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:126: error: `workspace' was 
>> not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:126: error: there are no 
>> arguments to `gsl_eigen_symmv_alloc' that depend on a template 
>> parameter, so a declaration of `gsl_eigen_symmv_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:127: error: there are no 
>> arguments to `gsl_eigen_symmv' that depend on a template parameter, 
>> so a declaration of `gsl_eigen_symmv' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:128: error: 
>> `GSL_EIGEN_SORT_VAL_DESC' was not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:128: error: there are no 
>> arguments to `gsl_eigen_symmv_sort' that depend on a template 
>> parameter, so a declaration of `gsl_eigen_symmv_sort' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:134: error: there are no 
>> arguments to `gsl_vector_get' that depend on a template parameter, so 
>> a declaration of `gsl_vector_get' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:137: error: there are no 
>> arguments to `gsl_matrix_get' that depend on a template parameter, so 
>> a declaration of `gsl_matrix_get' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:144: error: there are no 
>> arguments to `gsl_eigen_symmv_free' that depend on a template 
>> parameter, so a declaration of `gsl_eigen_symmv_free' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:145: error: there are no 
>> arguments to `gsl_matrix_free' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_free' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:146: error: there are no 
>> arguments to `gsl_matrix_free' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_free' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:147: error: there are no 
>> arguments to `gsl_vector_free' that depend on a template parameter, 
>> so a declaration of `gsl_vector_free' must be available
>>
>> In file included from ../../include/Rivet/Math/Matrices.hh:8,
>>
>>                  from ../../include/Rivet/Rivet.hh:78,
>>
>>                  from Beam.cc:2:
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:7:28: gsl/gsl_vector.h: No 
>> such file or directory
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:8:28: gsl/gsl_matrix.h: No 
>> such file or directory
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:9:27: gsl/gsl_eigen.h: No such 
>> file or directory
>>
>> In file included from ../../include/Rivet/Math/Matrices.hh:8,
>>
>>                  from ../../include/Rivet/Rivet.hh:78,
>>
>>                  from Beam.cc:2:
>>
>> ../../include/Rivet/Math/MatrixDiag.hh: In function 
>> `Rivet::EigenSystem<N> Rivet::diagonalize(const Rivet::Matrix<N>&)':
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: `gsl_matrix' was 
>> not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: `A' was not 
>> declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: there are no 
>> arguments to `gsl_matrix_alloc' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:116: error: (if you use 
>> `-fpermissive', G++ will accept your code, but allowing the use of an 
>> undeclared name is deprecated)
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:124: error: `vecs' was not 
>> declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:124: error: there are no 
>> arguments to `gsl_matrix_alloc' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:125: error: `gsl_vector' was 
>> not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:125: error: `vals' was not 
>> declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:125: error: there are no 
>> arguments to `gsl_vector_alloc' that depend on a template parameter, 
>> so a declaration of `gsl_vector_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:126: error: 
>> `gsl_eigen_symmv_workspace' was not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:126: error: `workspace' was 
>> not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:126: error: there are no 
>> arguments to `gsl_eigen_symmv_alloc' that depend on a template 
>> parameter, so a declaration of `gsl_eigen_symmv_alloc' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:127: error: there are no 
>> arguments to `gsl_eigen_symmv' that depend on a template parameter, 
>> so a declaration of `gsl_eigen_symmv' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:128: error: 
>> `GSL_EIGEN_SORT_VAL_DESC' was not declared in this scope
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:128: error: there are no 
>> arguments to `gsl_eigen_symmv_sort' that depend on a template 
>> parameter, so a declaration of `gsl_eigen_symmv_sort' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:134: error: there are no 
>> arguments to `gsl_vector_get' that depend on a template parameter, so 
>> a declaration of `gsl_vector_get' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:137: error: there are no 
>> arguments to `gsl_matrix_get' that depend on a template parameter, so 
>> a declaration of `gsl_matrix_get' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:144: error: there are no 
>> arguments to `gsl_eigen_symmv_free' that depend on a template 
>> parameter, so a declaration of `gsl_eigen_symmv_free' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:145: error: there are no 
>> arguments to `gsl_matrix_free' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_free' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:146: error: there are no 
>> arguments to `gsl_matrix_free' that depend on a template parameter, 
>> so a declaration of `gsl_matrix_free' must be available
>>
>> ../../include/Rivet/Math/MatrixDiag.hh:147: error: there are no 
>> arguments to `gsl_vector_free' that depend on a template parameter, 
>> so a declaration of `gsl_vector_free' must be available
>>
>> make[2]: *** [Beam.lo] Error 1
>>
>> make[2]: *** Waiting for unfinished jobs....
>>
>> make[2]: *** [ChargedFinalState.lo] Error 1
>>
>> make[2]: Leaving directory 
>> `/media/usbdisk/francav/swt/Rivet/build/Rivet-1.1.3/src/Projections'
>>
>> make[1]: *** [all-recursive] Error 1
>>
>> make[1]: Leaving directory 
>> `/media/usbdisk/francav/swt/Rivet/build/Rivet-1.1.3/src'
>>
>> make: *** [all-recursive] Error 1
>>
>> Traceback (most recent call last):
>>
>>   File "./rivet-bootstrap", line 425, in ?
>>
>>     conf_mk_mkinst(os.path.join(BUILDDIR, "rivet"), 
>> RIVET_CONFIGURE_FLAGS)
>>
>>   File "./rivet-bootstrap", line 181, in conf_mk_mkinst
>>
>>     sys.exit(1)
>>
>> SystemExit: 1
>>
>>
>> Thanks for your help.
>> Paolo Francavilla
>>
>


-- 
________________________________________
Paolo Francavilla
INFN Sez. di Pisa
Edificio C Polo Fibonacci, 
Largo B. Pontecorvo 3 56100 Pisa

Tel.: +390502214414
e-mail: paolo.francavilla at pi.infn.it   
	paolo.francavilla at cern.ch 




More information about the Rivet mailing list