Skip to content

Please consider the following formatting changes to #4240 #13

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
fkrizek merged 1 commit into from
Jan 5, 2024
Merged

Please consider the following formatting changes to #4240 #13

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
104 changes: 51 additions & 53 deletions PWGJE/Tasks/ChJetTriggerQATask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -106,40 +106,39 @@ struct ChJetTriggerQATask {
eventSelection = JetDerivedDataUtilities::initialiseEventSelection(static_cast<std::string>(evSel));
trackSelection = JetDerivedDataUtilities::initialiseTrackSelection(static_cast<std::string>(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
Expand All @@ -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;
Expand All @@ -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(),
Expand All @@ -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(),
Expand All @@ -223,8 +222,7 @@ struct ChJetTriggerQATask {
}
}


if(bAddSupplementHistosToOutput){
if (bAddSupplementHistosToOutput) {
if (leadingTrackPt > -1.) {
spectra.fill(HIST("ptphiLeadingTrack"), leadingTrackPt,
leadingTrackPhi);
Expand All @@ -244,21 +242,21 @@ struct ChJetTriggerQATask {
}

// access jet constituents as tracks
if(bAddBigHistosToOutput){
if (bAddBigHistosToOutput) {
for (auto& jct : jet.tracks_as<TrackCandidates>()) {
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()));
Expand All @@ -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
Expand All @@ -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<float>(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());
}
}
}
}
}
Expand Down