/*************************************************************************** adash@cern.ch - last modified on 03/04/2013 // // General macro to configure the RSN analysis task. // It calls all configs desired by the user, by means // of the boolean switches defined in the first lines. // --- // Inputs: // 1) flag to know if running on MC or data // 2) path where all configs are stored // --- // Returns: // kTRUE --> initialization successful // kFALSE --> initialization failed (some config gave errors) // ****************************************************************************/ AliRsnMiniAnalysisTask * AddTaskPhiPPb_TPC ( Bool_t isMC, Bool_t isPP, Int_t nmix = 0, AliPID::EParticleType type1, Float_t nsigmaKa = 2.0, Bool_t enableMonitor = kTRUE, Float_t maxDiffVzMix = 1.0, Float_t maxDiffMultMix = 10.0, TString outNameSuffix = "", TString optSy="Default" ) { // // -- INITIALIZATION ---------------------------------------------------------------------------- // retrieve analysis manager // AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); if (!mgr) { ::Error("AddTaskPhiPPb_TPC", "No analysis manager to connect to."); return NULL; } // create the task and configure TString taskName = Form("TPCPhiMeson%s%s_%i", (isPP? "pp" : "pPb"), (isMC ? "MC" : "Data"), (Int_t)type1 ); AliRsnMiniAnalysisTask *task = new AliRsnMiniAnalysisTask(taskName.Data(), isMC); //task->SelectCollisionCandidates(AliVEvent::kINT7); task->UseESDTriggerMask(AliVEvent::kINT7); if (isPP) task->UseMultiplicity("QUALITY"); else task->UseCentrality("V0A"); // set event mixing options task->UseContinuousMix(); //task->UseBinnedMix(); task->SetNMix(nmix); task->SetMaxDiffVz(maxDiffVzMix); task->SetMaxDiffMult(maxDiffMultMix); task->UseMC(isMC); ::Info("AddTaskPhiPPb_TPC", Form("Event mixing configuration: \n events to mix = %i \n max diff. vtxZ = cm %5.3f \n max diff multi = %5.3f \n", nmix, maxDiffVzMix, maxDiffMultMix)); mgr->AddTask(task); // // -- EVENT CUTS (same for all configs) --------------------------------------------------------- // // cut on primary vertex: // - 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); if (isPP) cutVertex->SetCheckPileUp(kTRUE); // set the check for pileup //set check for pileup in 2013 Bool_t rmFirstEvtChunk = kTRUE; Bool_t rejectPileUp = kTRUE; Int_t MinPlpContribSPD = 1; Bool_t useMVPileUpSelection = kTRUE; Int_t MinPlpContribMV = 5; AliRsnCutEventUtils *cutEventUtils = new AliRsnCutEventUtils("cutEventUtils", rmFirstEvtChunk, rejectPileUp); if (useMVPileUpSelection){ cutEventUtils->SetUseMVPlpSelection(useMVPileUpSelection); cutEventUtils->SetMinPlpContribMV(MinPlpContribMV); cutEventUtils->SetMinPlpContribSPD(MinPlpContribSPD); cutEventUtils->SetUseVertexSelection2013pA(kTRUE); ::Info("AddTaskPhiPPb_TPC", Form("Multiple-vtx Pile-up rejection settings: MinPlpContribMV = %i, MinPlpContribSPD = %i", MinPlpContribMV, MinPlpContribSPD)); } else { cutEventUtils->SetMinPlpContribSPD(MinPlpContribSPD); ::Info("AddTaskPhiPPb_TPC", Form("SPD Pile-up rejection settings: MinPlpContribSPD = %i", MinPlpContribSPD)); } ::Info("AddTaskPhiPPb_TPC", Form(":::::::::::::::::: Pile-up rejection mode: %s", (rejectPileUp?"ON":"OFF"))); ::Info("AddTaskPhiPPb_TPC", Form("::::::::::::: Remove first event in chunk: %s", (rmFirstEvtChunk?"ON":"OFF"))); // define and fill cut set for event cut AliRsnCutSet *eventCuts = new AliRsnCutSet("eventCuts", AliRsnTarget::kEvent); eventCuts->AddCut(cutVertex); eventCuts->AddCut(cutEventUtils); eventCuts->SetCutScheme(Form("%s&%s", cutEventUtils->GetName(), cutVertex->GetName())); // set cuts in task task->SetEventCuts(eventCuts); // // -- EVENT-ONLY COMPUTATIONS ------------------------------------------------------------------- // //vertex Int_t vtxID = task->CreateValue(AliRsnMiniValue::kVz, kFALSE); AliRsnMiniOutput *outVtx = task->CreateOutput("eventVtx", "HIST", "EVENT"); outVtx->AddAxis(vtxID, 400, -20.0, 20.0); //multiplicity or centrality Int_t multID = task->CreateValue(AliRsnMiniValue::kMult, kFALSE); AliRsnMiniOutput *outMult = task->CreateOutput("eventMult", "HIST", "EVENT"); if (isPP) outMult->AddAxis(multID, 400, 0.0, 400.0); else outMult->AddAxis(multID, 100, 0.0, 100.0); // // -- PAIR CUTS (common to all resonances) ------------------------------------------------------ // AliRsnCutMiniPair *cutY = new AliRsnCutMiniPair("cutRapidity", AliRsnCutMiniPair::kRapidityRange); cutY->SetRangeD(-0.465, 0.035);// 0 < y_cm < 0.5; y_cm = y_lab + 0.465 AliRsnCutSet *cutsPair = new AliRsnCutSet("pairCuts", AliRsnTarget::kMother); cutsPair->AddCut(cutY); cutsPair->SetCutScheme(cutY->GetName()); // // -- CONFIG ANALYSIS -------------------------------------------------------------------------- // gROOT->LoadMacro("$ALICE_ROOT/PWGLF/RESONANCES/macros/mini/ConfigPhiPPb_TPC.C"); if (!ConfigPhiPPb_TPC(task, isMC, isPP, "", cutsPair, type1, nsigmaKa, enableMonitor,optSy)) return 0x0; // // -- CONTAINERS -------------------------------------------------------------------------------- // TString outputFileName = AliAnalysisManager::GetCommonFileName(); // outputFileName += ":Rsn"; Printf("AddTaskPhiPPb_TPC - Set OutputFileName : \n %s\n", outputFileName.Data() ); AliAnalysisDataContainer *output = mgr->CreateContainer(Form("RsnOut_%s",outNameSuffix.Data()), TList::Class(), AliAnalysisManager::kOutputContainer, outputFileName); mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer()); mgr->ConnectOutput(task, 1, output); return task; }