Skip to content

Commit

Permalink
Merge pull request #11 from alibuild/alibot-cleanup-8158
Browse files Browse the repository at this point in the history 
Please consider the following formatting changes to AliceO2Group#8158
  • Loading branch information
@ljoergen
ljoergen authored Oct 29, 2024
2 parents a25f256 + d27b13a commit c6c360b
Showing 1 changed file with 21 additions and 22 deletions.
43 changes: 21 additions & 22 deletions PWGLF/Tasks/Nuspex/AngularCorrelationsInJets.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -169,7 +169,7 @@ struct AngularCorrelationsInJets {
registryData.add("hEventProtocol", "Event protocol", HistType::kTH1I, {{20, 0, 20}});
registryData.add("hTrackProtocol", "Track protocol", HistType::kTH1I, {{20, 0, 20}});
registryData.add("hNumPartInJet", "Number of particles in a jet", HistType::kTH1I, {{200, 0, 200}});
registryData.add("hNumJetsInEvent", "Number of jets selected", HistType::kTH1I, {{10,0,10}});
registryData.add("hNumJetsInEvent", "Number of jets selected", HistType::kTH1I, {{10, 0, 10}});

// (Pseudo)Rapidity
registryData.add("hJetRapidity", "Jet rapidity;#it{y}", HistType::kTH1F, {{200, -1, 1}});
Expand Down Expand Up @@ -217,14 +217,14 @@ struct AngularCorrelationsInJets {
registryData.add("hDCAzJetAntinuclei", "DCA_{z} of high purity antinuclei", HistType::kTH2F, {axisSpecs.ptAxisPos, axisSpecs.dcazAxis});

// Angular Distributions
registryQA.add("hPhiFullEvent", "#varphi in full event", HistType::kTH1F, {{1000,-6.3,6.3}});
registryQA.add("hPhiPtFullEvent", "#varphi vs. p_{T} in full event", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000,-6.3,6.3}});
registryQA.add("hPhiJet", "#varphi in jet", HistType::kTH1F, {{1000,-6.3,6.3}});
registryQA.add("hPhiPtJet", "#varphi vs. p_{T} in jet", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000,-6.3,6.3}});
registryQA.add("hEtaFullEvent", "#eta in full event", HistType::kTH1F, {{1000,-1,1}});
registryQA.add("hEtaPtFullEvent", "#eta vs. p_{T} in full event", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000,-1,1}});
registryQA.add("hEtaJet", "#eta in jet", HistType::kTH1F, {{1000,-1,1}});
registryQA.add("hEtaPtJet", "#eta vs. p_{T} in jet", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000,-1,1}});
registryQA.add("hPhiFullEvent", "#varphi in full event", HistType::kTH1F, {{1000, -6.3, 6.3}});
registryQA.add("hPhiPtFullEvent", "#varphi vs. p_{T} in full event", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000, -6.3, 6.3}});
registryQA.add("hPhiJet", "#varphi in jet", HistType::kTH1F, {{1000, -6.3, 6.3}});
registryQA.add("hPhiPtJet", "#varphi vs. p_{T} in jet", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000, -6.3, 6.3}});
registryQA.add("hEtaFullEvent", "#eta in full event", HistType::kTH1F, {{1000, -1, 1}});
registryQA.add("hEtaPtFullEvent", "#eta vs. p_{T} in full event", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000, -1, 1}});
registryQA.add("hEtaJet", "#eta in jet", HistType::kTH1F, {{1000, -1, 1}});
registryQA.add("hEtaPtJet", "#eta vs. p_{T} in jet", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000, -1, 1}});

registryData.add("hDeltaPhiSEFull", "#Delta#varphi of particles in single event", HistType::kTH1D, {axisSpecs.angDistPhiAxis});
registryData.add("hDeltaPhiSEJet", "#Delta#varphi of jet particles in single event", HistType::kTH1D, {axisSpecs.angDistPhiAxis});
Expand Down Expand Up @@ -254,12 +254,12 @@ struct AngularCorrelationsInJets {
// QA
registryQA.add("hPtDiff", "p_{T} difference PseudoJet/original track;#it{p}_{T} [GeV/#it{c}]", HistType::kTH1D, {{100, -5, 5}});
registryQA.add("hJetConeRadius", "Jet Radius;#it{R}", HistType::kTH1F, {{100, 0, 1}});
registryQA.add("hMaxRadiusVsPt", "Max Cone Radius vs p_{T}", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {100,0,1}});
registryQA.add("hRhoEstimatePerp", "Background #rho (perp)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200,0,20}});
registryQA.add("hRhoMEstimatePerp", "Background #rho_{m} (perp)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200,0,20}});
registryQA.add("hRhoEstimateArea", "Background #rho (area)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200,0,20}});
registryQA.add("hRhoMEstimateArea", "Background #rho_{m} (area)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200,0,20}});
registryQA.add("hJetBkgDeltaPt", "#Delta p_{T} Clustered Cone - Pure Jet", HistType::kTH1F, {{200,0,10}});
registryQA.add("hMaxRadiusVsPt", "Max Cone Radius vs p_{T}", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {100, 0, 1}});
registryQA.add("hRhoEstimatePerp", "Background #rho (perp)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200, 0, 20}});
registryQA.add("hRhoMEstimatePerp", "Background #rho_{m} (perp)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200, 0, 20}});
registryQA.add("hRhoEstimateArea", "Background #rho (area)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200, 0, 20}});
registryQA.add("hRhoMEstimateArea", "Background #rho_{m} (area)", HistType::kTH2F, {{axisSpecs.ptAxisPos}, {200, 0, 20}});
registryQA.add("hJetBkgDeltaPt", "#Delta p_{T} Clustered Cone - Pure Jet", HistType::kTH1F, {{200, 0, 10}});

registryQA.add("hTOFmass", "TOF mass", HistType::kTH2F, {axisSpecs.ptAxisPos, {1000, 0, 5, "#it{m} [GeV/#it{c}^{2}]"}});
registryQA.get<TH2>(HIST("hTOFmass"))->Sumw2();
Expand Down Expand Up @@ -796,7 +796,7 @@ struct AngularCorrelationsInJets {
continue;
double DeltaPhi = TVector2::Phi_0_2pi(constituent.phi() - jet.phi());
if (DeltaPhi > TMath::Pi())
DeltaPhi = DeltaPhi - 2*TMath::Pi();
DeltaPhi = DeltaPhi - 2 * TMath::Pi();
double DeltaEta = constituent.eta() - jet.eta();
double Delta = TMath::Sqrt(DeltaPhi * DeltaPhi + DeltaEta * DeltaEta);
registryQA.fill(HIST("hJetConeRadius"), Delta);
Expand All @@ -806,7 +806,7 @@ struct AngularCorrelationsInJets {
registryQA.fill(HIST("hMaxRadiusVsPt"), jet.pt(), maxRadius); // no entries - weird!

// QA for comparison with nuclei_in_jets
TVector3 pJet(0.,0.,0.);
TVector3 pJet(0., 0., 0.);
pJet.SetXYZ(jet.px(), jet.py(), jet.pz());
TVector3 UEAxis1(0.0, 0.0, 0.0);
TVector3 UEAxis2(0.0, 0.0, 0.0);
Expand Down Expand Up @@ -1032,7 +1032,7 @@ struct AngularCorrelationsInJets {
}
return jetCounter;
}

template <typename U>
void fillHistograms(U const& tracks)
{
Expand All @@ -1043,8 +1043,8 @@ struct AngularCorrelationsInJets {
std::vector<typename U::iterator> jetAll;
std::vector<std::pair<double, double>> fTempBufferFull;
fTempBufferFull.clear();
std::vector<fastjet::PseudoJet> jetInput; // input for jet finder
std::map<int, typename U::iterator> particles; // all selected particles in event
std::vector<fastjet::PseudoJet> jetInput; // input for jet finder
std::map<int, typename U::iterator> particles; // all selected particles in event
std::vector<typename U::iterator> particlesForCF; // particles for full event angular correlations
jetInput.clear();
particles.clear();
Expand Down Expand Up @@ -1091,7 +1091,6 @@ struct AngularCorrelationsInJets {
registryQA.fill(HIST("hPhiPtFullEvent"), track.pt(), track.phi());
registryQA.fill(HIST("hEtaFullEvent"), track.eta());
registryQA.fill(HIST("hEtaPtFullEvent"), track.pt(), track.eta());


fastjet::PseudoJet inputPseudoJet(track.px(), track.py(), track.pz(), track.energy(mass));
inputPseudoJet.set_user_index(index);
Expand All @@ -1101,7 +1100,7 @@ struct AngularCorrelationsInJets {

index++;
} // for (const auto& track : tracks)

if (jetInput.size() < 2)
return;
registryData.fill(HIST("hEventProtocol"), 2);
Expand Down

0 comments on commit c6c360b

Please sign in to comment.