SPINUP(I) = 9 for all particles in all events for calchep_batch generated .lhe files
I am facing a problem while using calchep generated .lhe files in pythia8 interface. The beam polarization option is not working and in the generated .lhe files, the last column of all the events (I guess the spin info for the corresponding particles, the standrad .lhe format SPINUP(I)) is set to value 9.0 for all particles. This signifies lack of spin info in .lhe files, Unknown or unpolarized particles generally be given SPINUP(I)=9. For that reason Pythia8 is crashing.
Please find my batch file attached. Is there a way to resolve this issue?
#######
Model: SM
Model changed: False
Gauge: Feynman
Process: e%,E%->n,N,A
Polarization1:-0.5
Polarization2:0.5
#Composite: le=e,m
#Composite: Le=E,M
Composite: n=ne,nm
Composite: N=Ne,Nm
Kinematics : 12 -> 34 , 5
Kinematics : 34 -> 3 , 4
pdf1: ISR & Beamstrahlung
pdf2: ISR & Beamstrahlung
Bunch x+y sizes (nm) : 480
Bunch length (mm) : 0.3
Number of particles : 2E+10
p1: 500
p2: 500
Cut parameter: T(A)
Cut invert: False
Cut min: 10
Cut max: 500
Cut parameter: E(A)
Cut invert: False
Cut min: 10
Cut max: 500
Cut parameter: Y(A)
Cut invert: False
Cut min: -2.5
Cut max: +2.5
Regularisation momentum: 12
Regularization mass: MZ
Regularization width: wZ
Regularization power: 2
Regularisation momentum: 34
Regularization mass: MZ
Regularization width: wZ
Regularization power: 2
Regularisation momentum: 13
Regularization mass: MW
Regularization width: wW
Regularization power: 2
Regularisation momentum: 24
Regularization mass: MW
Regularization width: wW
Regularization power: 2
Regularisation momentum: 24
Regularization mass: MW
Regularization width: wW
Regularization power: 2
Regularisation momentum: 14
Regularization mass: MW
Regularization width: wW
Regularization power: 2
Dist parameter: M(n,N)
Dist min: 0
Dist max: 600
Dist n bins: 100
Dist title: e,E->n,N,A
Dist x-title: M(n,N) (GeV)
Dist parameter: N(A)
Dist min: -5
Dist max: 5
Dist n bins: 300
Dist title: e,E->n,N,A
Dist x-title: N(A)
Dist parameter: T(A)
Dist min: 0
Dist max: 600
Dist n bins: 100
Dist title: e,E->n,N,A
Dist x-title: T(A) (GeV)
Dist parameter: E(A)
Dist min: 0
Dist max: 600
Dist n bins: 100
Dist title: e,E->n,N,A
Dist x-title: E(A) (GeV)
Dist parameter: Y(A)
Dist min: -4.0
Dist max: +4.0
Dist n bins: 100
Dist title: e,E->n,N,A
Dist x-title: Y(A)
Dist parameter: J(n,N)
Dist min: 0
Dist max: 7
Dist n bins: 300
Dist title: e,E->n,N,A
Dist x-title: J(n,N)
Number of events (per run step): 11000
Filename: ee_nNA1
NTuple: False
Cleanup: False
Max number of cpus: 8
nSess_1: 5
nCalls_1: 100000
nSess_2: 5
nCalls_2: 100000
#######
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- CalcHEP Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Alexander Belyaev
- Solved:
- Last query:
- Last reply: