Skip to content

DPG: add the flag for the specific track selection #1022

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
zhangbiao-phy wants to merge 2 commits into from
Closed

DPG: add the flag for the specific track selection #1022

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
618e8fc
DPG: add the flag for the specific track selection
zhangbiao-phy Jul 19, 2022
acb6049
DPG: add the flag for the specific track selection
zhangbiao-phy Jul 19, 2022
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
23 changes: 16 additions & 7 deletions DPG/Tasks/qaEfficiency.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -31,7 +31,9 @@
#include "TEfficiency.h"
#include "TList.h"

using namespace o2;
using namespace o2::framework;
using namespace o2::framework::expressions;

struct QaEfficiency {
// Particle information
Expand Down Expand Up @@ -63,8 +65,8 @@ struct QaEfficiency {
Configurable<bool> doTr{"do-tr", false, "Flag to run with the PDG code of tritons"};
Configurable<bool> doHe{"do-he", false, "Flag to run with the PDG code of helium 3"};
Configurable<bool> doAl{"do-al", false, "Flag to run with the PDG code of helium 4"};
// Track only selection
Configurable<bool> applyTrackSelection{"applyTrackSelection", false, "Flag to use the standard track selection when selecting numerator and denominator"};
// options to select only specific tracks
Configurable<int> trackSelection{"trackSelection", 1, "Track selection: 0 -> No Cut, 1 -> kGlobalTrack, 2 -> kGlobalTrackWoPtEta, 3 -> kGlobalTrackWoDCA, 4 -> kQualityTracks, 5 -> kInAcceptanceTracks"};
// Event selection
Configurable<int> nMinNumberOfContributors{"nMinNumberOfContributors", 2, "Minimum required number of contributors to the primary vertex"};
Configurable<float> vertexZMin{"vertex-z-min", -10.f, "Minimum position of the generated vertez in Z (cm)"};
Expand Down Expand Up @@ -1094,10 +1096,17 @@ struct QaEfficiency {
isCollisionSelected<true>(collision);
}

Filter trackFilter = (trackSelection.node() == 0) ||
((trackSelection.node() == 1) && requireGlobalTrackInFilter()) ||
((trackSelection.node() == 2) && requireGlobalTrackWoPtEtaInFilter()) ||
((trackSelection.node() == 3) && requireGlobalTrackWoDCAInFilter()) ||
((trackSelection.node() == 4) && requireQualityTracksInFilter()) ||
((trackSelection.node() == 5) && requireTrackCutInFilter(TrackSelectionFlags::kInAcceptanceTracks));
// MC process
using TrackTableMC = soa::Join<o2::aod::Tracks, o2::aod::TracksExtra, o2::aod::McTrackLabels, o2::aod::TrackSelection>;
void processMC(const o2::aod::McParticles& mcParticles,
const o2::soa::Join<o2::aod::Collisions, o2::aod::McCollisionLabels, o2::aod::EvSels>& collisions,
const o2::soa::Join<o2::aod::Tracks, o2::aod::TracksExtra, o2::aod::McTrackLabels, o2::aod::TrackSelection>& tracks,
soa::Filtered<TrackTableMC> const& tracks,
const o2::aod::McCollisions&)
{

Expand Down Expand Up @@ -1165,7 +1174,7 @@ struct QaEfficiency {
}

histos.fill(HIST("MC/trackSelection"), 10);
if (applyTrackSelection && !track.isGlobalTrack()) { // Check general cuts
if (trackSelection < 0) { // Check general cuts
continue;
}
histos.fill(HIST("MC/trackSelection"), 11);
Expand Down Expand Up @@ -1322,9 +1331,9 @@ struct QaEfficiency {
}
}
PROCESS_SWITCH(QaEfficiency, processMC, "process MC", false);

using TrackTableData = soa::Join<o2::aod::Tracks, o2::aod::TracksExtra, o2::aod::TrackSelection>;
void processData(o2::soa::Join<o2::aod::Collisions, o2::aod::EvSels>::iterator const& collision,
const o2::soa::Join<o2::aod::Tracks, o2::aod::TracksExtra, o2::aod::TrackSelection>& tracks)
soa::Filtered<TrackTableData> const& tracks)
{

if (!isCollisionSelected<false>(collision)) {
Expand All @@ -1347,7 +1356,7 @@ struct QaEfficiency {
continue;
}
histos.fill(HIST("Data/trackSelection"), 4);
if (applyTrackSelection) { // Check general cuts
if (trackSelection > 0) { // Check general cuts
if (!track.passedTrackType()) {
continue;
}
Expand Down