Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
41 changes: 22 additions & 19 deletions PWGCF/TableProducer/dptDptFilter.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -13,40 +13,43 @@
/// \brief Filters collisions and tracks according to selection criteria
/// \author victor.gonzalez.sebastian@gmail.com

#include <cmath>
#include <algorithm>
#include <string>
#include <vector>
#include "PWGCF/TableProducer/dptDptFilter.h"

#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"

#include "Framework/AnalysisTask.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/ASoAHelpers.h"
#include "CommonConstants/PhysicsConstants.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/Centrality.h"
#include "Common/Core/TableHelper.h"
#include "Common/Core/TrackSelection.h"
#include "Common/Core/TrackSelectionDefaults.h"
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/CollisionAssociationTables.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/PIDResponse.h"
#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "Common/DataModel/TrackSelectionTables.h"
#include "Common/DataModel/CollisionAssociationTables.h"

#include "CommonConstants/PhysicsConstants.h"
#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/O2DatabasePDGPlugin.h"
#include "Framework/runDataProcessing.h"
#include "Framework/RunningWorkflowInfo.h"
#include <TROOT.h>
#include <TPDGCode.h>
#include <TParameter.h>
#include <TList.h>
#include "Framework/runDataProcessing.h"

#include <TDirectory.h>
#include <TFolder.h>
#include <TH1.h>
#include <TH2.h>
#include <TH3.h>
#include <TList.h>
#include <TPDGCode.h>
#include <TParameter.h>
#include <TProfile3D.h>
#include <TROOT.h>

#include "PWGCF/TableProducer/dptDptFilter.h"
#include <algorithm>
#include <cmath>
#include <string>
#include <vector>

using namespace o2;
using namespace o2::framework;
Expand Down
43 changes: 23 additions & 20 deletions PWGCF/TableProducer/dptDptFilter.h
Original file line number Diff line number Diff line change
Expand Up @@ -16,33 +16,36 @@
#ifndef PWGCF_TABLEPRODUCER_DPTDPTFILTER_H_
#define PWGCF_TABLEPRODUCER_DPTDPTFILTER_H_

#include "PWGCF/Core/AnalysisConfigurableCuts.h"

#include "Common/Core/RecoDecay.h"
#include "Common/Core/TrackSelection.h"
#include "Common/Core/TrackSelectionDefaults.h"
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/Multiplicity.h"
#include "Common/DataModel/TrackSelectionTables.h"

#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "ReconstructionDataFormats/PID.h"
#include <CCDB/BasicCCDBManager.h>
#include <TPDGCode.h>
#include <TMCProcess.h>

#include <TF1.h>
#include <TList.h>
#include <vector>
#include <TMCProcess.h>
#include <TPDGCode.h>

#include <bitset>
#include <string>
#include <iomanip>
#include <fstream>
#include <locale>
#include <sstream>
#include <ranges>
#include <functional>
#include <iomanip>
#include <locale>
#include <map>

#include "ReconstructionDataFormats/PID.h"
#include "Framework/AnalysisTask.h"
#include "Framework/AnalysisDataModel.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/Multiplicity.h"
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/TrackSelectionTables.h"
#include "Common/Core/RecoDecay.h"
#include "Common/Core/TrackSelection.h"
#include "Common/Core/TrackSelectionDefaults.h"
#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include <ranges>
#include <sstream>
#include <string>
#include <vector>

namespace o2
{
Expand Down
51 changes: 28 additions & 23 deletions PWGCF/Tasks/dptDptCorrelations.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,26 @@
/// \brief implements two-particle correlations base data collection
/// \author victor.gonzalez.sebastian@gmail.com

#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/Core/PairCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

#include "Common/Core/RecoDecay.h"
#include "Common/Core/TableHelper.h"
#include "Common/Core/TrackSelection.h"
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/TrackSelectionTables.h"

#include "DataFormatsParameters/GRPObject.h"
#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/RunningWorkflowInfo.h"
#include "Framework/runDataProcessing.h"
#include <CCDB/BasicCCDBManager.h>

#include <TDirectory.h>
#include <TFolder.h>
#include <TH1.h>
Expand All @@ -24,28 +43,12 @@
#include <TProfile3D.h>
#include <TROOT.h>
#include <TVector2.h>

#include <cmath>
#include <cstdio>
#include <ctime>
#include <string>
#include <vector>
#include <cmath>
#include <ctime>

#include "Common/Core/TrackSelection.h"
#include "Common/Core/TableHelper.h"
#include "Common/Core/RecoDecay.h"
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/TrackSelectionTables.h"
#include "DataFormatsParameters/GRPObject.h"
#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/runDataProcessing.h"
#include "Framework/RunningWorkflowInfo.h"
#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/Core/PairCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

using namespace o2;
using namespace o2::framework;
Expand Down Expand Up @@ -279,15 +282,17 @@ struct DptDptCorrelations {
photon = p1+p2;
photon.M()*/

float tantheta1 = 1e10;
constexpr float kLARGETANTHETA = 1e10;
constexpr float kVERYSMALLETA = 1e-10;
float tantheta1 = kLARGETANTHETA;

if (track1.eta() < -1e-10 || track1.eta() > 1e-10) {
if (track1.eta() < -kVERYSMALLETA || track1.eta() > kVERYSMALLETA) {
float expTmp = std::exp(-track1.eta());
tantheta1 = 2.0 * expTmp / (1.0 - expTmp * expTmp);
}

float tantheta2 = 1e10;
if (track2.eta() < -1e-10 || track2.eta() > 1e-10) {
float tantheta2 = kLARGETANTHETA;
if (track2.eta() < -kVERYSMALLETA || track2.eta() > kVERYSMALLETA) {
float expTmp = std::exp(-track2.eta());
tantheta2 = 2.0 * expTmp / (1.0 - expTmp * expTmp);
}
Expand Down
9 changes: 5 additions & 4 deletions PWGCF/Tasks/dptDptFilterQa.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -13,15 +13,16 @@
/// \brief basic checks for the behavior of the filter task
/// \author victor.gonzalez.sebastian@gmail.com

#include <cmath>
#include <string>
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/runDataProcessing.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

#include <cmath>
#include <string>

using namespace o2;
using namespace o2::framework;
Expand Down
17 changes: 10 additions & 7 deletions PWGCF/Tasks/matchRecoGen.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -13,25 +13,24 @@
/// \brief basic check for the matching between generator level and detector level
/// \author victor.gonzalez.sebastian@gmail.com

#include <cmath>
#include <string>
#include <vector>
#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

#include "Common/Core/RecoDecay.h"
#include "Common/Core/TrackSelection.h"
#include "Common/Core/TrackSelectionDefaults.h"
#include "Common/Core/RecoDecay.h"
#include "Common/DataModel/Centrality.h"
#include "Common/DataModel/EventSelection.h"
#include "Common/DataModel/PIDResponse.h"
#include "Common/DataModel/TrackSelectionTables.h"

#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/O2DatabasePDGPlugin.h"
#include "Framework/runDataProcessing.h"
#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

#include <TDirectory.h>
#include <TFolder.h>
#include <TH1.h>
Expand All @@ -42,6 +41,10 @@
#include <TProfile3D.h>
#include <TROOT.h>

#include <cmath>
#include <string>
#include <vector>

using namespace o2;
using namespace o2::framework;
using namespace o2::soa;
Expand Down
61 changes: 32 additions & 29 deletions PWGCF/TwoParticleCorrelations/Tasks/dptDptEfficiencyAndQc.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -13,35 +13,38 @@
/// \brief Provides efficiency extraction and QC for track cuts and PID
/// \author victor.gonzalez.sebastian@gmail.com

#include <TH2F.h>
#include <TProfile2D.h>
#include <THnSparse.h>
#include <TPDGCode.h>
#include <TMCProcess.h>
#include <CCDB/BasicCCDBManager.h>
#include <vector>
#include <algorithm>
#include <cstdio>
#include <memory>
#include <string>
#include "ReconstructionDataFormats/PID.h"
#include "Common/Core/TrackSelection.h"
#include "Common/Core/TableHelper.h"
#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

#include "Common/Core/RecoDecay.h"
#include "Common/DataModel/TrackSelectionTables.h"
#include "Common/Core/TableHelper.h"
#include "Common/Core/TrackSelection.h"
#include "Common/DataModel/PIDResponse.h"
#include "Common/DataModel/TrackSelectionTables.h"

#include "Framework/ASoAHelpers.h"
#include "Framework/Expressions.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/runDataProcessing.h"
#include "Framework/Expressions.h"
#include "Framework/RunningWorkflowInfo.h"
#include "Framework/runDataProcessing.h"
#include "ReconstructionDataFormats/PID.h"
#include <CCDB/BasicCCDBManager.h>

#include "Math/MatrixFunctions.h"
#include "Math/SMatrix.h"
#include <TH2F.h>
#include <THnSparse.h>
#include <TMCProcess.h>
#include <TPDGCode.h>
#include <TProfile2D.h>

#include "PWGCF/Core/AnalysisConfigurableCuts.h"
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"
#include <algorithm>
#include <cstdio>
#include <memory>
#include <string>
#include <vector>

using namespace o2;
using namespace o2::framework;
Expand Down Expand Up @@ -741,41 +744,41 @@ struct PidDataCollectingEngine {
if constexpr (kindOfData == kReco) {
/* PID histograms */
std::vector<std::string> whenname{"Before", "After"};
constexpr char whenprefix[kNoOfSteps]{'B', 'A'};
constexpr char kWhenPrefix[kNoOfSteps]{'B', 'A'};
std::vector<std::string> whentitle{"before", ""};
for (uint ix = 0; ix < kNoOfSteps; ++ix) {
fhTPCdEdxSignalVsP[ix] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("tpcSignalVsP%c", whenprefix[ix]),
HNAMESTRING("tpcSignalVsP%c", kWhenPrefix[ix]),
HTITLESTRING("TPC dE/dx signal %s", whentitle[ix].c_str()), kTH2F, {pidPAxis, dEdxAxis});
fhTOFSignalVsP[ix] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("tofSignalVsP%c", whenprefix[ix]),
HNAMESTRING("tofSignalVsP%c", kWhenPrefix[ix]),
HTITLESTRING("TOF signal %s", whentitle[ix].c_str()),
kTH2F, {pidPAxis, {200, 0.0, 1.1, "#beta"}});
fhPvsTOFSqMass[ix] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("tofPvsMassSq%c", whenprefix[ix]),
HNAMESTRING("tofPvsMassSq%c", kWhenPrefix[ix]),
HTITLESTRING("Momentum versus #it{m}^{2} %s", whentitle[ix].c_str()),
kTH2F, {{140, 0.0, 1.4, "#it{m}^{2} ((GeV/c^{2})^{2})"}, pidPAxis});
for (uint isp = 0; isp < nmainsp; ++isp) {
fhTPCdEdxSignalDiffVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("tpcSignalDiffVsP%c_%s", whenprefix[ix], mainspnames[isp].c_str()),
HNAMESTRING("tpcSignalDiffVsP%c_%s", kWhenPrefix[ix], mainspnames[isp].c_str()),
HTITLESTRING("TPC dE/dx to the %s line %s", mainsptitles[isp].c_str(), whentitle[ix].c_str()),
kTH2F, {pidPAxis, {400, -200.0, 200.0, FORMATSTRING("dE/dx - <dE/dx>_{%s}", mainsptitles[isp].c_str())}});
fhTOFSignalDiffVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("tofSignalDiffVsP%c_%s", whenprefix[ix], mainspnames[isp].c_str()),
HNAMESTRING("tofSignalDiffVsP%c_%s", kWhenPrefix[ix], mainspnames[isp].c_str()),
HTITLESTRING("#Delta^{TOF_{%s}} %s", mainsptitles[isp].c_str(), whentitle[ix].c_str()),
kTH2F, {pidPAxis, {200, -1000.0, 1000.0, FORMATSTRING("t-t_{ev}-t_{exp_{%s}} (ps)", mainsptitles[isp].c_str())}});
fhTPCTOFSigmaVsP[ix][isp] = ADDHISTOGRAM(TH3, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("toftpcNSigmasVsP%c_%s", whenprefix[ix], mainspnames[isp].c_str()),
HNAMESTRING("toftpcNSigmasVsP%c_%s", kWhenPrefix[ix], mainspnames[isp].c_str()),
HTITLESTRING("n#sigma to the %s line %s", mainsptitles[isp].c_str(), whentitle[ix].c_str()),
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())}});
}
for (uint isp = 0; isp < nallmainsp; ++isp) {
fhTPCnSigmasVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("tpcNSigmasVsP%c_%s", whenprefix[ix], allmainspnames[isp].c_str()),
HNAMESTRING("tpcNSigmasVsP%c_%s", kWhenPrefix[ix], allmainspnames[isp].c_str()),
HTITLESTRING("TPC n#sigma to the %s line %s", allmainsptitles[isp].c_str(), whentitle[ix].c_str()),
kTH2F, {pidPAxis, {noOfNSigmaBins, minNSigma, maxNSigma, FORMATSTRING("n#sigma_{TPC}^{%s}", allmainsptitles[isp].c_str())}});
fhTOFnSigmasVsP[ix][isp] = ADDHISTOGRAM(TH2, DIRECTORYSTRING("%s/%s/%s", dirname, "PID", whenname[ix].c_str()),
HNAMESTRING("tofNSigmasVsP%c_%s", whenprefix[ix], allmainspnames[isp].c_str()),
HNAMESTRING("tofNSigmasVsP%c_%s", kWhenPrefix[ix], allmainspnames[isp].c_str()),
HTITLESTRING("TOF n#sigma to the %s line %s", allmainsptitles[isp].c_str(), whentitle[ix].c_str()),
kTH2F, {pidPAxis, {noOfNSigmaBins, minNSigma, maxNSigma, FORMATSTRING("n#sigma_{TOF}^{%s}", allmainsptitles[isp].c_str())}});
}
Expand Down
12 changes: 6 additions & 6 deletions PWGCF/TwoParticleCorrelations/Tasks/dptDptPerRunExtraQc.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -13,18 +13,18 @@
/// \brief basic per run check of the per analyzed species p vs TPC IW momentum
/// \author victor.gonzalez.sebastian@gmail.com

#include <array>
#include <cmath>
#include <unordered_map>
#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"

#include "Framework/ASoAHelpers.h"
#include "Framework/AnalysisDataModel.h"
#include "Framework/AnalysisTask.h"
#include "Framework/ASoAHelpers.h"
#include "Framework/HistogramRegistry.h"
#include "Framework/runDataProcessing.h"

#include "PWGCF/DataModel/DptDptFiltered.h"
#include "PWGCF/TableProducer/dptDptFilter.h"
#include <array>
#include <cmath>
#include <unordered_map>

using namespace o2;
using namespace o2::framework;
Expand Down
Loading
Loading