3.1.1 Beam Energies

The beam energies are spread from 1 GeV to 1 TeV (1000-1000000 MeV) Their values can be found here: PionBeamEnergies.txt

However, the single beam energy values are not concentrated in single files, but were randomly generated one after the other, thus every file contains every possible beam energy. Here the deposited energy (from the CalibHits , y-Axis) is plotted against the beam energy (Variable spcl_sume, x-Axis). This was done for negative (red) as well as for positive pions (blue).

Fig 1.: AllCluster _SumTotalEnergy plotted against the beam energy (for minus=red and plus=blue), no major differences visible

As a check, the number of events per beam energy was plotted in a chain, which tested the weights for different beam energies. The loop was supposed to lopp over all events, but stopped after 452200 events with the message

Event 452000

ientry < 0, breaking

WHY? As the number of max. events is set to 3000000, this is roughly only one sixth of the total number of entries of the entire chain.... Bug in TChain/Root??

# Number of events to be processed (default is 10)

theApp.EvtMax = 3000000

However, here is the plott of the beam energies, they are logarithmically distributed and each beam energy contains an equal amount of entries.

Fig 2.: Nevents per Beam energy

3.1.2 Eta und Phi characteristics of Single Pion Files

The files were generated with Pions uniformly distributed in eta:

Fig 3: Eta of the generated Pions - uniform distribution

Fig 4: Eta of ALL generated particles

3.1.3 Number of clusters

Fig. 7: Number of clusters per event (consistent with what has been seen in the ApplyWeights macro)

3.1.4 Cellnoise

Fig 8: Level of noise (units - MeV)

3.1.5 Energy of the Cluster vs. eta of the cluster (COG)

(here the z-axis is the number of entries)

Get Eta of the Cluster_COG:
-TMath::Log(TMath::Tan((TMath::ACos(cl_center_z_topo/TMath::Sqrt(cl_center_z_topo*cl_center_z_topo + cl_center_x_topo*cl_center_x_topo + cl_center_y_topo*cl_center_y_topo))/2))))

Scalarproduct a*b = |a| |b| cos (a,b) to get theta with b being the z-axis vector (0,0,1)
(cl_center_z_topo/TMath::Sqrt(cl_center_z_topo*cl_center_z_topo + cl_center_x_topo*cl_center_x_topo + cl_center_y_topo*cl_center_y_topo)

cl_e_topo is the variable for the cluster energy.

-- Main.kristin - 05 Jan 2007