diff --git a/PWGJE/Tasks/ChJetTriggerQATask.cxx b/PWGJE/Tasks/ChJetTriggerQATask.cxx index ec618db63d2..a9f15ccf85b 100644 --- a/PWGJE/Tasks/ChJetTriggerQATask.cxx +++ b/PWGJE/Tasks/ChJetTriggerQATask.cxx @@ -106,40 +106,39 @@ struct ChJetTriggerQATask { eventSelection = JetDerivedDataUtilities::initialiseEventSelection(static_cast(evSel)); trackSelection = JetDerivedDataUtilities::initialiseTrackSelection(static_cast(trackSelections)); - //Basic histos - spectra.add("vertexZ", "z vertex", {HistType::kTH1F, {{400, -20., +20.}}}); - spectra.add("ptphiTrackInclGood","pT vs phi inclusive good tracks",{HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("ptetaTrackInclGood","pT vs eta inclusive good tracks",{HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); - spectra.add("ptJetChInclFidVol","inclusive charged jet pT in fiducial volume",{HistType::kTH1F, {{200, 0., +200.}}}); - spectra.add("ptphiJetChInclFidVol","inclusive charged jet pT vs phi in fiducial volume",{HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("ptphiJetChInclFullVol","inclusive charged jet pT vs phi in full TPC volume",{HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("ptetaJetChInclFidVol","inclusive charged jet pT vs eta in fiducial volume",{HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); - spectra.add("ptetaJetChInclFullVol","inclusive charged jet pT vs eta in full TPC volume",{HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); - spectra.add("ptetaLeadingJetFullVol","pT vs eta leading jet",{HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); - spectra.add("ptphiLeadingJetFullVol","pT vs phi leading jet",{HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); - - //Supplementary plots - if(bAddSupplementHistosToOutput){ - spectra.add("ptJetChInclFullVol","inclusive charged jet pT in full volume",{HistType::kTH1F, {{200, 0., +200.}}}); - spectra.add("phietaTrackAllInclGood","phi vs eta all inclusive good tracks",{HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("phietaTrackHighPtInclGood","phi vs eta inclusive good tracks with pT > 10 GeV",{HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("phietaJetChInclFidVol","inclusive charged jet phi vs eta in fiducial volume",{HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("phietaJetChInclFullVol","inclusive charged jet phi vs eta in full TPC volume",{HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("phietaJetChInclHighPtFidVol","inclusive charged jet phi vs eta in fiducial volume",{HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("phietaJetChInclHighPtFullVol","inclusive charged jet phi vs eta in full TPC volume",{HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("ptetaLeadingTrack","pT vs eta leading tracks",{HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); - spectra.add("ptphiLeadingTrack","pT vs phi leading tracks",{HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); - spectra.add("jetAreaFullVol","area of all jets in full TPC volume",{HistType::kTH2F, {{100, 0., +100.}, {50, 0., 2.}}}); - spectra.add("jetAreaFidVol","area of all jets in fiducial volume",{HistType::kTH2F, {{100, 0., +100.}, {50, 0., 2.}}}); + // Basic histos + spectra.add("vertexZ", "z vertex", {HistType::kTH1F, {{400, -20., +20.}}}); + spectra.add("ptphiTrackInclGood", "pT vs phi inclusive good tracks", {HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("ptetaTrackInclGood", "pT vs eta inclusive good tracks", {HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); + spectra.add("ptJetChInclFidVol", "inclusive charged jet pT in fiducial volume", {HistType::kTH1F, {{200, 0., +200.}}}); + spectra.add("ptphiJetChInclFidVol", "inclusive charged jet pT vs phi in fiducial volume", {HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("ptphiJetChInclFullVol", "inclusive charged jet pT vs phi in full TPC volume", {HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("ptetaJetChInclFidVol", "inclusive charged jet pT vs eta in fiducial volume", {HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); + spectra.add("ptetaJetChInclFullVol", "inclusive charged jet pT vs eta in full TPC volume", {HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); + spectra.add("ptetaLeadingJetFullVol", "pT vs eta leading jet", {HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); + spectra.add("ptphiLeadingJetFullVol", "pT vs phi leading jet", {HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); + + // Supplementary plots + if (bAddSupplementHistosToOutput) { + spectra.add("ptJetChInclFullVol", "inclusive charged jet pT in full volume", {HistType::kTH1F, {{200, 0., +200.}}}); + spectra.add("phietaTrackAllInclGood", "phi vs eta all inclusive good tracks", {HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("phietaTrackHighPtInclGood", "phi vs eta inclusive good tracks with pT > 10 GeV", {HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("phietaJetChInclFidVol", "inclusive charged jet phi vs eta in fiducial volume", {HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("phietaJetChInclFullVol", "inclusive charged jet phi vs eta in full TPC volume", {HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("phietaJetChInclHighPtFidVol", "inclusive charged jet phi vs eta in fiducial volume", {HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("phietaJetChInclHighPtFullVol", "inclusive charged jet phi vs eta in full TPC volume", {HistType::kTH2F, {{80, -1., 1.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("ptetaLeadingTrack", "pT vs eta leading tracks", {HistType::kTH2F, {{100, 0., +100.}, {80, -1., 1.}}}); + spectra.add("ptphiLeadingTrack", "pT vs phi leading tracks", {HistType::kTH2F, {{100, 0., +100.}, {60, 0, TMath::TwoPi()}}}); + spectra.add("jetAreaFullVol", "area of all jets in full TPC volume", {HistType::kTH2F, {{100, 0., +100.}, {50, 0., 2.}}}); + spectra.add("jetAreaFidVol", "area of all jets in fiducial volume", {HistType::kTH2F, {{100, 0., +100.}, {50, 0., 2.}}}); } - //3D histogram - if(bAddBigHistosToOutput){ - spectra.add("fLeadJetChPtVsLeadingTrack","inclusive charged jet pT in TPC volume",{HistType::kTH2F, {{200, 0., +200.}, {200, 0., +200.}}}); - spectra.add("tracksThatWereNotJetConstituentsPtEtaPhi","PtEtaPhi of tracksThatWereNotjetConsituents in full volume", - {HistType::kTH3F, {{100, 0., +100.}, {40, -1., 1.}, {60, 0., TMath::TwoPi()}}}); + // 3D histogram + if (bAddBigHistosToOutput) { + spectra.add("fLeadJetChPtVsLeadingTrack", "inclusive charged jet pT in TPC volume", {HistType::kTH2F, {{200, 0., +200.}, {200, 0., +200.}}}); + spectra.add("tracksThatWereNotJetConstituentsPtEtaPhi", "PtEtaPhi of tracksThatWereNotjetConsituents in full volume", + {HistType::kTH3F, {{100, 0., +100.}, {40, -1., 1.}, {60, 0., TMath::TwoPi()}}}); } - } // declare filters on collisions @@ -164,11 +163,11 @@ struct ChJetTriggerQATask { return; } - if ((bLowPtTrigger && JetDerivedDataUtilities::selectChargedTrigger(collision, JetDerivedDataUtilities::JTrigSelCh::chargedLow)) || (bHighPtTrigger && JetDerivedDataUtilities::selectChargedTrigger(collision, JetDerivedDataUtilities::JTrigSelCh::chargedHigh)) || ((!bLowPtTrigger)&&(!bHighPtTrigger))) { - //bLowPtTrigger=1 and bHighPtTrigger=0 --> fill histos with low trigger only - //bLowPtTrigger=0 and bHighPtTrigger=1 --> fill histos with high trigger only - //bLowPtTrigger=1 and bHighPtTrigger=1 --> fill histos with mixture of low and high trigger - //bLowPtTrigger=0 and bHighPtTrigger=0 --> fill histos with minimum bias ie. ignore trigger decision + if ((bLowPtTrigger && JetDerivedDataUtilities::selectChargedTrigger(collision, JetDerivedDataUtilities::JTrigSelCh::chargedLow)) || (bHighPtTrigger && JetDerivedDataUtilities::selectChargedTrigger(collision, JetDerivedDataUtilities::JTrigSelCh::chargedHigh)) || ((!bLowPtTrigger) && (!bHighPtTrigger))) { + // bLowPtTrigger=1 and bHighPtTrigger=0 --> fill histos with low trigger only + // bLowPtTrigger=0 and bHighPtTrigger=1 --> fill histos with high trigger only + // bLowPtTrigger=1 and bHighPtTrigger=1 --> fill histos with mixture of low and high trigger + // bLowPtTrigger=0 and bHighPtTrigger=0 --> fill histos with minimum bias ie. ignore trigger decision float leadingJetPt = -1.0; float leadingJetEta = -2.0; @@ -191,10 +190,10 @@ struct ChJetTriggerQATask { continue; } - if(bAddBigHistosToOutput){ + if (bAddBigHistosToOutput) { v.SetPtEtaPhiM(trk.pt(), trk.eta(), trk.phi(), 0.139); acceptedTracks.push_back(GoodTrack(v, false, trk.globalIndex())); - } + } spectra.fill( HIST("ptphiTrackInclGood"), trk.pt(), @@ -203,11 +202,11 @@ struct ChJetTriggerQATask { HIST("ptetaTrackInclGood"), trk.pt(), trk.eta()); // Inclusive Track pT vs eta spectrum in TPC volume - if(bAddSupplementHistosToOutput){ + if (bAddSupplementHistosToOutput) { spectra.fill( HIST("phietaTrackAllInclGood"), trk.eta(), trk.phi()); // Inclusive Track pT vs eta spectrum in TPC volume - + if (trk.pt() > 5.0) { spectra.fill( HIST("phietaTrackHighPtInclGood"), trk.eta(), @@ -223,8 +222,7 @@ struct ChJetTriggerQATask { } } - - if(bAddSupplementHistosToOutput){ + if (bAddSupplementHistosToOutput) { if (leadingTrackPt > -1.) { spectra.fill(HIST("ptphiLeadingTrack"), leadingTrackPt, leadingTrackPhi); @@ -244,21 +242,21 @@ struct ChJetTriggerQATask { } // access jet constituents as tracks - if(bAddBigHistosToOutput){ + if (bAddBigHistosToOutput) { for (auto& jct : jet.tracks_as()) { for (UInt_t itr = 0; itr < acceptedTracks.size(); itr++) { if (acceptedTracks[itr].globalIndex == jct.globalIndex()) { - + acceptedTracks[itr].isJetConstituent = true; // initialization break; } } } - } + } } } - if(bAddBigHistosToOutput){ + if (bAddBigHistosToOutput) { for (UInt_t itr = 0; itr < acceptedTracks.size(); itr++) { if (!acceptedTracks[itr].isJetConstituent) { spectra.fill(HIST("tracksThatWereNotJetConstituentsPtEtaPhi"), acceptedTracks[itr].lv.Pt(), acceptedTracks[itr].lv.Eta(), TVector2::Phi_0_2pi(acceptedTracks[itr].lv.Phi())); @@ -271,7 +269,7 @@ struct ChJetTriggerQATask { spectra.fill(HIST("ptetaLeadingJetFullVol"), leadingJetPt, leadingJetEta); } - if(bAddBigHistosToOutput){ + if (bAddBigHistosToOutput) { if (leadingJetPt > -1. && leadingTrackPt > -1.) { spectra.fill(HIST("fLeadJetChPtVsLeadingTrack"), leadingTrackPt, leadingJetPt); // leading jet pT versus leading track pT @@ -285,28 +283,28 @@ struct ChJetTriggerQATask { spectra.fill(HIST("ptphiJetChInclFidVol"), jet.pt(), jet.phi()); spectra.fill(HIST("ptetaJetChInclFidVol"), jet.pt(), jet.eta()); - if(bAddSupplementHistosToOutput){ + if (bAddSupplementHistosToOutput) { spectra.fill(HIST("phietaJetChInclFidVol"), jet.eta(), jet.phi()); if (jet.pt() > 10.0) { spectra.fill(HIST("phietaJetChInclHighPtFidVol"), jet.eta(), jet.phi()); } spectra.fill(HIST("jetAreaFidVol"), jet.pt(), jet.area()); - } + } } if (fabs(jet.eta()) < static_cast(cfgTPCVolume)) { spectra.fill(HIST("ptphiJetChInclFullVol"), jet.pt(), jet.phi()); spectra.fill(HIST("ptetaJetChInclFullVol"), jet.pt(), jet.eta()); - if(bAddSupplementHistosToOutput){ - spectra.fill(HIST("ptJetChInclFullVol"), jet.pt()); + if (bAddSupplementHistosToOutput) { + spectra.fill(HIST("ptJetChInclFullVol"), jet.pt()); - spectra.fill(HIST("phietaJetChInclFullVol"), jet.eta(), jet.phi()); + spectra.fill(HIST("phietaJetChInclFullVol"), jet.eta(), jet.phi()); if (jet.pt() > 10.0) { spectra.fill(HIST("phietaJetChInclHighPtFullVol"), jet.eta(), jet.phi()); } spectra.fill(HIST("jetAreaFullVol"), jet.pt(), jet.area()); - } + } } } }