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
4 changes: 2 additions & 2 deletions PWGLF/DataModel/LFHStrangeCorrelationTables.h
Original file line number Diff line number Diff line change
Expand Up @@ -41,9 +41,9 @@ DECLARE_SOA_TABLE(TriggerTracks, "AOD", "TRIGGERTRACKS", o2::soa::Index<>, trigg
/// Table for storing assoc track indices
namespace assocPions
{
DECLARE_SOA_INDEX_COLUMN(Collision, collision); //!
DECLARE_SOA_INDEX_COLUMN(Collision, collision); //!
DECLARE_SOA_INDEX_COLUMN_FULL(Track, track, int, Tracks, "_Assoc"); //!
} // namespace triggerTracks
} // namespace assocPions
DECLARE_SOA_TABLE(AssocPions, "AOD", "ASSOCPIONS", o2::soa::Index<>, assocPions::CollisionId, assocPions::TrackId);
/// _________________________________________
/// Table for storing associated V0 indices
Expand Down
53 changes: 29 additions & 24 deletions PWGLF/TableProducer/hStrangeCorrelationFilter.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -203,34 +203,39 @@ struct hstrangecorrelationfilter {
if (!track.hasITS() && triggerRequireITS) {
continue; // skip, doesn't have ITS signal (skips lots of TPC-only!)
}
//prepared for Bayesian PID
// if (!track.bayesPi() > pionMinBayesProb) {
// continue;
// }
// if (track.bayesPi() < track.bayesPr() || track.bayesPi() < track.bayesKa()){
// continue;
// }
// if (track.tpcNSigmaPi() < assocPionNSigmaTPCFOF){
// continue;
// }
// if (track.tofSignal() > 0 && track.tofNSigmaPi() < assocPionNSigmaTPCFOF){
// continue;
// }
if (track.tofSignal() > 0){
if(std::sqrt(track.tofNSigmaPi()*track.tofNSigmaPi()+track.tpcNSigmaPi()*track.tpcNSigmaPi()) > assocPionNSigmaTPCFOF)
// prepared for Bayesian PID
// if (!track.bayesPi() > pionMinBayesProb) {
// continue;
// }
// if (track.bayesPi() < track.bayesPr() || track.bayesPi() < track.bayesKa()){
// continue;
// }
// if (track.tpcNSigmaPi() < assocPionNSigmaTPCFOF){
// continue;
// }
// if (track.tofSignal() > 0 && track.tofNSigmaPi() < assocPionNSigmaTPCFOF){
// continue;
// }
if (track.tofSignal() > 0) {
if (std::sqrt(track.tofNSigmaPi() * track.tofNSigmaPi() + track.tpcNSigmaPi() * track.tpcNSigmaPi()) > assocPionNSigmaTPCFOF)
continue;
if (track.tofNSigmaPr() < rejectSigma) continue;
if (track.tpcNSigmaPr() < rejectSigma) continue;
if (track.tofNSigmaKa() < rejectSigma) continue;
if (track.tpcNSigmaKa() < rejectSigma) continue;
}else{
if(track.tpcNSigmaPi() > assocPionNSigmaTPCFOF)
if (track.tofNSigmaPr() < rejectSigma)
continue;
if (track.tpcNSigmaPr() < rejectSigma)
continue;
if (track.tofNSigmaKa() < rejectSigma)
continue;
if (track.tpcNSigmaKa() < rejectSigma)
continue;
} else {
if (track.tpcNSigmaPi() > assocPionNSigmaTPCFOF)
continue;
if (track.tpcNSigmaPr() < rejectSigma)
continue;
if (track.tpcNSigmaKa() < rejectSigma)
continue;
if (track.tpcNSigmaPr() < rejectSigma) continue;
if (track.tpcNSigmaKa() < rejectSigma) continue;
}


assocPion(
track.collisionId(),
track.globalIndex());
Expand Down
16 changes: 8 additions & 8 deletions PWGLF/Tasks/hStrangeCorrelation.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -157,7 +157,7 @@ struct correlateStrangeness {
float deltaphi = ComputeDeltaPhi(trigg.phi(), assoc.phi());
float deltaeta = trigg.eta() - assoc.eta();
float ptassoc = assoc.pt();
if(!mixing)
if (!mixing)
histos.fill(HIST("sameEvent/Pion"), deltaphi, deltaeta, ptassoc, pvz, mult);
else
histos.fill(HIST("mixedEvent/Pion"), deltaphi, deltaeta, ptassoc, pvz, mult);
Expand Down Expand Up @@ -255,7 +255,7 @@ struct correlateStrangeness {
if (doCorrelationPion)
histos.add("sameEvent/Pion", "Pion", kTHnF, {axisDeltaPhi, axisDeltaEta, axisPtAssoc, ConfVtxBins, ConfMultBins});

if(doCorrelationK0Short || doCorrelationLambda || doCorrelationAntiLambda || doCorrelationXiMinus || doCorrelationXiPlus || doCorrelationOmegaMinus || doCorrelationOmegaPlus) {
if (doCorrelationK0Short || doCorrelationLambda || doCorrelationAntiLambda || doCorrelationXiMinus || doCorrelationXiPlus || doCorrelationOmegaMinus || doCorrelationOmegaPlus) {
histos.addClone("sameEvent/Signal/", "sameEvent/LeftBg/");
histos.addClone("sameEvent/Signal/", "sameEvent/RightBg/");
}
Expand Down Expand Up @@ -369,8 +369,8 @@ struct correlateStrangeness {
fillCorrelationsCascade(triggerTracks, associatedCascades, false, collision.posZ(), collision.centFT0M());
}
void processSameEventHPions(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms>::iterator const& collision,
aod::AssocPions const& associatedPions, aod::TriggerTracks const& triggerTracks,
TracksComplete const&)
aod::AssocPions const& associatedPions, aod::TriggerTracks const& triggerTracks,
TracksComplete const&)
{
// ________________________________________________
// Perform basic event selection
Expand All @@ -381,7 +381,7 @@ struct correlateStrangeness {
return;
}
// ________________________________________________
if (!doprocessSameEventHCascades && ! doprocessSameEventHV0s) {
if (!doprocessSameEventHCascades && !doprocessSameEventHV0s) {
histos.fill(HIST("MixingQA/hSECollisionBins"), colBinning.getBin({collision.posZ(), collision.centFT0M()}));
histos.fill(HIST("EventQA/hMult"), collision.centFT0M());
histos.fill(HIST("EventQA/hPvz"), collision.posZ());
Expand All @@ -391,7 +391,7 @@ struct correlateStrangeness {
auto pionTrack = pion.track_as<TracksComplete>();
histos.fill(HIST("hPionEtaVsPtVsPhi"), pionTrack.pt(), pionTrack.eta(), pionTrack.phi());
}
if (!doprocessSameEventHCascades && ! doprocessSameEventHV0s) {
if (!doprocessSameEventHCascades && !doprocessSameEventHV0s) {
for (auto const& triggerTrack : triggerTracks) {
auto track = triggerTrack.track_as<TracksComplete>();
histos.fill(HIST("hTrackEtaVsPtVsPhi"), track.pt(), track.eta(), track.phi());
Expand Down Expand Up @@ -455,8 +455,8 @@ struct correlateStrangeness {
}
}
void processMixedEventHPions(soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms> const& collisions,
aod::AssocPions const& assocPions, aod::TriggerTracks const& triggerTracks,
TracksComplete const&)
aod::AssocPions const& assocPions, aod::TriggerTracks const& triggerTracks,
TracksComplete const&)
{
for (auto& [collision1, collision2] : soa::selfCombinations(colBinning, 5, -1, collisions, collisions)) {
// ________________________________________________
Expand Down