|
| 1 | +The main purpose of this code is to fit the TPOL waveforms and create a |
| 2 | +tree to be analyzed in the next stage. |
| 3 | + |
| 4 | +An example fit of a TPOL waveform can be found in Fig. 9 of the TPOL |
| 5 | +analysis note ( http://meson.hldsite.com/tpol/analysisNoteV2.pdf ) |
| 6 | + |
| 7 | +The input trees from the collaboration are typically in a directory that |
| 8 | +looks like tree_TPOL/merged/ |
| 9 | + |
| 10 | +Example on how to run: For run 90063 with input file tree_TPOL_090063.root |
| 11 | +and output file stage1Out_090063.root you issue the command |
| 12 | + |
| 13 | + ./stage1 tree_TPOL_090063.root -o./stage1Out_090063.root -r90063 |
| 14 | + |
| 15 | +Some histograms that might be of interest: |
| 16 | + |
| 17 | +root -l stage1Out_090063.root |
| 18 | + hLT->Draw() |
| 19 | + |
| 20 | +The hLT histogram is a signal derived from a spice model (LT spice) of the |
| 21 | +electronics. |
| 22 | + |
| 23 | + hADC1->Draw() |
| 24 | + hADC2->Draw() |
| 25 | + . |
| 26 | + . |
| 27 | + . |
| 28 | + hADC100->Draw() |
| 29 | + |
| 30 | +The hADC1, hADC2, ... hADC100 histograms show the waveform fitting for the |
| 31 | +first 100 events. These histograms are saved to serve as a check on the |
| 32 | +fitting routine. The important information taken from these fits is energy |
| 33 | +deposition and timing. See line 315 of stage1.h to see how the fit function |
| 34 | +is defined. Essentially, par0 is a multiplicative factor for the hLT |
| 35 | +histogram, par1 is the timing shift of the hLT histogram, par2 is a constant |
| 36 | +additive factor and par 3,4,5 are used to represent a small sinusoidal |
| 37 | +ripple of noise. |
| 38 | + |
| 39 | + h2Hit->Draw("colz") |
| 40 | + |
| 41 | +The h2Hit histogram is occupied when there are two hits on the TPOL. The plot |
| 42 | +shows sector of first hit versus the sector of the second hit. Often the |
| 43 | +occurrence of a 2-hit event on TPOL is actually a single hit that shows up on |
| 44 | +neighboring sectors. |
| 45 | + |
| 46 | + hEDep->Draw() |
| 47 | + |
| 48 | +The hEDep histogram represents the energy deposition in the detector. |
| 49 | + |
| 50 | + hTDiff->Draw() |
| 51 | + |
| 52 | +The hTDiff histogram gives the timing difference between the PS and the tagger |
| 53 | + |
| 54 | + hEDiffVsESum->Draw("colz") |
| 55 | + |
| 56 | +The hEDiffVsESum histogram is for the PS and shows the energy difference between |
| 57 | +electron-positron and the energy sum. See Fig. 5 of the analysis note for an |
| 58 | +example plot |
| 59 | + |
0 commit comments