Warnings during pythia-pgs installation

Asked by Jose David Ruiz Alvarez on 2013-03-19

Dear MG5 team,

I've installed pythia in MG5 using the command "install pythia-pgs". However the installation ends with a satisfactory message I'm having a lot of warnings:

make[2]: Entering directory `/afs/cern.ch/work/j/jruizalv/public/MadGraph5_v1_5_8/pythia-pgs/libraries/PGS4'
test -d lib || mkdir -p lib

----- Now compiling PGSLIB -----

gfortran -fno-automatic -c -o ./lib/pgslib.o ./src/pgslib.f
./src/pgslib.f:2727.17:

      pt1 = v4pt(p1)
                 1
Warning: Rank mismatch in argument 'p' at (1) (rank-1 and scalar)
./src/pgslib.f:2728.17:

      pt2 = v4pt(p2)
                 1
Warning: Rank mismatch in argument 'p' at (1) (rank-1 and scalar)
./src/pgslib.f:2780.38:

          pgs_kt_beam_distance = v4pt(p1)**2
                                      1
Warning: Rank mismatch in argument 'p' at (1) (rank-1 and scalar)
./src/pgslib.f:2784.42:

          pgs_kt_beam_distance = 1.0/v4pt(p1)**2
                                          1
Warning: Rank mismatch in argument 'p' at (1) (rank-1 and scalar)
./src/pgslib.f:5267.33:

      del_phi = v4phi(p) - v4phi(q)
                                 1
Warning: Actual argument contains too few elements for dummy argument 'p' (3/4) at (1)
./src/pgslib.f:5267.22:

      del_phi = v4phi(p) - v4phi(q)
                      1
Warning: Actual argument contains too few elements for dummy argument 'p' (3/4) at (1)
ar -cr ./lib/libpgslib.a ./lib/pgslib.o
ranlib ./lib/libpgslib.a
rm ./lib/pgslib.o

----- Now compiling TAUOLA -----

gfortran -fno-automatic -c -o ./lib/tauola.o ./src/tauola.f
./src/tauola.f:50.27:

               call tdrand(mult,1)
                           1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:732.18:

      call tdrand(mult,1)
                  1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:199.21:

         call tdrand(rnuni,1)
                     1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:3201.20:

        CALL TDRAND(RN,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:3333.20:

        CALL TDRAND(RN,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:3419.20:

        CALL TDRAND(RN,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:4644.18:

      CALL TDRAND(RR1,1)
                  1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:3463.20:

        CALL TDRAND(RN,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:2958.20:

        CALL TDRAND(RN,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:3239.20:

        CALL TDRAND(RN,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:4311.18:

      CALL TDRAND(RR1,1)
                  1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:1657.20:

        CALL TDRAND(RMOD,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:3288.20:

        CALL TDRAND(RN,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:4084.18:

      CALL TDRAND(RR1,1)
                  1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:4172.18:

      CALL TDRAND(RN,1)
                  1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:4537.18:

      CALL TDRAND(RR1,1)
                  1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:3381.22:

          CALL TDRAND(RN,1)
                      1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:5714.18:

      CALL TDRAND(RRR,1)
                  1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:5282.20:

        CALL TDRAND(XIO,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:7161.20:

        CALL TDRAND(XIO,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:7170.20:

        CALL TDRAND(XIO,1)
                    1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:913.21:

         call tdrand(test,1)
                     1
Warning: Rank mismatch in argument 'rvec' at (1) (rank-1 and scalar)
./src/tauola.f:395.20:

      call dexay(-1,dummy)
                    1
Warning: Rank mismatch in argument 'pol' at (1) (rank-1 and scalar)
./src/tauola.f:2611.30:

          CALL DADMPI(-1,IDUM,PDUM,PDUM1,PDUM2)
                              1
Warning: Rank mismatch in argument 'hv' at (1) (rank-1 and scalar)
./src/tauola.f:2614.30:

          CALL DADMKK(-1,IDUM,PDUM,PDUM1,PDUM2)
                              1
Warning: Rank mismatch in argument 'hv' at (1) (rank-1 and scalar)
./src/tauola.f:2653.30:

          CALL DADMPI( 1,IDUM,PDUM,PDUM1,PDUM2)
                              1
Warning: Rank mismatch in argument 'hv' at (1) (rank-1 and scalar)
./src/tauola.f:2656.30:

          CALL DADMKK( 1,IDUM,PDUM,PDUM1,PDUM2)

The problem is after runnning on good lhe sample with the commands:
./bin/madevent
pythia-pgs run_01

I'm obtaining damaged .hep files. I know they are damaaged because when I try to do some checks with MadAnalysis I obtain the following error message:
ERROR: Block with the ID=1 is not managed by SampleAnalyzer
    *******************************************************
   Checking SampleAnalyzer output...
 ** ERROR: File called '/afs/cern.ch/work/j/jruizalv/public/Clean/CMSSW_6_1_0/src/MadAnalysis/madanalysis-current/W_4jets_matching_check_xqcut20/Output/_defaultset.saf' is not found.
 ** ERROR: errors have occured during the analysis.

This is because the warnings?

Thanks,
Jose

P.D.: I've checked with another samples, that I know are good, in MadAnalysis and I don't obtain the error.

Question information

Language:
English Edit question
Status:
Open
For:
MadGraph5_aMC@NLO Edit question
Assignee:
MadAnalysisTeam Edit question
Last query:
2013-04-01
Last reply:
2013-03-29

Hi Jose,

I've those warning on my computer too, and everything runs fine for pythia/pgs on my side.
So I don't think that this is the problem.

Otherwise the command:
/bin/madevent
pythia-pgs run_01

doesn't exists. This is directly crashing on my computer. I suppose that this is a typo, if not the code is expected to crash directly on such command.

Otherwise did you use MadAnalysis 4 or MadAnalysis5? In both case, could you send me by email (<email address hidden>) the (parton) lhe and the hep file?

Thanks,

Olivier

Not sure at 100%.

But in Any case, I have assign this to MA5 team for the moment (and they have access to lxplus for sure)
such that they can at least tell us what this error means.

Cheers,

Olivier

Benjamin Fuks (fuks) said : #3

Hi Jose, Olivier,

This indeed means that the sample is corrupted in some way. Could you please send me the mg5 scripts you are using, with the appropriate cards, to generate such a sample? I would like to try to reproduce the issue.

Cheers,

Benjamin

Hi Benjamin,

I tried what you told me to do. I took the folder for the process I have produced on lxplus and copied to IPNL serv to my madgraph package, then I corrected some stuffs as the path for madgraph that is found in me5_configuration.txt in the Cards directory and additionally I recompiled all the Source directory as pythia was complaining about some compilation issues of this directory. After that I ran pythia over the parton level sample, and it worked without any trouble. But then I go to MadAnalysis and I obtain the following errors

ma5>import /afs/cern.ch/user/j/jruizalv/public/TTbar_pythia_events.hep.gz
   -> Storing the file 'TTbar_pythia_events.hep.gz' in the dataset 'defaultset'.
ma5>set main.matching.check = true
ma5>set main.matching.njets = 6
ma5>submit TTbar_matching_check_xqcut20
   Creating folder 'TTbar_matching_check_xqcut20'...
   Copying 'SampleAnalyzer' source files...
   Inserting your selection into 'SampleAnalyzer'...
   Writing the list of datasets...
   Writing the command line history...
   Creating Makefiles...
   Compiling 'SampleAnalyzer'...
   Linking 'SampleAnalyzer'...
   Running 'SampleAnalyzer' over dataset 'defaultset'...
    *******************************************************
    * SampleAnalyzer 2.0 for MadAnalysis 5 - Welcome.
    * Configuration status:
     -> employing event weight information
     -> Applying the analysis = 'MadAnalysis5job'
     -> checking the merging procedure
    * Extracting event samples...
    * 1/1 /afs/cern.ch/user/j/jruizalv/public/TTbar_pythia_events.hep.gz
        => sample produced by an unknown generator (xsection assumed in pb).
ERROR: Corrupted STDHEP block: missing 4-momentum information. Event ignored.
#--------------------------------------------------------------------------
# FastJet release 3.0.3
# M. Cacciari, G.P. Salam and G. Soyez
# A software package for jet finding and analysis at colliders
# http://fastjet.fr
#
# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
# for scientific work and optionally PLB641(2006)57 [hep-ph/0512210].
#
# FastJet is provided without warranty under the terms of the GNU GPLv2.
# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
# and 3rd party plugin jet algorithms. See COPYING file for details.
#--------------------------------------------------------------------------
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
WARNING: Skipping a Infinity value in an histogram.
ERROR: Block with the ID=0 is not managed by SampleAnalyzer
    *******************************************************
   Checking SampleAnalyzer output...
 ** ERROR: File called '/afs/cern.ch/work/j/jruizalv/public/Clean/CMSSW_6_1_0/src/MadAnalysis/madanalysis-current/TTbar_matching_check_xqcut20/Output/_defaultset.saf' is not found.
 ** ERROR: errors have occured during the analysis.

In this space at IPNL serv is where I produced good hadronized samples that where read without any trouble by MadAnalysis.

Cheers,
Jose
________________________________________
From: <email address hidden> [<email address hidden>] on behalf of Benjamin Fuks [<email address hidden>]
Sent: 20 March 2013 20:01
To: Jose David Ruiz Alvarez
Subject: Re: [Question #224637]: Warnings during pythia-pgs installation

Your question #224637 on MadGraph5 changed:
https://answers.launchpad.net/madgraph5/+question/224637

Benjamin Fuks requested more information:
Hi Jose, Olivier,

This indeed means that the sample is corrupted in some way. Could you
please send me the mg5 scripts you are using, with the appropriate
cards, to generate such a sample? I would like to try to reproduce the
issue.

Cheers,

Benjamin

--
To answer this request for more information, you can either reply to
this email or enter your reply at the following page:
https://answers.launchpad.net/madgraph5/+question/224637

You received this question notification because you asked the question.

Benjamin Fuks (fuks) said : #5

Hi Jose,

Am I summarizing the problem correcting

(a) you generate everything on the IPNL machine, from the parton level events to the hadron level events. You then process the stdhep file with MadAnalysis 5 and there is no problem
(b) you generate the parton level event file on lxplus. Then you copy paste your working directory on the IPNL machine and use the latter to run pythia. The hadron level event file that is generated is found to be corrupted.

This is odd. However, the different system configurations on the IPNL machine and LXPLUS may explain why this does not work. I am afraid that I cannot help you here. Maybe Olivier, Johan or Eric could. Why can't you run everything on the lxplus machine? I have tried and this work well (I have tested simple processes).

Cheers,

Benjamin

Hi Benjamin,

Indeed is option (b). What happened is that I did a full, parton + pythia, production on lxplus and I found a corrupted file, that's why I posted this question. A full production only in IPNL machine works but then in order to save the job I did producing the parton level files, that are good, in lxplus you suggested to me to copy them to IPNL machine to hadronize them there with pythia. So, in order to avoid this technicalities I'll do everything again from the start in IPNL machine that I know works. However, this is kind of bug in the chain between MG and pythia.

Thanks,
Jose
________________________________________
From: <email address hidden> [<email address hidden>] on behalf of Benjamin Fuks [<email address hidden>]
Sent: 29 March 2013 19:06
To: Jose David Ruiz Alvarez
Subject: Re: [Question #224637]: Warnings during pythia-pgs installation

Your question #224637 on MadGraph5 changed:
https://answers.launchpad.net/madgraph5/+question/224637

    Status: Open => Needs information

Benjamin Fuks requested more information:
Hi Jose,

Am I summarizing the problem correcting

(a) you generate everything on the IPNL machine, from the parton level events to the hadron level events. You then process the stdhep file with MadAnalysis 5 and there is no problem
(b) you generate the parton level event file on lxplus. Then you copy paste your working directory on the IPNL machine and use the latter to run pythia. The hadron level event file that is generated is found to be corrupted.

This is odd. However, the different system configurations on the IPNL
machine and LXPLUS may explain why this does not work. I am afraid that
I cannot help you here. Maybe Olivier, Johan or Eric could. Why can't
you run everything on the lxplus machine? I have tried and this work
well (I have tested simple processes).

Cheers,

Benjamin

--
To answer this request for more information, you can either reply to
this email or enter your reply at the following page:
https://answers.launchpad.net/madgraph5/+question/224637

You received this question notification because you asked the question.

Can you help with this problem?

Provide an answer of your own, or ask Jose David Ruiz Alvarez for more information if necessary.

To post a message you must log in.