Skip to content

Commit 3fe8909

Browse files
jkinner2jkinner2
committed
use_double_precision_for_histograms_expected_with_large_content
1 parent bdcf463 commit 3fe8909

File tree

1 file changed

+31
-31
lines changed

1 file changed

+31
-31
lines changed

‎PWGJE/Tasks/photonChargedTriggerCorrelation.cxx‎

Lines changed: 31 additions & 31 deletions
Original file line numberDiff line numberDiff line change
@@ -684,40 +684,40 @@ struct PhotonChargedTriggerCorrelation {
684684
histos.get<TH1>(HIST("reco/info/h1_nEvents"))->GetXaxis()->SetBinLabel(2, "#it{N}_{ev}");
685685
histos.get<TH1>(HIST("reco/info/h1_nEvents"))->GetXaxis()->SetBinLabel(3, "#it{N}_{ev}^{trig}");
686686

687-
histos.add("reco/info/h2_zPvMult", "h2_zPvMult", kTHnSparseF, {axisZPv, axisMult}, true);
688-
histos.add("reco/info/h1_occupancy", "h1_occupancy", kTH1F, {axisOccupancy}, true);
687+
histos.add("reco/info/h2_zPvMult", "h2_zPvMult", kTHnSparseD, {axisZPv, axisMult}, true);
688+
histos.add("reco/info/h1_occupancy", "h1_occupancy", kTH1D, {axisOccupancy}, true);
689689

690690
// reco (correlation) analysis
691-
histos.add("reco/info/h2_zPvMult_trigEv", "h2_zPvMult_trigEv", kTHnSparseF, {axisZPv, axisMult}, true);
692-
histos.add("reco/info/h1_occupancy_trigEv", "h1_occupancy_trigEv", kTH1F, {axisOccupancy}, true);
693-
histos.add("reco/corr/h3_ptPhiEta_trig", "h3_ptPhiEta_trig", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
691+
histos.add("reco/info/h2_zPvMult_trigEv", "h2_zPvMult_trigEv", kTHnSparseD, {axisZPv, axisMult}, true);
692+
histos.add("reco/info/h1_occupancy_trigEv", "h1_occupancy_trigEv", kTH1D, {axisOccupancy}, true);
693+
histos.add("reco/corr/h3_ptPhiEta_trig", "h3_ptPhiEta_trig", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
694694

695695
// hadron
696-
histos.add("reco/plain/h3_ptPhiEta_hadron", "h3_ptPhiEta_hadron", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
697-
histos.add("reco/corr/h3_ptPhiEta_assoc_hadron", "h3_ptPhiEta_assoc_hadron", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
696+
histos.add("reco/plain/h3_ptPhiEta_hadron", "h3_ptPhiEta_hadron", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
697+
histos.add("reco/corr/h3_ptPhiEta_assoc_hadron", "h3_ptPhiEta_assoc_hadron", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
698698
histos.add("reco/corr/h6_corr_hadron", "h6_corr_hadron",
699699
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
700700
histos.add("reco/corr/h6_mix_hadron", "h6_mix_hadron",
701701
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
702702
// pipm
703-
histos.add("reco/plain/h3_ptPhiEta_pipm", "h3_ptPhiEta_pipm", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
704-
histos.add("reco/corr/h3_ptPhiEta_assoc_pipm", "h3_ptPhiEta_assoc_pipm", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
703+
histos.add("reco/plain/h3_ptPhiEta_pipm", "h3_ptPhiEta_pipm", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
704+
histos.add("reco/corr/h3_ptPhiEta_assoc_pipm", "h3_ptPhiEta_assoc_pipm", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
705705
histos.add("reco/corr/h6_corr_pipm", "h6_corr_pipm",
706706
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
707707
histos.add("reco/corr/h6_mix_pipm", "h6_mix_pipm",
708708
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
709709
// photonPCM
710-
histos.add("reco/plain/h3_ptPhiEta_photonPCM", "h3_ptPhiEta_photonPCM", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
711-
histos.add("reco/corr/h3_ptPhiEta_assoc_photonPCM", "h3_ptPhiEta_assoc_photonPCM", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
710+
histos.add("reco/plain/h3_ptPhiEta_photonPCM", "h3_ptPhiEta_photonPCM", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
711+
histos.add("reco/corr/h3_ptPhiEta_assoc_photonPCM", "h3_ptPhiEta_assoc_photonPCM", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
712712
histos.add("reco/corr/h6_corr_photonPCM", "h6_corr_photonPCM",
713713
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
714714
histos.add("reco/corr/h6_mix_photonPCM", "h6_mix_photonPCM",
715715
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
716716
// pi0PCM
717-
histos.add("reco/plain/h3_ptPhiEta_pi0PCM", "h3_ptPhiEta_pi0PCM", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
718-
histos.add("reco/plain/h4_ptMggZPvMult_pi0PCM", "h4_ptMggZPvMult_pi0PCM", kTHnSparseF, {axisPtAssoc, axisMgg, axisZPvBinning, axisMultBinning}, true);
719-
histos.add("reco/corr/h3_ptPhiEta_assoc_pi0PCM", "h3_ptPhiEta_assoc_pi0PCM", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
720-
histos.add("reco/corr/h4_ptMggZPvMult_assoc_pi0PCM", "h4_ptMggZPvMult_assoc_pi0PCM", kTHnSparseF, {axisPtAssoc, axisMgg, axisZPvBinning, axisMultBinning}, true);
717+
histos.add("reco/plain/h3_ptPhiEta_pi0PCM", "h3_ptPhiEta_pi0PCM", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
718+
histos.add("reco/plain/h4_ptMggZPvMult_pi0PCM", "h4_ptMggZPvMult_pi0PCM", kTHnSparseD, {axisPtAssoc, axisMgg, axisZPvBinning, axisMultBinning}, true);
719+
histos.add("reco/corr/h3_ptPhiEta_assoc_pi0PCM", "h3_ptPhiEta_assoc_pi0PCM", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
720+
histos.add("reco/corr/h4_ptMggZPvMult_assoc_pi0PCM", "h4_ptMggZPvMult_assoc_pi0PCM", kTHnSparseD, {axisPtAssoc, axisMgg, axisZPvBinning, axisMultBinning}, true);
721721
histos.add("reco/corr/h6_corr_pi0PCM", "h6_corr_pi0PCM",
722722
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
723723
histos.add("reco/corr/h6_mix_pi0PCM", "h6_mix_pi0PCM",
@@ -728,49 +728,49 @@ struct PhotonChargedTriggerCorrelation {
728728
histos.add("reco/corr/h6_mix_pi0PCMSide", "h6_mix_pi0PCMSide",
729729
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc, axisZPvBinning, axisMultBinning}, true);
730730
// event mixing for photon pairs
731-
histos.add("reco/plain/h2_zPvMult_pi0PCMMix", "h2_zPvMult_pi0PCMMix", kTHnSparseF, {axisZPv, axisMult}, true);
732-
histos.add("reco/plain/h3_ptPhiEta_pi0PCMMix", "h3_ptPhiEta_pi0PCMMix", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
733-
histos.add("reco/plain/h4_ptMggZPvMult_pi0PCMMix", "h4_ptMggZPvMult_pi0PCMMix", kTHnSparseF, {axisPtAssoc, axisMgg, axisZPvBinning, axisMultBinning}, true);
731+
histos.add("reco/plain/h2_zPvMult_pi0PCMMix", "h2_zPvMult_pi0PCMMix", kTHnSparseD, {axisZPv, axisMult}, true);
732+
histos.add("reco/plain/h3_ptPhiEta_pi0PCMMix", "h3_ptPhiEta_pi0PCMMix", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
733+
histos.add("reco/plain/h4_ptMggZPvMult_pi0PCMMix", "h4_ptMggZPvMult_pi0PCMMix", kTHnSparseD, {axisPtAssoc, axisMgg, axisZPvBinning, axisMultBinning}, true);
734734

735735
// mc info
736736
histos.add("mc/info/h1_nEvents_mcTrue", "h1_nEvents_mcTrue", kTH1D, {axisN});
737737
histos.add("mc/info/h1_nTriggerEvents_mcTrue", "h1_nTriggerEvents_mcTrue", kTH1D, {axisN});
738738

739-
histos.add("mc/info/h1_zPv_mcTrue", "h1_zPv_mcTrue", kTH1F, {axisZPv}, true);
740-
histos.add("mc/info/h1_mult_mcTrue", "h1_mult_mcTrue", kTH1F, {axisMult}, true);
739+
histos.add("mc/info/h1_zPv_mcTrue", "h1_zPv_mcTrue", kTH1D, {axisZPv}, true);
740+
histos.add("mc/info/h1_mult_mcTrue", "h1_mult_mcTrue", kTH1D, {axisMult}, true);
741741

742742
// reco and true collision correlations
743743
for (auto const& collision_type : {"true", "true_reco"}) {
744-
histos.add(std::format("mc/{}/corr/h3_ptPhiEta_trig", collision_type).data(), "h3_ptPhiEta_trig", kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
744+
histos.add(std::format("mc/{}/corr/h3_ptPhiEta_trig", collision_type).data(), "h3_ptPhiEta_trig", kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
745745
// hadron
746746
histos.add(std::format("mc/{}/corr/h3_ptPhiEta_assoc_hadron", collision_type).data(), "h3_ptPhiEta_assoc_hadron",
747-
kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
747+
kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
748748
histos.add(std::format("mc/{}/corr/h4_corr_hadron", collision_type).data(), "h4_corr_hadron",
749-
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
749+
kTHnSparseD, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
750750
// pipm
751751
histos.add(std::format("mc/{}/corr/h3_ptPhiEta_assoc_pipm", collision_type).data(), "h3_ptPhiEta_assoc_pipm",
752-
kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
752+
kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
753753
histos.add(std::format("mc/{}/corr/h4_corr_pipm", collision_type).data(), "h4_corr_pipm",
754-
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
754+
kTHnSparseD, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
755755
// photon
756756
histos.add(std::format("mc/{}/corr/h3_ptPhiEta_assoc_photon", collision_type).data(), "h3_ptPhiEta_assoc_photon",
757-
kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
757+
kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
758758
histos.add(std::format("mc/{}/corr/h4_corr_photon", collision_type).data(), "h4_corr_photon",
759-
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
759+
kTHnSparseD, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
760760
// pi0
761761
histos.add(std::format("mc/{}/corr/h3_ptPhiEta_assoc_pi0", collision_type).data(), "h3_ptPhiEta_assoc_pi0",
762-
kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
762+
kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
763763
histos.add(std::format("mc/{}/corr/h4_corr_pi0", collision_type).data(), "h4_corr_pi0",
764-
kTHnSparseF, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
764+
kTHnSparseD, {axisDPhi, axisDEta, axisPtTrig, axisPtAssoc}, true);
765765
}
766766

767767
// mc efficiency/purity
768768
std::function<void(std::string)> add_effHists =
769769
[&](std::string name_id) {
770770
histos.add(std::format("mc/eff/h3_ptPhiEta_{}", name_id).data(), "h3_ptPhiEta_mcReco_hadron",
771-
kTHnSparseF, {axisPtAssoc, axisPhi, axisEta}, true);
771+
kTHnSparseD, {axisPtAssoc, axisPhi, axisEta}, true);
772772
histos.add(std::format("mc/eff/h3_ptZPvMult_{}", name_id).data(), "h3_ptZPvMult_mcReco_hadron",
773-
kTHnSparseF, {axisPtAssoc, axisZPvSmol, axisMultSmol}, true);
773+
kTHnSparseD, {axisPtAssoc, axisZPvSmol, axisMultSmol}, true);
774774
};
775775
// mc tracks
776776
add_effHists("mcReco_hadron");

0 commit comments

Comments
 (0)