compile cuts.f

Asked by Asma Hadef

Hi Olivier,

After editing "cuts.f" file in subprocess directory, I would like to execute it to check if there is any debug in my code.

for a particular subprocess, for example /Subprocess/P0_gg_gttx_t_blvl_tx_bxlvl, I am trying to compile cuts.f using gfortran by running:
> gfortran -o cuts.exe cuts.f

Unfortunately, I got a long list of errors (See below please) and this happens even without adding my code (by using the original cuts.f).

I tried to use cuts.f90 instead cuts.f but this doesn't solve the problem
I also added a simple main program at the begining of the file in order to be executable as the following:

      program cuts
      write(*,*) pass_point(p)
      end program cuts

I am stuck!
How can I compile cuts.f file please?

Your help and assistance are greatly appreciated.

Best,
Asma

LIST OF ERRORS:

 In file genps.inc:9

     Included at cuts.f:43

      parameter (ng = 96, maxdim = 3*(max_particles-2)-1, maxinvar= 4*ma
                                                                       1
Error: Symbol 'ma' at (1) has no IMPLICIT type
 In file genps.inc:25

     Included at cuts.f:43

      parameter (maxprb = maxconfigs*maxplace*maxpoints)
                                   1
Error: Parameter 'maxconfigs' at (1) has not been declared or is a variable, which does not reduce to a constant expression
 In file genps.inc:27

     Included at cuts.f:43

      parameter (maxfprb = maxinvar*maxplace*maxpoints)
                                  1
Error: Parameter 'maxinvar' at (1) has not been declared or is a variable, which does not reduce to a constant expression
 In file cuts.f:81

      double precision etmin(nincoming+1:nexternal),etamax(nincoming+1:n
                                                                       1
Error: Expected another dimension in array declaration at (1)
 In file cuts.f:83

      double precision r2min(nincoming+1:nexternal,ni
                                                                       1
Error: Expected another dimension in array declaration at (1)
 In file cuts.f:85

      double precision etmax(nincoming+1:nexternal),etamin(nincoming+1:n
                                                                       1
Error: Expected another dimension in array declaration at (1)
 In file cuts.f:87

      double precision r2max(nincoming+1:nexternal,nincoming+1:nexternal
                                                                       1
Error: Expected another dimension in array declaration at (1)
 In file cuts.f:89

      double precision ptll_min(nexternal,nexternal),ptll_max(nexternal,
                                                                       1
Error: Bad array specification for an explicitly shaped array at (1)
 In file cuts.f:166

         write(*,formstr) 'Et >',(etmin(i),i=nincoming+1,nexternal)
                                      1
Error: Expected a right parenthesis in expression at (1)
 In file cuts.f:168

         write(*,formstr) 'Eta <',(etamax(i),i=nincoming+1,nexternal)
                                        1
Error: Expected a right parenthesis in expression at (1)
 In file cuts.f:173

     & (r2min(i,j),i=j+1,nexternal)
                      1
Error: Expected a right parenthesis in expression at (1)
 In file cuts.f:175

               r2min(i,j)=r2min(i,j)*dabs(r2min(i,j)) !Since r2 retur
              1
Error: Unclassifiable statement at (1)
 In file cuts.f:176

               r2max(i,j)=r2max(i,j)*dabs(r2max(i,j))
              1
Error: Unclassifiable statement at (1)
 In file cuts.f:320

           IF(mij.ge.500 .and. detaij.le.4 .and. isValid(i) .and. isVali
                                                                       1
Error: Syntax error in IF-expression at (1)
 In file cuts.f:323

           ENDIF
             1
Error: Expecting END DO statement at (1)
 In file cuts.f:373

         if(debug) write (*,*) 'pt(',i,')=',pt(p(0,i)),' ',etmin(i),
                                                                 1
Error: Syntax error in WRITE statement at (1)
 In file cuts.f:375

        notgood=(pt(p(0,i)) .lt. etmin(i)).or.
                                     1
Error: Expected a right parenthesis in expression at (1)
 In file cuts.f:415

      if(debug.and.ptemp(0).gt.0d0) write (*,*) 'Energy of leptons =',pt
                                                                       1
Error: Function 'pt' requires an argument list at (1)
 In file cuts.f:426

         if(dsqrt(SumDot(ptemp,ptemp2,1d0)).lt.mmnl.or.mmnlmax.ge.0d0.an
                                                                    1
Error: Syntax error in IF-expression at (1)
 In file cuts.f:431

      endif
        1
Error: Expecting END FUNCTION statement at (1)
 In file cuts.f:473

         notgood=(etamax(i).ge.0.and.abs(rap(p(0,i))) .gt. etamax(i)).or
                       1
Error: Expected a right parenthesis in expression at (1)
 In file cuts.f:489

            if(debug) write (*,*) 'r2(',i, ',' ,j,')=',dsqrt(r2(p(0,i),p
                                                                       1
Error: Syntax error in argument list at (1)
 In file cuts.f:490

            if(debug) write (*,*) dsqrt(r2min(j,i)),dsqrt(r2max(j,i))
                                            1
Error: Syntax error in argument list at (1)
 In file cuts.f:491

            if(r2min(j,i).gt.0.or.r2max(j,i).ge.0d0) then
                   1
Error: Syntax error in IF-expression at (1)
 In file cuts.f:493

               notgood=(tmp .lt. r2min(j,i)).or.
                                     1
Error: Expected a right parenthesis in expression at (1)
 In file cuts.f:500

            endif
              1
Error: Expecting END DO statement at (1)
 In file cuts.f:509

            if(debug)write (*,*) 'ptll(',i,',',j,')=',PtDot(p(0,i),p(0,j
                                                                       1
Error: Invalid form of array reference at (1)
 In file cuts.f:510

            if(debug)write (*,*) ptll_min(j,i),ptll_max(j,i)
                                        1
Error: Syntax error in WRITE statement at (1)
 In file cuts.f:511

            if(ptll_min(j,i).gt.0.or.ptll_max(j,i).ge.0d0) then
                      1
Error: Syntax error in IF-expression at (1)
 In file cuts.f:513

               notgood=(tmp .lt. ptll_min(j,i).or.
                                        1
Error: Expected a right parenthesis in expression at (1)
 In file cuts.f:520

            endif
              1
Error: Expecting END DO statement at (1)
 In file cuts.f:533

            if(debug) write (*,*) 's(',i,',',j,')=',Sumdot(p(0,i),p(0,j)
                                                                       1
Error: Syntax error in argument list at (1)
 In file cuts.f:693

     $ ptl1max.ge.0d0.or.ptl2max.ge.0d0.or.
                  1
Error: Bad character '1' in OPERATOR name at (1)
 In file cuts.f:733

               notgood=notgood.or.
              1
Error: Unclassifiable statement at (1)
 In file cuts.f:746

      endif
        1
Error: Expecting END FUNCTION statement at (1)
 In file cuts.f:844

        write(6,*) 'alpha_s for scale ',scale,' is ', G**2/(16d0*atan(1d
                                                                       1
Error: Missing exponent in real number at (1)
 In file cuts.f:458

         if(debug) write (*,*) 'p(0,',i,')=',p(0,i),' ',emin(i),':',em
                                                                       1
Error: Symbol 'em' at (1) has no IMPLICIT type
 In file run.inc:11

     Included at cuts.f:73

     $ fixed_couplings,ickkw,nhmult,hmult,asrwgtflavor
                                    1
Error: Symbol 'fixed_facfixed_couplings' at (1) has no IMPLICIT type
 In file cuts.f:91

      common/to_cuts/ etmin, emin, etamax, r2min, s_min,
                                         1
Error: Symbol 'etamax' at (1) has no IMPLICIT type
 In file cuts.f:91

      common/to_cuts/ etmin, emin, etamax, r2min, s_min,
                           1
Error: Symbol 'etmin' at (1) has no IMPLICIT type
 In file cuts.f:92

     $ etmax, emax, etamin, r2max, s_max, ptll_min, ptll_max, inclHt
                             1
Error: Symbol 'etamin' at (1) has no IMPLICIT type
 In file cuts.f:92

     $ etmax, emax, etamin, r2max, s_max, ptll_min, ptll_max, inclHt
               1
Error: Symbol 'etmax' at (1) has no IMPLICIT type
 In file run.inc:11

     Included at cuts.f:73

     $ fixed_couplings,ickkw,nhmult,hmult,asrwgtflavor
                                                       1
Error: Symbol 'hmult' at (1) has no IMPLICIT type
 In file cuts.f:92

     $ etmax, emax, etamin, r2max, s_max, ptll_min, ptll_max, inclHt
                                                               1
Error: Symbol 'ptll_max' at (1) has no IMPLICIT type
 In file cuts.f:92

     $ etmax, emax, etamin, r2max, s_max, ptll_min, ptll_max, inclHt
                                                     1
Error: Symbol 'ptll_min' at (1) has no IMPLICIT type
 In file cuts.f:92

     $ etmax, emax, etamin, r2max, s_max, ptll_min, ptll_max, inclHt
                                    1
Error: Symbol 'r2max' at (1) has no IMPLICIT type
 In file cuts.f:91

      common/to_cuts/ etmin, emin, etamax, r2min, s_min,
                                                1
Error: Symbol 'r2min' at (1) has no IMPLICIT type

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

Compiling only this file is clearly not enough since it has a couple of dependencies like (for the trivial ones) in the kinematical functions. But they are much more dependencies than that.

The way I debug this file is the following:

1) edit cuts.f in order to set the debug data statement to True.
2) compile the full directory
3) cd G1
4) ../madevent < input_app.txt

and then ctrl-c and read the log to see if everything is fine.

Cheers,

Olivier

Revision history for this message
Asma Hadef (asma-hadef-2012) said :
#2

Thanks Olivier Mattelaer, that solved my question.