|
35 | 35 | #include "Common/CCDB/RCTSelectionFlags.h" |
36 | 36 | #include "Common/CCDB/ctpRateFetcher.h" |
37 | 37 | #include "Common/Core/RecoDecay.h" |
38 | | -#include "Common/DataModel/Multiplicity.h" |
39 | 38 | #include "Common/DataModel/Centrality.h" |
40 | 39 | #include "Common/DataModel/EventSelection.h" |
| 40 | +#include "Common/DataModel/Multiplicity.h" |
41 | 41 |
|
42 | 42 | #include <CCDB/BasicCCDBManager.h> |
43 | 43 | #include <CCDB/CcdbApi.h> |
44 | | -#include <DetectorsBase/MatLayerCylSet.h> |
45 | | -#include <DetectorsBase/Propagator.h> |
46 | 44 | #include <CommonConstants/MathConstants.h> |
47 | 45 | #include <CommonConstants/PhysicsConstants.h> |
48 | 46 | #include <DCAFitter/FwdDCAFitterN.h> |
49 | 47 | #include <DataFormatsParameters/GRPMagField.h> |
| 48 | +#include <DetectorsBase/MatLayerCylSet.h> |
| 49 | +#include <DetectorsBase/Propagator.h> |
50 | 50 | #include <Framework/AnalysisDataModel.h> |
51 | 51 | #include <Framework/AnalysisHelpers.h> |
52 | 52 | #include <Framework/AnalysisTask.h> |
|
59 | 59 | #include <Framework/runDataProcessing.h> |
60 | 60 | #include <ReconstructionDataFormats/TrackFwd.h> |
61 | 61 |
|
62 | | -#include <Math/SVector.h> |
63 | | -#include <Math/SMatrix.h> |
64 | 62 | #include <Math/MatrixRepresentationsStatic.h> |
| 63 | +#include <Math/SMatrix.h> |
| 64 | +#include <Math/SVector.h> |
65 | 65 | #include <TH1.h> |
66 | 66 | #include <TH2.h> |
| 67 | +#include <TObject.h> |
67 | 68 | #include <TPDGCode.h> |
68 | 69 | #include <TVector3.h> |
69 | | -#include <TObject.h> |
70 | 70 |
|
71 | 71 | #include <array> |
72 | 72 | #include <cmath> |
@@ -230,7 +230,7 @@ struct forwardlambdakzeroanalysis { |
230 | 230 | } ccdbConfigurations; |
231 | 231 |
|
232 | 232 | o2::ccdb::CcdbApi ccdbApi; |
233 | | - Service<o2::ccdb::BasicCCDBManager> ccdb {}; |
| 233 | + Service<o2::ccdb::BasicCCDBManager> ccdb{}; |
234 | 234 | ctpRateFetcher rateFetcher; |
235 | 235 | int mRunNumber = 0; |
236 | 236 | float magField = 0.; |
@@ -1218,20 +1218,20 @@ struct forwardlambdakzeroanalysis { |
1218 | 1218 | BITSET(bitMap, selZmax); |
1219 | 1219 | } |
1220 | 1220 | // DCA proton and pion to PV for Lambda and AntiLambda decay hypotheses |
1221 | | - if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaPosToPVxy) < o2::track::DefaultDCA) && |
| 1221 | + if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaPosToPVxy) < o2::track::DefaultDCA) && |
1222 | 1222 | std::fabs(v0.dcaPosToPVxy) > v0Selections.dcaPosToPVxy) { |
1223 | 1223 | BITSET(bitMap, selDCAPosToPVxy); |
1224 | 1224 | } |
1225 | | - if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaNegToPVxy) < o2::track::DefaultDCA) && |
| 1225 | + if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaNegToPVxy) < o2::track::DefaultDCA) && |
1226 | 1226 | std::fabs(v0.dcaNegToPVxy) > v0Selections.dcaNegToPVxy) { |
1227 | 1227 | BITSET(bitMap, selDCANegToPVxy); |
1228 | 1228 | } |
1229 | 1229 | // DCA proton and pion to PV for Lambda and AntiLambda decay hypotheses |
1230 | | - if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaPosToPVz) < o2::track::DefaultDCA) && |
| 1230 | + if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaPosToPVz) < o2::track::DefaultDCA) && |
1231 | 1231 | std::fabs(v0.dcaPosToPVz) > v0Selections.dcaPosToPVz) { |
1232 | 1232 | BITSET(bitMap, selDCAPosToPVz); |
1233 | 1233 | } |
1234 | | - if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaNegToPVz) < o2::track::DefaultDCA) && |
| 1234 | + if ((!v0Selections.rejectFailedPropagation || std::fabs(v0.dcaNegToPVz) < o2::track::DefaultDCA) && |
1235 | 1235 | std::fabs(v0.dcaNegToPVz) > v0Selections.dcaNegToPVz) { |
1236 | 1236 | BITSET(bitMap, selDCANegToPVz); |
1237 | 1237 | } |
@@ -1461,18 +1461,18 @@ struct forwardlambdakzeroanalysis { |
1461 | 1461 | histos.fill(HIST("hPositiveMFTchi2"), v0.posChi2); |
1462 | 1462 | histos.fill(HIST("hNegativeMFTchi2"), v0.negChi2); |
1463 | 1463 | if (doprocessMonteCarlo) { |
1464 | | - histos.fill(HIST("hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt())/v0.positivePtMc()); |
1465 | | - histos.fill(HIST("hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt())/v0.negativePtMc()); |
1466 | | - histos.fill(HIST("hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz())/v0.positivePzMc()); |
1467 | | - histos.fill(HIST("hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz())/v0.negativePzMc()); |
| 1464 | + histos.fill(HIST("hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt()) / v0.positivePtMc()); |
| 1465 | + histos.fill(HIST("hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt()) / v0.negativePtMc()); |
| 1466 | + histos.fill(HIST("hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz()) / v0.positivePzMc()); |
| 1467 | + histos.fill(HIST("hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz()) / v0.negativePzMc()); |
1468 | 1468 |
|
1469 | 1469 | histos.fill(HIST("hV0XDiff"), v0.xMc, (v0.xMc - v0.X)); |
1470 | 1470 | histos.fill(HIST("hV0YDiff"), v0.yMc, (v0.yMc - v0.Y)); |
1471 | 1471 | histos.fill(HIST("hV0ZDiff"), v0.zMc, (v0.zMc - v0.Z)); |
1472 | | - histos.fill(HIST("hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius)/v0.RadiusMc); |
1473 | | - histos.fill(HIST("hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist)/v0.ZdistMc); |
1474 | | - histos.fill(HIST("hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT)/v0.pTMc); |
1475 | | - histos.fill(HIST("hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ)/v0.pZMc); |
| 1472 | + histos.fill(HIST("hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius) / v0.RadiusMc); |
| 1473 | + histos.fill(HIST("hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist) / v0.ZdistMc); |
| 1474 | + histos.fill(HIST("hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT) / v0.pTMc); |
| 1475 | + histos.fill(HIST("hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ) / v0.pZMc); |
1476 | 1476 | } |
1477 | 1477 | } |
1478 | 1478 |
|
@@ -1538,18 +1538,18 @@ struct forwardlambdakzeroanalysis { |
1538 | 1538 | histos.fill(HIST("K0Short/hPositiveMFTchi2"), v0.posChi2); |
1539 | 1539 | histos.fill(HIST("K0Short/hNegativeMFTchi2"), v0.negChi2); |
1540 | 1540 | if (doprocessMonteCarlo) { |
1541 | | - histos.fill(HIST("K0Short/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt())/v0.positivePtMc()); |
1542 | | - histos.fill(HIST("K0Short/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt())/v0.negativePtMc()); |
1543 | | - histos.fill(HIST("K0Short/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz())/v0.positivePzMc()); |
1544 | | - histos.fill(HIST("K0Short/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz())/v0.negativePzMc()); |
| 1541 | + histos.fill(HIST("K0Short/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt()) / v0.positivePtMc()); |
| 1542 | + histos.fill(HIST("K0Short/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt()) / v0.negativePtMc()); |
| 1543 | + histos.fill(HIST("K0Short/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz()) / v0.positivePzMc()); |
| 1544 | + histos.fill(HIST("K0Short/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz()) / v0.negativePzMc()); |
1545 | 1545 |
|
1546 | 1546 | histos.fill(HIST("K0Short/hV0XDiff"), v0.xMc, (v0.xMc - v0.X)); |
1547 | 1547 | histos.fill(HIST("K0Short/hV0YDiff"), v0.yMc, (v0.yMc - v0.Y)); |
1548 | 1548 | histos.fill(HIST("K0Short/hV0ZDiff"), v0.zMc, (v0.zMc - v0.Z)); |
1549 | | - histos.fill(HIST("K0Short/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius)/v0.RadiusMc); |
1550 | | - histos.fill(HIST("K0Short/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist)/v0.ZdistMc); |
1551 | | - histos.fill(HIST("K0Short/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT)/v0.pTMc); |
1552 | | - histos.fill(HIST("K0Short/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ)/v0.pZMc); |
| 1549 | + histos.fill(HIST("K0Short/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius) / v0.RadiusMc); |
| 1550 | + histos.fill(HIST("K0Short/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist) / v0.ZdistMc); |
| 1551 | + histos.fill(HIST("K0Short/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT) / v0.pTMc); |
| 1552 | + histos.fill(HIST("K0Short/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ) / v0.pZMc); |
1553 | 1553 | } |
1554 | 1554 | } |
1555 | 1555 | nK0Shorts++; |
@@ -1604,18 +1604,18 @@ struct forwardlambdakzeroanalysis { |
1604 | 1604 | histos.fill(HIST("Lambda/hPositiveMFTchi2"), v0.posChi2); |
1605 | 1605 | histos.fill(HIST("Lambda/hNegativeMFTchi2"), v0.negChi2); |
1606 | 1606 | if (doprocessMonteCarlo) { |
1607 | | - histos.fill(HIST("Lambda/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt())/v0.positivePtMc()); |
1608 | | - histos.fill(HIST("Lambda/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt())/v0.negativePtMc()); |
1609 | | - histos.fill(HIST("Lambda/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz())/v0.positivePzMc()); |
1610 | | - histos.fill(HIST("Lambda/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz())/v0.negativePzMc()); |
| 1607 | + histos.fill(HIST("Lambda/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt()) / v0.positivePtMc()); |
| 1608 | + histos.fill(HIST("Lambda/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt()) / v0.negativePtMc()); |
| 1609 | + histos.fill(HIST("Lambda/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz()) / v0.positivePzMc()); |
| 1610 | + histos.fill(HIST("Lambda/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz()) / v0.negativePzMc()); |
1611 | 1611 |
|
1612 | 1612 | histos.fill(HIST("Lambda/hV0XDiff"), v0.xMc, (v0.xMc - v0.X)); |
1613 | 1613 | histos.fill(HIST("Lambda/hV0YDiff"), v0.yMc, (v0.yMc - v0.Y)); |
1614 | 1614 | histos.fill(HIST("Lambda/hV0ZDiff"), v0.zMc, (v0.zMc - v0.Z)); |
1615 | | - histos.fill(HIST("Lambda/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius)/v0.RadiusMc); |
1616 | | - histos.fill(HIST("Lambda/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist)/v0.ZdistMc); |
1617 | | - histos.fill(HIST("Lambda/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT)/v0.pTMc); |
1618 | | - histos.fill(HIST("Lambda/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ)/v0.pZMc); |
| 1615 | + histos.fill(HIST("Lambda/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius) / v0.RadiusMc); |
| 1616 | + histos.fill(HIST("Lambda/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist) / v0.ZdistMc); |
| 1617 | + histos.fill(HIST("Lambda/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT) / v0.pTMc); |
| 1618 | + histos.fill(HIST("Lambda/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ) / v0.pZMc); |
1619 | 1619 | } |
1620 | 1620 | } |
1621 | 1621 | nLambdas++; |
@@ -1670,18 +1670,18 @@ struct forwardlambdakzeroanalysis { |
1670 | 1670 | histos.fill(HIST("AntiLambda/hPositiveMFTchi2"), v0.posChi2); |
1671 | 1671 | histos.fill(HIST("AntiLambda/hNegativeMFTchi2"), v0.negChi2); |
1672 | 1672 | if (doprocessMonteCarlo) { |
1673 | | - histos.fill(HIST("AntiLambda/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt())/v0.positivePtMc()); |
1674 | | - histos.fill(HIST("AntiLambda/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt())/v0.negativePtMc()); |
1675 | | - histos.fill(HIST("AntiLambda/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz())/v0.positivePzMc()); |
1676 | | - histos.fill(HIST("AntiLambda/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz())/v0.negativePzMc()); |
| 1673 | + histos.fill(HIST("AntiLambda/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt()) / v0.positivePtMc()); |
| 1674 | + histos.fill(HIST("AntiLambda/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt()) / v0.negativePtMc()); |
| 1675 | + histos.fill(HIST("AntiLambda/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz()) / v0.positivePzMc()); |
| 1676 | + histos.fill(HIST("AntiLambda/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz()) / v0.negativePzMc()); |
1677 | 1677 |
|
1678 | 1678 | histos.fill(HIST("AntiLambda/hV0XDiff"), v0.xMc, (v0.xMc - v0.X)); |
1679 | 1679 | histos.fill(HIST("AntiLambda/hV0YDiff"), v0.yMc, (v0.yMc - v0.Y)); |
1680 | 1680 | histos.fill(HIST("AntiLambda/hV0ZDiff"), v0.zMc, (v0.zMc - v0.Z)); |
1681 | | - histos.fill(HIST("AntiLambda/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius)/v0.RadiusMc); |
1682 | | - histos.fill(HIST("AntiLambda/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist)/v0.Zdist); |
1683 | | - histos.fill(HIST("AntiLambda/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT)/v0.pTMc); |
1684 | | - histos.fill(HIST("AntiLambda/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ)/v0.pZMc); |
| 1681 | + histos.fill(HIST("AntiLambda/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius) / v0.RadiusMc); |
| 1682 | + histos.fill(HIST("AntiLambda/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist) / v0.Zdist); |
| 1683 | + histos.fill(HIST("AntiLambda/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT) / v0.pTMc); |
| 1684 | + histos.fill(HIST("AntiLambda/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ) / v0.pZMc); |
1685 | 1685 | } |
1686 | 1686 | } |
1687 | 1687 | nAntiLambdas++; |
@@ -1736,18 +1736,18 @@ struct forwardlambdakzeroanalysis { |
1736 | 1736 | histos.fill(HIST("D0/hPositiveMFTchi2"), v0.posChi2); |
1737 | 1737 | histos.fill(HIST("D0/hNegativeMFTchi2"), v0.negChi2); |
1738 | 1738 | if (doprocessMonteCarlo) { |
1739 | | - histos.fill(HIST("D0/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt())/v0.positivePtMc()); |
1740 | | - histos.fill(HIST("D0/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt())/v0.negativePtMc()); |
1741 | | - histos.fill(HIST("D0/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz())/v0.positivePzMc()); |
1742 | | - histos.fill(HIST("D0/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz())/v0.negativePzMc()); |
| 1739 | + histos.fill(HIST("D0/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt()) / v0.positivePtMc()); |
| 1740 | + histos.fill(HIST("D0/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt()) / v0.negativePtMc()); |
| 1741 | + histos.fill(HIST("D0/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz()) / v0.positivePzMc()); |
| 1742 | + histos.fill(HIST("D0/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz()) / v0.negativePzMc()); |
1743 | 1743 |
|
1744 | 1744 | histos.fill(HIST("D0/hV0XDiff"), v0.xMc, (v0.xMc - v0.X)); |
1745 | 1745 | histos.fill(HIST("D0/hV0YDiff"), v0.yMc, (v0.yMc - v0.Y)); |
1746 | 1746 | histos.fill(HIST("D0/hV0ZDiff"), v0.zMc, (v0.zMc - v0.Z)); |
1747 | | - histos.fill(HIST("D0/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius)/v0.RadiusMc); |
1748 | | - histos.fill(HIST("D0/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist)/v0.ZdistMc); |
1749 | | - histos.fill(HIST("D0/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT)/v0.pTMc); |
1750 | | - histos.fill(HIST("D0/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ)/v0.pZMc); |
| 1747 | + histos.fill(HIST("D0/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius) / v0.RadiusMc); |
| 1748 | + histos.fill(HIST("D0/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist) / v0.ZdistMc); |
| 1749 | + histos.fill(HIST("D0/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT) / v0.pTMc); |
| 1750 | + histos.fill(HIST("D0/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ) / v0.pZMc); |
1751 | 1751 | } |
1752 | 1752 | } |
1753 | 1753 | nD0s++; |
@@ -1802,18 +1802,18 @@ struct forwardlambdakzeroanalysis { |
1802 | 1802 | histos.fill(HIST("AntiD0/hPositiveMFTchi2"), v0.posChi2); |
1803 | 1803 | histos.fill(HIST("AntiD0/hNegativeMFTchi2"), v0.negChi2); |
1804 | 1804 | if (doprocessMonteCarlo) { |
1805 | | - histos.fill(HIST("AntiD0/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt())/v0.positivePtMc()); |
1806 | | - histos.fill(HIST("AntiD0/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt())/v0.negativePtMc()); |
1807 | | - histos.fill(HIST("AntiD0/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz())/v0.positivePzMc()); |
1808 | | - histos.fill(HIST("AntiD0/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz())/v0.negativePzMc()); |
| 1805 | + histos.fill(HIST("AntiD0/hPositivePtResolution"), v0.positivePtMc(), (v0.positivePtMc() - v0.positivePt()) / v0.positivePtMc()); |
| 1806 | + histos.fill(HIST("AntiD0/hNegativePtResolution"), v0.negativePtMc(), (v0.negativePtMc() - v0.negativePt()) / v0.negativePtMc()); |
| 1807 | + histos.fill(HIST("AntiD0/hPositivePzResolution"), v0.positivePzMc(), (v0.positivePzMc() - v0.positivePz()) / v0.positivePzMc()); |
| 1808 | + histos.fill(HIST("AntiD0/hNegativePzResolution"), v0.negativePzMc(), (v0.negativePzMc() - v0.negativePz()) / v0.negativePzMc()); |
1809 | 1809 |
|
1810 | 1810 | histos.fill(HIST("AntiD0/hV0XDiff"), v0.xMc, (v0.xMc - v0.X)); |
1811 | 1811 | histos.fill(HIST("AntiD0/hV0YDiff"), v0.yMc, (v0.yMc - v0.Y)); |
1812 | 1812 | histos.fill(HIST("AntiD0/hV0ZDiff"), v0.zMc, (v0.zMc - v0.Z)); |
1813 | | - histos.fill(HIST("AntiD0/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius)/v0.RadiusMc); |
1814 | | - histos.fill(HIST("AntiD0/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist)/v0.ZdistMc); |
1815 | | - histos.fill(HIST("AntiD0/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT)/v0.pTMc); |
1816 | | - histos.fill(HIST("AntiD0/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ)/v0.pZMc); |
| 1813 | + histos.fill(HIST("AntiD0/hV0RadiusDiff"), v0.RadiusMc, (v0.RadiusMc - v0.Radius) / v0.RadiusMc); |
| 1814 | + histos.fill(HIST("AntiD0/hV0DistZDiff"), v0.ZdistMc, (v0.ZdistMc - v0.Zdist) / v0.ZdistMc); |
| 1815 | + histos.fill(HIST("AntiD0/hV0PtResolution"), v0.pTMc, (v0.pTMc - v0.pT) / v0.pTMc); |
| 1816 | + histos.fill(HIST("AntiD0/hV0PzResolution"), v0.pZMc, (v0.pZMc - v0.pZ) / v0.pZMc); |
1817 | 1817 | } |
1818 | 1818 | } |
1819 | 1819 | nAntiD0s++; |
@@ -2215,7 +2215,7 @@ struct forwardlambdakzeroanalysis { |
2215 | 2215 | dca[0] = o2::track::DefaultDCA; |
2216 | 2216 | dca[1] = o2::track::DefaultDCA; |
2217 | 2217 | dca[2] = o2::track::DefaultDCA; |
2218 | | - |
| 2218 | + |
2219 | 2219 | // propagate track to vertex including MCS effects if material budget included, simple propagation to Z otherwise |
2220 | 2220 | float x2x0 = 0; |
2221 | 2221 | auto mb = lut->getMatBudget(t.getX(), t.getY(), t.getZ(), collision.posX(), collision.posY(), collision.posZ()); |
@@ -2562,7 +2562,7 @@ struct forwardlambdakzeroanalysis { |
2562 | 2562 | // ______________________________________________________ |
2563 | 2563 | // Simulated processing (subscribes to MC information too) |
2564 | 2564 | template <typename TCollision, typename TMCCollisions, typename TMFTTracks, typename TTracks, typename TBCs, typename TMCParticles> |
2565 | | - void analyzeRecoedV0sInMonteCarlo(TCollision const& collision, TMCCollisions const& mcCollisions,TMFTTracks const& mftTracks, TTracks const& besttracks, TBCs const& bcs, TMCParticles const& mcParticles) |
| 2565 | + void analyzeRecoedV0sInMonteCarlo(TCollision const& collision, TMCCollisions const& mcCollisions, TMFTTracks const& mftTracks, TTracks const& besttracks, TBCs const& bcs, TMCParticles const& mcParticles) |
2566 | 2566 | { |
2567 | 2567 | // Fire up CCDB |
2568 | 2568 | initCCDB(bcs, collision); |
|
0 commit comments