@@ -715,56 +715,6 @@ struct FlowGfwOmegaXi {
715715 return ;
716716 }
717717
718- // // remove auto-corr
719- // template <char... chars>
720- // void fillProfilepTMass(const GFW::CorrConfig& corrconf, const GFW::CorrConfig& corrconfol, const ConstStr<chars...>& tarName, const int& ptbin, const int& PDGCode, const float& cent)
721- // {
722- // int nMassBins = 0;
723- // int nptbins = 0;
724- // TAxis* fMass = nullptr;
725- // TAxis* fpt = nullptr;
726- // if (PDGCode == kXiMinus) {
727- // nMassBins = cfgmassbins[2];
728- // nptbins = nXiPtBins;
729- // fpt = fXiPtAxis;
730- // fMass = fXiMass;
731- // } else if (PDGCode == kOmegaMinus) {
732- // nMassBins = cfgmassbins[3];
733- // nptbins = nOmegaPtBins;
734- // fpt = fOmegaPtAxis;
735- // fMass = fOmegaMass;
736- // } else if (PDGCode == kK0Short) {
737- // nMassBins = cfgmassbins[0];
738- // nptbins = nK0sPtBins;
739- // fpt = fK0sPtAxis;
740- // fMass = fK0sMass;
741- // } else if (PDGCode == kLambda0) {
742- // nMassBins = cfgmassbins[1];
743- // nptbins = nLambdaPtBins;
744- // fpt = fLambdaPtAxis;
745- // fMass = fLambdaMass;
746- // } else {
747- // LOGF(error, "Error, please put in correct PDGCode of K0s, Lambda, Xi or Omega");
748- // return;
749- // }
750- // for (int massbin = 1; massbin <= nMassBins; massbin++) {
751- // float dnx = 0;
752- // float val = 0;
753- // float dnxol = 0;
754- // dnx = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kTRUE).real();
755- // dnxol = fGFW->Calculate(corrconfol, (ptbin - 1) + ((massbin - 1) * nptbins), kTRUE).real();
756- // dnx = dnx - dnxol;
757- // if (dnx == 0)
758- // continue;
759- // val = (fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kFALSE).real() - fGFW->Calculate(corrconfol, (ptbin - 1) + ((massbin - 1) * nptbins), kFALSE).real()) / dnx;
760-
761- // if (std::fabs(val) < 1) {
762- // registry.fill(tarName, fpt->GetBinCenter(ptbin), fMass->GetBinCenter(massbin), cent, val, dnx);
763- // }
764- // }
765- // return;
766- // }
767-
768718 // input shared_ptr<TProfile3D>
769719 void fillProfilepTMass (const GFW::CorrConfig& corrconf, std::shared_ptr<TProfile3D> TProfile3D, const int & ptbin, const int & PDGCode, const float & cent)
770720 {
@@ -810,54 +760,6 @@ struct FlowGfwOmegaXi {
810760 return ;
811761 }
812762
813- // remove auto-corr
814- // void fillProfilepTMass(const GFW::CorrConfig& corrconf, const GFW::CorrConfig& corrconfol, std::shared_ptr<TProfile3D> TProfile3D, const int& ptbin, const int& PDGCode, const float& cent)
815- // {
816- // int nMassBins = 0;
817- // int nptbins = 0;
818- // TAxis* fMass = nullptr;
819- // TAxis* fpt = nullptr;
820- // if (PDGCode == kXiMinus) {
821- // nMassBins = cfgmassbins[2];
822- // nptbins = nXiPtBins;
823- // fpt = fXiPtAxis;
824- // fMass = fXiMass;
825- // } else if (PDGCode == kOmegaMinus) {
826- // nMassBins = cfgmassbins[3];
827- // nptbins = nOmegaPtBins;
828- // fpt = fOmegaPtAxis;
829- // fMass = fOmegaMass;
830- // } else if (PDGCode == kK0Short) {
831- // nMassBins = cfgmassbins[0];
832- // nptbins = nK0sPtBins;
833- // fpt = fK0sPtAxis;
834- // fMass = fK0sMass;
835- // } else if (PDGCode == kLambda0) {
836- // nMassBins = cfgmassbins[1];
837- // nptbins = nLambdaPtBins;
838- // fpt = fLambdaPtAxis;
839- // fMass = fLambdaMass;
840- // } else {
841- // LOGF(error, "Error, please put in correct PDGCode of K0s, Lambda, Xi or Omega");
842- // return;
843- // }
844- // for (int massbin = 1; massbin <= nMassBins; massbin++) {
845- // float dnx = 0;
846- // float val = 0;
847- // float dnxol = 0;
848- // dnx = fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kTRUE).real();
849- // dnxol = fGFW->Calculate(corrconfol, (ptbin - 1) + ((massbin - 1) * nptbins), kTRUE).real();
850- // dnx = dnx - dnxol;
851- // if (dnx == 0)
852- // continue;
853- // val = (fGFW->Calculate(corrconf, (ptbin - 1) + ((massbin - 1) * nptbins), kFALSE).real() - fGFW->Calculate(corrconfol, (ptbin - 1) + ((massbin - 1) * nptbins), kFALSE).real()) / dnx;
854- // if (std::fabs(val) < 1) {
855- // TProfile3D->Fill(fpt->GetBinCenter(ptbin), fMass->GetBinCenter(massbin), cent, val, dnx);
856- // }
857- // }
858- // return;
859- // }
860-
861763 void loadCorrections (uint64_t timestamp)
862764 {
863765 if (correctionsLoaded)
0 commit comments