Extracting cross-section for a process

[michael-pitt]

Dear experts,

I'm running a few examples, and I am interested in extracting the cross-section of a given process. Is there a simple way to get the cross-section?

Thank you in advance,
Michael

[sjdkay]

Hello Michael,

This depends a little bit upon what you want. If you just want the total cross section, it's fairly straightforward. In DEMPReaction.cc, you could just add a line to cout r_fSig after it is assigned on line 559 - https://github.com/JeffersonLab/DEMPgen/blob/v1p0_release/src/eic_evgen/process_routine/DEMP_Reaction.cc#L559C3-L559C38

Note that this would only give you the cross section for events that have passed other relevant cuts before this point.

If you want the separated components of the cross section, this is a bit more complex but I can try and outline/explain that if that's what you're after.

Also, I believe Love has a short piece of code to calculate the cross section (and L/T components) for the KLambda and KSigma parameterisation which can be run independently of the generator. Let me know if this is something you might want to use.

Thanks, Stephen