data variables - E1039-Collaboration/e1039-wiki GitHub Wiki

List and Meaning of Data Variables

This page is to accumulate the knowledge on the data variables of E906 and E1039. The similar information accumulated during E906 can be found on the UMich E906 TWiki page; http://twiki.npl.illinois.edu/bin/view/E906/TableStructure . That should be imported to this page time-by-time.

The data structure is definded as C++ classes in the following repositories. You are expected to look into them to find a complete list of the variables. The table below shows the specific class/function for each type of data.

Repository Event Detector Hit Beam Scaler
E1039 e1039-core/interface_main SQEvent SQHitVector n/a n/a
E906 seaquest-ktracker SRawEvent SRawEvent n/a n/a

Spectrometer Info (Main DAQ)

Event Info: SQEvent & SRawEvent

  • Turn ID and RF ID
    t.b.w.
  • QIE presum: SQEvent::get_qie_presum(i)
    Never used. Explained in DocDB ????.
  • QIE trigger count: SQEvent::get_qie_trigger_count()
    Never used. Explained in DocDB ????.
  • RF intensity: SQEvent::get_qie_rf_intensity(i) in E1039 & SRawEvent::fIntensity[16] in E906
    It is proportional to the number of protons. The coefficient varies with time/run, but you can use "1 QIE count = 30 protons" in rough estimate.qie_rf_intensity stores the raw ADC counts of RF-bucket intensities measured by the Cherenkov beam counter. It has 33 elements (RF-16...RF+16), where RF+00 corresponds to the triggered RF bucket. RF+01/RF-01 corresponds to the next/previous RF bucket.
  • Trigger:
    • An event can be fired by many triggers. To get the triggers fired : SQEvent::get_trigger()
    • Return is an integer.
    • Convert the decimal to binary to get the bit stream. Each resultant bit is assigned to a trigger.
    • Use get_trigger(SQEvent::TrigerMask) to see if a particular trigger is fired in the event considered.
    • Check TrigerMask for the trigger assignment.

MC Event Info: SQMCEvent

  • Event weight: SQMCEvent::get_weight()
    • The PYTHIA8 generator (PHPythia8) doesn't make use of it and thus it is always 1.
    • The E906 generator (SQPrimaryParticleGen) sets it to "cross section * luminosity". You have to count up the event yield with the event weight, e.g. the contents of a histogram bin have to be the sum of event weights.
    • When you are interested in the absolute yield, you normalize the (weighted) event yield by the simulated luminosity, where the simulated luminosity per job (=DST) is provided by PHGenIntegral. The normalized (weighted) event yield means the differential cross section of the bin. You could multiply it by a value of anticipated luminosity (like per one week, one month, etc.) to obtain the anticipated yield in the bin.

Hit Info: SQHit

Beam Info (Beam DAQ)

  • Detector and readout: DocDB 537
  • Live PoT: DocDB1212
  • Intensity:Chamber (Drift?) intensity interms of QIE unit. (Question from Abi: How it is different from RF+00? Is it collection of all RF buckets QIE that falls into the drift time (600 ns) ?
  • Intensity_p:
  • PotPerQie:

Scaler Info (Scaler DAQ)

ACNET Info

ACNET = Accelerator Control NETwork managed by Fermilab.

Name Unit Meaning
F:NS2FLO gpm Flow rate of LCW at NS2?
F:NS2SUP degF Temperature of LCW supply at NS2?
F:NS2RET degF Temperature of LCW return at NS2?
F:NM2V amps
F:NM2D1 amps
F:NM2H amps
F:NM3S amps FMag current.
F:NM4AN amps KMag current.
I:BEAM21 E12
I:VFOUT Hz
E:M3TGHM mm Mean of horizontal beam profile measured at NM3.
E:M3TGHS mm Sigma of horizontal beam profile measured at NM3.
E:M3TGVM mm Mean of vertical beam profile measured at NM3.
E:M3TGVS mm Sigma of vertical beam profile measured at NM3.
E:M3TGHI n/a Fitted center of horizontal beam profile?? Obsolete?
E:M3TGHF n/a Fitted width of horizontal beam profile?? Obsolete?
E:M3TGVI n/a Fitted center of vertical beam profile?? Obsolete?
E:M3TGVF n/a Fitted width of vertical beam profile?? Obsolete?
E:M2C2HM mm At NM2.
E:M2C2HS mm At NM2.
E:M2C2VM mm At NM2.
E:M2C2VS mm At NM2.
E:M2C2HI n/a At NM2.
E:M2C2HF n/a At NM2.
E:M2C2VI n/a At NM2.
E:M2C2VF n/a At NM2.
L:CBAR MB At NM2.
S:KTEVTC Torr At NM2.
F:NM2ION ppp Number of protons measured with Ion Chamber at NM2.
F:NM3ION ppp Number of protons measured with Ion Chamber at NM3.
F:NM3SEM ppp Number of protons measured with SEM at NM3.
F:NS7DFP psi
M:OUTTMP degF Outside temperature?
F:E906BM cnts Scaler counts of E906 Beam Monitor (which no longer exists).
S:G2SEM ppp Number of protons measured with SEM at G2.
F:NM3RST degF
F:NM3RRT degF
G:TURN13 Number of turns where a set of protons are injected into single RF bucket.
G:BNCH13 bnch Number of filled RF buckets per cycle (= 84 RF buckets). 84 at max.
G:NBSYD Number of filled RF buckets every 7 RF buckets. 7 at max.
S:MSEP1U Mils
S:F1SEM Ptns Number of protons measured with SEM at F1?
F:NM2Q1 Amps
F:NM2Q2 Amps
U:TODB25 Date
G:RD3161 mrem
G:RD3162 mrem
I:FTSDF 53k% Duty factor sampled at 53 kHz.
F:MT6SC1 Cnts
F:MC7SC1 Cnts
F:MC1D amps
F:MW1W Amps
F:NM4LCWFLOW gpm Flow rate of LCW at NM4.
F:NM4LCWP1 psig Pressure of LCW return at NM4.
F:NM4LCWP2 psig Pressure of LCW supply before regulator(?) at NM4.
F:NM4LCWP3 psig Pressure of LCW supply at NM4.
F:NM4LCWT1 degF Temperature of LCW return at NM4.
F:NM4LCWT2 degF Temperature of LCW supply at NM4.
F:NM4LCWT3 degF Invalid.
G:OUTTMP degF Outside temperature.
G:HUMID % Outside humidity.
G:BPRESS mBar Barometric pressure.