Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

This page holds a few example code-snippets for use in myana analysis. This analysis is written in C++ and the examples here use root for plotting. Compare with Pyana user examples to see how the same things can be done using the pyana analysis framework.

Beamline data (Bld)

To read out energy, charge and position of the beam from the beamline data, use getEBeam(). The function returns 0 if data is available and assignes values to its arguments from the current event. Call from within event().

Code Block

  double ebcharge; double ebenergy; double posx; double posy;
  double angx; double angy;
  int fail =  getEBeam(ebcharge, ebenergy, posx, posy, angx, angy);
  if (!fail) printf("ebeam: %f, %f, %f, %f, %f, %f \n", 
                     ebcharge, ebenergy, posx, posy, angx, angy);

To read out the energy from the front end enclosure (FEE) gas detector, use getFeeGasDet(). This assigns this event's value to an array of 4 numbers:

Code Block

  double shotEnergy[4];
  int fail = getFeeGasDet( shotEnergy );
  if (!fail) printf("GasDet energy: %f, %f, %f, %f \n ", shotEnergy[0], shotEnergy[1], shotEnergy[2], shotEnergy[3]);

To read out fit time and charge of the phase cavity, use getPhaseCavity() which assigns values to four arguments:

Code Block

  double phaseCavityTime1, phaseCavityTime2, phaseCavityCharge1, phaseCavityCharge2;
  int fail = getPhaseCavity( phaseCavityTime1, phaseCavityTime2, phaseCavityCharge1, phaseCavityCharge2 );
  if (!fail) printf("PhaseCavity: %f, %f, %f, %f \n", 
                     phaseCavityTime1,  phaseCavityTime2, phaseCavityCharge1, phaseCavityCharge2);

Display images from Princeton camera

...