Draft: Add momentum- and position-based pseudorapidity calculations to CompositeParticle

The CompositeParticle eta calculation uses the secondary vertex and PV positions to calculate pseudorapidity. This choice was not clear to users and doesn't give sensible results for non-displaced vertices. It could also provide worse resolution than a momentum based calculation for fully reconstructed particles.

This MR:

• Removes the CompositeParticle::eta
• Adds CompositeParticle::eta_position and CompositeParticle::eta_momentum
• Replaces eta with eta_position and eta_momentum in all selections. Unless the line authors' intentions were clear, I used eta_position for lines selecting partially reconstructed decays and eta_momentum for lines selecting fully reconstructed decays.

FYI @gligorov

device/combiners/src/FilterTwoSvs.cu

@@ -112,8 +112,8 @@ __global__ void FilterTwoSvs::filter_two_svs(
     if (dec) {
       // Kinematic cuts.
       dec &= vertex.minp() > minTrackPV1;
-      dec &= vertex.eta() > minEtaV1;
-      dec &= vertex.eta() < maxEtaV1;
+      dec &= vertex.eta_position() > minEtaV1;
+      dec &= vertex.eta_position() < maxEtaV1;
       dec &= vertex.dira() > minCosDiraV1;
       dec &= vertex.minpt() > minTrackPtV1;
       dec &= vertex.vertex().pt() > minPtV1;
@@ -135,8 +135,8 @@ __global__ void FilterTwoSvs::filter_two_svs(
     if (dec) {
       // Kinematic cuts.
       dec &= vertex.minp() > minTrackPV2;
-      dec &= vertex.eta() > minEtaV2;
-      dec &= vertex.eta() < maxEtaV2;
+      dec &= vertex.eta_position() > minEtaV2;
+      dec &= vertex.eta_position() < maxEtaV2;
       dec &= vertex.dira() > minCosDiraV2;
       dec &= vertex.minpt() > minTrackPtV2;
       dec &= vertex.vertex().pt() > minPtV2;