Int_t nmix = 5,
Double_t minYlab = -0.5,
Double_t maxYlab = 0.5,
+ Float_t mineta = -0.8,
+ Float_t maxeta = 0.8,
+ Float_t min_inv_mass = 0.6,
+ Float_t max_inv_mass = 2.2,
+ Int_t bins = 320,
Double_t dcaProduct = -1E-4,
Float_t maxDiffVzMix = 1.0,
Float_t maxDiffMultMix = 10.0,
Printf("========================== MC analysis - PID cuts used");
} else
Printf("========================== DATA analysis - PID cuts used");
- if (!ConfigD0(task, isPP, isMC, nsigmaTPCPi, nsigmaTPCKa, nsigmaTOFPi, nsigmaTOFKa, aodFilterBit, trackDCAcutMax, trackDCAcutMin, trackDCAZcutMax, NTPCcluster, minpt, maxSisters, checkP, minDCAcutFixed, maxDCAcutFixed, ptdepPIDcut, checkFeedDown, checkQuark, originDselection, "", cutsPairY, cutsPair)) return 0x0;
+ if (!ConfigD0(task, isPP, isMC, nsigmaTPCPi, nsigmaTPCKa, nsigmaTOFPi, nsigmaTOFKa, aodFilterBit, trackDCAcutMax, trackDCAcutMin, trackDCAZcutMax, NTPCcluster, minpt, maxSisters, checkP, minDCAcutFixed, maxDCAcutFixed, ptdepPIDcut, checkFeedDown, checkQuark, originDselection, mineta, maxeta, min_inv_mass, max_inv_mass, bins, "", cutsPairY, cutsPair)) return 0x0;
//
// -- CONTAINERS --------------------------------------------------------------------------------
-//
+/
// *** Configuration script for phi->KK analysis with 2010 runs ***
//
// A configuration script for RSN package needs to define the followings:
Bool_t checkFeedDown = kTRUE,
Bool_t checkQuark = kTRUE,
UShort_t originDselection = 0,
+ Float_t mineta = -0.8,
+ Float_t maxeta = 0.8,
+ Float_t min_inv_mass = 0.6,
+ Float_t max_inv_mass = 2.2,
+ Int_t bins = 320,
const char *suffix,
AliRsnCutSet *cutsPairY,
AliRsnCutSet *cutsPair
else cutQuality->SetDCARPtFormulaMin(formulaMin);
cutQuality->SetTPCminNClusters(NTPCcluster);
cutQuality->SetPtRange(minpt,1E20);
- cutQuality->SetEtaRange(-0.8, 0.8);
+ cutQuality->SetEtaRange(mineta, maxeta);
cutQuality->SetDCAZmax(trackDCAZcutMax);
cutQuality->SetSPDminNClusters(1);
cutQuality->SetITSminNClusters(0);
else cutQuality->SetDCARPtFormulaMin(formulaMin);
cutQuality->SetTPCminNClusters(NTPCcluster);
cutQuality->SetPtRange(minpt,1E20);
- cutQuality->SetEtaRange(-0.8, 0.8);
+ cutQuality->SetEtaRange(mineta, maxeta);
cutQuality->SetDCAZmax(trackDCAZcutMax);
cutQuality->SetSPDminNClusters(1);
cutQuality->SetITSminNClusters(0);
// axis X: invmass (or resolution)
//if (useIM[i])
- out->AddAxis(imID, 320, 0.6, 2.2);
+ out->AddAxis(imID, bins, min_inv_mass, max_inv_mass);
//else
// out->AddAxis(resID, 200, -0.02, 0.02);
// axis Y: transverse momentum
out->SetRejectCandidateIfNotFromQuark(checkQuark);
out->SetDselection(originDselection);
// binnings
- out->AddAxis(imID, 320, 0.6, 2.2);
+ out->AddAxis(imID, bins, min_inv_mass, max_inv_mass);
out->AddAxis(ptID, 200, 0.0, 20.0);
//out->AddAxis(yID, 100, -1, 1);
//out->AddAxis(dcapID, 100, -0.001, 0.001);
out->SetRejectCandidateIfNotFromQuark(checkQuark);
out->SetDselection(originDselection);
// binnings
- out->AddAxis(imID, 320, 0.6, 2.2);
+ out->AddAxis(imID, bins, min_inv_mass, max_inv_mass);
out->AddAxis(ptID, 200, 0.0, 20.0);
//out->AddAxis(yID, 100, -1, 1);
//out->AddAxis(dcapID, 100, -0.001, 0.001);
// pair cuts
out->SetPairCuts(cutsPairY);
// binnings
- out->AddAxis(imID, 320, 0.6, 2.2);
+ out->AddAxis(imID, bins, min_inv_mass, max_inv_mass);
out->AddAxis(ptID, 200, 0.0, 20.0);
//out->AddAxis(yID, 100, -1, 1);
// pair cuts
out->SetPairCuts(cutsPairY);
// binnings
- out->AddAxis(imID, 320, 0.6, 2.2);
+ out->AddAxis(imID, bins, min_inv_mass, max_inv_mass);
out->AddAxis(ptID, 200, 0.0, 20.0);
//out->AddAxis(yID, 100, -1, 1);