L2jet Investigation: Energy Discrepency between Real and Simulated data
One of the problems which came to light when comparing the diagnostic histos from L2jet for real data vs simulation was that the jet energy for the real data was much higher than that for simulation and much higher than makes sense. To start investigating this problem I have made text files listing the values of several parameters used in the projectAdc() function, which gets the adc values from the raw data blocks. The code for projectAdc() can be found here. This is the L2jetAlgo2008 code which uses the 3x3 JP. This does not use the new low noise ET that Ross has implemented for L2VirtualAlgo2009. I run over 5 events from the real data set: /star/institutions/mit/balewski/2008-L2events/R9067013.eve.bin-60k
and 5 events from the simu data set: /star/institutions/iucf/balewski/star/institutions/mit/balewski/StarTrigSimuEvents/pythia_pT7-9_rcf1225.eve.bin-348k
The output from the real data can be found here.
The output from the simu data can be found here.
Lines showing only the threshold and rawADC values are for towers who's adc did not pass threshold while lines showing all values are for towers who's adc's did pass threshold. For each event, all 4800 barrel towers are listed then the number of towers above threshold is listed, then all 720 endcap towers are listed and the number of those above threshold is given.
More analysis needs to be done, but the most striking feature is that there are ~10 times more towers passing threshold for the real data than for the simu.
I have also listed the values which appear in the arrays eve_patchEne and eve_phiEne. These arrays hold the adc4 values calculated earlier in projectAdc and are used by the scanPhi and scanEta functions to find the two highest jet patches. The values for the real data can be found here and the values for the simu data can be found here.
Finally, I list the output from the scanPhi and scanEta functions. The values for the real data can be found here and the values for the simu data can be found here.
The code I used can be found here: /star/u/pagebs/L2Trig/run9_2008-2009/myWork/initilizationTest/fromV1.0/StRoot/StTriggerUtilities/L2Emulator/L2jetAlgo/L2jetAlgo2008.cxx
The output from multiTest.c for the real data is here and the simu data is here.
- pagebs's blog
- Login or register to post comments