process p,p->q,q

Asked by pytnik89

I declared in gui seession process p,p->q,q; where p(u,U,d,D,c,C,s,S,b,B,G), q(u,U,d,D,c,C,s,S,b,B,t,T,G) and excluded (W+,W-,Z,h), It made diagrams and square it but when i start make and launch i get the error: Error in `/home/pytnik/Programs/calchep_3.6.14/bin/s_calchep': double free or corruption (top): 0x0000000000e595a0
*long log*
CalcHEP has finished with error code 134

Question information

Language:
English Edit question
Status:
Solved
For:
CalcHEP Edit question
Assignee:
No assignee Edit question
Solved by:
Alexander Pukhov
Solved:
Last query:
Last reply:
Revision history for this message
Best Alexander Pukhov (pukhov) said :
#1

I guess you have this problem in a model with hgg effective vertex.
Then I propose to download new version 3.6.15 from WEB where this bug was fixed.

Best
   Alexander Pukhov

Revision history for this message
pytnik89 (pytnik89) said :
#2

Thanks Alexander Pukhov, that solved my question.

Revision history for this message
Fabrizio Nesti (fnesti) said :
#3

Hi - still on this issue.

I also have a model with effective ggh vertex, and calchep just crashes after choosing "Enter Process":

./calchep: line 52: 17768 Abort trap $CALCHEP/bin/s_calchep $*
CalcHEP has finished with error code 134

Any clue? Thanks in advance!

Fabrizio

Revision history for this message
Fabrizio Nesti (fnesti) said :
#4

Sorry - addition, using latest version calchep_3.6.22 and feynrules 2.1 ( 16 May 2014).

f

Revision history for this message
Alexander Pukhov (pukhov) said :
#5

Use "Check models" first. May be it will detect bugs in the model. If not, send me your model.
Best
   Alexander Pukhov

Revision history for this message
Neil Christensen (neil-christensen-qft) said :
#6

Can you tell us more about this model. Did you create it yourself using FeynRules? Did you get it from the FeynRules model database? What steps have you taken to ensure it is correctly implemented? If the problem is on the FeynRules side, you can contact me directly.

Best wishes,
Neil

Revision history for this message
Fabrizio Nesti (fnesti) said :
#7

Hi - thanks for the quick reply.

Indeed CHECK MODEL returns an error in the constraints table "Unknown variable EE". The line is the following:

   AH |(47*pow(EE,2)*(1-(2*pow(MH,4)*pow(MT,-4))/987.-(14*pow(MH,2)*pow(MT,-2))/705.+(213*pow(MH,12)*pow(.........

but EE is defined shortly below (the electric charge :))

  EE |2*pow(aEW,0.5)*pow(Pi,0.5)

I generated the model with feynrules from our own .fr model file.
(The relative UFO and in particular this vertex has been working correctly in madgraph5.)

Basically it contains the standard ggh vertex (hence the functions AH, GH, AD, GD).

Revision history for this message
Fabrizio Nesti (fnesti) said :
#8

addition: the higgs part was taken from the feynrules page
http://feynrules.irmp.ucl.ac.be/wiki/HiggsEffectiveTheory

Revision history for this message
Neil Christensen (neil-christensen-qft) said :
#9

You need to define the parameters in order. So, you need to define EE before defining AH. Any parameters that depend on other parameters must have those other parameters defined earlier. You must do this at the FeynRules level. This is discussed in the FeynRules manual.

Best wishes,
Neil

Revision history for this message
Fabrizio Nesti (fnesti) said :
#10

HI, thanks a lot - I have checked, regenerated, and now the parameters are in order.

However now I stumble into a similar problem: in the generated func1.mdl there is no trace of the definition of gs (strong coupling), while I thought it was properly defined in the FR file.

I.e. from the FR file

  ...
   g1 == {
    ParameterType -> Internal,
    Value -> ee/Sqrt[1-2*sw^2],
    InteractionOrder -> {QED,1},
    TeX -> Subscript[g,1],
    Description -> "U(1)Y coupling constant at the Z pole"
  },
  gs == {
    ParameterType -> Internal,
    Value -> Sqrt[4 Pi aS],
    InteractionOrder -> {QCD,1},
    TeX -> Subscript[g,s],
    Description -> "Strong coupling constant at the Z pole"
  },
  vev == {
    ParameterType -> Internal,
    Value -> 1/Sqrt[Sqrt[2]*Gf],
    TeX -> Vev,
    InteractionOrder -> {EWL,-1},
    Description -> "Higgs vacuum expectation value"
  },
  ...

I only get in func1.mdl :

...
g1 |EE*pow(1-2*pow(sw,2),-0.5) % U(1)Y coupling constant at the Z pole
vev |pow(2,-0.25)*pow(Gf,-0.5) % Higgs vacuum expectation value
...

...quite weird in my understanding. (If I insert the gs coupling by hand in the func1.mdl, everything now works. )

Anyway thanks again!
cheers,
F