calchep_batch not changing parameter

Dear Calchep Experts

In my BSM model I am trying to generate .lhe files for different values of mass parameters (called MX and metap). I am using the following batch _file:

Model: I(1+2)HDMtypeI
Model changed: False
Gauge: Feynman
Process: p,p->j,~X,~X
Alias: p=u,u~,d,d~,s,s~,c,c~,b,b~,g
Alias: j=u,u~,d,d~,s,s~,c,c~,b,b~,g

pdf1: NNPDF23_lo_as_0130_qed(proton)
pdf2: NNPDF23_lo_as_0130_qed(proton)

p1: 6500
p2: 6500

Run parameter: MX
Run begin: 50
Run step size: 5
Run n steps: 5

#labc=-0.2
Parameter: metap=pow(pow(MX,2)-246*246/2*(-0.2),0.5)

parton dits. alpha : ON
alpha(MZ) : 0.118
alpha nf : 5
alpha order : NLO
mb(mb) : 4
Mtop(pole) : 173
alpha Q : T45
Factorization1 Qren
Factorization2 Qren
Showering Qren

#Cuts
cut parameter : T(~X,~X)
cut invert: False
cut min: 100
cut max:

Number of events (per run step): 10000
Filename: labc__022

#(default) Parallelization mode: local

nSess_1 : 5
nCalls_1: 100000
nSess_2 : 5
nCalls_2: 100000

as you can see from above, the I want to vary the parameter metap as a function of the parameter MX but when I look into the generated session.dat files I notice that metap is being held fixed at zero. I also noticed that the cut on the missing tranverse momentum is not being applied. Here as en example of a session.dat file that I am getting:

"#Subprocess 1 ( u~, g -> u~, ~X, ~X )
#Session_number 1
#Initial_state inP1=6.500000E+03 inP2=6.500000E+03
 Polarizations= { 0.000000E+00 0.000000E+00 }
  StrFun1="PDT:NNPDF23_lo_as_0130_qed(proton)" 2212
  StrFun2="PDT:NNPDF23_lo_as_0130_qed(proton)" 2212

#Physical_Parameters
      beta = 1.331565000000000E+00
      mu3p = 1.637345720000000E+02
      mu3m = 0.000000000000000E+00
    alpha1 = -2.400000000000000E-01
     metap = 0.000000000000000E+00
     metam = 0.000000000000000E+00
      leta = 1.490000000000000E+00
     aEWM1 = 1.279400000000000E+02
        Gf = 1.166370000000000E-05
        aS = 1.184000000000000E-01
         Q = 1.000000000000000E+02
      ymdo = 1.000000000000000E-02
      ymup = 1.000000000000000E-02
       yms = 2.000000000000000E-01
      MbMb = 4.230000000000000E+00
      McMc = 1.270000000000000E+00
       Mtp = 1.735000000000000E+02
       yme = 5.110000000000000E-04
       ymm = 1.057000000000000E-01
     ymtau = 1.777000000000000E+00
      cabi = 2.277360000000000E-01
        MZ = 9.118760000000000E+01
        Me = 5.110000000000000E-04
       MMU = 1.057000000000000E-01
       MTA = 1.777000000000000E+00
        MU = 1.000000000000000E-02
        MD = 1.000000000000000E-02
        MS = 2.000000000000000E-01
        Mh = 1.250000000000000E+02
        MH = 2.934600000000000E+02
        MA = 3.652000000000000E+02
        MX = 5.000000000000000E+01
       MXa = 8.962600000000000E+02
       MHP = 4.027300000000000E+02
       MXP = 9.000000000000000E+02
         E = 2.718281828450000E+00
        Pi = 3.141592653580000E+00
----

#Breit-Wigner BW range 2.700000
t-channel widths 0
GI trick in s- 0
GI trick in t- 1
#VVdecays 1

#alphaQCD alphaPDF=1 alpha(MZ)=1.180000E-01 NF=5 Order=2 MbMb=4.000000E+00 Mtp=1.730000E+02
#QCDscales
 Renormalization T45
 Factorization1 T45
 Factorization2 T45
 Showering T45

#Composites
*** Table ***
 Composites
  Name |> Comma separated list of particles <|
p |u,u~,d,d~,s,s~,c,c~,b,b~,g
j |u,u~,d,d~,s,s~,c,c~,b,b~,g
==========================================================================

#Cuts
*** Table ***
 Cuts
!| Parameter |> Min bound <|> Max bound <|
============================================

#Parallelization 1

#Distributions
*** Table ***
Distributions
Parameter_1|> Min_1 <|> Max_1 <|Parameter_2|> Min_2 <|> Max_2 <|
====================================================================

#Kinematical_scheme
12 -> 3 , 45
45 -> 4 , 5

#Regularization
*** Table ***
 Regularization
 Momentum |> Mass <|> Width <| Power|
=========================================

#Vegas_calls 100000x5 100000x5
#Vegas_integral 6.84200762631369669E-02 1.75161334626699161E+00 6.87012449837068662E-01 5 1 466560 1 -1
#Events 10000 1 10000
#Random B6AE37FF5CEE

#VEGAS_Grid Vegas_grid: dim=6 size=50
....
"

Could you please tell me why metap is not being scanned and why the cut I want to impose seems not being applied by the program?

Thanks in advance
Marco