|
13 | 13 | /// \brief Provides efficiency extraction and QC for track cuts and PID |
14 | 14 | /// \author victor.gonzalez.sebastian@gmail.com |
15 | 15 |
|
16 | | -#include <TH2F.h> |
17 | | -#include <TProfile2D.h> |
18 | | -#include <THnSparse.h> |
19 | | -#include <TPDGCode.h> |
20 | | -#include <TMCProcess.h> |
21 | | -#include <CCDB/BasicCCDBManager.h> |
22 | | -#include <vector> |
23 | | -#include <algorithm> |
24 | | -#include <cstdio> |
25 | | -#include <memory> |
26 | | -#include <string> |
27 | | -#include "ReconstructionDataFormats/PID.h" |
28 | | -#include "Common/Core/TrackSelection.h" |
29 | | -#include "Common/Core/TableHelper.h" |
| 16 | +#include "PWGCF/Core/AnalysisConfigurableCuts.h" |
| 17 | +#include "PWGCF/DataModel/DptDptFiltered.h" |
| 18 | +#include "PWGCF/TableProducer/dptDptFilter.h" |
| 19 | + |
30 | 20 | #include "Common/Core/RecoDecay.h" |
31 | | -#include "Common/DataModel/TrackSelectionTables.h" |
| 21 | +#include "Common/Core/TableHelper.h" |
| 22 | +#include "Common/Core/TrackSelection.h" |
32 | 23 | #include "Common/DataModel/PIDResponse.h" |
| 24 | +#include "Common/DataModel/TrackSelectionTables.h" |
| 25 | + |
33 | 26 | #include "Framework/ASoAHelpers.h" |
34 | | -#include "Framework/Expressions.h" |
35 | 27 | #include "Framework/AnalysisDataModel.h" |
36 | 28 | #include "Framework/AnalysisTask.h" |
37 | | -#include "Framework/runDataProcessing.h" |
| 29 | +#include "Framework/Expressions.h" |
38 | 30 | #include "Framework/RunningWorkflowInfo.h" |
| 31 | +#include "Framework/runDataProcessing.h" |
| 32 | +#include "ReconstructionDataFormats/PID.h" |
| 33 | +#include <CCDB/BasicCCDBManager.h> |
| 34 | + |
39 | 35 | #include "Math/MatrixFunctions.h" |
40 | 36 | #include "Math/SMatrix.h" |
| 37 | +#include <TH2F.h> |
| 38 | +#include <THnSparse.h> |
| 39 | +#include <TMCProcess.h> |
| 40 | +#include <TPDGCode.h> |
| 41 | +#include <TProfile2D.h> |
41 | 42 |
|
42 | | -#include "PWGCF/Core/AnalysisConfigurableCuts.h" |
43 | | -#include "PWGCF/DataModel/DptDptFiltered.h" |
44 | | -#include "PWGCF/TableProducer/dptDptFilter.h" |
| 43 | +#include <algorithm> |
| 44 | +#include <cstdio> |
| 45 | +#include <memory> |
| 46 | +#include <string> |
| 47 | +#include <vector> |
45 | 48 |
|
46 | 49 | using namespace o2; |
47 | 50 | using namespace o2::framework; |
@@ -741,41 +744,41 @@ struct PidDataCollectingEngine { |
741 | 744 | if constexpr (kindOfData == kReco) { |
742 | 745 | /* PID histograms */ |
743 | 746 | std::vector<std::string> whenname{"Before", "After"}; |
744 | | - constexpr char whenprefix[kNoOfSteps]{'B', 'A'}; |
| 747 | + constexpr char kWhenPrefix[kNoOfSteps]{'B', 'A'}; |
745 | 748 | std::vector<std::string> whentitle{"before", ""}; |
746 | 749 | for (uint ix = 0; ix < kNoOfSteps; ++ix) { |
747 | 750 | fhTPCdEdxSignalVsP[ix] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
748 | | - HNAMESTRING("tpcSignalVsP%c", whenprefix[ix]), |
| 751 | + HNAMESTRING("tpcSignalVsP%c", kWhenPrefix[ix]), |
749 | 752 | HTITLESTRING("TPC dE/dx signal %s", whentitle[ix].c_str()), kTH2F, {pidPAxis, dEdxAxis}); |
750 | 753 | fhTOFSignalVsP[ix] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
751 | | - HNAMESTRING("tofSignalVsP%c", whenprefix[ix]), |
| 754 | + HNAMESTRING("tofSignalVsP%c", kWhenPrefix[ix]), |
752 | 755 | HTITLESTRING("TOF signal %s", whentitle[ix].c_str()), |
753 | 756 | kTH2F, {pidPAxis, {200, 0.0, 1.1, "#beta"}}); |
754 | 757 | fhPvsTOFSqMass[ix] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
755 | | - HNAMESTRING("tofPvsMassSq%c", whenprefix[ix]), |
| 758 | + HNAMESTRING("tofPvsMassSq%c", kWhenPrefix[ix]), |
756 | 759 | HTITLESTRING("Momentum versus #it{m}^{2} %s", whentitle[ix].c_str()), |
757 | 760 | kTH2F, {{140, 0.0, 1.4, "#it{m}^{2} ((GeV/c^{2})^{2})"}, pidPAxis}); |
758 | 761 | for (uint isp = 0; isp < nmainsp; ++isp) { |
759 | 762 | fhTPCdEdxSignalDiffVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
760 | | - HNAMESTRING("tpcSignalDiffVsP%c_%s", whenprefix[ix], mainspnames[isp].c_str()), |
| 763 | + HNAMESTRING("tpcSignalDiffVsP%c_%s", kWhenPrefix[ix], mainspnames[isp].c_str()), |
761 | 764 | HTITLESTRING("TPC dE/dx to the %s line %s", mainsptitles[isp].c_str(), whentitle[ix].c_str()), |
762 | 765 | kTH2F, {pidPAxis, {400, -200.0, 200.0, FORMATSTRING("dE/dx - <dE/dx>_{%s}", mainsptitles[isp].c_str())}}); |
763 | 766 | fhTOFSignalDiffVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
764 | | - HNAMESTRING("tofSignalDiffVsP%c_%s", whenprefix[ix], mainspnames[isp].c_str()), |
| 767 | + HNAMESTRING("tofSignalDiffVsP%c_%s", kWhenPrefix[ix], mainspnames[isp].c_str()), |
765 | 768 | HTITLESTRING("#Delta^{TOF_{%s}} %s", mainsptitles[isp].c_str(), whentitle[ix].c_str()), |
766 | 769 | kTH2F, {pidPAxis, {200, -1000.0, 1000.0, FORMATSTRING("t-t_{ev}-t_{exp_{%s}} (ps)", mainsptitles[isp].c_str())}}); |
767 | 770 | fhTPCTOFSigmaVsP[ix][isp] = ADDHISTOGRAM(TH3, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
768 | | - HNAMESTRING("toftpcNSigmasVsP%c_%s", whenprefix[ix], mainspnames[isp].c_str()), |
| 771 | + HNAMESTRING("toftpcNSigmasVsP%c_%s", kWhenPrefix[ix], mainspnames[isp].c_str()), |
769 | 772 | HTITLESTRING("n#sigma to the %s line %s", mainsptitles[isp].c_str(), whentitle[ix].c_str()), |
770 | 773 | kTH3F, {pidPAxis, {noOfNSigmaBins, minNSigma, maxNSigma, FORMATSTRING("n#sigma_{TPC}^{%s}", mainsptitles[isp].c_str())}, {120, -6.0, 6.0, FORMATSTRING("n#sigma_{TOF}^{%s}", mainsptitles[isp].c_str())}}); |
771 | 774 | } |
772 | 775 | for (uint isp = 0; isp < nallmainsp; ++isp) { |
773 | 776 | fhTPCnSigmasVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
774 | | - HNAMESTRING("tpcNSigmasVsP%c_%s", whenprefix[ix], allmainspnames[isp].c_str()), |
| 777 | + HNAMESTRING("tpcNSigmasVsP%c_%s", kWhenPrefix[ix], allmainspnames[isp].c_str()), |
775 | 778 | HTITLESTRING("TPC n#sigma to the %s line %s", allmainsptitles[isp].c_str(), whentitle[ix].c_str()), |
776 | 779 | kTH2F, {pidPAxis, {noOfNSigmaBins, minNSigma, maxNSigma, FORMATSTRING("n#sigma_{TPC}^{%s}", allmainsptitles[isp].c_str())}}); |
777 | 780 | fhTOFnSigmasVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()), |
778 | | - HNAMESTRING("tofNSigmasVsP%c_%s", whenprefix[ix], allmainspnames[isp].c_str()), |
| 781 | + HNAMESTRING("tofNSigmasVsP%c_%s", kWhenPrefix[ix], allmainspnames[isp].c_str()), |
779 | 782 | HTITLESTRING("TOF n#sigma to the %s line %s", allmainsptitles[isp].c_str(), whentitle[ix].c_str()), |
780 | 783 | kTH2F, {pidPAxis, {noOfNSigmaBins, minNSigma, maxNSigma, FORMATSTRING("n#sigma_{TOF}^{%s}", allmainsptitles[isp].c_str())}}); |
781 | 784 | } |
|
0 commit comments