Skip to content

Commit ebf2081

Browse files
RrantuRrantu
committed
Simplify variable initialization in candidate loop
1 parent 13be68f commit ebf2081

File tree

1 file changed

+62
-193
lines changed

1 file changed

+62
-193
lines changed

‎PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx‎

Lines changed: 62 additions & 193 deletions
Original file line numberDiff line numberDiff line change
@@ -2653,138 +2653,7 @@ struct HfCandidateCreatorXic0Omegac0Mc {
26532653
int8_t flag = 0;
26542654
McMatchFlag debug{McMatchFlag::None};
26552655

2656-
float v0MassPull;
2657-
float xiMassPull;
2658-
float xic0MassPull;
2659-
2660-
float v0DauPosXDelta;
2661-
float v0DauPosYDelta;
2662-
float v0DauPosZDelta;
2663-
float v0DauPosPt;
2664-
float v0DauPosXPull;
2665-
float v0DauPosYPull;
2666-
float v0DauPosZPull;
2667-
2668-
float v0DauNegXDelta;
2669-
float v0DauNegYDelta;
2670-
float v0DauNegZDelta;
2671-
float v0DauNegPt;
2672-
float v0DauNegXPull;
2673-
float v0DauNegYPull;
2674-
float v0DauNegZPull;
2675-
2676-
float v0XDelta;
2677-
float v0YDelta;
2678-
float v0ZDelta;
2679-
float v0Pt;
2680-
float v0XPull;
2681-
float v0YPull;
2682-
float v0ZPull;
2683-
2684-
float lambdaXDelta;
2685-
float lambdaYDelta;
2686-
float lambdaZDelta;
2687-
2688-
float xiBachelorXDelta;
2689-
float xiBachelorYDelta;
2690-
float xiBachelorZDelta;
2691-
float xiBachelorPt;
2692-
float xiBachelorXPull;
2693-
float xiBachelorYPull;
2694-
float xiBachelorZPull;
2695-
2696-
float xiXDelta;
2697-
float xiYDelta;
2698-
float xiZDelta;
2699-
float xiPt;
2700-
float xiXPull;
2701-
float xiYPull;
2702-
float xiZPull;
2703-
2704-
float xic0BachelorXDelta;
2705-
float xic0BachelorYDelta;
2706-
float xic0BachelorZDelta;
2707-
float xic0BachelorPt;
2708-
float xic0BachelorXPull;
2709-
float xic0BachelorYPull;
2710-
float xic0BachelorZPull;
2711-
2712-
float xic0XDelta;
2713-
float xic0YDelta;
2714-
float xic0ZDelta;
2715-
float xic0Pt;
2716-
float xic0XPull;
2717-
float xic0YPull;
2718-
float xic0ZPull;
2719-
2720-
float cascXDelta;
2721-
float cascYDelta;
2722-
float cascZDelta;
2723-
27242656
for (const auto& candidate : candidates) {
2725-
v0MassPull = 0;
2726-
xiMassPull = 0;
2727-
xic0MassPull = 0;
2728-
2729-
v0DauPosXDelta = 0;
2730-
v0DauPosYDelta = 0;
2731-
v0DauPosZDelta = 0;
2732-
v0DauPosPt = 0;
2733-
v0DauPosXPull = 0;
2734-
v0DauPosYPull = 0;
2735-
v0DauPosZPull = 0;
2736-
2737-
v0DauNegXDelta = 0;
2738-
v0DauNegYDelta = 0;
2739-
v0DauNegZDelta = 0;
2740-
v0DauNegPt = 0;
2741-
v0DauNegXPull = 0;
2742-
v0DauNegYPull = 0;
2743-
v0DauNegZPull = 0;
2744-
2745-
v0XDelta = 0;
2746-
v0YDelta = 0;
2747-
v0ZDelta = 0;
2748-
v0Pt = 0;
2749-
v0XPull = 0;
2750-
v0YPull = 0;
2751-
v0ZPull = 0;
2752-
2753-
lambdaXDelta = 0;
2754-
lambdaYDelta = 0;
2755-
lambdaZDelta = 0;
2756-
2757-
xiBachelorXDelta = 0;
2758-
xiBachelorYDelta = 0;
2759-
xiBachelorZDelta = 0;
2760-
xiBachelorPt = 0;
2761-
xiBachelorXPull = 0;
2762-
xiBachelorYPull = 0;
2763-
xiBachelorZPull = 0;
2764-
2765-
xiXDelta = 0;
2766-
xiYDelta = 0;
2767-
xiZDelta = 0;
2768-
xiPt = 0;
2769-
xiXPull = 0;
2770-
xiYPull = 0;
2771-
xiZPull = 0;
2772-
2773-
xic0BachelorXDelta = 0;
2774-
xic0BachelorYDelta = 0;
2775-
xic0BachelorZDelta = 0;
2776-
xic0BachelorPt = 0;
2777-
xic0BachelorXPull = 0;
2778-
xic0BachelorYPull = 0;
2779-
xic0BachelorZPull = 0;
2780-
2781-
xic0XDelta = 0;
2782-
xic0YDelta = 0;
2783-
xic0ZDelta = 0;
2784-
xic0Pt = 0;
2785-
xic0XPull = 0;
2786-
xic0YPull = 0;
2787-
xic0ZPull = 0;
27882657

27892658
flag = 0;
27902659
debug = McMatchFlag::None;
@@ -2819,36 +2688,36 @@ struct HfCandidateCreatorXic0Omegac0Mc {
28192688
if (signV0 == 1) {
28202689
auto mcV0 = mcParticles.rawIteratorAt(indexRec - mcParticles.offset());
28212690

2822-
v0MassPull = (candidate.invMassLambda() - MassLambda0) / candidate.invMassV0Sig();
2691+
float v0MassPull = (candidate.invMassLambda() - MassLambda0) / candidate.invMassV0Sig();
28232692
registry.fill(HIST("hV0MassPullVsPt"), candidate.v0Pt(), v0MassPull);
28242693

2825-
v0DauPosXDelta = candidate.v0DauPosX() - mcV0DauPos.vx();
2826-
v0DauPosYDelta = candidate.v0DauPosY() - mcV0DauPos.vy();
2827-
v0DauPosZDelta = candidate.v0DauPosZ() - mcV0DauPos.vz();
2828-
v0DauPosPt = mcV0DauPos.pt();
2829-
v0DauPosXPull = v0DauPosXDelta / candidate.v0DauPosXError();
2830-
v0DauPosYPull = v0DauPosYDelta / candidate.v0DauPosYError();
2831-
v0DauPosZPull = v0DauPosZDelta / candidate.v0DauPosZError();
2832-
2833-
v0DauNegXDelta = candidate.v0DauNegX() - mcV0DauNeg.vx();
2834-
v0DauNegYDelta = candidate.v0DauNegY() - mcV0DauNeg.vy();
2835-
v0DauNegZDelta = candidate.v0DauNegZ() - mcV0DauNeg.vz();
2836-
v0DauNegPt = mcV0DauNeg.pt();
2837-
v0DauNegXPull = v0DauNegXDelta / candidate.v0DauNegXError();
2838-
v0DauNegYPull = v0DauNegYDelta / candidate.v0DauNegYError();
2839-
v0DauNegZPull = v0DauNegZDelta / candidate.v0DauNegZError();
2840-
2841-
v0XDelta = candidate.v0VtxX() - mcV0DauNeg.vx();
2842-
v0YDelta = candidate.v0VtxY() - mcV0DauNeg.vy();
2843-
v0ZDelta = candidate.v0VtxZ() - mcV0DauNeg.vz();
2844-
v0Pt = mcV0.pt();
2845-
v0XPull = v0XDelta / candidate.v0XError();
2846-
v0YPull = v0YDelta / candidate.v0YError();
2847-
v0ZPull = v0ZDelta / candidate.v0ZError();
2848-
2849-
lambdaXDelta = candidate.v0X() - mcV0DauNeg.vx();
2850-
lambdaYDelta = candidate.v0Y() - mcV0DauNeg.vy();
2851-
lambdaZDelta = candidate.v0Z() - mcV0DauNeg.vz();
2694+
float v0DauPosXDelta = candidate.v0DauPosX() - mcV0DauPos.vx();
2695+
float v0DauPosYDelta = candidate.v0DauPosY() - mcV0DauPos.vy();
2696+
float v0DauPosZDelta = candidate.v0DauPosZ() - mcV0DauPos.vz();
2697+
float v0DauPosPt = mcV0DauPos.pt();
2698+
float v0DauPosXPull = v0DauPosXDelta / candidate.v0DauPosXError();
2699+
float v0DauPosYPull = v0DauPosYDelta / candidate.v0DauPosYError();
2700+
float v0DauPosZPull = v0DauPosZDelta / candidate.v0DauPosZError();
2701+
2702+
float v0DauNegXDelta = candidate.v0DauNegX() - mcV0DauNeg.vx();
2703+
float v0DauNegYDelta = candidate.v0DauNegY() - mcV0DauNeg.vy();
2704+
float v0DauNegZDelta = candidate.v0DauNegZ() - mcV0DauNeg.vz();
2705+
float v0DauNegPt = mcV0DauNeg.pt();
2706+
float v0DauNegXPull = v0DauNegXDelta / candidate.v0DauNegXError();
2707+
float v0DauNegYPull = v0DauNegYDelta / candidate.v0DauNegYError();
2708+
float v0DauNegZPull = v0DauNegZDelta / candidate.v0DauNegZError();
2709+
2710+
float v0XDelta = candidate.v0VtxX() - mcV0DauNeg.vx();
2711+
float v0YDelta = candidate.v0VtxY() - mcV0DauNeg.vy();
2712+
float v0ZDelta = candidate.v0VtxZ() - mcV0DauNeg.vz();
2713+
float v0Pt = mcV0.pt();
2714+
float v0XPull = v0XDelta / candidate.v0XError();
2715+
float v0YPull = v0YDelta / candidate.v0YError();
2716+
float v0ZPull = v0ZDelta / candidate.v0ZError();
2717+
2718+
float lambdaXDelta = candidate.v0X() - mcV0DauNeg.vx();
2719+
float lambdaYDelta = candidate.v0Y() - mcV0DauNeg.vy();
2720+
float lambdaZDelta = candidate.v0Z() - mcV0DauNeg.vz();
28522721
registry.fill(HIST("hV0DauPosXDelta"), v0DauPosXDelta);
28532722
registry.fill(HIST("hV0DauPosYDelta"), v0DauPosYDelta);
28542723
registry.fill(HIST("hV0DauPosZDelta"), v0DauPosZDelta);
@@ -2900,30 +2769,30 @@ struct HfCandidateCreatorXic0Omegac0Mc {
29002769
if (indexRec > -1) {
29012770
// QA
29022771
if (signCasc == 1) {
2903-
xiMassPull = (candidate.invMassCascade() - MassXiMinus) / candidate.invMassXiSig();
2772+
float xiMassPull = (candidate.invMassCascade() - MassXiMinus) / candidate.invMassXiSig();
29042773
registry.fill(HIST("hXiMassPullVsPt"), candidate.xiPt(), xiMassPull);
29052774

2906-
xiBachelorXDelta = candidate.xiBachelorX() - mcXiBachelor.vx();
2907-
xiBachelorYDelta = candidate.xiBachelorY() - mcXiBachelor.vy();
2908-
xiBachelorZDelta = candidate.xiBachelorZ() - mcXiBachelor.vz();
2909-
xiBachelorPt = mcXiBachelor.pt();
2910-
xiBachelorXPull = xiBachelorXDelta / candidate.xiBachelorXError();
2911-
xiBachelorYPull = xiBachelorYDelta / candidate.xiBachelorYError();
2912-
xiBachelorZPull = xiBachelorZDelta / candidate.xiBachelorZError();
2775+
float xiBachelorXDelta = candidate.xiBachelorX() - mcXiBachelor.vx();
2776+
float xiBachelorYDelta = candidate.xiBachelorY() - mcXiBachelor.vy();
2777+
float xiBachelorZDelta = candidate.xiBachelorZ() - mcXiBachelor.vz();
2778+
float xiBachelorPt = mcXiBachelor.pt();
2779+
float xiBachelorXPull = xiBachelorXDelta / candidate.xiBachelorXError();
2780+
float xiBachelorYPull = xiBachelorYDelta / candidate.xiBachelorYError();
2781+
float xiBachelorZPull = xiBachelorZDelta / candidate.xiBachelorZError();
29132782

29142783
auto mcXi = mcParticles.rawIteratorAt(indexRec - mcParticles.offset());
29152784

2916-
xiXDelta = candidate.xiX() - mcXiBachelor.vx();
2917-
xiYDelta = candidate.xiY() - mcXiBachelor.vy();
2918-
xiZDelta = candidate.xiZ() - mcXiBachelor.vz();
2919-
xiPt = mcXi.pt();
2920-
xiXPull = xiXDelta / candidate.xiXError();
2921-
xiYPull = xiYDelta / candidate.xiYError();
2922-
xiZPull = xiZDelta / candidate.xiZError();
2785+
float xiXDelta = candidate.xiX() - mcXiBachelor.vx();
2786+
float xiYDelta = candidate.xiY() - mcXiBachelor.vy();
2787+
float xiZDelta = candidate.xiZ() - mcXiBachelor.vz();
2788+
float xiPt = mcXi.pt();
2789+
float xiXPull = xiXDelta / candidate.xiXError();
2790+
float xiYPull = xiYDelta / candidate.xiYError();
2791+
float xiZPull = xiZDelta / candidate.xiZError();
29232792

2924-
cascXDelta = candidate.xDecayVtxCascade() - mcXiBachelor.vx();
2925-
cascYDelta = candidate.yDecayVtxCascade() - mcXiBachelor.vy();
2926-
cascZDelta = candidate.zDecayVtxCascade() - mcXiBachelor.vz();
2793+
float cascXDelta = candidate.xDecayVtxCascade() - mcXiBachelor.vx();
2794+
float cascYDelta = candidate.yDecayVtxCascade() - mcXiBachelor.vy();
2795+
float cascZDelta = candidate.zDecayVtxCascade() - mcXiBachelor.vz();
29272796

29282797
registry.fill(HIST("hXiBachelorXDelta"), xiBachelorXDelta);
29292798
registry.fill(HIST("hXiBachelorYDelta"), xiBachelorYDelta);
@@ -2967,24 +2836,24 @@ struct HfCandidateCreatorXic0Omegac0Mc {
29672836

29682837
if (sign == 1) {
29692838
auto mcXic0 = mcParticles.rawIteratorAt(indexRec - mcParticles.offset());
2970-
xic0MassPull = (candidate.invMassCharmBaryon() - MassXiC0) / candidate.invMassXic0Sig();
2839+
float xic0MassPull = (candidate.invMassCharmBaryon() - MassXiC0) / candidate.invMassXic0Sig();
29712840
registry.fill(HIST("hXic0MassPullVsPt"), candidate.xic0Pt(), xic0MassPull);
29722841

2973-
xic0BachelorXDelta = candidate.xic0BachelorX() - mcXic0Bachelor.vx();
2974-
xic0BachelorYDelta = candidate.xic0BachelorY() - mcXic0Bachelor.vy();
2975-
xic0BachelorZDelta = candidate.xic0BachelorZ() - mcXic0Bachelor.vz();
2976-
xic0BachelorPt = mcXic0Bachelor.pt();
2977-
xic0BachelorXPull = xic0BachelorXDelta / candidate.xic0BachelorXError();
2978-
xic0BachelorYPull = xic0BachelorYDelta / candidate.xic0BachelorYError();
2979-
xic0BachelorZPull = xic0BachelorZDelta / candidate.xic0BachelorZError();
2980-
2981-
xic0XDelta = candidate.xDecayVtxCharmBaryon() - mcXic0Bachelor.vx();
2982-
xic0YDelta = candidate.yDecayVtxCharmBaryon() - mcXic0Bachelor.vy();
2983-
xic0ZDelta = candidate.zDecayVtxCharmBaryon() - mcXic0Bachelor.vz();
2984-
xic0Pt = mcXic0.pt();
2985-
xic0XPull = xic0XDelta / candidate.xic0XError();
2986-
xic0YPull = xic0YDelta / candidate.xic0YError();
2987-
xic0ZPull = xic0ZDelta / candidate.xic0ZError();
2842+
float xic0BachelorXDelta = candidate.xic0BachelorX() - mcXic0Bachelor.vx();
2843+
float xic0BachelorYDelta = candidate.xic0BachelorY() - mcXic0Bachelor.vy();
2844+
float xic0BachelorZDelta = candidate.xic0BachelorZ() - mcXic0Bachelor.vz();
2845+
float xic0BachelorPt = mcXic0Bachelor.pt();
2846+
float xic0BachelorXPull = xic0BachelorXDelta / candidate.xic0BachelorXError();
2847+
float xic0BachelorYPull = xic0BachelorYDelta / candidate.xic0BachelorYError();
2848+
float xic0BachelorZPull = xic0BachelorZDelta / candidate.xic0BachelorZError();
2849+
2850+
float xic0XDelta = candidate.xDecayVtxCharmBaryon() - mcXic0Bachelor.vx();
2851+
float xic0YDelta = candidate.yDecayVtxCharmBaryon() - mcXic0Bachelor.vy();
2852+
float xic0ZDelta = candidate.zDecayVtxCharmBaryon() - mcXic0Bachelor.vz();
2853+
float xic0Pt = mcXic0.pt();
2854+
float xic0XPull = xic0XDelta / candidate.xic0XError();
2855+
float xic0YPull = xic0YDelta / candidate.xic0YError();
2856+
float xic0ZPull = xic0ZDelta / candidate.xic0ZError();
29882857
registry.fill(HIST("hXic0BachelorXDelta"), xic0BachelorXDelta);
29892858
registry.fill(HIST("hXic0BachelorYDelta"), xic0BachelorYDelta);
29902859
registry.fill(HIST("hXic0BachelorZDelta"), xic0BachelorZDelta);

0 commit comments

Comments
 (0)