JeffersonLab / DEMPgen

A Deep Exclusive Meson Production (DEMP) event generator for SoLID and EIC simulations.
GNU General Public License v3.0
1 stars 6 forks source link

0 recorded events in ep -> e+p+pi0 prouction #12

Open michael-pitt opened 1 year ago

michael-pitt commented 1 year ago

Dear experts,

I'm running the code with the following settings, but unfortunately no events are recorded.

    "particle": "pi0",         // Choices: omega, pi+, pi0, K+
    "hadron": "proton",                // Choices: Neutron, proton, Lambda or Sigma0
    "ebeam": 5,                         // Electron beam energy in GeV
    "hbeam": 100,                        // Hadron beam energy in GeV
    "hbeam_part":"Proton",           // Hadron beam particle, proton, deut or helium3, work in progress, will need to be added to shell script later on as an argument
    "OutputType": "HEPMC3",        // choices: LUND (Docker), Pythia6 (ECCE Fun4All) or HEPMC3 (Athena)   
    "det_location": "ip6",         // choices: ip6 for STAR, ip8 for PHENIX
    "Ee_Low": 0.0, // The minimum scattered electron energy that will be generated as a fraction of the electron beam energy
    "Ee_High": 2.5, // The maximum scattered electron energy that will be generated as a fraction of the electron beam energy
    "e_Theta_Low": 0.0, // The minimum scattered electron angle (theta) that will be generated in degrees
    "e_Theta_High": 175.0, // The maximum scattered electron angle (theta) that will be generated in degrees

PS: During the production, this line: https://github.com/JeffersonLab/DEMPgen/blob/master/src/eic_evgen/process_routine/Pi0_prod.cc#L62 is executed. Should it be called only if LundFiles are written? The LundFiles folder is not created anyway.

Thank you in advance, Michael

sjdkay commented 1 year ago

Hi Michael,

Bill will know more about the pi0 section of the code, I've tagged him on this. Hopefully he can help out.

Thanks, Stephen

huberg commented 1 year ago

Dear Michael,

We have been concentrating our efforts on the t-channel pi+ and K+ generators for the initial release. The u-channel pi0 cross section model is not very sophisticated, as it is simply a port of a model used for JLab kinematics.

I'm guessing that you are getting zero events right now because the kinematics are outside the region of applicability of the model. Sorry that it's not very helpful to you yet, this is on our list of things to work on, but we haven't gotten there yet.

Regards, Garth Huber

Michael Pitt wrote:

Dear experts,

I'm running the code with the following settings, but unfortunately no events are recorded.

|"particle": "pi0", // Choices: omega, pi+, pi0, K+ "hadron": "proton", // Choices: Neutron, proton, Lambda or Sigma0 "ebeam": 5, // Electron beam energy in GeV "hbeam": 100, // Hadron beam energy in GeV "hbeam_part":"Proton", // Hadron beam particle, proton, deut or helium3, work in progress, will need to be added to shell script later on as an argument "OutputType": "HEPMC3", // choices: LUND (Docker), Pythia6 (ECCE Fun4All) or HEPMC3 (Athena) "det_location": "ip6", // choices: ip6 for STAR, ip8 for PHENIX "Ee_Low": 0.0, // The minimum scattered electron energy that will be generated as a fraction of the electron beam energy "Ee_High": 2.5, // The maximum scattered electron energy that will be generated as a fraction of the electron beam energy "e_Theta_Low": 0.0, // The minimum scattered electron angle (theta) that will be generated in degrees "e_Theta_High": 175.0, // The maximum scattered electron angle (theta) that will be generated in degrees |

PS: During the production, this line: https://github.com/JeffersonLab/DEMPgen/blob/master/src/eic_evgen/process_routine/Pi0_prod.cc#L62 is executed. Should it be called only if LundFiles are written? The |LundFiles| folder is not created anyway.

Thank you in advance, Michael

— Reply to this email directly, view it on GitHub https://github.com/JeffersonLab/DEMPgen/issues/12, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABXIPC67Q522FHTN6BCXRZLXPUO3JANCNFSM6AAAAAA2FXY4OE. You are receiving this because you are subscribed to this thread.Message ID: @.***>

-- Dr. Garth Huber, Dept. of Physics, Univ. of Regina, Regina, SK S4S0A2, Canada. tel: 1-306-585-4240. fax: 1-306-585-5659.