diff --git a/PWGEM/Dilepton/Core/DielectronCut.cxx b/PWGEM/Dilepton/Core/DielectronCut.cxx index 006a0fa5252..51d8d52fce5 100644 --- a/PWGEM/Dilepton/Core/DielectronCut.cxx +++ b/PWGEM/Dilepton/Core/DielectronCut.cxx @@ -206,7 +206,7 @@ void DielectronCut::ApplyPhiV(bool flag) void DielectronCut::SetPIDScheme(int scheme) { mPIDScheme = scheme; - LOG(info) << "Dielectron Cut, PID scheme: " << static_cast(mPIDScheme); + LOG(info) << "Dielectron Cut, PID scheme: " << mPIDScheme; } void DielectronCut::SetMinPinTOF(float min) { diff --git a/PWGEM/Dilepton/Core/Dilepton.h b/PWGEM/Dilepton/Core/Dilepton.h index 738634eeb3f..b39ba4da2f6 100644 --- a/PWGEM/Dilepton/Core/Dilepton.h +++ b/PWGEM/Dilepton/Core/Dilepton.h @@ -567,17 +567,17 @@ struct Dilepton { if (cfgApplySPresolution) { auto list = ccdb->getForTimeStamp(spresoPath, collision.timestamp()); - h2sp_resolution = reinterpret_cast(list->FindObject(spresoHistName.value.data())); + h2sp_resolution = dynamic_cast(list->FindObject(spresoHistName.value.data())); LOGF(info, "h2sp_resolution.GetBinContent(40, 1) = %f", h2sp_resolution->GetBinContent(40, 1)); } if (flowcorrectionGroup.cfgApplyWeightNUA) { auto list_nua = ccdb->getForTimeStamp(flowcorrectionGroup.nuaPath, collision.timestamp()); - h3nua_weights.push_back(reinterpret_cast(list_nua->FindObject("weights_uls_NUA"))); + h3nua_weights.push_back(dynamic_cast(list_nua->FindObject("weights_uls_NUA"))); LOGF(info, "h3nua_weights.at(0)->GetBinContent(5, 5, 5) = %f", h3nua_weights.at(0)->GetBinContent(5, 5, 5)); - h3nua_weights.push_back(reinterpret_cast(list_nua->FindObject("weights_lspp_NUA"))); + h3nua_weights.push_back(dynamic_cast(list_nua->FindObject("weights_lspp_NUA"))); LOGF(info, "h3nua_weights.at(1)->GetBinContent(5, 5, 5) = %f", h3nua_weights.at(1)->GetBinContent(5, 5, 5)); - h3nua_weights.push_back(reinterpret_cast(list_nua->FindObject("weights_lsmm_NUA"))); + h3nua_weights.push_back(dynamic_cast(list_nua->FindObject("weights_lsmm_NUA"))); LOGF(info, "h3nua_weights.at(2)->GetBinContent(5, 5, 5) = %f", h3nua_weights.at(2)->GetBinContent(5, 5, 5)); } } @@ -1309,33 +1309,33 @@ struct Dilepton { o2::framework::expressions::Filter pidFilter_electron = dielectroncuts.cfg_min_TPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < dielectroncuts.cfg_max_TPCNsigmaEl; o2::framework::expressions::Filter ttcaFilter_electron = ifnode(dielectroncuts.enableTTCA.node(), true, o2::aod::emprimaryelectron::isAssociatedToMPC == true); o2::framework::expressions::Filter prefilter_derived_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter_derived.node() && dielectroncuts.cfg_prefilter_bits_derived.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfbderived >= static_cast(0)); o2::framework::expressions::Filter prefilter_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter.node() && dielectroncuts.cfg_prefilter_bits.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfb >= static_cast(0)); - o2::framework::Partition positive_electrons = o2::aod::emprimaryelectron::sign > int8_t(0); - o2::framework::Partition negative_electrons = o2::aod::emprimaryelectron::sign < int8_t(0); + o2::framework::Partition positive_electrons = o2::aod::emprimaryelectron::sign > static_cast(0); + o2::framework::Partition negative_electrons = o2::aod::emprimaryelectron::sign < static_cast(0); // o2::framework::Preslice perCollision_muon = o2::aod::emprimarymuon::emeventId; o2::framework::Preslice perCollision_muon = o2::aod::emprimarymuon::emeventId; o2::framework::expressions::Filter trackFilter_muon = o2::aod::fwdtrack::trackType == dimuoncuts.cfg_track_type && dimuoncuts.cfg_min_pt_track < o2::aod::fwdtrack::pt && o2::aod::fwdtrack::pt < dimuoncuts.cfg_max_pt_track && dimuoncuts.cfg_min_eta_track < o2::aod::fwdtrack::eta && o2::aod::fwdtrack::eta < dimuoncuts.cfg_max_eta_track; o2::framework::expressions::Filter ttcaFilter_muon = ifnode(dimuoncuts.enableTTCA.node(), true, o2::aod::emprimarymuon::isAssociatedToMPC == true); - o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > int8_t(0); - o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < int8_t(0); + o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > static_cast(0); + o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < static_cast(0); TEMH* emh_pos = nullptr; TEMH* emh_neg = nullptr; diff --git a/PWGEM/Dilepton/Core/DileptonHadronMPC.h b/PWGEM/Dilepton/Core/DileptonHadronMPC.h index 0825c1d2f40..c9f4f87cccd 100644 --- a/PWGEM/Dilepton/Core/DileptonHadronMPC.h +++ b/PWGEM/Dilepton/Core/DileptonHadronMPC.h @@ -985,32 +985,32 @@ struct DileptonHadronMPC { o2::framework::expressions::Filter trackFilter_electron = dielectroncuts.cfg_min_pt_track < o2::aod::track::pt && dielectroncuts.cfg_min_eta_track < o2::aod::track::eta && o2::aod::track::eta < dielectroncuts.cfg_max_eta_track && nabs(o2::aod::track::dcaXY) < dielectroncuts.cfg_max_dcaxy && nabs(o2::aod::track::dcaZ) < dielectroncuts.cfg_max_dcaz; o2::framework::expressions::Filter ttcaFilter_electron = ifnode(dielectroncuts.enableTTCA.node(), o2::aod::emprimaryelectron::isAssociatedToMPC == true || o2::aod::emprimaryelectron::isAssociatedToMPC == false, o2::aod::emprimaryelectron::isAssociatedToMPC == true); o2::framework::expressions::Filter prefilter_derived_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter_derived.node() && dielectroncuts.cfg_prefilter_bits_derived.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfbderived >= static_cast(0)); o2::framework::expressions::Filter prefilter_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter.node() && dielectroncuts.cfg_prefilter_bits.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfb >= static_cast(0)); - o2::framework::Partition positive_electrons = o2::aod::emprimaryelectron::sign > int8_t(0); - o2::framework::Partition negative_electrons = o2::aod::emprimaryelectron::sign < int8_t(0); + o2::framework::Partition positive_electrons = o2::aod::emprimaryelectron::sign > static_cast(0); + o2::framework::Partition negative_electrons = o2::aod::emprimaryelectron::sign < static_cast(0); o2::framework::Preslice perCollision_muon = o2::aod::emprimarymuon::emeventId; o2::framework::expressions::Filter trackFilter_muon = o2::aod::fwdtrack::trackType == dimuoncuts.cfg_track_type && dimuoncuts.cfg_min_pt_track < o2::aod::fwdtrack::pt && o2::aod::fwdtrack::pt < dimuoncuts.cfg_max_pt_track && dimuoncuts.cfg_min_eta_track < o2::aod::fwdtrack::eta && o2::aod::fwdtrack::eta < dimuoncuts.cfg_max_eta_track && dimuoncuts.cfg_min_phi_track < o2::aod::fwdtrack::phi && o2::aod::fwdtrack::phi < dimuoncuts.cfg_max_phi_track; o2::framework::expressions::Filter ttcaFilter_muon = ifnode(dimuoncuts.enableTTCA.node(), true, o2::aod::emprimarymuon::isAssociatedToMPC == true); - o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > int8_t(0); - o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < int8_t(0); + o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > static_cast(0); + o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < static_cast(0); using RefTracks = o2::soa::Join; using RefTrack = RefTracks::iterator; diff --git a/PWGEM/Dilepton/Core/DileptonMC.h b/PWGEM/Dilepton/Core/DileptonMC.h index 219d6c0dec3..28a323f2345 100644 --- a/PWGEM/Dilepton/Core/DileptonMC.h +++ b/PWGEM/Dilepton/Core/DileptonMC.h @@ -2527,21 +2527,21 @@ struct DileptonMC { o2::framework::expressions::Filter pidFilter_electron = dielectroncuts.cfg_min_TPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < dielectroncuts.cfg_max_TPCNsigmaEl; o2::framework::expressions::Filter ttcaFilter_electron = ifnode(dielectroncuts.enableTTCA.node(), true, o2::aod::emprimaryelectron::isAssociatedToMPC == true); o2::framework::expressions::Filter prefilter_derived_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter_derived.node() && dielectroncuts.cfg_prefilter_bits_derived.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfbderived >= static_cast(0)); o2::framework::expressions::Filter prefilter_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter.node() && dielectroncuts.cfg_prefilter_bits.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfb >= static_cast(0)); o2::framework::Preslice perCollision_muon = o2::aod::emprimarymuon::emeventId; @@ -2554,10 +2554,10 @@ struct DileptonMC { o2::framework::expressions::Filter collisionFilter_occupancy_ft0c = eventcuts.cfgFT0COccupancyMin <= o2::aod::evsel::ft0cOccupancyInTimeRange && o2::aod::evsel::ft0cOccupancyInTimeRange < eventcuts.cfgFT0COccupancyMax; using FilteredMyCollisions = o2::soa::Filtered; - o2::framework::Partition positive_electrons = o2::aod::emprimaryelectron::sign > int8_t(0); // reconstructed tracks - o2::framework::Partition negative_electrons = o2::aod::emprimaryelectron::sign < int8_t(0); // reconstructed tracks - o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > int8_t(0); // reconstructed tracks - o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < int8_t(0); // reconstructed tracks + o2::framework::Partition positive_electrons = o2::aod::emprimaryelectron::sign > static_cast(0); // reconstructed tracks + o2::framework::Partition negative_electrons = o2::aod::emprimaryelectron::sign < static_cast(0); // reconstructed tracks + o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > static_cast(0); // reconstructed tracks + o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < static_cast(0); // reconstructed tracks o2::framework::Partition positive_electronsMC = o2::aod::mcparticle::pdgCode == -11; // e+ o2::framework::Partition negative_electronsMC = o2::aod::mcparticle::pdgCode == 11; // e- diff --git a/PWGEM/Dilepton/Core/SingleTrackQC.h b/PWGEM/Dilepton/Core/SingleTrackQC.h index 0530a64e6e2..7a23c3dfd82 100644 --- a/PWGEM/Dilepton/Core/SingleTrackQC.h +++ b/PWGEM/Dilepton/Core/SingleTrackQC.h @@ -896,21 +896,21 @@ struct SingleTrackQC { o2::framework::expressions::Filter ttcaFilter_electron = ifnode(dielectroncuts.enableTTCA.node(), true, o2::aod::emprimaryelectron::isAssociatedToMPC == true); o2::framework::expressions::Filter prefilter_derived_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter_derived.node() && dielectroncuts.cfg_prefilter_bits_derived.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfbderived >= static_cast(0)); o2::framework::expressions::Filter prefilter_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter.node() && dielectroncuts.cfg_prefilter_bits.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfb >= static_cast(0)); o2::framework::Preslice perCollision_muon = o2::aod::emprimarymuon::emeventId; diff --git a/PWGEM/Dilepton/Core/SingleTrackQCMC.h b/PWGEM/Dilepton/Core/SingleTrackQCMC.h index 314265b4e6b..9ad6c889cbc 100644 --- a/PWGEM/Dilepton/Core/SingleTrackQCMC.h +++ b/PWGEM/Dilepton/Core/SingleTrackQCMC.h @@ -1162,21 +1162,21 @@ struct SingleTrackQCMC { o2::framework::expressions::Filter ttcaFilter_electron = ifnode(dielectroncuts.enableTTCA.node(), true, o2::aod::emprimaryelectron::isAssociatedToMPC == true); o2::framework::expressions::Filter prefilter_derived_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter_derived.node() && dielectroncuts.cfg_prefilter_bits_derived.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfbderived >= static_cast(0)); o2::framework::expressions::Filter prefilter_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter.node() && dielectroncuts.cfg_prefilter_bits.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfb >= static_cast(0)); o2::framework::Preslice perCollision_muon = o2::aod::emprimarymuon::emeventId; diff --git a/PWGEM/Dilepton/TableProducer/prefilterDielectron.cxx b/PWGEM/Dilepton/TableProducer/prefilterDielectron.cxx index 8460b0caf5c..7e6c444740c 100644 --- a/PWGEM/Dilepton/TableProducer/prefilterDielectron.cxx +++ b/PWGEM/Dilepton/TableProducer/prefilterDielectron.cxx @@ -318,8 +318,8 @@ struct prefilterDielectron { SliceCache cache; Preslice perCollision_track = aod::emprimaryelectron::emeventId; - Partition posTracks = o2::aod::emprimaryelectron::sign > int8_t(0); - Partition negTracks = o2::aod::emprimaryelectron::sign < int8_t(0); + Partition posTracks = o2::aod::emprimaryelectron::sign > static_cast(0); + Partition negTracks = o2::aod::emprimaryelectron::sign < static_cast(0); Filter collisionFilter_centrality = (cfgCentMin < o2::aod::cent::centFT0M && o2::aod::cent::centFT0M < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0A && o2::aod::cent::centFT0A < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0C && o2::aod::cent::centFT0C < cfgCentMax); Filter collisionFilter_occupancy_track = eventcuts.cfgTrackOccupancyMin <= o2::aod::evsel::trackOccupancyInTimeRange && o2::aod::evsel::trackOccupancyInTimeRange < eventcuts.cfgTrackOccupancyMax; diff --git a/PWGEM/Dilepton/TableProducer/prefilterDimuon.cxx b/PWGEM/Dilepton/TableProducer/prefilterDimuon.cxx index aa8ab7bdc85..f9ddcc8d1e4 100644 --- a/PWGEM/Dilepton/TableProducer/prefilterDimuon.cxx +++ b/PWGEM/Dilepton/TableProducer/prefilterDimuon.cxx @@ -256,8 +256,8 @@ struct prefilterDimuon { SliceCache cache; Preslice perCollision_track = aod::emprimarymuon::emeventId; - Partition posTracks = o2::aod::emprimarymuon::sign > int8_t(0); - Partition negTracks = o2::aod::emprimarymuon::sign < int8_t(0); + Partition posTracks = o2::aod::emprimarymuon::sign > static_cast(0); + Partition negTracks = o2::aod::emprimarymuon::sign < static_cast(0); Filter collisionFilter_centrality = (cfgCentMin < o2::aod::cent::centFT0M && o2::aod::cent::centFT0M < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0A && o2::aod::cent::centFT0A < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0C && o2::aod::cent::centFT0C < cfgCentMax); Filter collisionFilter_occupancy_track = eventcuts.cfgTrackOccupancyMin <= o2::aod::evsel::trackOccupancyInTimeRange && o2::aod::evsel::trackOccupancyInTimeRange < eventcuts.cfgTrackOccupancyMax; diff --git a/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx b/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx index 78a104995b3..aaea7f087c3 100644 --- a/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx +++ b/PWGEM/Dilepton/TableProducer/skimmerPrimaryElectron.cxx @@ -1547,8 +1547,8 @@ struct prefilterPrimaryElectron { Partition posTracks = o2::aod::track::signed1Pt > 0.f; Partition negTracks = o2::aod::track::signed1Pt < 0.f; - Partition positrons = o2::aod::emprimaryelectron::sign > int8_t(0); - Partition electrons = o2::aod::emprimaryelectron::sign < int8_t(0); + Partition positrons = o2::aod::emprimaryelectron::sign > static_cast(0); + Partition electrons = o2::aod::emprimaryelectron::sign < static_cast(0); void processPrefilter_TTCA(MyCollisions const& collisions, aod::BCsWithTimestamps const&, MyTracks const&, aod::EMPrimaryElectrons const& primaryelectrons, aod::TrackAssoc const& trackIndices, aod::EMTOFNSigmas const& emtofs, aod::EMMLPIDs const& emmlpids) { @@ -1632,12 +1632,12 @@ struct prefilterPrimaryElectron { } for (int i = 0; i < static_cast(max_mee_vec.size()); i++) { if (v12.M() < max_mee_vec.at(i)) { - pfb_map[empos.globalIndex()] |= (uint8_t(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); + pfb_map[empos.globalIndex()] |= (static_cast(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); } } if (v12.M() < slope * phiv + intercept) { - pfb_map[empos.globalIndex()] |= (uint8_t(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); + pfb_map[empos.globalIndex()] |= (static_cast(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); } } // end of signal positon loop @@ -1674,12 +1674,12 @@ struct prefilterPrimaryElectron { } for (int i = 0; i < static_cast(max_mee_vec.size()); i++) { if (v12.M() < max_mee_vec.at(i)) { - pfb_map[emele.globalIndex()] |= (uint8_t(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); + pfb_map[emele.globalIndex()] |= (static_cast(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); } } if (v12.M() < slope * phiv + intercept) { - pfb_map[emele.globalIndex()] |= (uint8_t(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); + pfb_map[emele.globalIndex()] |= (static_cast(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); } } // end of signal electron loop } // end of loose positon loop @@ -1866,12 +1866,12 @@ struct prefilterPrimaryElectron { } for (int i = 0; i < static_cast(max_mee_vec.size()); i++) { if (v12.M() < max_mee_vec.at(i)) { - pfb_map[empos.globalIndex()] |= (uint8_t(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); + pfb_map[empos.globalIndex()] |= (static_cast(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); } } if (v12.M() < slope * phiv + intercept) { - pfb_map[empos.globalIndex()] |= (uint8_t(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); + pfb_map[empos.globalIndex()] |= (static_cast(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); } } // end of ULS pairing @@ -1898,12 +1898,12 @@ struct prefilterPrimaryElectron { } for (int i = 0; i < static_cast(max_mee_vec.size()); i++) { if (v12.M() < max_mee_vec.at(i)) { - pfb_map[emele.globalIndex()] |= (uint8_t(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); + pfb_map[emele.globalIndex()] |= (static_cast(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); } } if (v12.M() < slope * phiv + intercept) { - pfb_map[emele.globalIndex()] |= (uint8_t(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); + pfb_map[emele.globalIndex()] |= (static_cast(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); } } // end of ULS pairing diff --git a/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx b/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx index b3dbb015003..18b623375e2 100644 --- a/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx +++ b/PWGEM/Dilepton/TableProducer/skimmerPrimaryMFTTrack.cxx @@ -147,7 +147,7 @@ struct skimmerPrimaryMFTTrack { } o2::mch::TrackExtrap::setField(); const double centerMFT[3] = {0, 0, -61.4}; - o2::field::MagneticField* field = static_cast(TGeoGlobalMagField::Instance()->GetField()); + o2::field::MagneticField* field = dynamic_cast(TGeoGlobalMagField::Instance()->GetField()); mBz = field->getBz(centerMFT); // Get field at centre of MFT LOGF(info, "Bz at center of MFT = %f kZG", mBz); diff --git a/PWGEM/Dilepton/TableProducer/skimmerPrimaryMuon.cxx b/PWGEM/Dilepton/TableProducer/skimmerPrimaryMuon.cxx index 8ae3d9cd690..891e3bf941c 100644 --- a/PWGEM/Dilepton/TableProducer/skimmerPrimaryMuon.cxx +++ b/PWGEM/Dilepton/TableProducer/skimmerPrimaryMuon.cxx @@ -174,7 +174,7 @@ struct skimmerPrimaryMuon { } o2::mch::TrackExtrap::setField(); const double centerMFT[3] = {0, 0, -61.4}; - o2::field::MagneticField* field = static_cast(TGeoGlobalMagField::Instance()->GetField()); + o2::field::MagneticField* field = dynamic_cast(TGeoGlobalMagField::Instance()->GetField()); mBz = field->getBz(centerMFT); // Get field at centre of MFT LOGF(info, "Bz at center of MFT = %f kZG", mBz); diff --git a/PWGEM/Dilepton/TableProducer/treeCreatorElectronMLDDA.cxx b/PWGEM/Dilepton/TableProducer/treeCreatorElectronMLDDA.cxx index bd787d9846c..946a1e55e7c 100644 --- a/PWGEM/Dilepton/TableProducer/treeCreatorElectronMLDDA.cxx +++ b/PWGEM/Dilepton/TableProducer/treeCreatorElectronMLDDA.cxx @@ -1132,7 +1132,7 @@ struct TreeCreatorElectronMLDDA { std::map, float> mapTOFBeta; // map -> tof beta //! type of V0. 0: built solely for cascades (does not pass standard V0 cuts), 1: standard 2, 3: photon-like with TPC-only use. Regular analysis should always use type 1. - Filter v0Filter = o2::aod::v0data::v0Type == uint8_t(1) && o2::aod::v0data::v0cosPA > v0cuts.cfg_min_cospa&& o2::aod::v0data::dcaV0daughters v0cuts.cfg_min_dcaxy_v0leg&& nabs(o2::aod::v0data::dcanegtopv) > v0cuts.cfg_min_dcaxy_v0leg; + Filter v0Filter = o2::aod::v0data::v0Type == static_cast(1) && o2::aod::v0data::v0cosPA > v0cuts.cfg_min_cospa&& o2::aod::v0data::dcaV0daughters v0cuts.cfg_min_dcaxy_v0leg&& nabs(o2::aod::v0data::dcanegtopv) > v0cuts.cfg_min_dcaxy_v0leg; using filteredV0s = soa::Filtered; Filter cascadeFilter = o2::aod::cascdata::dcacascdaughters < cascadecuts.cfg_max_dcadau && nabs(o2::aod::cascdata::dcanegtopv) > cascadecuts.cfg_min_dcaxy_v0leg&& nabs(o2::aod::cascdata::dcanegtopv) > cascadecuts.cfg_min_dcaxy_v0leg&& nabs(o2::aod::cascdata::dcabachtopv) > cascadecuts.cfg_min_dcaxy_bachelor; diff --git a/PWGEM/Dilepton/TableProducer/treeCreatorMuonML.cxx b/PWGEM/Dilepton/TableProducer/treeCreatorMuonML.cxx index d9ff8eee3ba..3c249959985 100644 --- a/PWGEM/Dilepton/TableProducer/treeCreatorMuonML.cxx +++ b/PWGEM/Dilepton/TableProducer/treeCreatorMuonML.cxx @@ -183,7 +183,7 @@ struct TreeCreatorMuonML { } o2::mch::TrackExtrap::setField(); const double centerMFT[3] = {0, 0, -61.4}; - o2::field::MagneticField* field = static_cast(TGeoGlobalMagField::Instance()->GetField()); + o2::field::MagneticField* field = dynamic_cast(TGeoGlobalMagField::Instance()->GetField()); mBz = field->getBz(centerMFT); // Get field at centre of MFT LOGF(info, "Bz at center of MFT = %f kZG", mBz); diff --git a/PWGEM/Dilepton/Tasks/MCtemplates.cxx b/PWGEM/Dilepton/Tasks/MCtemplates.cxx index 9c69ee77dbd..417540a4f94 100644 --- a/PWGEM/Dilepton/Tasks/MCtemplates.cxx +++ b/PWGEM/Dilepton/Tasks/MCtemplates.cxx @@ -301,7 +301,7 @@ struct AnalysisTrackSelection { int i = 0; for (auto cut = fTrackCuts.begin(); cut != fTrackCuts.end(); cut++, i++) { if ((*cut).IsSelected(VarManager::fgValues)) { - filterMap |= (uint32_t(1) << i); + filterMap |= (static_cast(1) << i); if (fConfigQA && passEvFilter) { fHistMan->FillHistClass(fHistNamesReco[i].Data(), VarManager::fgValues); } @@ -323,23 +323,23 @@ struct AnalysisTrackSelection { for (auto sig = fMCSignals.begin(); sig != fMCSignals.end(); sig++, isig++) { if constexpr ((TTrackFillMap & VarManager::ObjTypes::ReducedTrack) > 0) { if ((*sig).CheckSignal(true, track.reducedMCTrack())) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } if constexpr ((TTrackFillMap & VarManager::ObjTypes::Track) > 0) { if ((*sig).CheckSignal(true, track.template mcParticle_as())) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } } // fill histograms for (unsigned int i = 0; i < fMCSignals.size(); i++) { - if (!(mcDecision & (uint32_t(1) << i))) { + if (!(mcDecision & (static_cast(1) << i))) { continue; } for (unsigned int j = 0; j < fTrackCuts.size(); j++) { - if (filterMap & (uint8_t(1) << j)) { + if (filterMap & (static_cast(1) << j)) { fHistMan->FillHistClass(fHistNamesMCMatched[j][i].Data(), VarManager::fgValues); } } // end loop over cuts @@ -531,7 +531,7 @@ struct AnalysisSameEventPairing { for (auto& [t1, t2] : combinations(tracks1, tracks2)) { if constexpr (TPairType == VarManager::kDecayToEE) { - twoTrackFilter = uint32_t(t1.isBarrelSelected()) & uint32_t(t2.isBarrelSelected()); + twoTrackFilter = static_cast(t1.isBarrelSelected()) & static_cast(t2.isBarrelSelected()); } if (!twoTrackFilter) { // the tracks must have at least one filter bit in common to continue continue; @@ -544,12 +544,12 @@ struct AnalysisSameEventPairing { for (auto sig = fRecMCSignals.begin(); sig != fRecMCSignals.end(); sig++, isig++) { if constexpr ((TTrackFillMap & VarManager::ObjTypes::ReducedTrack) > 0) { // for skimmed DQ model if ((*sig).CheckSignal(true, t1.reducedMCTrack(), t2.reducedMCTrack())) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } if constexpr ((TTrackFillMap & VarManager::ObjTypes::Track) > 0) { // for Framework data model if ((*sig).CheckSignal(true, t1.template mcParticle_as(), t2.template mcParticle_as())) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } } // end loop over MC signals @@ -559,11 +559,11 @@ struct AnalysisSameEventPairing { // Loop over all fulfilled cuts and fill pair histograms for (unsigned int icut = 0; icut < ncuts; icut++) { - if (twoTrackFilter & (uint8_t(1) << icut)) { + if (twoTrackFilter & (static_cast(1) << icut)) { if (t1.sign() * t2.sign() < 0) { fHistMan->FillHistClass(histNames[icut][0].Data(), VarManager::fgValues); for (unsigned int isig = 0; isig < fRecMCSignals.size(); isig++) { - if (mcDecision & (uint32_t(1) << isig)) { + if (mcDecision & (static_cast(1) << isig)) { fHistMan->FillHistClass(histNamesMCmatched[icut][isig][0].Data(), VarManager::fgValues); } } @@ -571,14 +571,14 @@ struct AnalysisSameEventPairing { if (t1.sign() > 0) { fHistMan->FillHistClass(histNames[icut][1].Data(), VarManager::fgValues); for (unsigned int isig = 0; isig < fRecMCSignals.size(); isig++) { - if (mcDecision & (uint32_t(1) << isig)) { + if (mcDecision & (static_cast(1) << isig)) { fHistMan->FillHistClass(histNamesMCmatched[icut][isig][1].Data(), VarManager::fgValues); } } } else { fHistMan->FillHistClass(histNames[icut][2].Data(), VarManager::fgValues); for (unsigned int isig = 0; isig < fRecMCSignals.size(); isig++) { - if (mcDecision & (uint32_t(1) << isig)) { + if (mcDecision & (static_cast(1) << isig)) { fHistMan->FillHistClass(histNamesMCmatched[icut][isig][2].Data(), VarManager::fgValues); } } diff --git a/PWGEM/Dilepton/Tasks/checkMCPairTemplate.cxx b/PWGEM/Dilepton/Tasks/checkMCPairTemplate.cxx index 62ef21abaab..bdf5f5a3621 100644 --- a/PWGEM/Dilepton/Tasks/checkMCPairTemplate.cxx +++ b/PWGEM/Dilepton/Tasks/checkMCPairTemplate.cxx @@ -3126,29 +3126,29 @@ struct checkMCPairTemplate { Filter trackFilter_electron = nabs(o2::aod::track::dcaXY) < dielectroncuts.cfg_max_dcaxy && nabs(o2::aod::track::dcaZ) < dielectroncuts.cfg_max_dcaz; Filter ttcaFilter_electron = ifnode(dielectroncuts.enableTTCA.node(), o2::aod::emprimaryelectron::isAssociatedToMPC == true || o2::aod::emprimaryelectron::isAssociatedToMPC == false, o2::aod::emprimaryelectron::isAssociatedToMPC == true); Filter prefilter_derived_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter_derived.node() && dielectroncuts.cfg_prefilter_bits_derived.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kMee))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kPhiV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimaryelectron::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfbderived >= static_cast(0)); Filter prefilter_electron = ifnode(dielectroncuts.cfg_apply_cuts_from_prefilter.node() && dielectroncuts.cfg_prefilter_bits.node() >= static_cast(1), - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && - ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_40MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_60MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_80MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_100MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_120MeV))) <= static_cast(0), true) && + ifnode((dielectroncuts.cfg_prefilter_bits.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) > static_cast(0), (o2::aod::emprimaryelectron::pfb & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_140MeV))) <= static_cast(0), true), o2::aod::emprimaryelectron::pfb >= static_cast(0)); Preslice perCollision_muon = aod::emprimarymuon::emeventId; Filter trackFilter_muon = o2::aod::fwdtrack::trackType == dimuoncuts.cfg_track_type; Filter ttcaFilter_muon = ifnode(dimuoncuts.enableTTCA.node(), o2::aod::emprimarymuon::isAssociatedToMPC == true || o2::aod::emprimarymuon::isAssociatedToMPC == false, o2::aod::emprimarymuon::isAssociatedToMPC == true); Filter prefilter_derived_muon = ifnode(dimuoncuts.cfg_apply_cuts_from_prefilter_derived.node() && dimuoncuts.cfg_prefilter_bits_derived.node() >= static_cast(1), - ifnode((dimuoncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimarymuon::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && - ifnode((dimuoncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimarymuon::pfbderived & static_cast(1 << int(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), + ifnode((dimuoncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) > static_cast(0), (o2::aod::emprimarymuon::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackLS))) <= static_cast(0), true) && + ifnode((dimuoncuts.cfg_prefilter_bits_derived.node() & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) > static_cast(0), (o2::aod::emprimarymuon::pfbderived & static_cast(1 << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBitDerived::kSplitOrMergedTrackULS))) <= static_cast(0), true), o2::aod::emprimarymuon::pfbderived >= static_cast(0)); Filter collisionFilter_centrality = (cfgCentMin < o2::aod::cent::centFT0M && o2::aod::cent::centFT0M < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0A && o2::aod::cent::centFT0A < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0C && o2::aod::cent::centFT0C < cfgCentMax); @@ -3157,10 +3157,10 @@ struct checkMCPairTemplate { Filter collisionFilter_occupancy_ft0c = eventcuts.cfgFT0COccupancyMin <= o2::aod::evsel::ft0cOccupancyInTimeRange && o2::aod::evsel::ft0cOccupancyInTimeRange < eventcuts.cfgFT0COccupancyMax; using FilteredMyCollisions = soa::Filtered; - Partition positive_electrons = o2::aod::emprimaryelectron::sign > int8_t(0); // reconstructed tracks - Partition negative_electrons = o2::aod::emprimaryelectron::sign < int8_t(0); // reconstructed tracks - Partition positive_muons = o2::aod::emprimarymuon::sign > int8_t(0); // reconstructed tracks - Partition negative_muons = o2::aod::emprimarymuon::sign < int8_t(0); // reconstructed tracks + Partition positive_electrons = o2::aod::emprimaryelectron::sign > static_cast(0); // reconstructed tracks + Partition negative_electrons = o2::aod::emprimaryelectron::sign < static_cast(0); // reconstructed tracks + Partition positive_muons = o2::aod::emprimarymuon::sign > static_cast(0); // reconstructed tracks + Partition negative_muons = o2::aod::emprimarymuon::sign < static_cast(0); // reconstructed tracks Partition positive_electronsMC = o2::aod::mcparticle::pdgCode == -11; // e+ Partition negative_electronsMC = o2::aod::mcparticle::pdgCode == 11; // e- diff --git a/PWGEM/Dilepton/Tasks/checkMCTemplate.cxx b/PWGEM/Dilepton/Tasks/checkMCTemplate.cxx index ad299df81f7..4144e07c43a 100644 --- a/PWGEM/Dilepton/Tasks/checkMCTemplate.cxx +++ b/PWGEM/Dilepton/Tasks/checkMCTemplate.cxx @@ -316,7 +316,7 @@ struct checkMCTemplate { } o2::mch::TrackExtrap::setField(); const double centerMFT[3] = {0, 0, -61.4}; - o2::field::MagneticField* field = static_cast(TGeoGlobalMagField::Instance()->GetField()); + o2::field::MagneticField* field = dynamic_cast(TGeoGlobalMagField::Instance()->GetField()); mBzMFT = field->getBz(centerMFT); // Get field at centre of MFT LOGF(info, "Bz at center of MFT = %f kZG", mBzMFT); mRunNumber = bc.runNumber(); diff --git a/PWGEM/Dilepton/Tasks/createResolutionMap.cxx b/PWGEM/Dilepton/Tasks/createResolutionMap.cxx index 71d1d91139b..4a1c5f485f3 100644 --- a/PWGEM/Dilepton/Tasks/createResolutionMap.cxx +++ b/PWGEM/Dilepton/Tasks/createResolutionMap.cxx @@ -374,7 +374,7 @@ struct CreateResolutionMap { } o2::mch::TrackExtrap::setField(); const double centerMFT[3] = {0, 0, -61.4}; - o2::field::MagneticField* field = static_cast(TGeoGlobalMagField::Instance()->GetField()); + o2::field::MagneticField* field = dynamic_cast(TGeoGlobalMagField::Instance()->GetField()); mBzMFT = field->getBz(centerMFT); // Get field at centre of MFT LOGF(info, "Bz at center of MFT = %f kZG", mBzMFT); mRunNumber = bc.runNumber(); diff --git a/PWGEM/Dilepton/Tasks/dimuonV1.cxx b/PWGEM/Dilepton/Tasks/dimuonV1.cxx index b069475d676..c61790532fc 100644 --- a/PWGEM/Dilepton/Tasks/dimuonV1.cxx +++ b/PWGEM/Dilepton/Tasks/dimuonV1.cxx @@ -529,8 +529,8 @@ struct dimuonV1 { using FilteredMyMuons = o2::soa::Filtered; using FilteredMyMuon = FilteredMyMuons::iterator; - o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > int8_t(0); - o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < int8_t(0); + o2::framework::Partition positive_muons = o2::aod::emprimarymuon::sign > static_cast(0); + o2::framework::Partition negative_muons = o2::aod::emprimarymuon::sign < static_cast(0); // MyEMH_muon* emh_pos = nullptr; // MyEMH_muon* emh_neg = nullptr; diff --git a/PWGEM/Dilepton/Tasks/emEfficiencyEE.cxx b/PWGEM/Dilepton/Tasks/emEfficiencyEE.cxx index a739ff9ca50..d42edda01d5 100644 --- a/PWGEM/Dilepton/Tasks/emEfficiencyEE.cxx +++ b/PWGEM/Dilepton/Tasks/emEfficiencyEE.cxx @@ -269,7 +269,7 @@ struct AnalysisEventQa { for (auto& event : events) { VarManager::ResetValues(0, VarManager::kNEventWiseVariables); VarManager::FillEvent(event); - eventFilter = uint32_t(event.isEventSelected()); + eventFilter = static_cast(event.isEventSelected()); if (!eventFilter) continue; // fEvCounters++; @@ -382,7 +382,7 @@ struct AnalysisMCEvent { for (auto& event : events) { Int_t globalindexmc = -1; - eventFilter = uint32_t(event.isEventSelected()); + eventFilter = static_cast(event.isEventSelected()); if (!eventFilter) continue; @@ -715,7 +715,7 @@ struct AnalysisTrackSelection { // fill event information which might be needed in histograms that combine track and event properties VarManager::FillEvent(event); // if(!fEventCut->IsSelected(VarManager::fgValues)) continue; - eventFilter = uint32_t(event.isEventSelected()); + eventFilter = static_cast(event.isEventSelected()); if (!eventFilter) { pass = false; } @@ -923,7 +923,7 @@ struct AnalysisTrackSelection { int i = 0; for (auto cut = fTrackCuts.begin(); cut != fTrackCuts.end(); cut++, i++) { if ((*cut).IsSelected(VarManager::fgValues)) { - filterMap |= (uint32_t(1) << i); + filterMap |= (static_cast(1) << i); if (fConfigQA && pass) { fHistManQA->FillHistClass(fHistNamesRecoQA[i].Data(), VarManager::fgValues); } @@ -942,14 +942,14 @@ struct AnalysisTrackSelection { if constexpr ((TTrackFillMap & VarManager::ObjTypes::ReducedTrack) > 0) { if ((*sig).CheckSignal(true, track.reducedMCTrack())) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } if constexpr ((TTrackFillMap & VarManager::ObjTypes::Track) > 0) { if (track.has_mcParticle()) { auto mctrack = track.template mcParticle_as(); if ((*sig).CheckSignal(true, mctrack)) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } } @@ -957,11 +957,11 @@ struct AnalysisTrackSelection { // fill histograms for (unsigned int i = 0; i < fMCSignals.size(); i++) { - if (!(mcDecision & (uint32_t(1) << i))) { + if (!(mcDecision & (static_cast(1) << i))) { continue; } for (unsigned int j = 0; j < fTrackCuts.size(); j++) { - if (filterMap & (uint8_t(1) << j)) { + if (filterMap & (static_cast(1) << j)) { if (track.sign() < 0) { fHistRecNegPart[j * fMCSignals.size() + i]->Fill(track.pt(), track.eta(), track.phi()); } else { @@ -1142,7 +1142,7 @@ struct AnalysisTrackSelection { int i = 0; for (auto cut = fTrackCuts.begin(); cut != fTrackCuts.end(); cut++, i++) { if ((*cut).IsSelected(VarManager::fgValues)) { - filterMap |= (uint32_t(1) << i); + filterMap |= (static_cast(1) << i); if (fConfigQA) { fHistManQA->FillHistClass(fHistNamesRecoQA[i].Data(), VarManager::fgValues); } @@ -1159,7 +1159,7 @@ struct AnalysisTrackSelection { if (track.has_mcParticle()) { auto mctrack = track.template mcParticle_as(); if ((*sig).CheckSignal(true, mctrack)) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); mctrackindex = mctrack.globalIndex(); } } @@ -1168,7 +1168,7 @@ struct AnalysisTrackSelection { // // Double reconstructed track only for the signal (they should not be redundant or crossing!!) // for (unsigned int i = 0; i < fMCSignals.size(); i++) { - // if (!(mcDecision & (uint32_t(1) << i))) { + // if (!(mcDecision & (static_cast(1) << i))) { // continue; // } @@ -1183,7 +1183,7 @@ struct AnalysisTrackSelection { } // track cuts for (unsigned int j = 0; j < fTrackCuts.size(); j++) { - if (filterMap & (uint8_t(1) << j)) { + if (filterMap & (static_cast(1) << j)) { if (!(fRecTrackLabels[j].find(mctrackindex) != fRecTrackLabels[j].end())) { fRecTrackLabels[j][mctrackindex] = fRecCounters[j]; @@ -1199,7 +1199,7 @@ struct AnalysisTrackSelection { // fill histograms for (unsigned int i = 0; i < fMCSignals.size(); i++) { - if (!(mcDecision & (uint32_t(1) << i))) { + if (!(mcDecision & (static_cast(1) << i))) { continue; } Double_t mmcpt = -10000.; @@ -1244,7 +1244,7 @@ struct AnalysisTrackSelection { } // track cuts for (unsigned int j = 0; j < fTrackCuts.size(); j++) { - if (filterMap & (uint8_t(1) << j)) { + if (filterMap & (static_cast(1) << j)) { if (track.sign() < 0) { fHistRecNegPart[j * fMCSignals.size() + i]->Fill(track.pt(), track.eta(), track.phi()); } else { @@ -1733,7 +1733,7 @@ struct AnalysisSameEventPairing { for (auto& [t1, t2] : combinations(soa::CombinationsStrictlyUpperIndexPolicy(tracks, tracks))) { // for (auto& [t1, t2] : combinations(tracks, tracks)) { - twoTrackFilter = uint32_t(t1.isBarrelSelected()) & uint32_t(t2.isBarrelSelected()); + twoTrackFilter = static_cast(t1.isBarrelSelected()) & static_cast(t2.isBarrelSelected()); if (!twoTrackFilter) { // the tracks must have at least one filter bit in common to continue continue; @@ -1759,7 +1759,7 @@ struct AnalysisSameEventPairing { // Fill the QA for all passing tracks if (fConfigQA) { for (unsigned int j = 0; j < fTrackCuts.size(); j++) { - if (twoTrackFilter & (uint8_t(1) << j)) { + if (twoTrackFilter & (static_cast(1) << j)) { if (!fConfigFillLS) { fHistManQA->FillHistClass(Form("PairsBarrelSEPM_%s", fTrackCuts.at(j).GetName()), VarManager::fgValues); } else { @@ -1788,7 +1788,7 @@ struct AnalysisSameEventPairing { continue; } if ((*sig).CheckSignal(true, t1.reducedMCTrack(), t2.reducedMCTrack())) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } @@ -1801,7 +1801,7 @@ struct AnalysisSameEventPairing { continue; } if ((*sig).CheckSignal(true, mct1, mct2)) { - mcDecision |= (uint32_t(1) << isig); + mcDecision |= (static_cast(1) << isig); } } } @@ -1809,12 +1809,12 @@ struct AnalysisSameEventPairing { } // end of loop MC signals for (unsigned int i = 0; i < fMCSignals.size(); i++) { - if (!(mcDecision & (uint32_t(1) << i))) { + if (!(mcDecision & (static_cast(1) << i))) { continue; } if (recfidcut) { for (unsigned int j = 0; j < fTrackCuts.size(); j++) { - if (twoTrackFilter & (uint8_t(1) << j)) { + if (twoTrackFilter & (static_cast(1) << j)) { if (!fConfigFillLS) { fHistRecPair[j * fMCSignals.size() + i]->Fill(VarManager::fgValues[VarManager::kMass], VarManager::fgValues[VarManager::kPt]); if (fConfigQA) @@ -1874,7 +1874,7 @@ struct AnalysisSameEventPairing { } if (genfidcut) { for (unsigned int j = 0; j < fTrackCuts.size(); j++) { - if (twoTrackFilter & (uint8_t(1) << j)) { + if (twoTrackFilter & (static_cast(1) << j)) { if (!fConfigFillLS) { fHistRecPairMC[j * fMCSignals.size() + i]->Fill(mass, pairpt); } else { diff --git a/PWGEM/Dilepton/Tasks/matchingMFT.cxx b/PWGEM/Dilepton/Tasks/matchingMFT.cxx index 6d95d3abaa2..08a1ddb5bde 100644 --- a/PWGEM/Dilepton/Tasks/matchingMFT.cxx +++ b/PWGEM/Dilepton/Tasks/matchingMFT.cxx @@ -257,7 +257,7 @@ struct matchingMFT { } o2::mch::TrackExtrap::setField(); const double centerMFT[3] = {0, 0, -61.4}; - o2::field::MagneticField* field = static_cast(TGeoGlobalMagField::Instance()->GetField()); + o2::field::MagneticField* field = dynamic_cast(TGeoGlobalMagField::Instance()->GetField()); mBz = field->getBz(centerMFT); // Get field at centre of MFT LOGF(info, "Bz at center of MFT = %f kZG", mBz); diff --git a/PWGEM/Dilepton/Tasks/tableReaderBarrel.cxx b/PWGEM/Dilepton/Tasks/tableReaderBarrel.cxx index 46cfe0d6e90..95a6dd3bc6b 100644 --- a/PWGEM/Dilepton/Tasks/tableReaderBarrel.cxx +++ b/PWGEM/Dilepton/Tasks/tableReaderBarrel.cxx @@ -385,37 +385,37 @@ struct AnalysisTrackSelection { bool CheckSize() { - auto veceta = (std::vector)trackcuts.cfgEtaMax; - auto vecptmin = (std::vector)trackcuts.cfgPtMin; - auto vecptmax = (std::vector)trackcuts.cfgPtMax; - auto vecDCAxymax = (std::vector)trackcuts.cfgDCAxyMax; - auto vecDCAzmax = (std::vector)trackcuts.cfgDCAzMax; - auto vecIsSPDfirst = (std::vector)trackcuts.cfgIsSPDfirst; - auto vecIsSPDany = (std::vector)trackcuts.cfgIsSPDany; - auto vecIsITSibAny = (std::vector)trackcuts.cfgIsITSibAny; - auto vecITSchi2max = (std::vector)trackcuts.cfgITSchi2Max; - auto vecITSnclsmin = (std::vector)trackcuts.cfgITSnclsMin; - auto vecITSnclsmax = (std::vector)trackcuts.cfgITSnclsMax; - auto vecTPCchi2max = (std::vector)trackcuts.cfgTPCchi2Max; - auto vecTPCnclsmin = (std::vector)trackcuts.cfgTPCnclsMin; - auto vecTPCnclsmax = (std::vector)trackcuts.cfgTPCnclsMax; - auto vecTPCnclsCRmin = (std::vector)trackcuts.cfgTPCnclsCRMin; - auto vecTPCnclsCRmax = (std::vector)trackcuts.cfgTPCnclsCRMax; - auto vecIsDalitzLeg = (std::vector)trackcuts.cfgIsDalitzLeg; - auto vecPIDmode = (std::vector)pidcuts.cfgPIDmode; - auto vecrejbadtof = (std::vector)pidcuts.cfgRejBadTOF; - auto vecTPCnsigmaelmin = (std::vector)pidcuts.cfgTPCNSigmaElMin; - auto vecTPCnsigmaelmax = (std::vector)pidcuts.cfgTPCNSigmaElMax; - auto vecTPCnsigmapimin = (std::vector)pidcuts.cfgTPCNSigmaPiMin; - auto vecTPCnsigmapimax = (std::vector)pidcuts.cfgTPCNSigmaPiMax; - auto vecTPCnsigmaprmin = (std::vector)pidcuts.cfgTPCNSigmaPrMin; - auto vecTPCnsigmaprmax = (std::vector)pidcuts.cfgTPCNSigmaPrMax; - auto vecTPCnsigmakamin = (std::vector)pidcuts.cfgTPCNSigmaKaMin; - auto vecTPCnsigmakamax = (std::vector)pidcuts.cfgTPCNSigmaKaMax; - auto vecTPCnsigmamumin = (std::vector)pidcuts.cfgTPCNSigmaMuMin; - auto vecTPCnsigmamumax = (std::vector)pidcuts.cfgTPCNSigmaMuMax; - auto vecTOFnsigmaelmin = (std::vector)pidcuts.cfgTOFNSigmaElMin; - auto vecTOFnsigmaelmax = (std::vector)pidcuts.cfgTOFNSigmaElMax; + auto veceta = std::vector(trackcuts.cfgEtaMax); + auto vecptmin = std::vector(trackcuts.cfgPtMin); + auto vecptmax = std::vector(trackcuts.cfgPtMax); + auto vecDCAxymax = std::vector(trackcuts.cfgDCAxyMax); + auto vecDCAzmax = std::vector(trackcuts.cfgDCAzMax); + auto vecIsSPDfirst = std::vector(trackcuts.cfgIsSPDfirst); + auto vecIsSPDany = std::vector(trackcuts.cfgIsSPDany); + auto vecIsITSibAny = std::vector(trackcuts.cfgIsITSibAny); + auto vecITSchi2max = std::vector(trackcuts.cfgITSchi2Max); + auto vecITSnclsmin = std::vector(trackcuts.cfgITSnclsMin); + auto vecITSnclsmax = std::vector(trackcuts.cfgITSnclsMax); + auto vecTPCchi2max = std::vector(trackcuts.cfgTPCchi2Max); + auto vecTPCnclsmin = std::vector(trackcuts.cfgTPCnclsMin); + auto vecTPCnclsmax = std::vector(trackcuts.cfgTPCnclsMax); + auto vecTPCnclsCRmin = std::vector(trackcuts.cfgTPCnclsCRMin); + auto vecTPCnclsCRmax = std::vector(trackcuts.cfgTPCnclsCRMax); + auto vecIsDalitzLeg = std::vector(trackcuts.cfgIsDalitzLeg); + auto vecPIDmode = std::vector(pidcuts.cfgPIDmode); + auto vecrejbadtof = std::vector(pidcuts.cfgRejBadTOF); + auto vecTPCnsigmaelmin = std::vector(pidcuts.cfgTPCNSigmaElMin); + auto vecTPCnsigmaelmax = std::vector(pidcuts.cfgTPCNSigmaElMax); + auto vecTPCnsigmapimin = std::vector(pidcuts.cfgTPCNSigmaPiMin); + auto vecTPCnsigmapimax = std::vector(pidcuts.cfgTPCNSigmaPiMax); + auto vecTPCnsigmaprmin = std::vector(pidcuts.cfgTPCNSigmaPrMin); + auto vecTPCnsigmaprmax = std::vector(pidcuts.cfgTPCNSigmaPrMax); + auto vecTPCnsigmakamin = std::vector(pidcuts.cfgTPCNSigmaKaMin); + auto vecTPCnsigmakamax = std::vector(pidcuts.cfgTPCNSigmaKaMax); + auto vecTPCnsigmamumin = std::vector(pidcuts.cfgTPCNSigmaMuMin); + auto vecTPCnsigmamumax = std::vector(pidcuts.cfgTPCNSigmaMuMax); + auto vecTOFnsigmaelmin = std::vector(pidcuts.cfgTOFNSigmaElMin); + auto vecTOFnsigmaelmax = std::vector(pidcuts.cfgTOFNSigmaElMax); if (veceta.size() != fConfigNbTrackCut) return false; @@ -484,20 +484,20 @@ struct AnalysisTrackSelection { AnalysisCompositeCut* GetPIDCut(unsigned int i) { - auto vecPIDmode = (std::vector)pidcuts.cfgPIDmode; - auto vecrejbadtof = (std::vector)pidcuts.cfgRejBadTOF; - auto vecTPCnsigmaelmin = (std::vector)pidcuts.cfgTPCNSigmaElMin; - auto vecTPCnsigmaelmax = (std::vector)pidcuts.cfgTPCNSigmaElMax; - auto vecTPCnsigmapimin = (std::vector)pidcuts.cfgTPCNSigmaPiMin; - auto vecTPCnsigmapimax = (std::vector)pidcuts.cfgTPCNSigmaPiMax; - auto vecTPCnsigmaprmin = (std::vector)pidcuts.cfgTPCNSigmaPrMin; - auto vecTPCnsigmaprmax = (std::vector)pidcuts.cfgTPCNSigmaPrMax; - auto vecTPCnsigmakamin = (std::vector)pidcuts.cfgTPCNSigmaKaMin; - auto vecTPCnsigmakamax = (std::vector)pidcuts.cfgTPCNSigmaKaMax; - auto vecTPCnsigmamumin = (std::vector)pidcuts.cfgTPCNSigmaMuMin; - auto vecTPCnsigmamumax = (std::vector)pidcuts.cfgTPCNSigmaMuMax; - auto vecTOFnsigmaelmin = (std::vector)pidcuts.cfgTOFNSigmaElMin; - auto vecTOFnsigmaelmax = (std::vector)pidcuts.cfgTOFNSigmaElMax; + auto vecPIDmode = std::vector(pidcuts.cfgPIDmode); + auto vecrejbadtof = std::vector(pidcuts.cfgRejBadTOF); + auto vecTPCnsigmaelmin = std::vector(pidcuts.cfgTPCNSigmaElMin); + auto vecTPCnsigmaelmax = std::vector(pidcuts.cfgTPCNSigmaElMax); + auto vecTPCnsigmapimin = std::vector(pidcuts.cfgTPCNSigmaPiMin); + auto vecTPCnsigmapimax = std::vector(pidcuts.cfgTPCNSigmaPiMax); + auto vecTPCnsigmaprmin = std::vector(pidcuts.cfgTPCNSigmaPrMin); + auto vecTPCnsigmaprmax = std::vector(pidcuts.cfgTPCNSigmaPrMax); + auto vecTPCnsigmakamin = std::vector(pidcuts.cfgTPCNSigmaKaMin); + auto vecTPCnsigmakamax = std::vector(pidcuts.cfgTPCNSigmaKaMax); + auto vecTPCnsigmamumin = std::vector(pidcuts.cfgTPCNSigmaMuMin); + auto vecTPCnsigmamumax = std::vector(pidcuts.cfgTPCNSigmaMuMax); + auto vecTOFnsigmaelmin = std::vector(pidcuts.cfgTOFNSigmaElMin); + auto vecTOFnsigmaelmax = std::vector(pidcuts.cfgTOFNSigmaElMax); AnalysisCompositeCut* cut_tpc_nSigma = new AnalysisCompositeCut(Form("pid_TPCnSigma_%d", i), Form("pid_TPCnSigma_%d", i), kTRUE); AnalysisCut* cuttpc = new AnalysisCut(Form("pidcuttpc_%d", i), Form("pidcuttpc_%d", i)); @@ -547,28 +547,28 @@ struct AnalysisTrackSelection { AnalysisCut* GetTrackCut(unsigned int i) { - auto veceta = (std::vector)trackcuts.cfgEtaMax; - auto vecptmin = (std::vector)trackcuts.cfgPtMin; - auto vecptmax = (std::vector)trackcuts.cfgPtMax; + auto veceta = std::vector(trackcuts.cfgEtaMax); + auto vecptmin = std::vector(trackcuts.cfgPtMin); + auto vecptmax = std::vector(trackcuts.cfgPtMax); - auto vecDCAxymax = (std::vector)trackcuts.cfgDCAxyMax; - auto vecDCAzmax = (std::vector)trackcuts.cfgDCAzMax; + auto vecDCAxymax = std::vector(trackcuts.cfgDCAxyMax); + auto vecDCAzmax = std::vector(trackcuts.cfgDCAzMax); - auto vecIsSPDfirst = (std::vector)trackcuts.cfgIsSPDfirst; - auto vecIsSPDany = (std::vector)trackcuts.cfgIsSPDany; - auto vecIsITSibAny = (std::vector)trackcuts.cfgIsITSibAny; + auto vecIsSPDfirst = std::vector(trackcuts.cfgIsSPDfirst); + auto vecIsSPDany = std::vector(trackcuts.cfgIsSPDany); + auto vecIsITSibAny = std::vector(trackcuts.cfgIsITSibAny); - auto vecITSchi2max = (std::vector)trackcuts.cfgITSchi2Max; - auto vecITSnclsmin = (std::vector)trackcuts.cfgITSnclsMin; - auto vecITSnclsmax = (std::vector)trackcuts.cfgITSnclsMax; + auto vecITSchi2max = std::vector(trackcuts.cfgITSchi2Max); + auto vecITSnclsmin = std::vector(trackcuts.cfgITSnclsMin); + auto vecITSnclsmax = std::vector(trackcuts.cfgITSnclsMax); - auto vecTPCchi2max = (std::vector)trackcuts.cfgTPCchi2Max; - auto vecTPCnclsmin = (std::vector)trackcuts.cfgTPCnclsMin; - auto vecTPCnclsmax = (std::vector)trackcuts.cfgTPCnclsMax; - auto vecTPCnclsCRmin = (std::vector)trackcuts.cfgTPCnclsCRMin; - auto vecTPCnclsCRmax = (std::vector)trackcuts.cfgTPCnclsCRMax; + auto vecTPCchi2max = std::vector(trackcuts.cfgTPCchi2Max); + auto vecTPCnclsmin = std::vector(trackcuts.cfgTPCnclsMin); + auto vecTPCnclsmax = std::vector(trackcuts.cfgTPCnclsMax); + auto vecTPCnclsCRmin = std::vector(trackcuts.cfgTPCnclsCRMin); + auto vecTPCnclsCRmax = std::vector(trackcuts.cfgTPCnclsCRMax); - auto vecIsDalitzLeg = (std::vector)trackcuts.cfgIsDalitzLeg; + auto vecIsDalitzLeg = std::vector(trackcuts.cfgIsDalitzLeg); AnalysisCut* cut = new AnalysisCut(Form("tracksel%d", i), Form("tracksel%d", i)); cut->AddCut(VarManager::kPt, vecptmin.at(i), vecptmax.at(i)); @@ -986,11 +986,11 @@ struct AnalysisSameEventPairing { bool CheckSize() { - auto vecRej = (std::vector)paircuts.cfgRej; - auto vecmassmin = (std::vector)paircuts.cfgMassMin; - auto vecmassmax = (std::vector)paircuts.cfgMassMax; - auto vecopamax = (std::vector)paircuts.cfgOpAngMax; - auto vecphivmin = (std::vector)paircuts.cfgPhiVMin; + auto vecRej = std::vector(paircuts.cfgRej); + auto vecmassmin = std::vector(paircuts.cfgMassMin); + auto vecmassmax = std::vector(paircuts.cfgMassMax); + auto vecopamax = std::vector(paircuts.cfgOpAngMax); + auto vecphivmin = std::vector(paircuts.cfgPhiVMin); if (vecRej.size() != fConfigNbPairCut) return false; @@ -1007,11 +1007,11 @@ struct AnalysisSameEventPairing { AnalysisCompositeCut* GetPairCut(unsigned int i) { - auto vecRej = (std::vector)paircuts.cfgRej; - auto vecmassmin = (std::vector)paircuts.cfgMassMin; - auto vecmassmax = (std::vector)paircuts.cfgMassMax; - auto vecopamax = (std::vector)paircuts.cfgOpAngMax; - auto vecphivmin = (std::vector)paircuts.cfgPhiVMin; + auto vecRej = std::vector(paircuts.cfgRej); + auto vecmassmin = std::vector(paircuts.cfgMassMin); + auto vecmassmax = std::vector(paircuts.cfgMassMax); + auto vecopamax = std::vector(paircuts.cfgOpAngMax); + auto vecphivmin = std::vector(paircuts.cfgPhiVMin); if (vecRej.at(i)) { AnalysisCompositeCut* cut_pair = new AnalysisCompositeCut(Form("cut_pair_%d", i), Form("cut_pair_%d", i), kFALSE); @@ -1099,7 +1099,7 @@ struct AnalysisSameEventPairing { } // end loop (track cuts) } // end if (track cuts) - VarManager::SetCollisionSystem((TString)fCollisionSystem, fCenterMassEnergy); // set collision system and center of mass energy + VarManager::SetCollisionSystem(TString(fCollisionSystem), fCenterMassEnergy); // set collision system and center of mass energy DefineHistograms(fHistMan, histNames.Data(), fConfigAddSEPHistogram); // define all histograms VarManager::SetUseVars(fHistMan->GetUsedVars()); // provide the list of required variables so that VarManager knows what to fill diff --git a/PWGEM/Dilepton/Tasks/taggingHFE.cxx b/PWGEM/Dilepton/Tasks/taggingHFE.cxx index 27bf830e4f5..1d051f43a58 100644 --- a/PWGEM/Dilepton/Tasks/taggingHFE.cxx +++ b/PWGEM/Dilepton/Tasks/taggingHFE.cxx @@ -2684,7 +2684,7 @@ struct taggingHFE { // std::vector> stored_trackIds; //! type of V0. 0: built solely for cascades (does not pass standard V0 cut), 1: standard 2, 3: photon-like with TPC-only use. Regular analysis should always use type 1. - Filter v0Filter = o2::aod::v0data::v0Type == uint8_t(1) && o2::aod::v0data::v0cosPA > v0Cut.cfg_min_cospa&& o2::aod::v0data::dcaV0daughters v0Cut.cfg_min_dcaxy&& nabs(o2::aod::v0data::dcanegtopv) > v0Cut.cfg_min_dcaxy; + Filter v0Filter = o2::aod::v0data::v0Type == static_cast(1) && o2::aod::v0data::v0cosPA > v0Cut.cfg_min_cospa&& o2::aod::v0data::dcaV0daughters v0Cut.cfg_min_dcaxy&& nabs(o2::aod::v0data::dcanegtopv) > v0Cut.cfg_min_dcaxy; using filteredV0s = soa::Filtered; Filter cascadeFilter = nabs(o2::aod::cascdata::dcanegtopv) > cascadeCut.cfg_min_dcaxy_v0leg&& nabs(o2::aod::cascdata::dcanegtopv) > cascadeCut.cfg_min_dcaxy_v0leg&& nabs(o2::aod::cascdata::dcabachtopv) > cascadeCut.cfg_min_dcaxy_bachelor; diff --git a/PWGEM/Dilepton/Utils/ElectronModule.h b/PWGEM/Dilepton/Utils/ElectronModule.h index 9f9d4a9131e..83453f055e6 100644 --- a/PWGEM/Dilepton/Utils/ElectronModule.h +++ b/PWGEM/Dilepton/Utils/ElectronModule.h @@ -1600,12 +1600,12 @@ class ElectronModule registry.fill(HIST("Prefilter/Pair/hMvsPhiV"), phiv, v12.M()); if (v12.M() < fElectronPFCut.slope * phiv + fElectronPFCut.intercept) { - pfb |= (uint8_t(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); + pfb |= (static_cast(1) << static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPC)); } for (int i = 0; i < static_cast(max_mee_vec.size()); i++) { if (v12.M() < max_mee_vec.at(i)) { - pfb |= (uint8_t(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); + pfb |= (static_cast(1) << (static_cast(o2::aod::pwgem::dilepton::utils::pairutil::DileptonPrefilterBit::kElFromPi0_20MeV) + i)); } } // end of mee loop diff --git a/PWGEM/Dilepton/Utils/MCUtilities.h b/PWGEM/Dilepton/Utils/MCUtilities.h index 51753c2da09..56203d1cebd 100644 --- a/PWGEM/Dilepton/Utils/MCUtilities.h +++ b/PWGEM/Dilepton/Utils/MCUtilities.h @@ -914,8 +914,8 @@ int searchMothers(T& p, U& mcParticles, int pdg, bool equal) } quark_id = i; //} else if ((static_cast(std::abs(mpdg) / 100) == pdg || static_cast(std::abs(mpdg) / 1000) == pdg) && mpdg * p.pdgCode() > 0) { // check for other mothers with flavour content - } else if ((static_cast(std::abs(mpdg) / 100) == pdg || static_cast(std::abs(mpdg) / 1000) == pdg)) { // check for other mothers with flavour content to allow for beauty and charm - if (quark_id > -1 || next_mother_id > -1) { // we already found a possible candidate in the list of mothers, so now we have (at least) two + } else if (((std::abs(mpdg) / 100) == pdg || (std::abs(mpdg) / 1000) == pdg)) { // check for other mothers with flavour content to allow for beauty and charm + if (quark_id > -1 || next_mother_id > -1) { // we already found a possible candidate in the list of mothers, so now we have (at least) two // LOG(warning) << "Flavour tracking is ambiguous. Stopping here."; return -1; } diff --git a/PWGEM/Dilepton/Utils/MomentumSmearer.h b/PWGEM/Dilepton/Utils/MomentumSmearer.h index 878f3403109..cfbebd00427 100644 --- a/PWGEM/Dilepton/Utils/MomentumSmearer.h +++ b/PWGEM/Dilepton/Utils/MomentumSmearer.h @@ -110,22 +110,22 @@ class MomentumSmearer void getResolutionHistos(TList* list) { - fResoPt = reinterpret_cast(list->FindObject(fResPtHistName)); + fResoPt = dynamic_cast(list->FindObject(fResPtHistName)); if (!fResoPt) { LOGP(fatal, "Could not open {} from file {}", fResPtHistName.Data(), fResFileName.Data()); } - fResoEta = reinterpret_cast(list->FindObject(fResEtaHistName)); + fResoEta = dynamic_cast(list->FindObject(fResEtaHistName)); if (!fResoEta) { LOGP(fatal, "Could not open {} from file {}", fResEtaHistName.Data(), fResFileName.Data()); } - fResoPhi_Pos = reinterpret_cast(list->FindObject(fResPhiPosHistName)); + fResoPhi_Pos = dynamic_cast(list->FindObject(fResPhiPosHistName)); if (!fResoPhi_Pos) { LOGP(fatal, "Could not open {} from file {}", fResPhiPosHistName.Data(), fResFileName.Data()); } - fResoPhi_Neg = reinterpret_cast(list->FindObject(fResPhiNegHistName)); + fResoPhi_Neg = dynamic_cast(list->FindObject(fResPhiNegHistName)); if (!fResoPhi_Neg) { LOGP(fatal, "Could not open {} from file {}", fResPhiNegHistName.Data(), fResFileName.Data()); } @@ -137,8 +137,8 @@ class MomentumSmearer int nBinsPt = axisPt->GetNbins(); fVecReso.resize(nBinsPt); for (int i = 0; i < nBinsPt; i++) { - auto h1 = reinterpret_cast(fReso->ProjectionY(Form("h1reso%s_pt%d", suffix, i), i + 1, i + 1)); - h1->Scale(1.f, "width"); // convert ntrack to probability density + auto h1 = dynamic_cast(fReso->ProjectionY(Form("h1reso%s_pt%d", suffix, i), i + 1, i + 1)); // FIXME: invalid cast + h1->Scale(1.f, "width"); // convert ntrack to probability density fVecReso[i] = h1; } } @@ -156,7 +156,7 @@ class MomentumSmearer fVecResoND.resize(fNCenBins, std::vector>>>(fNPtBins, std::vector>>(fNEtaBins, std::vector>(fNPhiBins, std::vector(fNChBins))))); // fVecResoND.resize(fNPtBins, std::vector>>(fNEtaBins, std::vector>(fNPhiBins, std::vector(fNChBins)))); - // auto h3 = reinterpret_cast(hs_reso->Projection(4, 5, 6)); + // auto h3 = reinterpret_cast(hs_reso->Projection(4, 5, 6)); // FIXME: invalid cast // h3->SetName(Form("h3reso_pt%d_eta%d_phi%d_ch%d", 0, 0, 0, 0)); // fVecResoND[0][0][0][0] = h3; @@ -173,7 +173,7 @@ class MomentumSmearer continue; } hs_reso->GetAxis(4)->SetRange(ich + 1, ich + 1); - auto h3 = reinterpret_cast(hs_reso->Projection(5, 6, 7)); + auto h3 = hs_reso->Projection(5, 6, 7); h3->SetName(Form("h3reso_cen%d_pt%d_eta%d_phi%d_ch%d", icen, ipt, ieta, iphi, ich)); h3->Scale(1.f, "width"); // convert ntrack to probability density fVecResoND[icen][ipt][ieta][iphi][ich] = h3; @@ -214,7 +214,7 @@ class MomentumSmearer return; } if (fFile->GetListOfKeys()->Contains("ccdb_object")) { - listRes = reinterpret_cast(fFile->Get("ccdb_object")); + listRes = dynamic_cast(fFile->Get("ccdb_object")); } else { for (TObject* keyAsObj : *(fFile->GetListOfKeys())) { auto key = dynamic_cast(keyAsObj); @@ -231,7 +231,7 @@ class MomentumSmearer if (fDoNDSmearing) { if (fResType != 0) { - fResoND = reinterpret_cast(listRes->FindObject(fResNDHistName)); + fResoND = dynamic_cast(listRes->FindObject(fResNDHistName)); if (!fResoND) { LOGP(fatal, "Could not open {} from file {}", fResNDHistName.Data(), fResFileName.Data()); } @@ -239,22 +239,22 @@ class MomentumSmearer } } else { if (fResType != 0) { - fResoPt = reinterpret_cast(listRes->FindObject(fResPtHistName)); + fResoPt = dynamic_cast(listRes->FindObject(fResPtHistName)); if (!fResoPt) { LOGP(fatal, "Could not open {} from file {}", fResPtHistName.Data(), fResFileName.Data()); } - fResoEta = reinterpret_cast(listRes->FindObject(fResEtaHistName)); + fResoEta = dynamic_cast(listRes->FindObject(fResEtaHistName)); if (!fResoEta) { LOGP(fatal, "Could not open {} from file {}", fResEtaHistName.Data(), fResFileName.Data()); } - fResoPhi_Pos = reinterpret_cast(listRes->FindObject(fResPhiPosHistName)); + fResoPhi_Pos = dynamic_cast(listRes->FindObject(fResPhiPosHistName)); if (!fResoPhi_Pos) { LOGP(fatal, "Could not open {} from file {}", fResPhiPosHistName.Data(), fResFileName.Data()); } - fResoPhi_Neg = reinterpret_cast(listRes->FindObject(fResPhiNegHistName)); + fResoPhi_Neg = dynamic_cast(listRes->FindObject(fResPhiNegHistName)); if (!fResoPhi_Neg) { LOGP(fatal, "Could not open {} from file {}", fResPhiNegHistName.Data(), fResFileName.Data()); } @@ -289,7 +289,7 @@ class MomentumSmearer return; } if (fFile->GetListOfKeys()->Contains("ccdb_object")) { - listEff = reinterpret_cast(fFile->Get("ccdb_object")); + listEff = dynamic_cast(fFile->Get("ccdb_object")); } else { for (TObject* keyAsObj : *(fFile->GetListOfKeys())) { auto key = dynamic_cast(keyAsObj); @@ -302,7 +302,7 @@ class MomentumSmearer } } if (fEffType != 0) { - fEff = reinterpret_cast(listEff->FindObject(fEffHistName)); + fEff = listEff->FindObject(fEffHistName); if (!fEff) { LOGP(fatal, "Could not open {} from file {}", fEffHistName.Data(), fEffFileName.Data()); } @@ -345,7 +345,7 @@ class MomentumSmearer return; } if (fFile->GetListOfKeys()->Contains("ccdb_object")) { - listDCA = reinterpret_cast(fFile->Get("ccdb_object")); + listDCA = dynamic_cast(fFile->Get("ccdb_object")); } else { for (TObject* keyAsObj : *(fFile->GetListOfKeys())) { auto key = dynamic_cast(keyAsObj); @@ -358,7 +358,7 @@ class MomentumSmearer } } if (fDCAType != 0) { - fDCA = reinterpret_cast(listDCA->FindObject(fDCAHistName)); + fDCA = dynamic_cast(listDCA->FindObject(fDCAHistName)); if (!fDCA) { LOGP(fatal, "Could not open {} from file {}", fDCAHistName.Data(), fDCAFileName.Data()); } @@ -481,7 +481,7 @@ class MomentumSmearer } if (fEffType == 1) { - TH1F* hist = reinterpret_cast(fEff); + TH1F* hist = dynamic_cast(fEff); int ptbin = hist->GetXaxis()->FindBin(pt); int ptbin_max = hist->GetXaxis()->GetNbins(); // make sure that no underflow or overflow bins are used @@ -493,7 +493,7 @@ class MomentumSmearer } if (fEffType == 2) { - TH2F* hist = reinterpret_cast(fEff); + TH2F* hist = dynamic_cast(fEff); int ptbin = hist->GetXaxis()->FindBin(pt); int ptbin_max = hist->GetXaxis()->GetNbins(); int etabin = hist->GetYaxis()->FindBin(eta); @@ -511,7 +511,7 @@ class MomentumSmearer } if (fEffType == 3) { - TH3F* hist = reinterpret_cast(fEff); + TH3F* hist = dynamic_cast(fEff); int ptbin = hist->GetXaxis()->FindBin(pt); int ptbin_max = hist->GetXaxis()->GetNbins(); int etabin = hist->GetYaxis()->FindBin(eta); diff --git a/PWGEM/PhotonMeson/Core/HistogramsLibrary.cxx b/PWGEM/PhotonMeson/Core/HistogramsLibrary.cxx index 64e5c5cbd36..44e770bce73 100644 --- a/PWGEM/PhotonMeson/Core/HistogramsLibrary.cxx +++ b/PWGEM/PhotonMeson/Core/HistogramsLibrary.cxx @@ -33,16 +33,16 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const { if (TString(histClass) == "Event") { list->Add(new TH1F("hCollisionCounter", "hCollisionCounter", 10, 0.5f, 10.5f)); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(1, "all"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(2, "No TF border"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(3, "No ITS ROF border"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(4, "No Same Bunch Pileup"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(5, "Is Vertex ITSTPC"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(6, "Is Good Zvtx FT0vsPV"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(7, "FT0AND"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(8, "sel8"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(9, "|Z_{vtx}| < 10 cm"); - reinterpret_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(10, "accepted"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(1, "all"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(2, "No TF border"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(3, "No ITS ROF border"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(4, "No Same Bunch Pileup"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(5, "Is Vertex ITSTPC"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(6, "Is Good Zvtx FT0vsPV"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(7, "FT0AND"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(8, "sel8"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(9, "|Z_{vtx}| < 10 cm"); + dynamic_cast(list->FindObject("hCollisionCounter"))->GetXaxis()->SetBinLabel(10, "accepted"); list->Add(new TH1F("hZvtx", "vertex z; Z_{vtx} (cm)", 100, -50, +50)); list->Add(new TH1F("hMultNTracksPV", "hMultNTracksPV; N_{track} to PV", 6001, -0.5, 6000.5)); @@ -254,16 +254,16 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const TH2F* hMvsPhiV_uls_same = new TH2F("hMvsPhiV_uls_same", "m_{ee} vs. #varphi_{V};#varphi_{V} (rad.);m_{ee} (GeV/c^{2})", 72, 0, M_PI, 100, 0.0f, 0.1f); hMvsPhiV_uls_same->Sumw2(); list->Add(hMvsPhiV_uls_same); - list->Add(reinterpret_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lspp_same"))); - list->Add(reinterpret_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lsmm_same"))); + list->Add(dynamic_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lspp_same"))); + list->Add(dynamic_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lsmm_same"))); if (TString(subGroup).Contains("mix")) { - list->Add(reinterpret_cast(hs_dilepton_uls_same->Clone("hs_dilepton_uls_mix"))); - list->Add(reinterpret_cast(hs_dilepton_lspp_same->Clone("hs_dilepton_lspp_mix"))); - list->Add(reinterpret_cast(hs_dilepton_lsmm_same->Clone("hs_dilepton_lsmm_mix"))); - list->Add(reinterpret_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_uls_mix"))); - list->Add(reinterpret_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lspp_mix"))); - list->Add(reinterpret_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lsmm_mix"))); + list->Add(dynamic_cast(hs_dilepton_uls_same->Clone("hs_dilepton_uls_mix"))); + list->Add(dynamic_cast(hs_dilepton_lspp_same->Clone("hs_dilepton_lspp_mix"))); + list->Add(dynamic_cast(hs_dilepton_lsmm_same->Clone("hs_dilepton_lsmm_mix"))); + list->Add(dynamic_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_uls_mix"))); + list->Add(dynamic_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lspp_mix"))); + list->Add(dynamic_cast(hMvsPhiV_uls_same->Clone("hMvsPhiV_lsmm_mix"))); } if (TString(subGroup).Contains("dca")) { @@ -280,14 +280,14 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const hs_dilepton_uls_dca_same->Sumw2(); list->Add(hs_dilepton_uls_dca_same); - hs_dilepton_lspp_dca_same = reinterpret_cast(hs_dilepton_uls_dca_same->Clone("hs_dilepton_lspp_dca_same")); - hs_dilepton_lsmm_dca_same = reinterpret_cast(hs_dilepton_uls_dca_same->Clone("hs_dilepton_lsmm_dca_same")); + hs_dilepton_lspp_dca_same = dynamic_cast(hs_dilepton_uls_dca_same->Clone("hs_dilepton_lspp_dca_same")); + hs_dilepton_lsmm_dca_same = dynamic_cast(hs_dilepton_uls_dca_same->Clone("hs_dilepton_lsmm_dca_same")); list->Add(hs_dilepton_lspp_dca_same); list->Add(hs_dilepton_lsmm_dca_same); - THnSparseF* hs_dilepton_uls_dca_mix = reinterpret_cast(hs_dilepton_uls_dca_same->Clone("hs_dilepton_uls_dca_mix")); - THnSparseF* hs_dilepton_lspp_dca_mix = reinterpret_cast(hs_dilepton_lspp_dca_same->Clone("hs_dilepton_lspp_dca_mix")); - THnSparseF* hs_dilepton_lsmm_dca_mix = reinterpret_cast(hs_dilepton_lsmm_dca_same->Clone("hs_dilepton_lsmm_dca_mix")); + THnSparseF* hs_dilepton_uls_dca_mix = dynamic_cast(hs_dilepton_uls_dca_same->Clone("hs_dilepton_uls_dca_mix")); + THnSparseF* hs_dilepton_lspp_dca_mix = dynamic_cast(hs_dilepton_lspp_dca_same->Clone("hs_dilepton_lspp_dca_mix")); + THnSparseF* hs_dilepton_lsmm_dca_mix = dynamic_cast(hs_dilepton_lsmm_dca_same->Clone("hs_dilepton_lsmm_dca_mix")); list->Add(hs_dilepton_uls_dca_mix); list->Add(hs_dilepton_lspp_dca_mix); list->Add(hs_dilepton_lsmm_dca_mix); @@ -336,9 +336,9 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const list->Add(hs_dilepton_lsmm_same); if (TString(subGroup).Contains("mix")) { - list->Add(reinterpret_cast(hs_dilepton_uls_same->Clone("hs_dilepton_uls_mix"))); - list->Add(reinterpret_cast(hs_dilepton_lspp_same->Clone("hs_dilepton_lspp_mix"))); - list->Add(reinterpret_cast(hs_dilepton_lsmm_same->Clone("hs_dilepton_lsmm_mix"))); + list->Add(dynamic_cast(hs_dilepton_uls_same->Clone("hs_dilepton_uls_mix"))); + list->Add(dynamic_cast(hs_dilepton_lspp_same->Clone("hs_dilepton_lspp_mix"))); + list->Add(dynamic_cast(hs_dilepton_lsmm_same->Clone("hs_dilepton_lsmm_mix"))); } } else { LOGF(info, "EE or MuMu are supported."); @@ -427,14 +427,14 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const list->Add(new TH1F("hDistToBC", "distance to bad channel;#it{d};#it{N}_{cluster}", 100, 0, 1500)); } list->Add(new TH2F("hClusterQualityCuts", "Energy at which clusters are removed by a given cut;;#it{E} (GeV)", 8, -0.5, 7.5, 500, 0, 50)); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(1, "In"); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(2, "Energy"); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(3, "NCell"); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(4, "M02"); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(5, "Timing"); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(6, "Track matching"); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(7, "Exotic"); - reinterpret_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(8, "Out"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(1, "In"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(2, "Energy"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(3, "NCell"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(4, "M02"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(5, "Timing"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(6, "Track matching"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(7, "Exotic"); + dynamic_cast(list->FindObject("hClusterQualityCuts"))->GetXaxis()->SetBinLabel(8, "Out"); } } @@ -444,18 +444,18 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const list->Add(new TH2F("hPt_SPQ2FT0A", "p_{T,#gamma} vs. SP;p_{T,#gamma} (GeV/c);u_{2}^{#gamma} #upoint Q_{2}^{FT0A}", 400, 0.0f, 20, 40, -10, +10)); list->Add(new TH2F("hPt_SPQ2FT0C", "p_{T,#gamma} vs. SP;p_{T,#gamma} (GeV/c);u_{2}^{#gamma} #upoint Q_{2}^{FT0C}", 400, 0.0f, 20, 40, -10, +10)); list->Add(new TH2F("hPt_SPQ2FV0A", "p_{T,#gamma} vs. SP;p_{T,#gamma} (GeV/c);u_{2}^{#gamma} #upoint Q_{2}^{FV0A}", 400, 0.0f, 20, 40, -10, +10)); - reinterpret_cast(list->FindObject("hPt_SPQ2FT0M"))->Sumw2(); - reinterpret_cast(list->FindObject("hPt_SPQ2FT0A"))->Sumw2(); - reinterpret_cast(list->FindObject("hPt_SPQ2FT0C"))->Sumw2(); - reinterpret_cast(list->FindObject("hPt_SPQ2FV0A"))->Sumw2(); + dynamic_cast(list->FindObject("hPt_SPQ2FT0M"))->Sumw2(); + dynamic_cast(list->FindObject("hPt_SPQ2FT0A"))->Sumw2(); + dynamic_cast(list->FindObject("hPt_SPQ2FT0C"))->Sumw2(); + dynamic_cast(list->FindObject("hPt_SPQ2FV0A"))->Sumw2(); } else { list->Add(new TH1F("hPt", "pT of photon candidates;p_{T,#gamma} (GeV/c)", 400, 0.0f, 20)); - reinterpret_cast(list->FindObject("hPt"))->Sumw2(); + dynamic_cast(list->FindObject("hPt"))->Sumw2(); } list->Add(new TH1F("hY", "rapidity of photon candidates;y", 40, -2.0f, 2.0f)); list->Add(new TH1F("hPhi", "azimuthal angle of photon candidates;#varphi (rad.)", 180, 0, 2 * M_PI)); - reinterpret_cast(list->FindObject("hY"))->Sumw2(); - reinterpret_cast(list->FindObject("hPhi"))->Sumw2(); + dynamic_cast(list->FindObject("hY"))->Sumw2(); + dynamic_cast(list->FindObject("hPhi"))->Sumw2(); if (TString(subGroup) == "mc") { list->Add(new TH1F("hPt_Photon_Primary", "pT;p_{T} (GeV/c)", 2000, 0.0f, 20)); // for MC efficiency @@ -500,18 +500,18 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const hs_same_tmp->Sumw2(); std::string_view sp_names[4] = {"FT0M", "FT0A", "FT0C", "FV0A"}; for (int i = 0; i < 4; i++) { - THnSparseF* hs = reinterpret_cast(hs_same_tmp->Clone(Form("hs_Same_SPQ2%s", sp_names[i].data()))); + THnSparseF* hs = dynamic_cast(hs_same_tmp->Clone(Form("hs_Same_SPQ2%s", sp_names[i].data()))); hs->SetTitle(Form("diphoton info %s", sp_names[i].data())); hs->GetAxis(2)->SetTitle(Form("u_{2}^{#gamma#gamma} #upoint Q_{2}^{%s}", sp_names[i].data())); list->Add(hs); } } else { list->Add(new TH2F("hMggPt_Same", "m_{#gamma#gamma} vs. p_{T};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma#gamma} (GeV/c)", nmgg - 1, mgg, npTgg - 1, pTgg)); - reinterpret_cast(list->FindObject("hMggPt_Same"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Same"))->Sumw2(); } list->Add(new TH2F("hMggPt_Mixed", "m_{#gamma#gamma} vs. p_{T};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma#gamma} (GeV/c)", nmgg - 1, mgg, npTgg - 1, pTgg)); - reinterpret_cast(list->FindObject("hMggPt_Mixed"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Mixed"))->Sumw2(); if (TString(subGroup).Contains("PCMPHOS") || TString(subGroup).Contains("PCMEMC")) { list->Add(new TH2F("hdEtadPhi", "#Delta#eta vs. #Delta#varphi;#Delta#varphi (rad.);#Delta#eta", 200, -1, +1, 200, -1, +1)); @@ -522,16 +522,16 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const if (TString(subGroup).Contains("EMCEMC")) { list->Add(new TH2F("hMggPt_Same_RotatedBkg", "m_{#gamma#gamma} vs. p_{T};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma#gamma} (GeV/c)", nmgg - 1, mgg, npTgg - 1, pTgg)); - reinterpret_cast(list->FindObject("hMggPt_Same_RotatedBkg"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Same_RotatedBkg"))->Sumw2(); } } if (TString(histClass) == "gammagamma_mass_pt_mc") { list->Add(new TH2F("hMggPt_Pi0_Primary", "m_{#gamma#gamma} vs. p_{T};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma#gamma} (GeV/c)", nmgg - 1, mgg, npTgg - 1, pTgg)); list->Add(new TH2F("hMggPt_Pi0_FromWD", "m_{#gamma#gamma} vs. p_{T};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma#gamma} (GeV/c)", nmgg - 1, mgg, npTgg - 1, pTgg)); list->Add(new TH2F("hMggPt_Eta_Primary", "m_{#gamma#gamma} vs. p_{T};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma#gamma} (GeV/c)", nmgg - 1, mgg, npTgg - 1, pTgg)); - reinterpret_cast(list->FindObject("hMggPt_Pi0_Primary"))->Sumw2(); - reinterpret_cast(list->FindObject("hMggPt_Pi0_FromWD"))->Sumw2(); - reinterpret_cast(list->FindObject("hMggPt_Eta_Primary"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Pi0_Primary"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Pi0_FromWD"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Eta_Primary"))->Sumw2(); } if (TString(histClass) == "material_budget_study") { @@ -640,36 +640,36 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const if (TString(histClass) == "tagging_pi0") { list->Add(new TH2F("hMggPt_Same", "m_{ee#gamma} vs. p_{T,#gamma};m_{ee#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c)", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); list->Add(new TH2F("hMggPt_Mixed", "m_{ee#gamma} vs. p_{T,#gamma};m_{ee#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c)", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); - reinterpret_cast(list->FindObject("hMggPt_Same"))->Sumw2(); - reinterpret_cast(list->FindObject("hMggPt_Mixed"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Same"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Mixed"))->Sumw2(); } if (TString(histClass) == "tagging_pi0_mc") { if (TString(subGroup) == "pcm") { list->Add(new TH1F("hPt_v0photon_Pi0_Primary", "reconstcuted v0 photon from primary #pi^{0};p_{T,#gamma} (GeV/c);N_{#gamma}^{#pi^{0}}", npTgg10 - 1, pTgg10)); // denominator for conditional probability - reinterpret_cast(list->FindObject("hPt_v0photon_Pi0_Primary"))->Sumw2(); + dynamic_cast(list->FindObject("hPt_v0photon_Pi0_Primary"))->Sumw2(); list->Add(new TH1F("hPt_v0photon_Pi0_FromWD", "reconstcuted v0 photon from #pi^{0} from WD;p_{T,#gamma} (GeV/c);N_{#gamma}^{#pi^{0}}", npTgg10 - 1, pTgg10)); // denominator for conditional probability - reinterpret_cast(list->FindObject("hPt_v0photon_Pi0_FromWD"))->Sumw2(); + dynamic_cast(list->FindObject("hPt_v0photon_Pi0_FromWD"))->Sumw2(); list->Add(new TH1F("hPt_v0photon_Pi0_hs", "reconstcuted v0 photon from #pi^{0} from hadronic shower in materials;p_{T,#gamma} (GeV/c);N_{#gamma}^{#pi^{0}}", npTgg10 - 1, pTgg10)); // denominator for conditional probability - reinterpret_cast(list->FindObject("hPt_v0photon_Pi0_hs"))->Sumw2(); + dynamic_cast(list->FindObject("hPt_v0photon_Pi0_hs"))->Sumw2(); } else if (TString(subGroup) == "pair") { list->Add(new TH2F("hMggPt_Pi0_Primary", "reconstructed m_{ee#gamma} vs. p_{T,#gamma} from primary #pi^{0};m_{ee#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c);N_{#gamma}^{tagged #pi^{0}}", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); // numerator for conditional probability - reinterpret_cast(list->FindObject("hMggPt_Pi0_Primary"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Pi0_Primary"))->Sumw2(); list->Add(new TH2F("hMggPt_Pi0_FromWD", "reconstructed m_{ee#gamma} vs. p_{T,#gamma} from #pi^{0} from WD;m_{ee#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c);N_{#gamma}^{tagged #pi^{0}}", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); // numerator for conditional probability - reinterpret_cast(list->FindObject("hMggPt_Pi0_FromWD"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Pi0_FromWD"))->Sumw2(); list->Add(new TH2F("hMggPt_Pi0_hs", "reconstructed m_{ee#gamma} vs. p_{T,#gamma} from #pi^{0} from hadronic shower in material;m_{ee#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c);N_{#gamma}^{tagged #pi^{0}}", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); // numerator for conditional probability - reinterpret_cast(list->FindObject("hMggPt_Pi0_hs"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Pi0_hs"))->Sumw2(); } } if (TString(histClass) == "tag_and_probe") { list->Add(new TH2F("hMggPt_Probe_Same", "m_{#gamma#gamma} vs. p_{T,#gamma};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c)", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); list->Add(new TH2F("hMggPt_Probe_Mixed", "m_{#gamma#gamma} vs. p_{T,#gamma};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c)", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); - reinterpret_cast(list->FindObject("hMggPt_Probe_Same"))->Sumw2(); - reinterpret_cast(list->FindObject("hMggPt_Probe_Mixed"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Probe_Same"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_Probe_Mixed"))->Sumw2(); list->Add(new TH2F("hMggPt_PassingProbe_Same", "m_{#gamma#gamma} vs. p_{T,#gamma};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c)", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); list->Add(new TH2F("hMggPt_PassingProbe_Mixed", "m_{#gamma#gamma} vs. p_{T,#gamma};m_{#gamma#gamma} (GeV/c^{2});p_{T,#gamma} (GeV/c)", nmgg04 - 1, mgg04, npTgg10 - 1, pTgg10)); - reinterpret_cast(list->FindObject("hMggPt_PassingProbe_Same"))->Sumw2(); - reinterpret_cast(list->FindObject("hMggPt_PassingProbe_Mixed"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_PassingProbe_Same"))->Sumw2(); + dynamic_cast(list->FindObject("hMggPt_PassingProbe_Mixed"))->Sumw2(); } if (TString(histClass) == "photon_hbt") { @@ -700,7 +700,7 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const hs_q_same->SetBinEdges(1, m_hbt); hs_q_same->SetBinEdges(2, dca_hbt); hs_q_same->SetBinEdges(3, dca_hbt); - hs_q_mix = reinterpret_cast(hs_q_same->Clone("hs_q_mix")); + hs_q_mix = dynamic_cast(hs_q_same->Clone("hs_q_mix")); list->Add(hs_q_same); list->Add(hs_q_mix); } else if (TString(subGroup) == "3d") { @@ -715,7 +715,7 @@ void o2::aod::pwgem::photon::histogram::DefineHistograms(THashList* list, const hs_q_same->SetBinEdges(1, m_hbt); hs_q_same->SetBinEdges(2, dca_hbt); hs_q_same->SetBinEdges(3, dca_hbt); - hs_q_mix = reinterpret_cast(hs_q_same->Clone("hs_q_mix")); + hs_q_mix = dynamic_cast(hs_q_same->Clone("hs_q_mix")); list->Add(hs_q_same); list->Add(hs_q_mix); } @@ -725,7 +725,7 @@ THashList* o2::aod::pwgem::photon::histogram::AddHistClass(THashList* list, cons { if (list->FindObject(histClass)) { LOGF(info, "HistogramsLibrary::AddHistClass(): Cannot add histogram class %s because it already exists.", histClass); - return static_cast(list->FindObject(histClass)); + return dynamic_cast(list->FindObject(histClass)); } auto* sublist = new THashList(); diff --git a/PWGEM/PhotonMeson/Core/HistogramsLibrary.h b/PWGEM/PhotonMeson/Core/HistogramsLibrary.h index 368acef59ce..335df4b9858 100644 --- a/PWGEM/PhotonMeson/Core/HistogramsLibrary.h +++ b/PWGEM/PhotonMeson/Core/HistogramsLibrary.h @@ -52,62 +52,62 @@ template void FillHistClass(THashList* list, const char* subGroup, T1 const& obj1 /*, const float weight = 1.f*/) { if constexpr (htype1 == EMHistType::kEvent || htype1 == EMHistType::kEvent_Cent || htype1 == EMHistType::kEvent_Cent_Qvec) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("all", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("all", 1.f); if (obj1.selection_bit(o2::aod::evsel::kNoTimeFrameBorder)) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("No TF border", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("No TF border", 1.f); } if (obj1.selection_bit(o2::aod::evsel::kNoITSROFrameBorder)) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("No ITS ROF border", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("No ITS ROF border", 1.f); } if (obj1.selection_bit(o2::aod::evsel::kNoSameBunchPileup)) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("No Same Bunch Pileup", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("No Same Bunch Pileup", 1.f); } if (obj1.selection_bit(o2::aod::evsel::kIsVertexITSTPC)) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("Is Vertex ITSTPC", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("Is Vertex ITSTPC", 1.f); } if (obj1.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV)) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("Is Good Zvtx FT0vsPV", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("Is Good Zvtx FT0vsPV", 1.f); } if (obj1.selection_bit(o2::aod::evsel::kIsTriggerTVX)) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("FT0AND", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("FT0AND", 1.f); } if (obj1.sel8()) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("sel8", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("sel8", 1.f); } if (std::abs(obj1.posZ()) < maxZ) { - reinterpret_cast(list->FindObject("hCollisionCounter"))->Fill("|Z_{vtx}| < 10 cm", 1.f); + dynamic_cast(list->FindObject("hCollisionCounter"))->Fill("|Z_{vtx}| < 10 cm", 1.f); } - reinterpret_cast(list->FindObject("hZvtx"))->Fill(obj1.posZ()); - reinterpret_cast(list->FindObject("hMultNTracksPV"))->Fill(obj1.multNTracksPV()); - reinterpret_cast(list->FindObject("hMultNTracksPVeta1"))->Fill(obj1.multNTracksPVeta1()); - reinterpret_cast(list->FindObject("hMultFT0"))->Fill(obj1.multFT0A(), obj1.multFT0C()); - reinterpret_cast(list->FindObject("hCentFT0A"))->Fill(obj1.centFT0A()); - reinterpret_cast(list->FindObject("hCentFT0C"))->Fill(obj1.centFT0C()); - reinterpret_cast(list->FindObject("hCentFT0M"))->Fill(obj1.centFT0M()); - reinterpret_cast(list->FindObject("hCentFT0MvsMultNTracksPV"))->Fill(obj1.centFT0M(), obj1.multNTracksPV()); - reinterpret_cast(list->FindObject("hMultFT0MvsMultNTracksPV"))->Fill(obj1.multFT0A() + obj1.multFT0C(), obj1.multNTracksPV()); + dynamic_cast(list->FindObject("hZvtx"))->Fill(obj1.posZ()); + dynamic_cast(list->FindObject("hMultNTracksPV"))->Fill(obj1.multNTracksPV()); + dynamic_cast(list->FindObject("hMultNTracksPVeta1"))->Fill(obj1.multNTracksPVeta1()); + dynamic_cast(list->FindObject("hMultFT0"))->Fill(obj1.multFT0A(), obj1.multFT0C()); + dynamic_cast(list->FindObject("hCentFT0A"))->Fill(obj1.centFT0A()); + dynamic_cast(list->FindObject("hCentFT0C"))->Fill(obj1.centFT0C()); + dynamic_cast(list->FindObject("hCentFT0M"))->Fill(obj1.centFT0M()); + dynamic_cast(list->FindObject("hCentFT0MvsMultNTracksPV"))->Fill(obj1.centFT0M(), obj1.multNTracksPV()); + dynamic_cast(list->FindObject("hMultFT0MvsMultNTracksPV"))->Fill(obj1.multFT0A() + obj1.multFT0C(), obj1.multNTracksPV()); if constexpr (htype1 == EMHistType::kEvent_Cent_Qvec) { - reinterpret_cast(list->FindObject("hQ2xFT0M_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xft0m()); - reinterpret_cast(list->FindObject("hQ2yFT0M_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yft0m()); - reinterpret_cast(list->FindObject("hQ2xFT0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xft0a()); - reinterpret_cast(list->FindObject("hQ2yFT0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yft0a()); - reinterpret_cast(list->FindObject("hQ2xFT0C_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xft0c()); - reinterpret_cast(list->FindObject("hQ2yFT0C_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yft0c()); - // reinterpret_cast(list->FindObject("hQ2xFV0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xfv0a()); - // reinterpret_cast(list->FindObject("hQ2yFV0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yfv0a()); - reinterpret_cast(list->FindObject("hQ2xBPos_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xbpos()); - reinterpret_cast(list->FindObject("hQ2yBPos_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2ybpos()); - reinterpret_cast(list->FindObject("hQ2xBNeg_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xbneg()); - reinterpret_cast(list->FindObject("hQ2yBNeg_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2ybneg()); + dynamic_cast(list->FindObject("hQ2xFT0M_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xft0m()); + dynamic_cast(list->FindObject("hQ2yFT0M_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yft0m()); + dynamic_cast(list->FindObject("hQ2xFT0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xft0a()); + dynamic_cast(list->FindObject("hQ2yFT0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yft0a()); + dynamic_cast(list->FindObject("hQ2xFT0C_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xft0c()); + dynamic_cast(list->FindObject("hQ2yFT0C_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yft0c()); + // static_cast(list->FindObject("hQ2xFV0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xfv0a()); + // static_cast(list->FindObject("hQ2yFV0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2yfv0a()); + dynamic_cast(list->FindObject("hQ2xBPos_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xbpos()); + dynamic_cast(list->FindObject("hQ2yBPos_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2ybpos()); + dynamic_cast(list->FindObject("hQ2xBNeg_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2xbneg()); + dynamic_cast(list->FindObject("hQ2yBNeg_CentFT0C"))->Fill(obj1.centFT0C(), obj1.q2ybneg()); - reinterpret_cast(list->FindObject("hEP2FT0M_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2ft0m()); - reinterpret_cast(list->FindObject("hEP2FT0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2ft0a()); - reinterpret_cast(list->FindObject("hEP2FT0C_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2ft0c()); - // reinterpret_cast(list->FindObject("hEP2FV0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2fv0a()); - reinterpret_cast(list->FindObject("hEP2BPos_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2bpos()); - reinterpret_cast(list->FindObject("hEP2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2bneg()); + dynamic_cast(list->FindObject("hEP2FT0M_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2ft0m()); + dynamic_cast(list->FindObject("hEP2FT0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2ft0a()); + dynamic_cast(list->FindObject("hEP2FT0C_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2ft0c()); + // static_cast(list->FindObject("hEP2FV0A_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2fv0a()); + dynamic_cast(list->FindObject("hEP2BPos_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2bpos()); + dynamic_cast(list->FindObject("hEP2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), obj1.ep2bneg()); std::array q2ft0m = {obj1.q2xft0m(), obj1.q2yft0m()}; std::array q2ft0a = {obj1.q2xft0a(), obj1.q2yft0a()}; @@ -116,110 +116,110 @@ void FillHistClass(THashList* list, const char* subGroup, T1 const& obj1 /*, con std::array q2bpos = {obj1.q2xbpos(), obj1.q2ybpos()}; std::array q2bneg = {obj1.q2xbneg(), obj1.q2ybneg()}; - reinterpret_cast(list->FindObject("hQ2FT0MQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0m, q2bpos)); - reinterpret_cast(list->FindObject("hQ2FT0MQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0m, q2bneg)); - reinterpret_cast(list->FindObject("hQ2FT0AQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0a, q2bpos)); - reinterpret_cast(list->FindObject("hQ2FT0AQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0a, q2bneg)); - reinterpret_cast(list->FindObject("hQ2FT0CQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0c, q2bpos)); - reinterpret_cast(list->FindObject("hQ2FT0CQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0c, q2bneg)); - // reinterpret_cast(list->FindObject("hQ2FV0AQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2fv0a, q2bpos)); - // reinterpret_cast(list->FindObject("hQ2FV0AQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2fv0a, q2bneg)); - reinterpret_cast(list->FindObject("hQ2BPosQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2bpos, q2bneg)); + dynamic_cast(list->FindObject("hQ2FT0MQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0m, q2bpos)); + dynamic_cast(list->FindObject("hQ2FT0MQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0m, q2bneg)); + dynamic_cast(list->FindObject("hQ2FT0AQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0a, q2bpos)); + dynamic_cast(list->FindObject("hQ2FT0AQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0a, q2bneg)); + dynamic_cast(list->FindObject("hQ2FT0CQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0c, q2bpos)); + dynamic_cast(list->FindObject("hQ2FT0CQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2ft0c, q2bneg)); + // static_cast(list->FindObject("hQ2FV0AQ2BPos_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2fv0a, q2bpos)); + // static_cast(list->FindObject("hQ2FV0AQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2fv0a, q2bneg)); + dynamic_cast(list->FindObject("hQ2BPosQ2BNeg_CentFT0C"))->Fill(obj1.centFT0C(), RecoDecay::dotProd(q2bpos, q2bneg)); } } else if constexpr (htype1 == EMHistType::kV0) { - reinterpret_cast(list->FindObject("hPt"))->Fill(obj1.pt()); - reinterpret_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); - reinterpret_cast(list->FindObject("hRadius"))->Fill(obj1.vz(), obj1.v0radius()); - reinterpret_cast(list->FindObject("hCosPA"))->Fill(obj1.cospa()); - reinterpret_cast(list->FindObject("hCosPA_Rxy"))->Fill(obj1.v0radius(), obj1.cospa()); - reinterpret_cast(list->FindObject("hPCA"))->Fill(obj1.pca()); - reinterpret_cast(list->FindObject("hPCA_CosPA"))->Fill(obj1.cospa(), obj1.pca()); - reinterpret_cast(list->FindObject("hPCA_Rxy"))->Fill(obj1.v0radius(), obj1.pca()); - reinterpret_cast(list->FindObject("hDCAxyz"))->Fill(obj1.dcaXYtopv(), obj1.dcaZtopv()); - reinterpret_cast(list->FindObject("hAPplot"))->Fill(obj1.alpha(), obj1.qtarm()); - reinterpret_cast(list->FindObject("hMassGamma"))->Fill(obj1.v0radius(), obj1.mGamma()); - reinterpret_cast(list->FindObject("hGammaRxy"))->Fill(obj1.vx(), obj1.vy()); - reinterpret_cast(list->FindObject("hKFChi2vsM"))->Fill(obj1.mGamma(), obj1.chiSquareNDF()); - reinterpret_cast(list->FindObject("hKFChi2vsR"))->Fill(obj1.v0radius(), obj1.chiSquareNDF()); - reinterpret_cast(list->FindObject("hKFChi2vsX"))->Fill(obj1.vx(), obj1.chiSquareNDF()); - reinterpret_cast(list->FindObject("hKFChi2vsY"))->Fill(obj1.vy(), obj1.chiSquareNDF()); - reinterpret_cast(list->FindObject("hKFChi2vsZ"))->Fill(obj1.vz(), obj1.chiSquareNDF()); + dynamic_cast(list->FindObject("hPt"))->Fill(obj1.pt()); + dynamic_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); + dynamic_cast(list->FindObject("hRadius"))->Fill(obj1.vz(), obj1.v0radius()); + dynamic_cast(list->FindObject("hCosPA"))->Fill(obj1.cospa()); + dynamic_cast(list->FindObject("hCosPA_Rxy"))->Fill(obj1.v0radius(), obj1.cospa()); + dynamic_cast(list->FindObject("hPCA"))->Fill(obj1.pca()); + dynamic_cast(list->FindObject("hPCA_CosPA"))->Fill(obj1.cospa(), obj1.pca()); + dynamic_cast(list->FindObject("hPCA_Rxy"))->Fill(obj1.v0radius(), obj1.pca()); + dynamic_cast(list->FindObject("hDCAxyz"))->Fill(obj1.dcaXYtopv(), obj1.dcaZtopv()); + dynamic_cast(list->FindObject("hAPplot"))->Fill(obj1.alpha(), obj1.qtarm()); + dynamic_cast(list->FindObject("hMassGamma"))->Fill(obj1.v0radius(), obj1.mGamma()); + dynamic_cast(list->FindObject("hGammaRxy"))->Fill(obj1.vx(), obj1.vy()); + dynamic_cast(list->FindObject("hKFChi2vsM"))->Fill(obj1.mGamma(), obj1.chiSquareNDF()); + dynamic_cast(list->FindObject("hKFChi2vsR"))->Fill(obj1.v0radius(), obj1.chiSquareNDF()); + dynamic_cast(list->FindObject("hKFChi2vsX"))->Fill(obj1.vx(), obj1.chiSquareNDF()); + dynamic_cast(list->FindObject("hKFChi2vsY"))->Fill(obj1.vy(), obj1.chiSquareNDF()); + dynamic_cast(list->FindObject("hKFChi2vsZ"))->Fill(obj1.vz(), obj1.chiSquareNDF()); } else if constexpr (htype1 == EMHistType::kV0Leg) { - reinterpret_cast(list->FindObject("hPt"))->Fill(obj1.pt()); - reinterpret_cast(list->FindObject("hQoverPt"))->Fill(obj1.sign() / obj1.pt()); - reinterpret_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); - reinterpret_cast(list->FindObject("hDCAxyz"))->Fill(obj1.dcaXY(), obj1.dcaZ()); - reinterpret_cast(list->FindObject("hNclsITS"))->Fill(obj1.itsNCls()); - reinterpret_cast(list->FindObject("hNclsTPC"))->Fill(obj1.tpcNClsFound()); - reinterpret_cast(list->FindObject("hNcrTPC"))->Fill(obj1.tpcNClsCrossedRows()); - reinterpret_cast(list->FindObject("hTPCNcr2Nf"))->Fill(obj1.tpcCrossedRowsOverFindableCls()); - reinterpret_cast(list->FindObject("hTPCNcls2Nf"))->Fill(obj1.tpcFoundOverFindableCls()); - reinterpret_cast(list->FindObject("hChi2TPC"))->Fill(obj1.tpcChi2NCl()); - reinterpret_cast(list->FindObject("hChi2ITS"))->Fill(obj1.itsChi2NCl()); - reinterpret_cast(list->FindObject("hITSClusterMap"))->Fill(obj1.itsClusterMap()); - reinterpret_cast(list->FindObject("hMeanClusterSizeITS"))->Fill(obj1.meanClusterSizeITS() * std::cos(std::atan(obj1.tgl()))); - reinterpret_cast(list->FindObject("hTPCdEdx"))->Fill(obj1.tpcInnerParam(), obj1.tpcSignal()); - reinterpret_cast(list->FindObject("hTPCNsigmaEl"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaEl()); - reinterpret_cast(list->FindObject("hTPCNsigmaPi"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaPi()); - reinterpret_cast(list->FindObject("hXY"))->Fill(obj1.x(), obj1.y()); - reinterpret_cast(list->FindObject("hZX"))->Fill(obj1.z(), obj1.x()); - reinterpret_cast(list->FindObject("hZY"))->Fill(obj1.z(), obj1.y()); + dynamic_cast(list->FindObject("hPt"))->Fill(obj1.pt()); + dynamic_cast(list->FindObject("hQoverPt"))->Fill(obj1.sign() / obj1.pt()); + dynamic_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); + dynamic_cast(list->FindObject("hDCAxyz"))->Fill(obj1.dcaXY(), obj1.dcaZ()); + dynamic_cast(list->FindObject("hNclsITS"))->Fill(obj1.itsNCls()); + dynamic_cast(list->FindObject("hNclsTPC"))->Fill(obj1.tpcNClsFound()); + dynamic_cast(list->FindObject("hNcrTPC"))->Fill(obj1.tpcNClsCrossedRows()); + dynamic_cast(list->FindObject("hTPCNcr2Nf"))->Fill(obj1.tpcCrossedRowsOverFindableCls()); + dynamic_cast(list->FindObject("hTPCNcls2Nf"))->Fill(obj1.tpcFoundOverFindableCls()); + dynamic_cast(list->FindObject("hChi2TPC"))->Fill(obj1.tpcChi2NCl()); + dynamic_cast(list->FindObject("hChi2ITS"))->Fill(obj1.itsChi2NCl()); + dynamic_cast(list->FindObject("hITSClusterMap"))->Fill(obj1.itsClusterMap()); + dynamic_cast(list->FindObject("hMeanClusterSizeITS"))->Fill(obj1.meanClusterSizeITS() * std::cos(std::atan(obj1.tgl()))); + dynamic_cast(list->FindObject("hTPCdEdx"))->Fill(obj1.tpcInnerParam(), obj1.tpcSignal()); + dynamic_cast(list->FindObject("hTPCNsigmaEl"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaEl()); + dynamic_cast(list->FindObject("hTPCNsigmaPi"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaPi()); + dynamic_cast(list->FindObject("hXY"))->Fill(obj1.x(), obj1.y()); + dynamic_cast(list->FindObject("hZX"))->Fill(obj1.z(), obj1.x()); + dynamic_cast(list->FindObject("hZY"))->Fill(obj1.z(), obj1.y()); } else if constexpr (htype1 == EMHistType::kTrack) { - reinterpret_cast(list->FindObject("hPt"))->Fill(obj1.pt()); - reinterpret_cast(list->FindObject("hQoverPt"))->Fill(obj1.sign() / obj1.pt()); - reinterpret_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); - reinterpret_cast(list->FindObject("hDCAxyz"))->Fill(obj1.dcaXY(), obj1.dcaZ()); - reinterpret_cast(list->FindObject("hDCAxyzSigma"))->Fill(obj1.dcaXY() / std::sqrt(obj1.cYY()), obj1.dcaZ() / std::sqrt(obj1.cZZ())); - reinterpret_cast(list->FindObject("hDCAxyRes_Pt"))->Fill(obj1.pt(), std::sqrt(obj1.cYY()) * 1e+4); // convert cm to um - reinterpret_cast(list->FindObject("hDCAzRes_Pt"))->Fill(obj1.pt(), std::sqrt(obj1.cZZ()) * 1e+4); // convert cm to um - reinterpret_cast(list->FindObject("hNclsITS"))->Fill(obj1.itsNCls()); - reinterpret_cast(list->FindObject("hNclsTPC"))->Fill(obj1.tpcNClsFound()); - reinterpret_cast(list->FindObject("hNcrTPC"))->Fill(obj1.tpcNClsCrossedRows()); - reinterpret_cast(list->FindObject("hTPCNcr2Nf"))->Fill(obj1.tpcCrossedRowsOverFindableCls()); - reinterpret_cast(list->FindObject("hTPCNcls2Nf"))->Fill(obj1.tpcFoundOverFindableCls()); - reinterpret_cast(list->FindObject("hChi2TPC"))->Fill(obj1.tpcChi2NCl()); - reinterpret_cast(list->FindObject("hChi2ITS"))->Fill(obj1.itsChi2NCl()); - reinterpret_cast(list->FindObject("hITSClusterMap"))->Fill(obj1.itsClusterMap()); - reinterpret_cast(list->FindObject("hMeanClusterSizeITS"))->Fill(obj1.meanClusterSizeITS() * std::cos(std::atan(obj1.tgl()))); - reinterpret_cast(list->FindObject("hTPCdEdx"))->Fill(obj1.tpcInnerParam(), obj1.tpcSignal()); - reinterpret_cast(list->FindObject("hTPCNsigmaEl"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaEl()); - reinterpret_cast(list->FindObject("hTPCNsigmaMu"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaMu()); - reinterpret_cast(list->FindObject("hTPCNsigmaPi"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaPi()); - reinterpret_cast(list->FindObject("hTPCNsigmaKa"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaKa()); - reinterpret_cast(list->FindObject("hTPCNsigmaPr"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaPr()); - reinterpret_cast(list->FindObject("hTOFbeta"))->Fill(obj1.tpcInnerParam(), obj1.beta()); - reinterpret_cast(list->FindObject("hTOFNsigmaEl"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaEl()); - reinterpret_cast(list->FindObject("hTOFNsigmaMu"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaMu()); - reinterpret_cast(list->FindObject("hTOFNsigmaPi"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaPi()); - reinterpret_cast(list->FindObject("hTOFNsigmaKa"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaKa()); - reinterpret_cast(list->FindObject("hTOFNsigmaPr"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaPr()); + dynamic_cast(list->FindObject("hPt"))->Fill(obj1.pt()); + dynamic_cast(list->FindObject("hQoverPt"))->Fill(obj1.sign() / obj1.pt()); + dynamic_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); + dynamic_cast(list->FindObject("hDCAxyz"))->Fill(obj1.dcaXY(), obj1.dcaZ()); + dynamic_cast(list->FindObject("hDCAxyzSigma"))->Fill(obj1.dcaXY() / std::sqrt(obj1.cYY()), obj1.dcaZ() / std::sqrt(obj1.cZZ())); + dynamic_cast(list->FindObject("hDCAxyRes_Pt"))->Fill(obj1.pt(), std::sqrt(obj1.cYY()) * 1e+4); // convert cm to um + dynamic_cast(list->FindObject("hDCAzRes_Pt"))->Fill(obj1.pt(), std::sqrt(obj1.cZZ()) * 1e+4); // convert cm to um + dynamic_cast(list->FindObject("hNclsITS"))->Fill(obj1.itsNCls()); + dynamic_cast(list->FindObject("hNclsTPC"))->Fill(obj1.tpcNClsFound()); + dynamic_cast(list->FindObject("hNcrTPC"))->Fill(obj1.tpcNClsCrossedRows()); + dynamic_cast(list->FindObject("hTPCNcr2Nf"))->Fill(obj1.tpcCrossedRowsOverFindableCls()); + dynamic_cast(list->FindObject("hTPCNcls2Nf"))->Fill(obj1.tpcFoundOverFindableCls()); + dynamic_cast(list->FindObject("hChi2TPC"))->Fill(obj1.tpcChi2NCl()); + dynamic_cast(list->FindObject("hChi2ITS"))->Fill(obj1.itsChi2NCl()); + dynamic_cast(list->FindObject("hITSClusterMap"))->Fill(obj1.itsClusterMap()); + dynamic_cast(list->FindObject("hMeanClusterSizeITS"))->Fill(obj1.meanClusterSizeITS() * std::cos(std::atan(obj1.tgl()))); + dynamic_cast(list->FindObject("hTPCdEdx"))->Fill(obj1.tpcInnerParam(), obj1.tpcSignal()); + dynamic_cast(list->FindObject("hTPCNsigmaEl"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaEl()); + dynamic_cast(list->FindObject("hTPCNsigmaMu"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaMu()); + dynamic_cast(list->FindObject("hTPCNsigmaPi"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaPi()); + dynamic_cast(list->FindObject("hTPCNsigmaKa"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaKa()); + dynamic_cast(list->FindObject("hTPCNsigmaPr"))->Fill(obj1.tpcInnerParam(), obj1.tpcNSigmaPr()); + dynamic_cast(list->FindObject("hTOFbeta"))->Fill(obj1.tpcInnerParam(), obj1.beta()); + dynamic_cast(list->FindObject("hTOFNsigmaEl"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaEl()); + dynamic_cast(list->FindObject("hTOFNsigmaMu"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaMu()); + dynamic_cast(list->FindObject("hTOFNsigmaPi"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaPi()); + dynamic_cast(list->FindObject("hTOFNsigmaKa"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaKa()); + dynamic_cast(list->FindObject("hTOFNsigmaPr"))->Fill(obj1.tpcInnerParam(), obj1.tofNSigmaPr()); } else if constexpr (htype1 == EMHistType::kPHOSCluster) { - reinterpret_cast(list->FindObject("hPt"))->Fill(obj1.pt()); - reinterpret_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); - reinterpret_cast(list->FindObject("hEvsNcell"))->Fill(obj1.e(), obj1.nCells()); - reinterpret_cast(list->FindObject("hEvsM02"))->Fill(obj1.e(), obj1.m02()); - reinterpret_cast(list->FindObject("hEvsM20"))->Fill(obj1.e(), obj1.m20()); - reinterpret_cast(list->FindObject("hDistToBC"))->Fill(obj1.distanceToBadChannel()); - reinterpret_cast(list->FindObject(Form("hClusterXZM%d", obj1.mod())))->Fill(obj1.cellx(), obj1.cellz()); + dynamic_cast(list->FindObject("hPt"))->Fill(obj1.pt()); + dynamic_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); + dynamic_cast(list->FindObject("hEvsNcell"))->Fill(obj1.e(), obj1.nCells()); + dynamic_cast(list->FindObject("hEvsM02"))->Fill(obj1.e(), obj1.m02()); + dynamic_cast(list->FindObject("hEvsM20"))->Fill(obj1.e(), obj1.m20()); + dynamic_cast(list->FindObject("hDistToBC"))->Fill(obj1.distanceToBadChannel()); + static_cast(list->FindObject(Form("hClusterXZM%d", obj1.mod())))->Fill(obj1.cellx(), obj1.cellz()); } else if constexpr (htype1 == EMHistType::kEMCCluster) { if (TString(subGroup) == "2D") { - reinterpret_cast(list->FindObject("hNCell"))->Fill(obj1.nCells(), obj1.e()); - reinterpret_cast(list->FindObject("hM02"))->Fill(obj1.m02(), obj1.e()); - reinterpret_cast(list->FindObject("hM20"))->Fill(obj1.m20(), obj1.e()); - reinterpret_cast(list->FindObject("hTime"))->Fill(obj1.time(), obj1.e()); - reinterpret_cast(list->FindObject("hDistToBC"))->Fill(obj1.distanceToBadChannel(), obj1.e()); + dynamic_cast(list->FindObject("hNCell"))->Fill(obj1.nCells(), obj1.e()); + dynamic_cast(list->FindObject("hM02"))->Fill(obj1.m02(), obj1.e()); + dynamic_cast(list->FindObject("hM20"))->Fill(obj1.m20(), obj1.e()); + dynamic_cast(list->FindObject("hTime"))->Fill(obj1.time(), obj1.e()); + dynamic_cast(list->FindObject("hDistToBC"))->Fill(obj1.distanceToBadChannel(), obj1.e()); } else { - reinterpret_cast(list->FindObject("hNCell"))->Fill(obj1.nCells()); - reinterpret_cast(list->FindObject("hM02"))->Fill(obj1.m02()); - reinterpret_cast(list->FindObject("hM20"))->Fill(obj1.m20()); - reinterpret_cast(list->FindObject("hTime"))->Fill(obj1.time()); - reinterpret_cast(list->FindObject("hDistToBC"))->Fill(obj1.distanceToBadChannel()); + dynamic_cast(list->FindObject("hNCell"))->Fill(obj1.nCells()); + dynamic_cast(list->FindObject("hM02"))->Fill(obj1.m02()); + dynamic_cast(list->FindObject("hM20"))->Fill(obj1.m20()); + dynamic_cast(list->FindObject("hTime"))->Fill(obj1.time()); + dynamic_cast(list->FindObject("hDistToBC"))->Fill(obj1.distanceToBadChannel()); } - reinterpret_cast(list->FindObject("hPt"))->Fill(obj1.pt()); - reinterpret_cast(list->FindObject("hE"))->Fill(obj1.e()); - reinterpret_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); + dynamic_cast(list->FindObject("hPt"))->Fill(obj1.pt()); + dynamic_cast(list->FindObject("hE"))->Fill(obj1.e()); + dynamic_cast(list->FindObject("hEtaPhi"))->Fill(obj1.phi(), obj1.eta()); for (size_t itrack = 0; itrack < obj1.deltaEta().size(); itrack++) { // Fill TrackEtaPhi histogram with delta phi and delta eta of all tracks saved in the vectors in skimmerGammaCalo.cxx - reinterpret_cast(list->FindObject("hTrackEtaPhi"))->Fill(obj1.deltaPhi()[itrack], obj1.deltaEta()[itrack]); + dynamic_cast(list->FindObject("hTrackEtaPhi"))->Fill(obj1.deltaPhi()[itrack], obj1.deltaEta()[itrack]); } } } diff --git a/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h b/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h index 22ee95e1bc3..5e8502a4e86 100644 --- a/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h +++ b/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGamma.h @@ -282,8 +282,8 @@ struct Pi0EtaToGammaGamma { o2::framework::PresliceOptional perEMCClusterMT = o2::aod::mintm::minClusterId; o2::framework::PresliceOptional perEMCClusterMS = o2::aod::mintm::minClusterId; - o2::framework::Partition>> positrons = o2::aod::emprimaryelectron::sign > int8_t(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt&& nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; - o2::framework::Partition>> electrons = o2::aod::emprimaryelectron::sign < int8_t(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; + o2::framework::Partition>> positrons = o2::aod::emprimaryelectron::sign > static_cast(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt&& nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; + o2::framework::Partition>> electrons = o2::aod::emprimaryelectron::sign < static_cast(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; o2::aod::pwgem::dilepton::utils::EventMixingHandler, std::pair, o2::aod::pwgem::photonmeson::utils::EMPhoton>* emh1 = nullptr; o2::aod::pwgem::dilepton::utils::EventMixingHandler, std::pair, o2::aod::pwgem::photonmeson::utils::EMPhoton>* emh2 = nullptr; diff --git a/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGammaMC.h b/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGammaMC.h index e8c2333ccd3..0d1325e7338 100644 --- a/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGammaMC.h +++ b/PWGEM/PhotonMeson/Core/Pi0EtaToGammaGammaMC.h @@ -482,8 +482,8 @@ struct Pi0EtaToGammaGammaMC { o2::framework::PresliceOptional perEMCClusterMT = o2::aod::mintm::minClusterId; o2::framework::PresliceOptional perEMCClusterMS = o2::aod::mintm::minClusterId; - o2::framework::Partition>> positrons = o2::aod::emprimaryelectron::sign > int8_t(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt&& nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; - o2::framework::Partition>> electrons = o2::aod::emprimaryelectron::sign < int8_t(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; + o2::framework::Partition>> positrons = o2::aod::emprimaryelectron::sign > static_cast(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt&& nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; + o2::framework::Partition>> electrons = o2::aod::emprimaryelectron::sign < static_cast(0) && dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track; //--------------------------------------------------------------------------- // In the following are tags defined which help to select the correct preslice and cuts diff --git a/PWGEM/PhotonMeson/Core/TaggingPi0.h b/PWGEM/PhotonMeson/Core/TaggingPi0.h index 64080cd2a2d..9abb932afbb 100644 --- a/PWGEM/PhotonMeson/Core/TaggingPi0.h +++ b/PWGEM/PhotonMeson/Core/TaggingPi0.h @@ -429,8 +429,8 @@ struct TaggingPi0 { // o2::framework::Preslice perCollision_phos = o2::aod::phoscluster::pmeventId; o2::framework::Preslice perCollision_electron = o2::aod::emprimaryelectronda::pmeventId; - o2::framework::Partition positrons = o2::aod::emprimaryelectron::sign > int8_t(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt&& nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); - o2::framework::Partition electrons = o2::aod::emprimaryelectron::sign < int8_t(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt && nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); + o2::framework::Partition positrons = o2::aod::emprimaryelectron::sign > static_cast(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt&& nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); + o2::framework::Partition electrons = o2::aod::emprimaryelectron::sign < static_cast(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt && nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); MyEMH* emh1 = nullptr; MyEMH* emh2 = nullptr; diff --git a/PWGEM/PhotonMeson/Core/TaggingPi0MC.h b/PWGEM/PhotonMeson/Core/TaggingPi0MC.h index 4588f70182d..128e8c99c58 100644 --- a/PWGEM/PhotonMeson/Core/TaggingPi0MC.h +++ b/PWGEM/PhotonMeson/Core/TaggingPi0MC.h @@ -404,8 +404,8 @@ struct TaggingPi0MC { // o2::framework::Preslice perCollision_phos = o2::aod::phoscluster::pmeventId; o2::framework::Preslice perCollision_electron = o2::aod::emprimaryelectronda::pmeventId; - o2::framework::Partition positrons = o2::aod::emprimaryelectron::sign > int8_t(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt&& nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); - o2::framework::Partition electrons = o2::aod::emprimaryelectron::sign < int8_t(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt && nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); + o2::framework::Partition positrons = o2::aod::emprimaryelectron::sign > static_cast(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt&& nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); + o2::framework::Partition electrons = o2::aod::emprimaryelectron::sign < static_cast(0) && static_cast(dileptoncuts.cfg_min_pt_track) < o2::aod::track::pt && nabs(o2::aod::track::eta) < static_cast(dileptoncuts.cfg_max_eta_track) && static_cast(dileptoncuts.cfg_min_TPCNsigmaEl) < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < static_cast(dileptoncuts.cfg_max_TPCNsigmaEl); template void runTruePairing(TCollisions const& collisions, diff --git a/PWGEM/PhotonMeson/Legacy/skimmerDalitzEE.cxx b/PWGEM/PhotonMeson/Legacy/skimmerDalitzEE.cxx index fc64bbba7fd..30f9133aa04 100644 --- a/PWGEM/PhotonMeson/Legacy/skimmerDalitzEE.cxx +++ b/PWGEM/PhotonMeson/Legacy/skimmerDalitzEE.cxx @@ -338,8 +338,8 @@ struct skimmerDalitzEE { return npair; } - Partition posTracks = o2::aod::emprimaryelectron::sign > int8_t(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; - Partition negTracks = o2::aod::emprimaryelectron::sign < int8_t(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; + Partition posTracks = o2::aod::emprimaryelectron::sign > static_cast(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; + Partition negTracks = o2::aod::emprimaryelectron::sign < static_cast(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; void processAnalysis(MyCollisions const& collisions, MyTracks const&) { for (auto& collision : collisions) { @@ -367,8 +367,8 @@ struct skimmerDalitzEE { } PROCESS_SWITCH(skimmerDalitzEE, processAnalysis, "Process dalitz ee for analysis", true); - Partition posTracks_cefp = o2::aod::emprimaryelectron::sign > int8_t(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; - Partition negTracks_cefp = o2::aod::emprimaryelectron::sign < int8_t(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; + Partition posTracks_cefp = o2::aod::emprimaryelectron::sign > static_cast(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; + Partition negTracks_cefp = o2::aod::emprimaryelectron::sign < static_cast(0) && o2::aod::track::pt > minpt&& nabs(o2::aod::track::eta) < maxeta&& minTPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl&& o2::aod::pidtpc::tpcNSigmaEl < maxTPCNsigmaEl; void processCEFP(aod::Collisions const& collisions, MyTracksCEFP const&, aod::BCsWithTimestamps const&) { for (auto& collision : collisions) { diff --git a/PWGEM/PhotonMeson/TableProducer/photonconversionbuilder.h b/PWGEM/PhotonMeson/TableProducer/photonconversionbuilder.h index 483f38098b7..4213c366045 100644 --- a/PWGEM/PhotonMeson/TableProducer/photonconversionbuilder.h +++ b/PWGEM/PhotonMeson/TableProducer/photonconversionbuilder.h @@ -1499,7 +1499,7 @@ struct PhotonConversionBuilder { vecV0Dedup.clear(); } // end of build - o2::framework::expressions::Filter v0Filter = o2::aod::v0::v0Type > (uint8_t)0; + o2::framework::expressions::Filter v0Filter = o2::aod::v0::v0Type > static_cast(0); using FilteredV0s = o2::soa::Filtered; enum PairFate { diff --git a/PWGEM/PhotonMeson/Tasks/dalitzEEQC.cxx b/PWGEM/PhotonMeson/Tasks/dalitzEEQC.cxx index 43f82e92041..44a744cceb5 100644 --- a/PWGEM/PhotonMeson/Tasks/dalitzEEQC.cxx +++ b/PWGEM/PhotonMeson/Tasks/dalitzEEQC.cxx @@ -356,8 +356,8 @@ struct DalitzEEQC { Filter trackFilter = dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track && o2::aod::track::tpcChi2NCl < dileptoncuts.cfg_max_chi2tpc && o2::aod::track::itsChi2NCl < dileptoncuts.cfg_max_chi2its && nabs(o2::aod::track::dcaXY) < dileptoncuts.cfg_max_dcaxy && nabs(o2::aod::track::dcaZ) < dileptoncuts.cfg_max_dcaz; Filter pidFilter = dileptoncuts.cfg_min_TPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < dileptoncuts.cfg_max_TPCNsigmaEl && (o2::aod::pidtpc::tpcNSigmaPi < dileptoncuts.cfg_min_TPCNsigmaPi || dileptoncuts.cfg_max_TPCNsigmaPi < o2::aod::pidtpc::tpcNSigmaPi); using FilteredMyTracks = soa::Filtered; - Partition posTracks = o2::aod::emprimaryelectron::sign > int8_t(0); - Partition negTracks = o2::aod::emprimaryelectron::sign < int8_t(0); + Partition posTracks = o2::aod::emprimaryelectron::sign > static_cast(0); + Partition negTracks = o2::aod::emprimaryelectron::sign < static_cast(0); std::vector used_trackIds; void processQC(FilteredMyCollisions const& collisions, FilteredMyTracks const& /*tracks*/) diff --git a/PWGEM/PhotonMeson/Tasks/dalitzEEQCMC.cxx b/PWGEM/PhotonMeson/Tasks/dalitzEEQCMC.cxx index c0dd75f3545..8ba01ad39be 100644 --- a/PWGEM/PhotonMeson/Tasks/dalitzEEQCMC.cxx +++ b/PWGEM/PhotonMeson/Tasks/dalitzEEQCMC.cxx @@ -506,8 +506,8 @@ struct DalitzEEQCMC { Filter trackFilter = dileptoncuts.cfg_min_pt_track < o2::aod::track::pt && nabs(o2::aod::track::eta) < dileptoncuts.cfg_max_eta_track && o2::aod::track::tpcChi2NCl < dileptoncuts.cfg_max_chi2tpc && o2::aod::track::itsChi2NCl < dileptoncuts.cfg_max_chi2its && nabs(o2::aod::track::dcaXY) < dileptoncuts.cfg_max_dcaxy && nabs(o2::aod::track::dcaZ) < dileptoncuts.cfg_max_dcaz; Filter pidFilter = dileptoncuts.cfg_min_TPCNsigmaEl < o2::aod::pidtpc::tpcNSigmaEl && o2::aod::pidtpc::tpcNSigmaEl < dileptoncuts.cfg_max_TPCNsigmaEl && (o2::aod::pidtpc::tpcNSigmaPi < dileptoncuts.cfg_min_TPCNsigmaPi || dileptoncuts.cfg_max_TPCNsigmaPi < o2::aod::pidtpc::tpcNSigmaPi); using FilteredMyMCTracks = soa::Filtered; - Partition posTracks = o2::aod::emprimaryelectron::sign > int8_t(0); - Partition negTracks = o2::aod::emprimaryelectron::sign < int8_t(0); + Partition posTracks = o2::aod::emprimaryelectron::sign > static_cast(0); + Partition negTracks = o2::aod::emprimaryelectron::sign < static_cast(0); Filter collisionFilter_centrality = (cfgCentMin < o2::aod::cent::centFT0M && o2::aod::cent::centFT0M < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0A && o2::aod::cent::centFT0A < cfgCentMax) || (cfgCentMin < o2::aod::cent::centFT0C && o2::aod::cent::centFT0C < cfgCentMax); using FilteredMyCollisions = soa::Filtered; diff --git a/PWGEM/PhotonMeson/Tasks/phosQC.cxx b/PWGEM/PhotonMeson/Tasks/phosQC.cxx index b4302803248..1d5c5ecd6ba 100644 --- a/PWGEM/PhotonMeson/Tasks/phosQC.cxx +++ b/PWGEM/PhotonMeson/Tasks/phosQC.cxx @@ -64,11 +64,11 @@ struct phosQC { // create sub lists first. o2::aod::pwgem::photon::histogram::AddHistClass(fMainList, "Event"); - THashList* list_ev = reinterpret_cast(fMainList->FindObject("Event")); + THashList* list_ev = dynamic_cast(fMainList->FindObject("Event")); o2::aod::pwgem::photon::histogram::DefineHistograms(list_ev, "Event"); o2::aod::pwgem::photon::histogram::AddHistClass(fMainList, "Cluster"); - THashList* list_cluster = reinterpret_cast(fMainList->FindObject("Cluster")); + THashList* list_cluster = dynamic_cast(fMainList->FindObject("Cluster")); for (const auto& cut : fPHOSCuts) { std::string cutname = cut.getName(); @@ -78,7 +78,7 @@ struct phosQC { // for Clusters for (auto& cut : fPHOSCuts) { const std::string& cutname = cut.getName(); - THashList* list = reinterpret_cast(fMainList->FindObject("Cluster")->FindObject(cutname.c_str())); + THashList* list = dynamic_cast(fMainList->FindObject("Cluster")->FindObject(cutname.c_str())); o2::aod::pwgem::photon::histogram::DefineHistograms(list, "Cluster", "PHOS"); } } @@ -105,16 +105,16 @@ struct phosQC { DefineCuts(); addhistograms(); // please call this after DefinCuts(); - fOutputEvent.setObject(reinterpret_cast(fMainList->FindObject("Event"))); - fOutputCluster.setObject(reinterpret_cast(fMainList->FindObject("Cluster"))); + fOutputEvent.setObject(dynamic_cast(fMainList->FindObject("Event"))); + fOutputCluster.setObject(dynamic_cast(fMainList->FindObject("Cluster"))); } Preslice perCollision = aod::skimmedcluster::collisionId; void processQC(MyCollisions const& collisions, aod::PHOSClusters const& clusters) { - THashList* list_ev = static_cast(fMainList->FindObject("Event")); - THashList* list_cluster = static_cast(fMainList->FindObject("Cluster")); + THashList* list_ev = dynamic_cast(fMainList->FindObject("Event")); + THashList* list_cluster = dynamic_cast(fMainList->FindObject("Cluster")); for (auto& collision : collisions) { @@ -122,27 +122,27 @@ struct phosQC { continue; } - reinterpret_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(1.0); + dynamic_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(1.0); if (!collision.sel8()) { continue; } - reinterpret_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(2.0); + dynamic_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(2.0); if (collision.numContrib() < 0.5) { continue; } - reinterpret_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(3.0); + dynamic_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(3.0); if (abs(collision.posZ()) > 10.0) { continue; } - reinterpret_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(4.0); - reinterpret_cast(fMainList->FindObject("Event")->FindObject("hZvtx"))->Fill(collision.posZ()); + dynamic_cast(fMainList->FindObject("Event")->FindObject("hCollisionCounter"))->Fill(4.0); + dynamic_cast(fMainList->FindObject("Event")->FindObject("hZvtx"))->Fill(collision.posZ()); o2::aod::pwgem::photon::histogram::FillHistClass(list_ev, "", collision); auto clusters_per_coll = clusters.sliceBy(perCollision, collision.collisionId()); for (const auto& cut : fPHOSCuts) { - THashList* list_cluster_cut = static_cast(list_cluster->FindObject(cut.getName().c_str())); + THashList* list_cluster_cut = dynamic_cast(list_cluster->FindObject(cut.getName().c_str())); int ng = 0; for (auto& cluster : clusters_per_coll) { @@ -151,7 +151,7 @@ struct phosQC { ng++; } } // end of v0 loop - reinterpret_cast(fMainList->FindObject("Cluster")->FindObject(cut.getName().c_str())->FindObject("hNgamma"))->Fill(ng); + dynamic_cast(fMainList->FindObject("Cluster")->FindObject(cut.getName().c_str())->FindObject("hNgamma"))->Fill(ng); } // end of cut loop } // end of collision loop } // end of process diff --git a/PWGEM/PhotonMeson/Tasks/prefilterPhoton.cxx b/PWGEM/PhotonMeson/Tasks/prefilterPhoton.cxx index 8c1b1ed5bcb..989aa53c7fc 100644 --- a/PWGEM/PhotonMeson/Tasks/prefilterPhoton.cxx +++ b/PWGEM/PhotonMeson/Tasks/prefilterPhoton.cxx @@ -593,8 +593,8 @@ struct prefilterPhoton { Filter collisionFilter_occupancy_ft0c = eventcuts.cfgFT0COccupancyMin <= o2::aod::evsel::ft0cOccupancyInTimeRange && o2::aod::evsel::ft0cOccupancyInTimeRange < eventcuts.cfgFT0COccupancyMax; using FilteredMyCollisions = soa::Filtered; - Partition posTracks = o2::aod::emprimaryelectron::sign > int8_t(0); - Partition negTracks = o2::aod::emprimaryelectron::sign < int8_t(0); + Partition posTracks = o2::aod::emprimaryelectron::sign > static_cast(0); + Partition negTracks = o2::aod::emprimaryelectron::sign < static_cast(0); void processPCMPCM(FilteredMyCollisions const& collisions, MyV0Photons const& v0s, aod::V0Legs const& v0legs) { diff --git a/PWGEM/Tasks/phosElId.cxx b/PWGEM/Tasks/phosElId.cxx index 6c753a2b704..a08310375e2 100644 --- a/PWGEM/Tasks/phosElId.cxx +++ b/PWGEM/Tasks/phosElId.cxx @@ -192,7 +192,7 @@ struct PhosElId { dcazfilter = nabs(aod::track::dcaZ) < DCAzMax, itschi2filter = aod::track::itsChi2NCl < ITSchi2Max, tpcchi2filter = aod::track::tpcChi2NCl < TPCchi2Max, - mapfilter = (aod::track::itsClusterMap & uint8_t(1)) > 0; + mapfilter = (aod::track::itsClusterMap & static_cast(1)) > 0; Service ccdb; std::unique_ptr geomPHOS; @@ -346,10 +346,10 @@ struct PhosElId { geomPHOS = std::make_unique("PHOS"); fSigma_dz = new TF1("fSigma_dz", "[0]/(x+[1])^[2]+pol1(3)", 0.3, 10); - fSigma_dz->SetParameters(((std::vector)pSigmadZ).at(0), ((std::vector)pSigmadZ).at(1), ((std::vector)pSigmadZ).at(2), ((std::vector)pSigmadZ).at(3), ((std::vector)pSigmadZ).at(4)); + fSigma_dz->SetParameters((std::vector(pSigmadZ)).at(0), (std::vector(pSigmadZ)).at(1), (std::vector(pSigmadZ)).at(2), (std::vector(pSigmadZ)).at(3), (std::vector(pSigmadZ)).at(4)); fSigma_dx = new TF1("fSigma_dx", "[0]/x^[1]+[2]", 0.1, 10); - fSigma_dx->SetParameters(((std::vector)pSigmadX).at(0), ((std::vector)pSigmadX).at(1), ((std::vector)pSigmadX).at(2)); + fSigma_dx->SetParameters((std::vector(pSigmadX)).at(0), (std::vector(pSigmadX)).at(1), (std::vector(pSigmadX)).at(2)); for (int i = 0; i < mAmountOfModules; i++) { fMeandXPosMod[i] = new TF1(Form("funcMeandx_pos_mod%i", i + 1), "[0]/(x+[1])^[2]", 0.1, 10); @@ -398,7 +398,7 @@ struct PhosElId { if (!track.has_collision() || !track.hasTPC()) continue; - if (track.itsNCls() < ITSnclsMin || track.itsNCls() > ITSnclsMax || !((track.itsClusterMap() & uint8_t(1)) > 0)) + if (track.itsNCls() < ITSnclsMin || track.itsNCls() > ITSnclsMax || !((track.itsClusterMap() & static_cast(1)) > 0)) continue; if (track.tpcNClsFound() < TPCnclsMin || track.tpcNClsFound() > TPCnclsMax) continue; @@ -648,7 +648,7 @@ struct PhosElId { if (!matchedTrack.has_collision() || !matchedTrack.hasTPC()) continue; - if (matchedTrack.itsNCls() < ITSnclsMin || matchedTrack.itsNCls() > ITSnclsMax || !((matchedTrack.itsClusterMap() & uint8_t(1)) > 0)) + if (matchedTrack.itsNCls() < ITSnclsMin || matchedTrack.itsNCls() > ITSnclsMax || !((matchedTrack.itsClusterMap() & static_cast(1)) > 0)) continue; if (matchedTrack.tpcNClsFound() < TPCnclsMin || matchedTrack.tpcNClsFound() > TPCnclsMax) continue; @@ -735,7 +735,7 @@ struct PhosElId { if (!track.has_collision() || !track.hasTPC()) continue; - if (track.itsNCls() < ITSnclsMin || track.itsNCls() > ITSnclsMax || !((track.itsClusterMap() & uint8_t(1)) > 0)) + if (track.itsNCls() < ITSnclsMin || track.itsNCls() > ITSnclsMax || !((track.itsClusterMap() & static_cast(1)) > 0)) continue; if (track.tpcNClsFound() < TPCnclsMin || track.tpcNClsFound() > TPCnclsMax) continue; @@ -1117,7 +1117,7 @@ struct PhosElId { if (!matchedTrack.has_collision() || !matchedTrack.hasTPC()) continue; - if (matchedTrack.itsNCls() < ITSnclsMin || matchedTrack.itsNCls() > ITSnclsMax || !((matchedTrack.itsClusterMap() & uint8_t(1)) > 0)) + if (matchedTrack.itsNCls() < ITSnclsMin || matchedTrack.itsNCls() > ITSnclsMax || !((matchedTrack.itsClusterMap() & static_cast(1)) > 0)) continue; if (matchedTrack.tpcNClsFound() < TPCnclsMin || matchedTrack.tpcNClsFound() > TPCnclsMax) continue; @@ -1170,13 +1170,13 @@ struct PhosElId { bool isWithinNSigma(int16_t& mod, float p, float deltaZ, float deltaX, bool positiveCharge) { int modMinus1 = mod - 1; - if (std::fabs(deltaZ - ((std::vector)pPhosShiftZ).at(modMinus1)) > NsigmaTrackMatch * fSigma_dz->Eval(p)) + if (std::fabs(deltaZ - (std::vector(pPhosShiftZ)).at(modMinus1)) > NsigmaTrackMatch * fSigma_dz->Eval(p)) return false; if (positiveCharge) { - if (std::fabs(deltaX - fMeandXPosMod[modMinus1]->Eval(p) - ((std::vector)pPhosShiftX).at(modMinus1)) > NsigmaTrackMatch * fSigma_dx->Eval(p)) + if (std::fabs(deltaX - fMeandXPosMod[modMinus1]->Eval(p) - (std::vector(pPhosShiftX)).at(modMinus1)) > NsigmaTrackMatch * fSigma_dx->Eval(p)) return false; } else { - if (std::fabs(deltaX - fMeandXNegMod[modMinus1]->Eval(p) - ((std::vector)pPhosShiftX).at(modMinus1)) > NsigmaTrackMatch * fSigma_dx->Eval(p)) + if (std::fabs(deltaX - fMeandXNegMod[modMinus1]->Eval(p) - (std::vector(pPhosShiftX)).at(modMinus1)) > NsigmaTrackMatch * fSigma_dx->Eval(p)) return false; } return true; @@ -1407,7 +1407,7 @@ struct MassSpectra { continue; if (TPCel.itsChi2NCl() > fITSchi2Max) continue; - if (TPCel.itsNCls() < fITSnclsMin || TPCel.itsNCls() > fITSnclsMax || !((TPCel.itsClusterMap() & uint8_t(1)) > 0)) + if (TPCel.itsNCls() < fITSnclsMin || TPCel.itsNCls() > fITSnclsMax || !((TPCel.itsClusterMap() & static_cast(1)) > 0)) continue; if (TPCel.tpcChi2NCl() > fTPCchi2Max) continue; @@ -1751,7 +1751,7 @@ struct TpcElIdMassSpectrum { return false; if (track.tpcChi2NCl() >= TPCchi2Max) return false; - if (!((track.itsClusterMap() & uint8_t(1)) > 0)) + if (!((track.itsClusterMap() & static_cast(1)) > 0)) return false; if (track.itsNCls() < ITSnclsMin || track.itsNCls() > ITSnclsMax) return false; @@ -1831,7 +1831,7 @@ struct TpcElIdMassSpectrum { return false; if (track.tpcChi2NCl() >= TPCchi2Max) return false; - if (!((track.itsClusterMap() & uint8_t(1)) > 0)) + if (!((track.itsClusterMap() & static_cast(1)) > 0)) return false; if (track.itsNCls() < ITSnclsMin || track.itsNCls() > ITSnclsMax) return false; @@ -1856,7 +1856,7 @@ struct TpcElIdMassSpectrum { auto isGoodTagElectron = [&](const MyTracks::iterator& track) -> bool { if (!track.has_collision() || !track.hasTPC()) return false; - if (!((track.itsClusterMap() & uint8_t(1)) > 0)) + if (!((track.itsClusterMap() & static_cast(1)) > 0)) return false; if (track.itsChi2NCl() > ITSchi2Max || track.tpcChi2NCl() > TPCchi2Max) return false; @@ -1889,7 +1889,7 @@ struct TpcElIdMassSpectrum { auto isGoodProbeBaseTrack = [&](const MyTracks::iterator& track) -> bool { if (!track.has_collision() || !track.hasTPC()) return false; - if (!((track.itsClusterMap() & uint8_t(1)) > 0)) + if (!((track.itsClusterMap() & static_cast(1)) > 0)) return false; if (track.itsChi2NCl() > ITSchi2Max || track.tpcChi2NCl() > TPCchi2Max) return false; @@ -2219,7 +2219,7 @@ struct TpcElIdMassSpectrum { return false; if (track.tpcChi2NCl() >= TPCchi2Max) return false; - if (!((track.itsClusterMap() & uint8_t(1)) > 0)) + if (!((track.itsClusterMap() & static_cast(1)) > 0)) return false; if (track.itsNCls() < ITSnclsMin || track.itsNCls() > ITSnclsMax) return false; @@ -2258,7 +2258,7 @@ struct TpcElIdMassSpectrum { auto isGoodTagElectron = [&](auto const& track) -> bool { if (!track.has_collision() || !track.hasTPC()) return false; - if (!((track.itsClusterMap() & uint8_t(1)) > 0)) + if (!((track.itsClusterMap() & static_cast(1)) > 0)) return false; if (track.itsChi2NCl() > ITSchi2Max || track.tpcChi2NCl() > TPCchi2Max) return false; @@ -2291,7 +2291,7 @@ struct TpcElIdMassSpectrum { auto isGoodProbeBaseTrack = [&](auto const& track) -> bool { if (!track.has_collision() || !track.hasTPC()) return false; - if (!((track.itsClusterMap() & uint8_t(1)) > 0)) + if (!((track.itsClusterMap() & static_cast(1)) > 0)) return false; if (track.itsChi2NCl() > ITSchi2Max || track.tpcChi2NCl() > TPCchi2Max) return false;