Generating MC for specific reactions with gen_amp

Current builds of the GlueX software are located at /sciclone/home10/jrstevens01/analysisGluexI/builds/ and can be set in your environment with this command

source /sciclone/home10/jrstevens01/analysisGluexI/builds/setup.csh

Example generator configuration

An example generation file for the K+K-pi0 reaction shows how amplitudes can be constructed for different resonance contributions:

reaction Pi0K+K- Beam Proton Pi0 K+ K-

sum Pi0K+K- xpol

amplitude Pi0K+K-::xpol::eta1295_a0pi_S BreitWigner eta1295 0 34 2 amplitude Pi0K+K-::xpol::eta1295_a0pi_S BreitWigner a0 0 3 4

amplitude Pi0K+K-::xpol::f11420_KStarPlusK_P BreitWigner f11420 1 23 4 amplitude Pi0K+K-::xpol::f11420_KStarPlusK_P BreitWigner KStar 0 2 3

amplitude Pi0K+K-::xpol::f11420_KStarMinusK_P BreitWigner f11420 1 24 3 amplitude Pi0K+K-::xpol::f11420_KStarMinusK_P BreitWigner KStar 0 2 4

initialize Pi0K+K-::xpol::eta1295_a0pi_S cartesian 1.0 0.0 initialize Pi0K+K-::xpol::f11420_KStarPlusK_P cartesian 1.0 0.0 initialize Pi0K+K-::xpol::f11420_KStarMinusK_P cartesian 1.0 0.0

For other final states a similar set of amplitude definitions for the BreitWigner amplitudes can be defined.

Generating events with MCWrapper

Once the generator config file is prepared you can create an MCWrapper configuration file which will generate events and run the full reconstruction to produce the Analysis TTrees for your reaction. See for example,

To generate events in the local interactive mode for testing you can use the command

gluex_MC.py MC.config 30496 100 batch=0

To submit jobs which generate many more events

gluex_MC.py MC.config 30496 25000 per_file=1000 batch=1