Task: find the optimal pairing of r-φ and stereo modules for the TIB layers.
Strategy followed in my macro:

1. Count number of modules in the two samples (r-φ and stereo).
2. The smallest sample is labeled "sample 1", the other "sample 2".
3. The V_depl average is computed for sample 1.
4. Modules in sample 2 are sorted according to the distance in V_depl from the average of sample 1.
5. Since sample 2 is larger than sample 1 (by definition), the extra modules too far from the V_depl average of sample 1 are thrown away, keeping the same number of modules as in sample 1.
6. Both sample 1 and sample 2 are sorted by V_depl.
7. Pair module 1 of the 1st sample with the closest (in V_depl) module of the 2nd sample, then pair module 2 with the closest of the remaining, and so on.
8. Apply "simulated annealing" to check if more optimal pairings are possible.

Files:

• ROOT macro
• Database (r-φ: 1.1.1.1, stereo left: 1.1.2.2, stereo right: 1.1.3.3)
• SQL query:

  select module_id,v_depl_v from info_module where type='1.1.x.x'

• Input files:
  • r-φ
  • stereo left
  • stereo right
• Output file

Usage:

root
.x ds.C(1,"inputfile1","inputfile2","outputfile",cut,0)
.q
more outputfile