Skip to content

Commit 12f4425

Browse files
committed
New selections on eta and rigidity
Adds new configurable parameters for the He3 Lambda analysis: - `etaMax`: Maximum He3 pseudorapidity - `minTPCrigidity`: Minimum He3 rigidity These parameters are used to further refine the He3 selection criteria in the analysis.
1 parent ed48a3c commit 12f4425

File tree

1 file changed

+4
-5
lines changed

1 file changed

+4
-5
lines changed

PWGLF/TableProducer/Nuspex/he3LambdaAnalysis.cxx

Lines changed: 4 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -162,6 +162,8 @@ struct he3LambdaAnalysis {
162162
std::string prefix = "cfgHe3";
163163
Configurable<float> ptMin{"ptMin", 1.0f, "Minimum He3 pT"};
164164
Configurable<float> ptMax{"ptMax", 10.0f, "Maximum He3 pT"};
165+
Configurable<float> etaMax{"etaMax", 0.8f, "Maximum He3 pseudorapidity"};
166+
Configurable<float> minTPCrigidity{"minTPCrigidity", 0.5f, "Minimum He3 rigidity"};
165167
Configurable<float> nSigmaTPCMax{"nSigmaTPCMax", 4.0f, "Maximum He3 TPC nSigma"};
166168
Configurable<float> dcaxyMax{"dcaxyMax", 0.5f, "Maximum He3 DCA xy"};
167169
Configurable<float> dcazMax{"dcazMax", 0.5f, "Maximum He3 DCA z"};
@@ -311,14 +313,11 @@ struct he3LambdaAnalysis {
311313
}
312314
hTPCsignalAll->Fill(track.tpcInnerParam() * track.sign(), track.tpcSignal());
313315
const float pt = track.pt() * 2.0f;
314-
if (pt < cfgHe3.ptMin || pt > cfgHe3.ptMax) {
315-
continue; // Skip tracks outside pT range
316-
}
317316
float expTPCSignal = o2::tpc::BetheBlochAleph(track.tpcInnerParam() * 2.0f / constants::physics::MassHelium3, mBBparamsHe[0], mBBparamsHe[1], mBBparamsHe[2], mBBparamsHe[3], mBBparamsHe[4]);
318317
double nSigmaTPC = (track.tpcSignal() - expTPCSignal) / (expTPCSignal * mBBparamsHe[5]);
319318
hTPCnSigmaAll->Fill(track.tpcInnerParam() * track.sign(), nSigmaTPC);
320-
if (std::abs(nSigmaTPC) > cfgHe3.nSigmaTPCMax) {
321-
continue; // Skip tracks with TPC nSigma outside range
319+
if (pt < cfgHe3.ptMin || pt > cfgHe3.ptMax || std::abs(track.eta()) > cfgHe3.etaMax || track.tpcInnerParam() < cfgHe3.minTPCrigidity || std::abs(nSigmaTPC) > cfgHe3.nSigmaTPCMax) {
320+
continue; // Skip tracks outside He3 PID+kinematics selection criteria
322321
}
323322
setTrackParCov(track, trackParCov);
324323
std::array<float, 2> dcaInfo;

0 commit comments

Comments
 (0)