From 8ad9631a5682bf23bbc851dec6c5f01b97c7f4f6 Mon Sep 17 00:00:00 2001 From: Raphaelle Bailhache Date: Wed, 20 Jul 2022 13:07:33 +0200 Subject: [PATCH 1/2] Some more cuts, signal definition and number of ITS clusters for skimmed data for Run 3 --- PWGDQ/Core/CutsLibrary.h | 13 +++++++++++++ PWGDQ/Core/MCSignalLibrary.h | 13 +++++++++++++ PWGDQ/Core/VarManager.h | 2 +- 3 files changed, 27 insertions(+), 1 deletion(-) diff --git a/PWGDQ/Core/CutsLibrary.h b/PWGDQ/Core/CutsLibrary.h index 0628d37104c..27746658441 100644 --- a/PWGDQ/Core/CutsLibrary.h +++ b/PWGDQ/Core/CutsLibrary.h @@ -284,6 +284,13 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName) return cut; } + if (!nameStr.compare("lmee_TPCTrackRun3_lowPt")) { + cut->AddCut(GetAnalysisCut("lmeeLowBKine")); + cut->AddCut(GetAnalysisCut("TightTPCTrackRun3")); + cut->AddCut(GetAnalysisCut("standardPrimaryTrack")); + return cut; + } + if (!nameStr.compare("lmee_GlobalTrackRun3_TPC_ePID_lowPt")) { cut->AddCut(GetAnalysisCut("lmeeLowBKine")); cut->AddCut(GetAnalysisCut("TightGlobalTrackRun3")); @@ -460,6 +467,12 @@ AnalysisCut* o2::aod::dqcuts::GetAnalysisCut(const char* cutName) return cut; } + if (!nameStr.compare("TightTPCTrackRun3")) { + cut->AddCut(VarManager::kTPCchi2, 0.0, 4.0); + cut->AddCut(VarManager::kTPCnclsCR, 80.0, 161.); + return cut; + } + if (!nameStr.compare("electronStandardQuality")) { cut->AddCut(VarManager::kIsSPDany, 0.5, 1.5); cut->AddCut(VarManager::kIsITSrefit, 0.5, 1.5); diff --git a/PWGDQ/Core/MCSignalLibrary.h b/PWGDQ/Core/MCSignalLibrary.h index 58994d06f2e..a5949543588 100644 --- a/PWGDQ/Core/MCSignalLibrary.h +++ b/PWGDQ/Core/MCSignalLibrary.h @@ -45,6 +45,12 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name) signal = new MCSignal(name, "Primary electrons", {prong}, {-1}); // define the signal using the full constructor return signal; } + if (!nameStr.compare("pionPrimary")) { + MCProng prong(1, {211}, {true}, {false}, {0}, {0}, {false}); // define 1-generation prong using the full constructor + prong.SetSourceBit(0, MCProng::kPhysicalPrimary); // set source to be ALICE primary particles + signal = new MCSignal(name, "Primary electrons", {prong}, {-1}); // define the signal using the full constructor + return signal; + } if (!nameStr.compare("photon")) { MCProng prong(1, {22}, {true}, {false}, {0}, {0}, {false}); // define 1-generation prong using the full constructor signal = new MCSignal(name, "Photon", {prong}, {-1}); // define the signal using the full constructor @@ -357,6 +363,13 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name) signal = new MCSignal(name, "Electrons from open charmed hadron decays from b hadron decays", {prong}, {-1}); return signal; } + if (!nameStr.compare("HFdecayToE")) { + MCProng prong(2, {902, 11}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}); + // prong.SetSourceBit(0, MCProng::kPhysicalPrimary, false); // set source to be ALICE primary particles + prong.SetSignalInTime(true); // set direction to check for daughters (true, in time) or for mothers (false, back in time) + signal = new MCSignal(name, "Open charm and beauty to electrons", {prong}, {-1}); + return signal; + } // if (!nameStr.compare("LFQtoPC")) { // MCProng prong(3, {900, 22, 11}, {true, true, true}, {false, false, false}, {0, 0, 0}, {0, 0, 0}, {false, false, false}); // prong.SetSignalInTime(true); // set direction to check for daughters (true, in time) or for mothers (false, back in time) diff --git a/PWGDQ/Core/VarManager.h b/PWGDQ/Core/VarManager.h index 3b46db921be..ccf75bb8e13 100644 --- a/PWGDQ/Core/VarManager.h +++ b/PWGDQ/Core/VarManager.h @@ -667,7 +667,7 @@ void VarManager::FillTrack(T const& track, float* values) if constexpr ((fillMap & ReducedTrackBarrel) > 0) { if (fgUsedVars[kITSncls]) { values[kITSncls] = 0.0; - for (int i = 0; i < 6; ++i) { + for (int i = 0; i < 7; ++i) { values[kITSncls] += ((track.itsClusterMap() & (1 << i)) ? 1 : 0); } } From 3a54519fd799cf50cf34badbdfedc18f7c0bf821 Mon Sep 17 00:00:00 2001 From: Raphaelle Bailhache Date: Wed, 20 Jul 2022 13:30:14 +0200 Subject: [PATCH 2/2] Fix --- PWGDQ/Core/MCSignalLibrary.h | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/PWGDQ/Core/MCSignalLibrary.h b/PWGDQ/Core/MCSignalLibrary.h index a5949543588..eaf21491dfd 100644 --- a/PWGDQ/Core/MCSignalLibrary.h +++ b/PWGDQ/Core/MCSignalLibrary.h @@ -46,7 +46,7 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name) return signal; } if (!nameStr.compare("pionPrimary")) { - MCProng prong(1, {211}, {true}, {false}, {0}, {0}, {false}); // define 1-generation prong using the full constructor + MCProng prong(1, {211}, {true}, {false}, {0}, {0}, {false}); // define 1-generation prong using the full constructor prong.SetSourceBit(0, MCProng::kPhysicalPrimary); // set source to be ALICE primary particles signal = new MCSignal(name, "Primary electrons", {prong}, {-1}); // define the signal using the full constructor return signal; @@ -366,8 +366,8 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name) if (!nameStr.compare("HFdecayToE")) { MCProng prong(2, {902, 11}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}); // prong.SetSourceBit(0, MCProng::kPhysicalPrimary, false); // set source to be ALICE primary particles - prong.SetSignalInTime(true); // set direction to check for daughters (true, in time) or for mothers (false, back in time) - signal = new MCSignal(name, "Open charm and beauty to electrons", {prong}, {-1}); + prong.SetSignalInTime(true); // set direction to check for daughters (true, in time) or for mothers (false, back in time) + signal = new MCSignal(name, "Open charm and beauty to electrons", {prong}, {-1}); return signal; } // if (!nameStr.compare("LFQtoPC")) {