Updated macros from Massimo

diff --git a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskD0.C b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskD0.C
index 8819f0031b0691806e08383b68a1f63911b697eb..1693ad4edbdd41c1dfe2994c9b41185c7a131aff 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskD0.C
+++ b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskD0.C
@@ -44,6 +44,11 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
    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,
@@ -196,7 +201,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskD0
        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 --------------------------------------------------------------------------------

diff --git a/PWGLF/RESONANCES/macros/mini/ConfigD0.C b/PWGLF/RESONANCES/macros/mini/ConfigD0.C
index a83e04cea7aa39536147de8fa893d76591702eaf..dc5281405b23e46b0ad772b6c306b3788b13e779 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/ConfigD0.C
+++ b/PWGLF/RESONANCES/macros/mini/ConfigD0.C
@@ -1,4 +1,4 @@
-//
+/
 // *** Configuration script for phi->KK analysis with 2010 runs ***
 // 
 // A configuration script for RSN package needs to define the followings:
@@ -31,6 +31,11 @@ Bool_t ConfigD0
    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
@@ -66,7 +71,7 @@ Bool_t ConfigD0
    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);
@@ -100,7 +105,7 @@ Bool_t ConfigD0
    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);
@@ -173,7 +178,7 @@ Bool_t ConfigD0
 
       // 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
@@ -235,7 +240,7 @@ Bool_t ConfigD0
    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);
@@ -261,7 +266,7 @@ Bool_t ConfigD0
    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);
@@ -342,7 +347,7 @@ Bool_t ConfigD0
    // 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);
 
@@ -359,7 +364,7 @@ Bool_t ConfigD0
    // 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);