Skip to content
Snippets Groups Projects
Commit f789a906 authored by Theodore Efremov's avatar Theodore Efremov :hibiscus:
Browse files

[AlPhaPha] ReadMe done

parent 9ca92f5c
No related branches found
No related tags found
No related merge requests found
Pipeline #394402 passed
# This file needs to be commented
# AoQ
The goal of this directory of macro is to generate and apply the proper calibration in order to retrieve the mass of the fission fragment detected in VAMOS.
To do so we need:
- To generate proper time offset for the velocity measurement for each MWPC section;
- To generate a set of IC segment coefficient for each MWPC section to recover a proper Z
---
## Toff
To generate the proper Toff we align the peak seen on the 1D histogram of the AoQ on an integer value.
All the macro are stored in the Toff subdirectry. There is currently 3 macros :
- ToffHighRes.C
- MergeToff.C
- ToffAllSec.sh
to use this subdir you just have to do :
```
. ToffAllSec.sh
```
The resulting toff will be in Output/ToffAoQ.cal and should be added to the toff already used in the analysis.
### ToffHighRes.C
This macro use the TSpectrum Search High Res method to find the peaks in the AoQ histogram and then use MIGRAD to find the optimal value of Toff to make one of the peak fall close to 3.
If there is no peak found near 3 it tries to align on 4.
On line 146 you can add or remove a limit in number of events.
It then store the result of the minimisation in Output/Toff/ but those output are only used by MergeToff.C
### MergeToff.C
Simply merge the result of the toff from Output/Toff/ into a single file named ToffAoQ.cal
### ToffAllSec.sh
Launch ToffHighres on all section on all your core expect 4.
Before using this make sure the single macro ToffHighRes work well on a section with a low amount of event.
---
## Q vs AoQ
To generate the proper coefficient for the linear combination of the IC Section we try to straighten a given Q on the QvsAoQ plot.
To do so we must proceed in multiple step :
1. Generate a cut on a given charge state e.g 39 on each of the MWPC section. To do so use DrawAoQ.C
'''
root
.L DrawAoQ.C
DrawAoQ(SectionNumber,0)
'''
2. Store them in a file named cut_QvsAoQ.root. (The cut are named Cut<SecNumber>)
3. Fill the root file named for the minimisation with FillSelectTree.C
```
root FillSelectTree.C
```
4. Try the minimisation on one of the section using Minimization.C further details on this step will be found in next section.
5. If all well and good use MiniParra.sh to launch the minimisation on all section
### Minimization.C
This macro minimize the distance of each point of the 2D histogram Q:AoQ to be as close as possible to the charge selected in the associated cut.
The degrees of freedom are the coefficient in the linear combination of the IC made to reconstruct the Etot.
The initial parameter are on line **26** and **29** and **60** also on line **66** and **69** you will find the up and low bounds for the minimisation.
You should play around with those value until you find a satisfying minimisation.
The number of entries to do the minimisation is on line **128 - 129** you should use all your entries when doing the full minimisation.
Otherwise section >16 will not converge.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment