Outcoming particles not sorted in energy

Asked by Andrea Paolo Puppin

hi , sorry to bother you
this is my code:
process ee = A, e1 => E1,e1,e1

beams = A, e1 => gaussian

beams_momentum = 5 keV, 5 GeV

gaussian_spread1= 0%

gaussian_spread2= 0%

plot E_Theta_Phi_E1

plot E_Theta_Phi_e1

plot E_Theta_Phi_e2

sample_format = mokka

beams_theta = 0,0.1

#?rebuild_events = true

analysis= record E_Theta_Phi_E1 (eval E["E1"], eval Theta["E1"], eval Phi["E1"]) and

          record E_Theta_Phi_e1 (eval E[extract index 1 [e1]], eval Theta[extract index 1 [e1]], eval Phi[extract index 1 [e1]])

          and

          record E_Theta_Phi_e2 (eval E[extract index 2 [e1]], eval Theta[extract index 2 [e1]], eval Phi[extract index 2 [e1]])

write_analysis

simulate (ee) {

luminosity = 0.63 / 1.e8 fbarn

 }

 compile_analysis
i just wanna ask how exactly whizard recognize the two electron after the process of Tpp.
I thought whizard called, for example, the more energetic electron(which should be the incoming electron) with index 1 and the less energetic with index 2.
But if i'm goingo to study theyr energy i obtain two identical profile, with some electron very energetic and some electron with energy near 20 Mev.

Question information

Language:
English Edit question
Status:
Solved
For:
WHIZARD Edit question
Assignee:
Juergen Reuter Edit question
Last query:
Last reply:
Revision history for this message
Juergen Reuter (j.r.reuter) said :
#1

Hi Andrea,
when you extract the index of a particle, it is really just according to the listing of the final state, and it depends on the phase-space point, which one is more or less energetic. If you want the electrons sorted according to energy, you have to sort them by the corresponding SINDARIN command. Please read the corresponding sections in the manual in Chap. 5, especially on the extract and sort keywords on p. 86 (for v3.1.2).
Cheers,
    JRR

Revision history for this message
Juergen Reuter (j.r.reuter) said :
#2

Marking as 'solved'.