@@ -224,17 +224,13 @@ struct V0PtInvMassPlots {
224224 const auto & posDaughterTrack = v0.posTrack_as <DaughterTracks>();
225225 const auto & negDaughterTrack = v0.negTrack_as <DaughterTracks>();
226226 // Armenteros-Podolandski Plot Values
227- Double_t pv0 = std::sqrt ((v0.px () * v0.px ()) + (v0.py () * v0.py ()) + (v0.pz () * v0.pz ()));
228- Double_t pposdauparallelv0 = ((v0.posTrack_as <DaughterTracks>().px () * v0.px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz ())) / pv0;
229- Double_t pnegdauparallelv0 = ((v0.negTrack_as <DaughterTracks>().px () * v0.px ()) + (v0.negTrack_as <DaughterTracks>().py () * v0.py ()) + (v0.negTrack_as <DaughterTracks>().pz () * v0.pz ())) / pv0;
230- Double_t qValueposdau = std::sqrt (((v0.posTrack_as <DaughterTracks>().px () * v0.posTrack_as <DaughterTracks>().px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.posTrack_as <DaughterTracks>().py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.posTrack_as <DaughterTracks>().pz ())) - (pposdauparallelv0 * pposdauparallelv0));
231- Double_t qValuenegdau = std::sqrt (((v0.negTrack_as <DaughterTracks>().px () * v0.negTrack_as <DaughterTracks>().px ()) + (v0.negTrack_as <DaughterTracks>().py () * v0.negTrack_as <DaughterTracks>().py ()) + (v0.negTrack_as <DaughterTracks>().pz () * v0.negTrack_as <DaughterTracks>().pz ())) - (pnegdauparallelv0 * pnegdauparallelv0));
232- // Double_t qValue = std::sqrt(qValueposdau + qValuenegdau);
233- Double_t plpos = (v0.posTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.posTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
234- Double_t plneg = (v0.negTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.negTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.negTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
235- Double_t aValue = (plpos - plneg) / (plpos + plneg);
236- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlot" ), aValue, qValueposdau);
237- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlot" ), aValue, qValuenegdau);
227+ double pv0 = std::sqrt ((v0.px () * v0.px ()) + (v0.py () * v0.py ()) + (v0.pz () * v0.pz ()));
228+ double pposdauparallelv0 = ((v0.posTrack_as <DaughterTracks>().px () * v0.px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz ())) / pv0;
229+ double qValue = std::sqrt (((v0.posTrack_as <DaughterTracks>().px () * v0.posTrack_as <DaughterTracks>().px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.posTrack_as <DaughterTracks>().py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.posTrack_as <DaughterTracks>().pz ())) - (pposdauparallelv0 * pposdauparallelv0));
230+ double plpos = (v0.posTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.posTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
231+ double plneg = (v0.negTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.negTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.negTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
232+ double aValue = (plpos - plneg) / (plpos + plneg);
233+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlot" ), aValue, qValue);
238234 rPtAnalysis.fill (HIST (" hVertexZRec" ), collision.posZ ());
239235 rPtAnalysis.fill (HIST (" hV0PtAll" ), v0.pt ());
240236 // Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)
@@ -257,8 +253,7 @@ struct V0PtInvMassPlots {
257253 rPtAnalysis.fill (HIST (" hK0shEtaNegDau" ), v0.negTrack_as <DaughterTracks>().eta ());
258254 rPtAnalysis.fill (HIST (" hK0shEtaDaughters" ), v0.negTrack_as <DaughterTracks>().eta ());
259255 rPtAnalysis.fill (HIST (" hK0shEtaDaughters" ), v0.posTrack_as <DaughterTracks>().eta ());
260- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotK0Short" ), aValue, qValueposdau);
261- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotK0Short" ), aValue, qValuenegdau);
256+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotK0Short" ), aValue, qValue);
262257 for (int i = 0 ; i < 20 ; i++) {
263258 // getting the pt value in #_# for and converting it to a number #.# for use, we get two values which correspond to the range of each bin
264259 std::string pt1 = pthistos::kaonptbins[i]; // getting the lower string-value of the bin
@@ -290,8 +285,7 @@ struct V0PtInvMassPlots {
290285 rPtAnalysis.fill (HIST (" hLambdaReconstructedPtSpectrum" ), v0.pt ());
291286 rPtAnalysis.fill (HIST (" hLambdaEtaDaughters" ), v0.negTrack_as <DaughterTracks>().eta ());
292287 rPtAnalysis.fill (HIST (" hLambdaEtaDaughters" ), v0.posTrack_as <DaughterTracks>().eta ());
293- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotLambda" ), aValue, qValueposdau);
294- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotLambda" ), aValue, qValuenegdau);
288+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotLambda" ), aValue, qValue);
295289 for (int i = 0 ; i < 20 ; i++) {
296290 // same as above with kzerosh we fill the 20 lambda namespace histograms within their Pt range
297291 std::string pt1 = pthistos::lambdaPtBins[i];
@@ -322,8 +316,7 @@ struct V0PtInvMassPlots {
322316 rPtAnalysis.fill (HIST (" hAntilambdaReconstructedPtSpectrum" ), v0.pt ());
323317 rPtAnalysis.fill (HIST (" hAntiLambdaEtaDaughters" ), v0.negTrack_as <DaughterTracks>().eta ());
324318 rPtAnalysis.fill (HIST (" hAntiLambdaEtaDaughters" ), v0.posTrack_as <DaughterTracks>().eta ());
325- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotAntiLambda" ), aValue, qValueposdau);
326- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotAntiLambda" ), aValue, qValuenegdau);
319+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotAntiLambda" ), aValue, qValue);
327320 for (int i = 0 ; i < 20 ; i++) {
328321 // same as above with kzerosh and lambda we fill the 20 anti-lambda namespace histograms within their Pt range
329322 std::string pt1 = pthistos::antiLambdaPtBins[i];
@@ -357,17 +350,13 @@ struct V0PtInvMassPlots {
357350 const auto & posDaughterTrack = v0.posTrack_as <DaughterTracks>();
358351 const auto & negDaughterTrack = v0.negTrack_as <DaughterTracks>();
359352 // Armenteros-Podolandski Plot Values
360- Double_t pv0 = std::sqrt ((v0.px () * v0.px ()) + (v0.py () * v0.py ()) + (v0.pz () * v0.pz ()));
361- Double_t pposdauparallelv0 = ((v0.posTrack_as <DaughterTracks>().px () * v0.px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz ())) / pv0;
362- Double_t pnegdauparallelv0 = ((v0.negTrack_as <DaughterTracks>().px () * v0.px ()) + (v0.negTrack_as <DaughterTracks>().py () * v0.py ()) + (v0.negTrack_as <DaughterTracks>().pz () * v0.pz ())) / pv0;
363- Double_t qValueposdau = std::sqrt (((v0.posTrack_as <DaughterTracks>().px () * v0.posTrack_as <DaughterTracks>().px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.posTrack_as <DaughterTracks>().py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.posTrack_as <DaughterTracks>().pz ())) - (pposdauparallelv0 * pposdauparallelv0));
364- Double_t qValuenegdau = std::sqrt (((v0.negTrack_as <DaughterTracks>().px () * v0.negTrack_as <DaughterTracks>().px ()) + (v0.negTrack_as <DaughterTracks>().py () * v0.negTrack_as <DaughterTracks>().py ()) + (v0.negTrack_as <DaughterTracks>().pz () * v0.negTrack_as <DaughterTracks>().pz ())) - (pnegdauparallelv0 * pnegdauparallelv0));
365- // Double_t qValue = std::sqrt(qValueposdau + qValuenegdau);
366- Double_t plpos = (v0.posTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.posTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
367- Double_t plneg = (v0.negTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.negTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.negTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
368- Double_t aValue = (plpos - plneg) / (plpos + plneg);
369- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlot" ), aValue, qValueposdau);
370- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlot" ), aValue, qValuenegdau);
353+ double pv0 = std::sqrt ((v0.px () * v0.px ()) + (v0.py () * v0.py ()) + (v0.pz () * v0.pz ()));
354+ double pposdauparallelv0 = ((v0.posTrack_as <DaughterTracks>().px () * v0.px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz ())) / pv0;
355+ double qValue = std::sqrt (((v0.posTrack_as <DaughterTracks>().px () * v0.posTrack_as <DaughterTracks>().px ()) + (v0.posTrack_as <DaughterTracks>().py () * v0.posTrack_as <DaughterTracks>().py ()) + (v0.posTrack_as <DaughterTracks>().pz () * v0.posTrack_as <DaughterTracks>().pz ())) - (pposdauparallelv0 * pposdauparallelv0));
356+ double plpos = (v0.posTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.posTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.posTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
357+ double plneg = (v0.negTrack_as <DaughterTracks>().px () * v0.px () / pv0) + (v0.negTrack_as <DaughterTracks>().py () * v0.py () / pv0) + (v0.negTrack_as <DaughterTracks>().pz () * v0.pz () / pv0);
358+ double aValue = (plpos - plneg) / (plpos + plneg);
359+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlot" ), aValue, qValue);
371360 rPtAnalysis.fill (HIST (" hVertexZRec" ), collision.posZ ());
372361 rPtAnalysis.fill (HIST (" hV0PtAll" ), v0.pt ());
373362 if (std::abs (v0.posTrack_as <DaughterTracks>().eta ()) < etadau && std::abs (v0.negTrack_as <DaughterTracks>().eta ()) < etadau) { // daughters pseudorapidity cut
@@ -392,8 +381,7 @@ struct V0PtInvMassPlots {
392381 rPtAnalysis.fill (HIST (" hK0shEtaNegDau" ), v0.negTrack_as <DaughterTracks>().eta ());
393382 rPtAnalysis.fill (HIST (" hK0shEtaDaughters" ), v0.negTrack_as <DaughterTracks>().eta ());
394383 rPtAnalysis.fill (HIST (" hK0shEtaDaughters" ), v0.posTrack_as <DaughterTracks>().eta ());
395- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotK0Short" ), aValue, qValueposdau);
396- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotK0Short" ), aValue, qValuenegdau);
384+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotK0Short" ), aValue, qValue);
397385 for (int i = 0 ; i < 20 ; i++) { // same as above MC-process we fill the namespace histos with the kaon invariant mass of the particle within the pt range of the histo
398386 std::string pt1 = pthistos::kaonptbins[i];
399387 std::string pt2 = pthistos::kaonptbins[i + 1 ];
@@ -426,8 +414,7 @@ struct V0PtInvMassPlots {
426414 rPtAnalysis.fill (HIST (" hMassLambdaAllAfterCuts" ), v0.mLambda ());
427415 rPtAnalysis.fill (HIST (" hLambdaEtaDaughters" ), v0.negTrack_as <DaughterTracks>().eta ());
428416 rPtAnalysis.fill (HIST (" hLambdaEtaDaughters" ), v0.posTrack_as <DaughterTracks>().eta ());
429- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotLambda" ), aValue, qValueposdau);
430- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotLambda" ), aValue, qValuenegdau);
417+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotLambda" ), aValue, qValue);
431418 for (int i = 0 ; i < 20 ; i++) { // same as above MC-process we fill the namespace histos with the lambda invariant mass of the particle within the pt range of the histo
432419 std::string pt1 = pthistos::lambdaPtBins[i];
433420 std::string pt2 = pthistos::lambdaPtBins[i + 1 ];
@@ -460,8 +447,7 @@ struct V0PtInvMassPlots {
460447 rPtAnalysis.fill (HIST (" hMassAntilambdaAllAfterCuts" ), v0.mAntiLambda ());
461448 rPtAnalysis.fill (HIST (" hAntiLambdaEtaDaughters" ), v0.negTrack_as <DaughterTracks>().eta ());
462449 rPtAnalysis.fill (HIST (" hAntiLambdaEtaDaughters" ), v0.posTrack_as <DaughterTracks>().eta ());
463- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotAntiLambda" ), aValue, qValueposdau);
464- rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotAntiLambda" ), aValue, qValuenegdau);
450+ rPtAnalysis.fill (HIST (" hArmenterosPodolanskiPlotAntiLambda" ), aValue, qValue);
465451 for (int i = 0 ; i < 20 ; i++) { // same as above MC-process we fill the namespace histos with the antilambda invariant mass of the particle within the pt range of the histo
466452 std::string pt1 = pthistos::antiLambdaPtBins[i];
467453 std::string pt2 = pthistos::antiLambdaPtBins[i + 1 ];
0 commit comments