Skip to content

Commit 1e36b14

Browse files
JimunLeejimun_lee
JimunLee
and
jimun_lee
authored
[PWGJE] [PWGJE]:Added switch options (#11014)
Co-authored-by: jimun_lee <jimun.lee@cern.ch>
1 parent 9cce6ac commit 1e36b14

File tree

1 file changed

+114
-87
lines changed

1 file changed

+114
-87
lines changed

‎PWGJE/Tasks/phiInJets.cxx‎

Lines changed: 114 additions & 87 deletions
Original file line numberDiff line numberDiff line change
@@ -84,6 +84,10 @@ struct phiInJets {
8484
Configurable<bool> cfgMCRecMBHists{"cfgMCRecMBHists", false, "Enables MCRec MB Hists"};
8585
Configurable<bool> cfgMCRecInsideHists{"cfgMCRecInsideHists", false, "Enables MCRec Inside Hists"};
8686
Configurable<bool> cfgMCRecRotationalHists{"cfgMCRecRotationalHists", false, "Enables MCRotational Hists"};
87+
Configurable<bool> cfgPIDQAHists{"cfgPIDQAHists", false, "Enables PIDQA Hists"};
88+
Configurable<bool> cfgDaughterQAHists{"cfgDaughterQAHists", false, "Enables DaughterQA Hists"};
89+
Configurable<bool> cfgJetQAHists{"cfgJetQAHists", false, "Enables JetQA Hists"};
90+
8791
Configurable<bool> cfgMCGenHists{"cfgMCGenHists", false, "Enables MCGenHists"};
8892
Configurable<bool> cfgMCGenMATCHEDHists{"cfgMCGenMATCHEDHists", false, "Enables MCGenMATCHEDHists"};
8993
Configurable<bool> cfgMCRecMATCHEDHists{"cfgMCRecMATCHEDHists", false, "Enables MCRecMATCHEDHists"};
@@ -152,56 +156,58 @@ struct phiInJets {
152156

153157
if (cfgMCRecHists) {
154158
JEhistos.add("nEvents_MCRec", "nEvents_MCRec", kTH1F, {{4, 0.0, 4.0}});
159+
if (cfgJetQAHists) {
160+
JEhistos.add("h_jet_pt", "jet pT;#it{p}_{T,jet} (GeV/#it{c});entries", {HistType::kTH1F, {{4000, 0., 200.}}});
161+
JEhistos.add("h_jet_eta", "jet #eta;#eta_{jet};entries", {HistType::kTH1F, {{100, -1.0, 1.0}}});
162+
JEhistos.add("h_jet_phi", "jet #phi;#phi_{jet};entries", {HistType::kTH1F, {{80, -1.0, 7.}}});
163+
}
164+
if (cfgPIDQAHists) {
165+
JEhistos.add("ptJEHistogramPion", "ptJEHistogramPion", kTH1F, {PtAxis});
166+
JEhistos.add("ptJEHistogramKaon", "ptJEHistogramKaon", kTH1F, {PtAxis});
167+
JEhistos.add("ptJEHistogramProton", "ptJEHistogramProton", kTH1F, {PtAxis});
168+
JEhistos.add("ptJEHistogramPhi", "ptJEHistogramPhi", kTH1F, {PtAxis});
169+
JEhistos.add("ptJEHistogramPhi_JetTrigger", "ptJEHistogramPhi_JetTrigger", kTH1F, {PtAxis});
170+
}
155171

156-
JEhistos.add("h_jet_pt", "jet pT;#it{p}_{T,jet} (GeV/#it{c});entries", {HistType::kTH1F, {{4000, 0., 200.}}});
157-
JEhistos.add("h_jet_eta", "jet #eta;#eta_{jet};entries", {HistType::kTH1F, {{100, -1.0, 1.0}}});
158-
JEhistos.add("h_jet_phi", "jet #phi;#phi_{jet};entries", {HistType::kTH1F, {{80, -1.0, 7.}}});
159-
160-
JEhistos.add("ptJEHistogramPion", "ptJEHistogramPion", kTH1F, {PtAxis});
161-
JEhistos.add("ptJEHistogramKaon", "ptJEHistogramKaon", kTH1F, {PtAxis});
162-
JEhistos.add("ptJEHistogramProton", "ptJEHistogramProton", kTH1F, {PtAxis});
163-
JEhistos.add("ptJEHistogramPhi", "ptJEHistogramPhi", kTH1F, {PtAxis});
164-
JEhistos.add("ptJEHistogramPhi_JetTrigger", "ptJEHistogramPhi_JetTrigger", kTH1F, {PtAxis});
165-
JEhistos.add("minvJEHistogramPhi", "minvJEHistogramPhi", kTH1F, {MinvAxis});
166-
JEhistos.add("hNRealPhiVPhiCand", "hNRealPhiVPhiCand", kTH2F, {{10, 0, 10}, {10, 0, 10}});
167-
JEhistos.add("hNRealPhiWithJetVPhiCand", "hNRealPhiWithJetVPhiCand", kTH2F, {{10, 0, 10}, {10, 0, 10}});
168-
JEhistos.add("hNRealPhiInJetVPhiCand", "hNRealPhiInJetVPhiCand", kTH2F, {{10, 0, 10}, {10, 0, 10}});
169-
170-
JEhistos.add("hMCRec_nonmatch_hUSS_KtoKangle_v_pt", "hMCRec_nonmatch_hUSS_KtoKangle_v_pt", kTH2F, {axisEta, PtAxis});
171-
JEhistos.add("hMCRec_nonmatch_hUSS_Kangle_v_pt", "hMCRec_nonmatch_hUSS_Kangle_v_pt", kTH2F, {axisEta, PtAxis});
172-
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE_pt_v_eta", "hMCRec_nonmatch_hUSS_INSIDE_pt_v_eta", kTH2F, {PtAxis, axisEta});
173-
JEhistos.add("JetVsPhi_REC", "JetVsPhi_REC", kTH2F, {{4000, 0., 200.}, {200, 0, 20.0}});
172+
if (cfgDaughterQAHists) {
173+
JEhistos.add("hNRealPhiVPhiCand", "hNRealPhiVPhiCand", kTH2F, {{10, 0, 10}, {10, 0, 10}});
174+
JEhistos.add("hNRealPhiWithJetVPhiCand", "hNRealPhiWithJetVPhiCand", kTH2F, {{10, 0, 10}, {10, 0, 10}});
175+
JEhistos.add("hNRealPhiInJetVPhiCand", "hNRealPhiInJetVPhiCand", kTH2F, {{10, 0, 10}, {10, 0, 10}});
176+
JEhistos.add("hMCRec_nonmatch_hUSS_KtoKangle_v_pt", "hMCRec_nonmatch_hUSS_KtoKangle_v_pt", kTH2F, {axisEta, PtAxis});
177+
JEhistos.add("hMCRec_nonmatch_hUSS_Kangle_v_pt", "hMCRec_nonmatch_hUSS_Kangle_v_pt", kTH2F, {axisEta, PtAxis});
178+
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE_pt_v_eta", "hMCRec_nonmatch_hUSS_INSIDE_pt_v_eta", kTH2F, {PtAxis, axisEta});
179+
}
174180

175181
// used for Minv closure tests
176182
// MB
177183
if (cfgMCRecMBHists) {
178184
JEhistos.add("hMCRec_hUSS", "hMCRec_hUSS", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
179185
JEhistos.add("hMCRec_hLSS", "hMCRec_hLSS", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
180186
JEhistos.add("hMCRecTrue_hUSS", "hMCRecTrue_hUSS", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
181-
JEhistos.add("hMCRecTrue_hLSS", "hMCRecTrue_hLSS", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
182187
JEhistos.add("hMCRec_R_distribution", "hMCRec_R_distribution", kTH1F, {{100, 0.0, 2 * TMath::Pi()}});
183188
JEhistos.add("hMCRec_dPhi_distribution", "hMCRec_dPhi_distribution", kTH1F, {{80, -5.0, 7.0}});
184189
JEhistos.add("hMCRec_dEta_distribution", "hMCRec_dEta_distribution", kTH1F, {{100, -2.0, 2.0}});
190+
191+
if (cfgMCRecRotationalHists) {
192+
JEhistos.add("hMCRec_hUSS_Rotational", "hMCRec_hUSS_Rotational", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
193+
JEhistos.add("hMCRec_R_Rotation_distribution", "hMCRec_R_Rotation_distribution", HistType::kTH1F, {{100, 0.0, 2 * TMath::Pi()}});
194+
JEhistos.add("hMCRec_dPhi_rot_distribution", "hMCRec_dPhi_rot_distribution", kTH1F, {{80, -5.0, 7.0}});
195+
JEhistos.add("hMCRec_dEta_rot_distribution", "hMCRec_dEta_rot_distribution", kTH1F, {{100, -2.0, 2.0}});
196+
JEhistos.add("hMCRec_dEta_qa_rot_distribution", "hMCRec_dEta_qa_rot_distribution", kTH1F, {{100, -4.0, 2.0}});
197+
}
185198
}
186-
if (cfgMCRecRotationalHists) {
187-
JEhistos.add("hMCRec_hUSS_Rotational", "hMCRec_hUSS_Rotational", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
188-
JEhistos.add("hMCRec_R_Rotation_distribution", "hMCRec_R_Rotation_distribution", HistType::kTH1F, {{100, 0.0, 2 * TMath::Pi()}});
189-
JEhistos.add("hMCRec_dPhi_rot_distribution", "hMCRec_dPhi_rot_distribution", kTH1F, {{80, -5.0, 7.0}});
190-
JEhistos.add("hMCRec_dEta_rot_distribution", "hMCRec_dEta_rot_distribution", kTH1F, {{100, -2.0, 2.0}});
191-
JEhistos.add("hMCRec_dEta_qa_rot_distribution", "hMCRec_dEta_qa_rot_distribution", kTH1F, {{100, -4.0, 2.0}});
199+
// INSIDE
200+
if (cfgMCRecInsideHists) {
201+
JEhistos.add("hMCRec_hUSS_INSIDE", "hMCRec_hUSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
202+
JEhistos.add("hMCRec_hLSS_INSIDE", "hMCRec_hLSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
203+
JEhistos.add("hMCRecTrue_hUSS_INSIDE", "hMCRecTrue_hUSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
204+
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE", "hMCRec_nonmatch_hUSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
205+
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE_1D", "hMCRec_nonmatch_hUSS_INSIDE_1D", kTH1F, {MinvAxis});
206+
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE_1D_2_3", "hMCRec_nonmatch_hUSS_INSIDE_1D_2_3", kTH1F, {MinvAxis});
207+
JEhistos.add("JetVsPhi_REC", "JetVsPhi_REC", kTH2F, {{4000, 0., 200.}, {200, 0, 20.0}});
208+
JEhistos.add("minvJEHistogramPhi", "minvJEHistogramPhi", kTH1F, {MinvAxis});
192209
}
193210
}
194-
// INSIDE
195-
if (cfgMCRecInsideHists) {
196-
JEhistos.add("hMCRec_hUSS_INSIDE", "hMCRec_hUSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
197-
JEhistos.add("hMCRec_hLSS_INSIDE", "hMCRec_hLSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
198-
JEhistos.add("hMCRecTrue_hUSS_INSIDE", "hMCRecTrue_hUSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
199-
JEhistos.add("hMCRecTrue_hLSS_INSIDE", "hMCRecTrue_hLSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
200-
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE", "hMCRec_nonmatch_hUSS_INSIDE", kTHnSparseF, {dRAxis, PtAxis, MinvAxis});
201-
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE_1D", "hMCRec_nonmatch_hUSS_INSIDE_1D", kTH1F, {MinvAxis});
202-
JEhistos.add("hMCRec_nonmatch_hUSS_INSIDE_1D_2_3", "hMCRec_nonmatch_hUSS_INSIDE_1D_2_3", kTH1F, {MinvAxis});
203-
}
204-
205211
if (cfgMCGenHists) {
206212
JEhistos.add("nEvents_MCGen", "nEvents_MCGen", kTH1F, {{4, 0.0, 4.0}});
207213

@@ -667,9 +673,11 @@ struct phiInJets {
667673
mcd_pt.push_back(mcdjet.pt());
668674
mcd_eta.push_back(mcdjet.eta());
669675
mcd_phi.push_back(mcdjet.phi());
670-
JEhistos.fill(HIST("h_jet_pt"), mcdjet.pt());
671-
JEhistos.fill(HIST("h_jet_eta"), mcdjet.eta());
672-
JEhistos.fill(HIST("h_jet_phi"), mcdjet.phi());
676+
if (cfgJetQAHists) {
677+
JEhistos.fill(HIST("h_jet_pt"), mcdjet.pt());
678+
JEhistos.fill(HIST("h_jet_eta"), mcdjet.eta());
679+
JEhistos.fill(HIST("h_jet_phi"), mcdjet.phi());
680+
}
673681
}
674682
if (hasJets)
675683
JEhistos.fill(HIST("nEvents_MCRec"), 2.5);
@@ -690,14 +698,15 @@ struct phiInJets {
690698

691699
if (track.has_mcParticle()) {
692700
auto mcParticle = track.mcParticle();
693-
694-
if (mcParticle.isPhysicalPrimary() && fabs(mcParticle.eta()) <= cfgtrkMaxEta) {
695-
if (abs(mcParticle.pdgCode()) == 211)
696-
JEhistos.fill(HIST("ptJEHistogramPion"), mcParticle.pt());
697-
if (abs(mcParticle.pdgCode()) == 321)
698-
JEhistos.fill(HIST("ptJEHistogramKaon"), mcParticle.pt());
699-
if (abs(mcParticle.pdgCode()) == 2212)
700-
JEhistos.fill(HIST("ptJEHistogramProton"), mcParticle.pt());
701+
if (cfgPIDQAHists) {
702+
if (mcParticle.isPhysicalPrimary() && fabs(mcParticle.eta()) <= cfgtrkMaxEta) {
703+
if (abs(mcParticle.pdgCode()) == 211)
704+
JEhistos.fill(HIST("ptJEHistogramPion"), mcParticle.pt());
705+
if (abs(mcParticle.pdgCode()) == 321)
706+
JEhistos.fill(HIST("ptJEHistogramKaon"), mcParticle.pt());
707+
if (abs(mcParticle.pdgCode()) == 2212)
708+
JEhistos.fill(HIST("ptJEHistogramProton"), mcParticle.pt());
709+
}
701710
}
702711
}
703712
for (const auto& track2 : tracks) {
@@ -719,9 +728,10 @@ struct phiInJets {
719728

720729
double dPhi = TVector2::Phi_mpi_pi(originalTrack.phi() - originalTrack2.phi());
721730
double dEta = originalTrack.eta() - originalTrack2.eta();
722-
JEhistos.fill(HIST("hMCRec_dPhi_distribution"), dPhi);
723-
JEhistos.fill(HIST("hMCRec_dEta_distribution"), dEta);
724-
731+
if (cfgMCRecMBHists) {
732+
JEhistos.fill(HIST("hMCRec_dPhi_distribution"), dPhi);
733+
JEhistos.fill(HIST("hMCRec_dEta_distribution"), dEta);
734+
}
725735
double dR = TMath::Sqrt(dPhi * dPhi + dEta * dEta);
726736
double dR_rot = 0;
727737
TLorentzVector lDecayDaughter1, lDecayDaughter2, lResonance, lRotationalTrack, lRotationalResonance;
@@ -760,21 +770,24 @@ struct phiInJets {
760770
double dEta_rot_qa = TMath::Abs(lDecayDaughter2.Eta()) - TMath::Abs(lRotationalTrack.Eta());
761771

762772
dR_rot = TMath::Sqrt(dPhi_rot * dPhi_rot + dEta_rot * dEta_rot);
763-
764-
JEhistos.fill(HIST("hMCRec_dPhi_rot_distribution"), dPhi_rot);
765-
JEhistos.fill(HIST("hMCRec_dEta_rot_distribution"), dEta_rot);
766-
JEhistos.fill(HIST("hMCRec_dEta_qa_rot_distribution"), dEta_rot_qa);
773+
if (cfgMCRecRotationalHists) {
774+
JEhistos.fill(HIST("hMCRec_dPhi_rot_distribution"), dPhi_rot);
775+
JEhistos.fill(HIST("hMCRec_dEta_rot_distribution"), dEta_rot);
776+
JEhistos.fill(HIST("hMCRec_dEta_qa_rot_distribution"), dEta_rot_qa);
777+
}
767778
}
779+
768780
//-----------------------------------------------------------------------
769781
lResonance = lDecayDaughter1 + lDecayDaughter2;
770782
if (cfgIsKstar)
771783
lRotationalResonance = lDecayDaughter1 + lRotationalTrack;
772784
if (fabs(lResonance.Eta()) > cfgtrkMaxEta)
773785
continue;
774786

775-
if (lResonance.M() > 1.005 && lResonance.M() < 1.035)
776-
PhiCand++;
777-
787+
if (cfgDaughterQAHists) {
788+
if (lResonance.M() > 1.005 && lResonance.M() < 1.035)
789+
PhiCand++;
790+
}
778791
//==================
779792
// 1.MB REC Closure
780793
//==================
@@ -802,9 +815,11 @@ struct phiInJets {
802815
double jetpt = 0;
803816

804817
for (std::size_t i = 0; i < mcd_pt.size(); i++) {
805-
if (i == 0) {
806-
if (lResonance.M() > 1.005 && lResonance.M() < 1.035) {
807-
RealPhiCandWithJet++;
818+
if (cfgDaughterQAHists) {
819+
if (i == 0) {
820+
if (lResonance.M() > 1.005 && lResonance.M() < 1.035) {
821+
RealPhiCandWithJet++;
822+
}
808823
}
809824
}
810825
double phidiff = TVector2::Phi_mpi_pi(mcd_phi[i] - lResonance.Phi());
@@ -818,7 +833,7 @@ struct phiInJets {
818833
double phidiff_K2 = TVector2::Phi_mpi_pi(mcd_phi[i] - lDecayDaughter2.Phi());
819834
double etadiff_K2 = mcd_eta[i] - lDecayDaughter2.Eta();
820835
double R_K2 = TMath::Sqrt((etadiff_K2 * etadiff_K2) + (phidiff_K2 * phidiff_K2));
821-
if (R < cfgjetR) {
836+
if (cfgDaughterQAHists && R < cfgjetR) {
822837
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_Kangle_v_pt"), R_K1, lResonance.Pt());
823838
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_Kangle_v_pt"), R_K2, lResonance.Pt());
824839
}
@@ -888,20 +903,22 @@ struct phiInJets {
888903
if (mothers1[0] != mothers2[0])
889904
continue; // Kaons not from the same phi
890905

891-
double phidiff_Kaons = TVector2::Phi_mpi_pi(lDecayDaughter2.Phi() - lDecayDaughter1.Phi());
892-
double etadiff_Kaons = lDecayDaughter2.Eta() - lDecayDaughter1.Eta();
893-
double R_Kaons = TMath::Sqrt((etadiff_Kaons * etadiff_Kaons) + (phidiff_Kaons * phidiff_Kaons));
894-
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_KtoKangle_v_pt"), R_Kaons, lResonance.Pt());
895-
JEhistos.fill(HIST("ptJEHistogramPhi"), lResonance.Pt());
906+
if (cfgDaughterQAHists) {
907+
double phidiff_Kaons = TVector2::Phi_mpi_pi(lDecayDaughter2.Phi() - lDecayDaughter1.Phi());
908+
double etadiff_Kaons = lDecayDaughter2.Eta() - lDecayDaughter1.Eta();
909+
double R_Kaons = TMath::Sqrt((etadiff_Kaons * etadiff_Kaons) + (phidiff_Kaons * phidiff_Kaons));
910+
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_KtoKangle_v_pt"), R_Kaons, lResonance.Pt());
911+
}
896912

913+
if (cfgPIDQAHists) {
914+
JEhistos.fill(HIST("ptJEHistogramPhi"), lResonance.Pt());
915+
}
897916
//=====================
898917
// 4.MB True Closure
899918
//=====================
900919
if (cfgMCRecMBHists) {
901920
if (originalTrack.sign() * originalTrack2.sign() < 0) {
902921
JEhistos.fill(HIST("hMCRecTrue_hUSS"), 1.0, lResonance.Pt(), lResonance.M());
903-
} else if (originalTrack.sign() * originalTrack2.sign() > 0) {
904-
JEhistos.fill(HIST("hMCRecTrue_hLSS"), 1.0, lResonance.Pt(), lResonance.M());
905922
}
906923
}
907924
//===========================
@@ -911,45 +928,55 @@ struct phiInJets {
911928
if (jetFlag) {
912929
if (originalTrack.sign() * originalTrack2.sign() < 0) {
913930
JEhistos.fill(HIST("hMCRecTrue_hUSS_INSIDE"), 1.0, lResonance.Pt(), lResonance.M());
914-
} else if (originalTrack.sign() * originalTrack2.sign() > 0) {
915-
JEhistos.fill(HIST("hMCRecTrue_hLSS_INSIDE"), 1.0, lResonance.Pt(), lResonance.M());
916931
}
917932
}
918933
}
919-
if (lResonance.M() > 1.005 && lResonance.M() < 1.035)
920-
RealPhiCand++;
921-
934+
if (cfgDaughterQAHists) {
935+
if (lResonance.M() > 1.005 && lResonance.M() < 1.035)
936+
RealPhiCand++;
937+
}
922938
// Now we do jets
923939
if (cfgMCRecInsideHists) {
924940
if (cfgSingleJet)
925941
if (goodjets > 1)
926942
jetpt = DistinguishJetsMC(mcd_pt, mcd_phi, mcd_eta, lResonance);
927943

928944
if (jetFlag) {
929-
if (lResonance.M() > 1.005 && lResonance.M() < 1.035)
930-
RealPhiCandInJet++;
931-
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE_pt_v_eta"), lResonance.Pt(), lResonance.Eta());
932-
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE_1D"), lResonance.M());
933-
if (lResonance.Pt() > 2.0 && lResonance.Pt() < 3)
934-
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE_1D_2_3"), lResonance.M());
935-
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE"), jetpt, lResonance.Pt(), lResonance.M());
945+
if (cfgDaughterQAHists) {
946+
if (lResonance.M() > 1.005 && lResonance.M() < 1.035)
947+
RealPhiCandInJet++;
948+
}
949+
if (cfgDaughterQAHists) {
950+
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE_pt_v_eta"), lResonance.Pt(), lResonance.Eta());
951+
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE_1D"), lResonance.M());
952+
if (lResonance.Pt() > 2.0 && lResonance.Pt() < 3)
953+
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE_1D_2_3"), lResonance.M());
954+
JEhistos.fill(HIST("hMCRec_nonmatch_hUSS_INSIDE"), jetpt, lResonance.Pt(), lResonance.M());
955+
}
936956
}
937957

938958
if (hasJets) {
939-
JEhistos.fill(HIST("ptJEHistogramPhi_JetTrigger"), lResonance.Pt());
959+
if (cfgPIDQAHists)
960+
JEhistos.fill(HIST("ptJEHistogramPhi_JetTrigger"), lResonance.Pt());
961+
940962
auto triggerjet = std::min_element(mcd_pt.begin(), mcd_pt.end());
941963
double triggerjet_pt = *triggerjet;
942-
JEhistos.fill(HIST("JetVsPhi_REC"), triggerjet_pt, lResonance.Pt());
964+
if (!cfgIsKstar) {
965+
JEhistos.fill(HIST("JetVsPhi_REC"), triggerjet_pt, lResonance.Pt());
966+
}
967+
}
968+
if (cfgDaughterQAHists) {
969+
JEhistos.fill(HIST("minvJEHistogramPhi"), lResonance.M());
943970
}
944-
JEhistos.fill(HIST("minvJEHistogramPhi"), lResonance.M());
945971
}
946972
} // mcpart check
947973
} // tracks2
948974
} // tracks1
949-
JEhistos.fill(HIST("hNRealPhiVPhiCand"), PhiCand, RealPhiCand);
950-
JEhistos.fill(HIST("hNRealPhiWithJetVPhiCand"), PhiCand, RealPhiCandWithJet);
951-
JEhistos.fill(HIST("hNRealPhiInJetVPhiCand"), PhiCand, RealPhiCandInJet);
952-
975+
if (cfgDaughterQAHists) {
976+
JEhistos.fill(HIST("hNRealPhiVPhiCand"), PhiCand, RealPhiCand);
977+
JEhistos.fill(HIST("hNRealPhiWithJetVPhiCand"), PhiCand, RealPhiCandWithJet);
978+
JEhistos.fill(HIST("hNRealPhiInJetVPhiCand"), PhiCand, RealPhiCandInJet);
979+
}
953980
// Jet Eff
954981
}
955982
PROCESS_SWITCH(phiInJets, processRec, "pikp detector level MC JE", true);

0 commit comments

Comments
 (0)