diff --git a/PWGUD/TableProducer/SGCandProducer.cxx b/PWGUD/TableProducer/SGCandProducer.cxx index 35f000b9378..55b6e4fef60 100644 --- a/PWGUD/TableProducer/SGCandProducer.cxx +++ b/PWGUD/TableProducer/SGCandProducer.cxx @@ -33,11 +33,11 @@ struct SGCandProducer { Configurable SGCuts{"SGCuts", {}, "SG event cuts"}; Configurable saveAllTracks{"saveAllTracks", true, "save only PV contributors or all tracks associated to a collision"}; Configurable savenonPVCITSOnlyTracks{"savenonPVCITSOnlyTracks", false, "save non PV contributors with ITS only information"}; - Configurable rejectAtTFBoundary{"rejectAtTFBoundary", true, "reject collisions at a TF boundary"}; - Configurable noITSROFrameBorder{"noITSROFrameBorder", true, "reject ITS RO Frame Border"}; - Configurable noSameBunchPileUp{"noSameBunchPileUp", true, "reject SameBunchPileUp"}; - Configurable IsGoodVertex{"IsGoodVertex", false, "Select FT0 PV vertex matching"}; - Configurable ITSTPCVertex{"ITSTPCVertex", true, "reject ITS-only vertex"}; //if one wants to look at Single Gap pp events + Configurable rejectAtTFBoundary{"rejectAtTFBoundary", true, "reject collisions at a TF boundary"}; + Configurable noITSROFrameBorder{"noITSROFrameBorder", true, "reject ITS RO Frame Border"}; + Configurable noSameBunchPileUp{"noSameBunchPileUp", true, "reject SameBunchPileUp"}; + Configurable IsGoodVertex{"IsGoodVertex", false, "Select FT0 PV vertex matching"}; + Configurable ITSTPCVertex{"ITSTPCVertex", true, "reject ITS-only vertex"}; // if one wants to look at Single Gap pp events // SG selector SGSelector sgSelector; @@ -153,26 +153,26 @@ struct SGCandProducer { { LOGF(debug, " collision %d", collision.globalIndex()); registry.get(HIST("reco/Stat"))->Fill(0., 1.); - //reject collisions at TF boundaries + // reject collisions at TF boundaries if (rejectAtTFBoundary && !collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) { return; } - //reject collisions at ITS RO TF boundaries + // reject collisions at ITS RO TF boundaries if (noITSROFrameBorder && !collision.selection_bit(aod::evsel::kNoITSROFrameBorder)) { return; } // registry.get(HIST("reco/Stat"))->Fill(1., 1.); - //reject Same Bunch PileUp + // reject Same Bunch PileUp if (noSameBunchPileUp && !collision.selection_bit(aod::evsel::kNoSameBunchPileup)) { return; } // registry.get(HIST("reco/Stat"))->Fill(1., 1.); - //check vertex matching to FT0 + // check vertex matching to FT0 if (IsGoodVertex && !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV)) { return; } // registry.get(HIST("reco/Stat"))->Fill(1., 1.); - //reject ITS Only vertices + // reject ITS Only vertices if (ITSTPCVertex && !collision.selection_bit(aod::evsel::kIsVertexITSTPC)) { return; } @@ -228,13 +228,15 @@ struct SGCandProducer { } // update SGTracks tables for (auto& track : tracks) { - if (track.pt() > sameCuts.minPt() && track.eta() > sameCuts.minEta() && track.eta() < sameCuts.maxEta()){ - if (track.isPVContributor()) updateUDTrackTables(outputCollisions.lastIndex(), track, bc.globalBC()); - else if (saveAllTracks) { - if (track.itsClusterSizes() && track.itsChi2NCl() > 0 && ((track.tpcNClsFindable() == 0 && savenonPVCITSOnlyTracks) || track.tpcNClsFindable() > 50)) updateUDTrackTables(outputCollisions.lastIndex(), track, bc.globalBC()); - // if (track.isPVContributor()) updateUDTrackTables(outputCollisions.lastIndex(), track, bc.globalBC()); - } - } + if (track.pt() > sameCuts.minPt() && track.eta() > sameCuts.minEta() && track.eta() < sameCuts.maxEta()) { + if (track.isPVContributor()) + updateUDTrackTables(outputCollisions.lastIndex(), track, bc.globalBC()); + else if (saveAllTracks) { + if (track.itsClusterSizes() && track.itsChi2NCl() > 0 && ((track.tpcNClsFindable() == 0 && savenonPVCITSOnlyTracks) || track.tpcNClsFindable() > 50)) + updateUDTrackTables(outputCollisions.lastIndex(), track, bc.globalBC()); + // if (track.isPVContributor()) updateUDTrackTables(outputCollisions.lastIndex(), track, bc.globalBC()); + } + } } // update SGFwdTracks tables if (sameCuts.withFwdTracks()) { diff --git a/PWGUD/Tasks/SG_FourPi_Analyzer.cxx b/PWGUD/Tasks/SG_FourPi_Analyzer.cxx index 8771d268479..21cbd9847f6 100644 --- a/PWGUD/Tasks/SG_FourPi_Analyzer.cxx +++ b/PWGUD/Tasks/SG_FourPi_Analyzer.cxx @@ -32,7 +32,7 @@ using namespace std; using namespace o2; using namespace o2::aod; -//using namespace o2::aod::track::v001; +// using namespace o2::aod::track::v001; using namespace o2::framework; using namespace o2::framework::expressions; #define mpion 0.1396 @@ -92,9 +92,9 @@ struct SG_FourPi_Analyzer { // Single gap either side TLorentzVector v01; // int truegapSide = sgSelector.trueGap(collision); - //int truegapSide = sgSelector.trueGap(collision, FV0_cut, ZDC_cut); - float FIT_cut[5]={FV0_cut, FT0A_cut, FT0C_cut, FDDA_cut, FDDC_cut}; - //int truegapSide = sgSelector.trueGap(collision, *FIT_cut, ZDC_cut); + // int truegapSide = sgSelector.trueGap(collision, FV0_cut, ZDC_cut); + float FIT_cut[5] = {FV0_cut, FT0A_cut, FT0C_cut, FDDA_cut, FDDC_cut}; + // int truegapSide = sgSelector.trueGap(collision, *FIT_cut, ZDC_cut); int truegapSide = sgSelector.trueGap(collision, FIT_cut[0], FIT_cut[1], FIT_cut[3], ZDC_cut); registry.fill(HIST("GapSide"), gapSide); registry.fill(HIST("TrueGapSide"), truegapSide); @@ -103,70 +103,78 @@ struct SG_FourPi_Analyzer { // Look for D0 and D0bar float sign = 0; for (auto t : tracks) { - int itsNCls = t.itsNCls(); -// int uditsNCls = t.uditsNCls(); -// int itsNCls = o2::aod::track::v001::ITSNCls; -// if (t.itsNCls()) {std::cout << t.itsClusterSizes()<; + // if (t.itsNCls()) {std::cout << t.itsClusterSizes()<