Skip to content

[PWGLF] added few configurables for calibration and event selections #15575

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
prottayCMT wants to merge 1 commit into
base: master
Choose a base branch
from
+26 −26
Draft

[PWGLF] added few configurables for calibration and event selections #15575

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
48 changes: 23 additions & 25 deletions PWGLF/TableProducer/Common/spvector.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -86,6 +86,7 @@ struct spvector {
Configurable<float> cfgCutCentralityMin{"cfgCutCentralityMin", 0.0f, "Centrality cut Min"};
Configurable<bool> additionalEvSel{"additionalEvSel", false, "additionalEvSel"};
Configurable<bool> usemem{"usemem", true, "usemem"};
Configurable<bool> usecfactor{"usecfactor", false, "use c factor"};

struct : ConfigurableGroup {
Configurable<int> QxyNbins{"QxyNbins", 100, "Number of bins in QxQy histograms"};
Expand Down Expand Up @@ -293,11 +294,6 @@ struct spvector {
int lastRunNumber = -999;
TH2D* gainprofile;
TProfile* gainprofilevxy;
/*THnF* hrecentereSp;
TH2F* hrecenterecentSp;
TH2F* hrecenterevxSp;
TH2F* hrecenterevySp;
TH2F* hrecenterevzSp;*/
std::array<THnF*, 6> hrecentereSpA; // Array of 6 histograms
std::array<TH2F*, 6> hrecenterecentSpA; // Array of 5 histograms
std::array<TH2F*, 6> hrecenterevxSpA; // Array of 5 histograms
Expand All @@ -306,16 +302,9 @@ struct spvector {
TProfile3D* shiftprofileA;
TProfile3D* shiftprofileC;

// Bool_t Correctcoarse(int64_t ts, Configurable<std::string>& ConfRecentereSpp, bool useRecentereSp, int currentRunNumber, int lastRunNumber, auto centrality, auto vx, auto vy, auto vz, auto& qxZDCA, auto& qyZDCA, auto& qxZDCC, auto& qyZDCC)
//{
Bool_t Correctcoarse(const THnF* hrecentereSp, auto centrality, auto vx, auto vy, auto vz, auto& qxZDCA, auto& qyZDCA, auto& qxZDCC, auto& qyZDCC)
{

/*
if (useRecentereSp && (currentRunNumber != lastRunNumber)) {
hrecentereSp = ccdb->getForTimeStamp<THnF>(ConfRecentereSpp.value, ts);
}*/

int binCoords[5];

// Get axes of the THnSparse
Expand Down Expand Up @@ -357,22 +346,13 @@ struct spvector {
return kTRUE;
}

// Bool_t Correctfine(int64_t ts, Configurable<std::string>& ConfRecenterecentSpp, Configurable<std::string>& ConfRecenterevxSpp, Configurable<std::string>& ConfRecenterevySpp, Configurable<std::string>& ConfRecenterevzSpp, bool useRecenterefineSp, int currentRunNumber, int lastRunNumber, auto centrality, auto vx, auto vy, auto vz, auto& qxZDCA, auto& qyZDCA, auto& qxZDCC, auto& qyZDCC)
//{
Bool_t Correctfine(TH2F* hrecenterecentSp, TH2F* hrecenterevxSp, TH2F* hrecenterevySp, TH2F* hrecenterevzSp, auto centrality, auto vx, auto vy, auto vz, auto& qxZDCA, auto& qyZDCA, auto& qxZDCC, auto& qyZDCC)
{

if (!hrecenterecentSp || !hrecenterevxSp || !hrecenterevySp || !hrecenterevzSp) {
std::cerr << "Error: One or more histograms are null." << std::endl;
return false;
}
/*
if (useRecenterefineSp && (currentRunNumber != lastRunNumber)) {
hrecenterecentSp = ccdb->getForTimeStamp<TH2F>(ConfRecenterecentSpp.value, ts);
hrecenterevxSp = ccdb->getForTimeStamp<TH2F>(ConfRecenterevxSpp.value, ts);
hrecenterevySp = ccdb->getForTimeStamp<TH2F>(ConfRecenterevySpp.value, ts);
hrecenterevzSp = ccdb->getForTimeStamp<TH2F>(ConfRecenterevzSpp.value, ts);
}*/

double meanxAcent = hrecenterecentSp->GetBinContent(hrecenterecentSp->FindBin(centrality + 0.00001, 0.5));
double meanyAcent = hrecenterecentSp->GetBinContent(hrecenterecentSp->FindBin(centrality + 0.00001, 1.5));
Expand Down Expand Up @@ -452,6 +432,7 @@ struct spvector {
auto znaEnergy = zdc.energySectorZNA();
auto zncEnergycommon = zdc.energyCommonZNC();
auto znaEnergycommon = zdc.energyCommonZNA();
auto beamEne = 5.36 * 0.5;

if (znaEnergycommon <= 0.0 || zncEnergycommon <= 0.0) {
triggerevent = false;
Expand Down Expand Up @@ -504,6 +485,8 @@ struct spvector {
auto alphaZDC = 0.395;
constexpr double x[4] = {-1.75, 1.75, -1.75, 1.75};
constexpr double y[4] = {-1.75, -1.75, 1.75, 1.75};
double zncEnergycommonsum = 0.0;
double znaEnergycommonsum = 0.0;

histos.fill(HIST("ZDCAmpCommon"), 0.5, vz, znaEnergycommon);
histos.fill(HIST("ZDCAmpCommon"), 1.5, vz, zncEnergycommon);
Expand All @@ -522,6 +505,7 @@ struct spvector {
return;
} else {
double ampl = gainequal * znaEnergy[iChA];
znaEnergycommonsum += ampl;
if (followpub) {
ampl = TMath::Power(ampl, alphaZDC);
}
Expand All @@ -537,6 +521,7 @@ struct spvector {
return;
} else {
double ampl = gainequal * zncEnergy[iChA - 4];
zncEnergycommonsum += ampl;
if (followpub) {
ampl = TMath::Power(ampl, alphaZDC);
}
Expand All @@ -548,13 +533,26 @@ struct spvector {
}
}

auto cZNC = 1.0;
auto cZNA = 1.0;

if (sumA > 0) {
qxZDCA = qxZDCA / sumA;
qyZDCA = qyZDCA / sumA;
float nSpecnA = znaEnergycommonsum / beamEne;
if (usecfactor)
cZNA = 1.89358 - 0.71262 / (nSpecnA + 0.71789);
else
cZNA = 1.0;
qxZDCA = cZNA * (qxZDCA / sumA);
qyZDCA = cZNA * (qyZDCA / sumA);
}
if (sumC > 0) {
qxZDCC = qxZDCC / sumC;
qyZDCC = qyZDCC / sumC;
float nSpecnC = zncEnergycommonsum / beamEne;
if (usecfactor)
cZNC = 1.89358 - 0.71262 / (nSpecnC + 0.71789);
else
cZNC = 1.0;
qxZDCC = cZNC * (qxZDCC / sumC);
qyZDCC = cZNC * (qyZDCC / sumC);
}

if (sumA <= 1e-4 || sumC <= 1e-4) {
Expand Down
4 changes: 3 additions & 1 deletion PWGLF/TableProducer/Strangeness/lambdaspincorrelation.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -73,6 +73,8 @@ struct lambdaspincorrelation {
} rctCut;
Configurable<bool> useNoCollInTimeRangeStandard{"useNoCollInTimeRangeStandard", false, "Apply kNoCollInTimeRangeStandard selection bit"};
Configurable<bool> useGoodITSLayersAll{"useGoodITSLayersAll", true, "Apply kIsGoodITSLayersAll selection bit"};
Configurable<bool> useNoSameBunchPileup{"useNoSameBunchPileup", true, "Apply NoSameBunchPileup selection bit"};
Configurable<bool> useGoodZvtxFT0vsPV{"useGoodZvtxFT0vsPV", true, "Apply GoodZvtxFT0vsPV selection bit"};
// mixing
// Produce derived tables
Configurable<int> cfgCutOccupancy{"cfgCutOccupancy", 2000, "Occupancy cut"};
Expand Down Expand Up @@ -270,7 +272,7 @@ struct lambdaspincorrelation {
int occupancy = collision.trackOccupancyInTimeRange();
histos.fill(HIST("hEvtSelInfo"), 0.5);
// if ((!rctCut.requireRCTFlagChecker || rctChecker(collision)) && collision.selection_bit(aod::evsel::kNoSameBunchPileup) && collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) && collision.selection_bit(aod::evsel::kNoTimeFrameBorder) && collision.selection_bit(aod::evsel::kNoITSROFrameBorder) && collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard) && collision.sel8() && collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll) && occupancy < cfgCutOccupancy) {
if ((!rctCut.requireRCTFlagChecker || rctChecker(collision)) && collision.selection_bit(aod::evsel::kNoSameBunchPileup) && collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) && collision.selection_bit(aod::evsel::kNoTimeFrameBorder) && collision.selection_bit(aod::evsel::kNoITSROFrameBorder) && (!useNoCollInTimeRangeStandard || collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) && collision.sel8() && (!useGoodITSLayersAll || collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) && occupancy < cfgCutOccupancy) {
if ((!rctCut.requireRCTFlagChecker || rctChecker(collision)) && (!useNoSameBunchPileup || collision.selection_bit(aod::evsel::kNoSameBunchPileup)) && (!useGoodZvtxFT0vsPV || collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV)) && collision.selection_bit(aod::evsel::kNoTimeFrameBorder) && collision.selection_bit(aod::evsel::kNoITSROFrameBorder) && (!useNoCollInTimeRangeStandard || collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) && collision.sel8() && (!useGoodITSLayersAll || collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) && occupancy < cfgCutOccupancy) {
histos.fill(HIST("hEvtSelInfo"), 1.5);
for (const auto& v0 : V0s) {
// LOGF(info, "v0 index 0 : (%d)", v0.index());
Expand Down
Loading