@@ -140,17 +140,19 @@ struct HfTaskElectronWeakBoson {
140140 Configurable<std::vector<double >> centralityBins{" centralityBins" , {0 , 20 , 60 , 100 }, " centrality bins" };
141141
142142 // QA for Z->ee
143- Configurable<bool > enableZeeRecoQA{" enableZeeRecoQA" , false , " Enable QA for Z->ee reconstruction" };
143+ Configurable<bool > enableZeeTrkRecoQA{" enableZeeTrkRecoQA" , false , " Enable QA for track info Z->ee reconstruction" };
144+ Configurable<bool > enableZeeEmcRecoQA{" enableZeeEmcRecoQA" , false , " Enable QA for Emc info Z->ee reconstruction" };
144145 Configurable<int > ZeeRecoQAtype{" ZeeRecoQAtype" , 1 , " QA for Z->ee reconstruction" };
145146 Configurable<float > massZMinQA{" massZMinQA" , 0.1 , " minimum mass cut for Zee Reco QA" };
146- // CCDB service object
147- Service<o2::ccdb::BasicCCDBManager> ccdb{};
148147 // UE
149148 Configurable<int > nRandomCones{" nRandomCones" , 100 , " number of random cones" };
150149 Configurable<float > rcHardE{" rcHardE" , 5.0 , " hard cluster veto energy" };
151150 Configurable<float > rcVetoR{" rcVetoR" , 0.4 , " veto radius" };
152151 Configurable<bool > useUEsub{" useUEsub" , true , " apply UE subtraction in isolation" };
153152
153+ // CCDB service object
154+ Service<o2::ccdb::BasicCCDBManager> ccdb{};
155+
154156 struct HfElectronCandidate {
155157 float pt, eta, phi, dcaxyTrk, dcazTrk, eop, energyIso, momIso, dedxTrk, m02Emc;
156158 int ntrackIso, nclusterTPC, nclusterITS;
@@ -812,7 +814,7 @@ struct HfTaskElectronWeakBoson {
812814 // LOG(info) << "R mim = " << rMin;
813815 registry.fill (HIST (" hTrMatch_mim" ), dPhiMin, dEtaMin);
814816 }
815- if (enableZeeRecoQA && track.pt () > ptZeeMin) {
817+ if ((enableZeeTrkRecoQA || enableZeeEmcRecoQA) && track.pt () > ptZeeMin) {
816818 if (track.sign () < 0 ) {
817819 selectedElectronsIso.emplace_back (
818820 track.pt (),
@@ -871,7 +873,7 @@ struct HfTaskElectronWeakBoson {
871873 }
872874 } // end of Z-hadron correlation
873875 // Z->ee QA
874- if (enableZeeRecoQA ) {
876+ if (enableZeeTrkRecoQA || enableZeeEmcRecoQA ) {
875877 if (!selectedElectronsIso.empty () && !selectedPositronsIso.empty ()) {
876878 // signal
877879 for (const auto & trackEle : selectedElectronsIso) {
@@ -882,10 +884,10 @@ struct HfTaskElectronWeakBoson {
882884 double const invMass = RecoDecay::m (std::array{child1, child2}, std::array{o2::constants::physics::MassElectron, o2::constants::physics::MassElectron});
883885 float const sectorpos = trackPos.phi / o2::constants::math::SectorSpanRad;
884886 if (invMass > massZMinQA) {
885- if (ZeeRecoQAtype == 1 ) {
887+ if (enableZeeTrkRecoQA ) {
886888 registry.fill (HIST (" hInvMassZeeTrkQA" ), invMass, trackEle.pt , trackPos.pt , trackEle.dcaxyTrk , trackPos.dcaxyTrk , trackPos.dcazTrk , trackEle.nclusterTPC , trackPos.nclusterTPC , trackEle.nclusterITS , trackPos.nclusterITS , sectorneg, sectorpos, trackEle.eop , trackPos.eop , trackEle.energyIso , trackPos.energyIso , trackEle.momIso , trackPos.momIso , trackEle.ntrackIso , trackPos.ntrackIso );
887889 }
888- if (ZeeRecoQAtype == 2 ) {
890+ if (enableZeeEmcRecoQA ) {
889891 registry.fill (HIST (" hInvMassZeeEmcQA" ), invMass, trackEle.pt , trackPos.pt , trackEle.eop , trackPos.eop , trackEle.m02Emc , trackPos.m02Emc , trackEle.dedxTrk , trackPos.dedxTrk , trackEle.energyIso , trackPos.energyIso , trackEle.momIso , trackPos.momIso , trackEle.ntrackIso , trackPos.ntrackIso );
890892 }
891893 }
@@ -902,10 +904,10 @@ struct HfTaskElectronWeakBoson {
902904 float const sectorpos = trackEle2.phi / o2::constants::math::SectorSpanRad;
903905 double const invMass = RecoDecay::m (std::array{child1, child2}, std::array{o2::constants::physics::MassElectron, o2::constants::physics::MassElectron});
904906 if (invMass > massZMinQA) {
905- if (ZeeRecoQAtype == 1 ) {
907+ if (enableZeeTrkRecoQA ) {
906908 registry.fill (HIST (" hInvMassZeeTrkQAbg" ), invMass, trackEle.pt , trackEle2.pt , trackEle.dcaxyTrk , trackEle2.dcaxyTrk , trackEle2.dcazTrk , trackEle.nclusterTPC , trackEle2.nclusterTPC , trackEle.nclusterITS , trackEle2.nclusterITS , sectorneg, sectorpos, trackEle.eop , trackEle2.eop , trackEle.energyIso , trackEle2.energyIso , trackEle.momIso , trackEle2.momIso , trackEle.ntrackIso , trackEle2.ntrackIso );
907909 }
908- if (ZeeRecoQAtype == 2 ) {
910+ if (enableZeeEmcRecoQA ) {
909911 registry.fill (HIST (" hInvMassZeeEmcQAbg" ), invMass, trackEle.pt , trackEle2.pt , trackEle.eop , trackEle2.eop , trackEle.m02Emc , trackEle2.m02Emc , trackEle.dedxTrk , trackEle2.dedxTrk , trackEle.energyIso , trackEle2.energyIso , trackEle.momIso , trackEle2.momIso , trackEle.ntrackIso , trackEle2.ntrackIso );
910912 }
911913 }
@@ -922,10 +924,10 @@ struct HfTaskElectronWeakBoson {
922924 float const sectorpos = trackPos2.phi / o2::constants::math::SectorSpanRad;
923925 double const invMass = RecoDecay::m (std::array{child1, child2}, std::array{o2::constants::physics::MassElectron, o2::constants::physics::MassElectron});
924926 if (invMass > massZMinQA) {
925- if (ZeeRecoQAtype == 1 ) {
927+ if (enableZeeTrkRecoQA ) {
926928 registry.fill (HIST (" hInvMassZeeTrkQAbg" ), invMass, trackPos.pt , trackPos2.pt , trackPos.dcaxyTrk , trackPos2.dcaxyTrk , trackPos2.dcazTrk , trackPos.nclusterTPC , trackPos2.nclusterTPC , trackPos.nclusterITS , trackPos2.nclusterITS , sectorneg, sectorpos, trackPos.eop , trackPos2.eop , trackPos.energyIso , trackPos2.energyIso , trackPos.momIso , trackPos2.momIso , trackPos.ntrackIso , trackPos2.ntrackIso );
927929 }
928- if (ZeeRecoQAtype == 2 ) {
930+ if (enableZeeEmcRecoQA ) {
929931 registry.fill (HIST (" hInvMassZeeEmcQAbg" ), invMass, trackPos.pt , trackPos2.pt , trackPos.eop , trackPos2.eop , trackPos.m02Emc , trackPos2.m02Emc , trackPos.dedxTrk , trackPos2.dedxTrk , trackPos.energyIso , trackPos2.energyIso , trackPos.momIso , trackPos2.momIso , trackPos.ntrackIso , trackPos2.ntrackIso );
930932 }
931933 }
0 commit comments