Merge branch 'feature-movesplit'

[u/mrichter/AliRoot.git] / PWGLF / RESONANCES / macros / mini / AddAnalysisTaskSigmaStarPP.C

diff --git a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
index 399d998ec46a6b5ba650b51f4c19b2acbef0c42d..781fff3b36d2905b2b0f4cfc182824c3d43696e8 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
+++ b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
@@ -15,19 +15,26 @@
 //
 ****************************************************************************/
 
-AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStar
+AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPP
 (
    Bool_t      isMC,
    Bool_t      isPP,
+   Float_t     cutV = 10.0,
    Int_t       aodFilterBit = 5,
-   Int_t       piPIDCut = 3.0,
-   Int_t       pPIDCut = 3.0,
-   Int_t       nmix = 0,
+   Float_t     piPIDCut = 3.0,
+   Float_t     pPIDCut = 3.0,
+   Float_t     trackDCAcut = 7.0,
+   Float_t     massTol = 0.01,
+   Float_t     lambdaDCA = 0.3,
+   Float_t     lambdaCosPoinAn = 0.99,
+   Float_t     lambdaDaughDCA = 0.5,
+   Int_t       NTPCcluster = 70,
+   Int_t       nmix = 5,
    Float_t     maxDiffVzMix = 1.0,
    Float_t     maxDiffMultMix = 10.0,
    Float_t     maxDiffAngleMixDeg = 20.0,
    Int_t       aodN = 0,
-   TString     outNameSuffix = ""
+   TString     outNameSuffix = "Sigma1385"
 )
 {  
   //
@@ -42,7 +49,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStar
    } 
 
    // create the task and configure 
-   TString taskName = Form("SigmaStarPP%s%s_%i%i", (isPP? "pp" : "PbPb"), (isMC ? "MC" : "Data"));
+   TString taskName = Form("SigmaStar%s%s_%.1f_%d_%.1f_%.1f_%.2f_%.2f_%.1f_%.2f_%.1f", (isPP? "pp" : "PbPb"), (isMC ? "MC" : "Data"),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA);
    AliRsnMiniAnalysisTask *task = new AliRsnMiniAnalysisTask(taskName.Data(), isMC);
    if (!isMC && !isPP){
      Printf(Form("========== SETTING USE CENTRALITY PATCH AOD049 : %s", (aodN==49)? "yes" : "no"));
@@ -59,7 +66,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStar
    task->SetMaxDiffVz(maxDiffVzMix);
    task->SetMaxDiffMult(maxDiffMultMix);
    if (!isPP) task->SetMaxDiffAngle(maxDiffAngleMixDeg*TMath::DegToRad()); //set angle diff in rad
-   ::Info("AddAnalysisTaskSigmaStar", Form("Event mixing configuration: \n events to mix = %i \n max diff. vtxZ = cm %5.3f \n max diff multi = %5.3f \n max diff EP angle = %5.3f deg", nmix, maxDiffVzMix, maxDiffMultMix, (isPP ? 0.0 : maxDiffAngleMixDeg)));
+   ::Info("AddAnalysisTaskSigmaStarPP", Form("Event mixing configuration: \n events to mix = %i \n max diff. vtxZ = cm %5.3f \n max diff multi = %5.3f \n max diff EP angle = %5.3f deg", nmix, maxDiffVzMix, maxDiffMultMix, (isPP ? 0.0 : maxDiffAngleMixDeg)));
    
    mgr->AddTask(task);
    
@@ -70,7 +77,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStar
    // - 2nd argument --> |Vz| range
    // - 3rd argument --> minimum required number of contributors
    // - 4th argument --> tells if TPC stand-alone vertexes must be accepted
-   AliRsnCutPrimaryVertex *cutVertex = new AliRsnCutPrimaryVertex("cutVertex", 10.0, 0, kFALSE);
+   AliRsnCutPrimaryVertex *cutVertex = new AliRsnCutPrimaryVertex("cutVertex", cutV, 0, kFALSE);
    if (isPP) cutVertex->SetCheckPileUp(kTRUE);   // set the check for pileup
    
    // define and fill cut set for event cut
@@ -98,15 +105,17 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStar
    
    //event plane (only for PbPb)
    Int_t planeID = task->CreateValue(AliRsnMiniValue::kPlaneAngle, kFALSE);
-   AliRsnMiniOutput *outPlane = task->CreateOutput("eventPlane", "HIST", "EVENT");
-   if (!isPP)
+   AliRsnMiniOutput *outPlane = 0x0; //task->CreateOutput("eventPlane", "HIST", "EVENT");
+   if (!isPP){
+     outPlane = task->CreateOutput("eventPlane", "HIST", "EVENT");
      outPlane->AddAxis(planeID, 180, 0.0, TMath::Pi());
+     }
    
    //
    // -- PAIR CUTS (common to all resonances) ------------------------------------------------------
    //
    AliRsnCutMiniPair *cutY = new AliRsnCutMiniPair("cutRapidity", AliRsnCutMiniPair::kRapidityRange);
-   cutY->SetRangeD(-0.8, 0.8);
+   cutY->SetRangeD(-0.5, 0.5);
    
    AliRsnCutSet *cutsPair = new AliRsnCutSet("pairCuts", AliRsnTarget::kMother);
    cutsPair->AddCut(cutY);
@@ -116,22 +125,19 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStar
    // -- CONFIG ANALYSIS --------------------------------------------------------------------------
    gROOT->LoadMacro("$ALICE_ROOT/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C");
    if (isMC) {
-       Printf("========================== MC analysis - PID cuts used");
-       piPIDCut = 1E20; pPIDCut = 1E20;
-       //ConfigSigmaStar(task, isPP, isMC, piPIDCut, pPIDCut, aodFilterBit, "", cutsPair);
+       Printf("========================== MC analysis - PID cuts not used");
    } else 
      Printf("========================== DATA analysis - PID cuts used");
-     //ConfigSigmaStar(task, isPP, isMC, piPIDCut, pPIDCut, aodFilterBit, "", cutsPair);
-   if (!ConfigSigmaStar(task, isMC, isPP, piPIDCuts, pPIDCuts, aodFilterBit, "", cutsPair) return 0x0;
+   if (!ConfigSigmaStar(task, isPP, isMC, piPIDCut, pPIDCut, aodFilterBit, trackDCAcut, massTol, lambdaDCA, lambdaCosPoinAn, lambdaDaughDCA, NTPCcluster, "", cutsPair)) return 0x0;
    
    //
    // -- CONTAINERS --------------------------------------------------------------------------------
    //
    TString outputFileName = AliAnalysisManager::GetCommonFileName();
    //  outputFileName += ":Rsn";
-   Printf("AddAnalysisTaskSigmaStar - Set OutputFileName : \n %s\n", outputFileName.Data() );
+   Printf("AddAnalysisTaskSigmaStarPP - Set OutputFileName : \n %s\n", outputFileName.Data() );
    
-   AliAnalysisDataContainer *output = mgr->CreateContainer(Form("RsnOut_%s",outNameSuffix.Data()), 
+   AliAnalysisDataContainer *output = mgr->CreateContainer(Form("RsnOut_%s_%.1f_%d_%.1f_%.1f_%.2f_%.2f_%.1f_%.2f_%.1f",outNameSuffix.Data(),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA), 
                                                           TList::Class(), 
                                                           AliAnalysisManager::kOutputContainer, 
                                                           outputFileName);