/** * @file AddTaskForwardQA.C * @author Christian Holm Christensen * @date Wed Mar 23 12:14:03 2011 * * @brief Include the Forward QA task in a train. * * @ingroup pwglf_forward_scripts_tasks */ /** * @defgroup pwglf_forward_qa Quality Assurance * * Code to deal with Quality Assurance * * @ingroup pwglf_forward_topical */ /** * This is the macro to include the Forward QA task in a train. * * @param mc Monte-carlo input * @param useCent Use centrality * * @return newly constructured task object * * @ingroup pwglf_forward_eloss */ AliAnalysisTask* AddTaskForwardQA(Bool_t mc=false, Bool_t useCent=false) { // --- Load libraries ---------------------------------------------- gROOT->LoadClass("AliAODForwardMult", "libPWGLFforward2"); // --- Get analysis manager ---------------------------------------- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); if (!mgr) { Error("AddTaskForwardQA", "No analysis manager to connect to."); return NULL; } // --- Make the task and add it to the manager --------------------- AliForwardQATask* task = new AliForwardQATask("forwardQA"); mgr->AddTask(task); // --- Cuts -------------------------------------------------------- // Old style cuts // Double_t nXi = mc ? 2 : 2; //HHD test // Bool_t includeSigma = false; //true; // High cuts for sharing filter AliFMDMultCuts cHigh; cHigh.SetMPVFraction(0.7); cHigh.SetMultCuts(-1); cHigh.SetNXi(2); // cHigh.SetMultCuts(100); // Low cuts for sharing and density calculator AliFMDMultCuts cLow; // cLow.SetMultCuts(0.1, 0.1, 0.12, 0.1, 0.12); cLow.SetMultCuts(0.3, 0.3, 0.3, 0.3, 0.3); // Density cuts AliFMDMultCuts cDensity; // cDensity.SetMPVFraction(0.7); // cDensity.SetMultCuts(100); cDensity.SetMultCuts(0.3, 0.3, 0.3, 0.3, 0.3); // --- Set parameters on the algorithms ---------------------------- // Set the number of SPD tracklets for which we consider the event a // low flux event task->GetEventInspector().SetLowFluxCut(1000); // Set the maximum error on v_z [cm] task->GetEventInspector().SetMaxVzErr(0.2); // Disable use of code from 1st Physics task->GetEventInspector().SetUseFirstPhysicsVtx(kFALSE); // --- Set parameters on energy loss fitter ------------------------ // Set the eta axis to use - note, this overrides whatever is used // by the rest of the algorithms - but only for the energy fitter // algorithm. task->GetEnergyFitter().SetEtaAxis(200, -4, 6); // Set maximum energy loss to consider task->GetEnergyFitter().SetMaxE(15); // Set number of energy loss bins task->GetEnergyFitter().SetNEbins(500); // Set whether to use increasing bin sizes task->GetEnergyFitter().SetUseIncreasingBins(true); // Set whether to do fit the energy distributions task->GetEnergyFitter().SetDoFits(kTRUE); // Set whether to make the correction object task->GetEnergyFitter().SetDoMakeObject(kFALSE); // Set the low cut used for energy task->GetEnergyFitter().SetLowCut(0.4); // Set the number of bins to subtract from maximum of distributions // to get the lower bound of the fit range task->GetEnergyFitter().SetFitRangeBinWidth(4); // Set the minimum number of entries in the distribution before // trying to fit to the data task->GetEnergyFitter().SetMinEntries(1000); // --- Sharing filter ---------------------------------------------- // Enable use of angle corrected signals in the algorithm task->GetSharingFilter().SetUseAngleCorrectedSignals(true); // Disable use of angle corrected signals in the algorithm task->GetSharingFilter().SetZeroSharedHitsBelowThreshold(false); // Enable use of angle corrected signals in the algorithm task->GetSharingFilter().SetHCuts(cHigh); // Multiplicity cut // task->GetSharingFilter().GetHCuts().SetMultCuts(-1); // Set the number of xi's (width of landau peak) to stop at // task->GetSharingFilter().GetHCuts().SetNXi(nXi); // Set whether or not to include sigma in cut // task->GetSharingFilter().GetHCuts().SetIncludeSigma(includeSigma); // Lower cuts from object task->GetSharingFilter().SetLCuts(cLow); // Use simplified sharing algorithm task->GetSharingFilter().SetUseSimpleSharing(kTRUE); // --- Density calculator ------------------------------------------ // Set the maximum number of particle to try to reconstruct task->GetDensityCalculator().SetMaxParticles(3); // Wet whether to use poisson statistics to estimate N_ch task->GetDensityCalculator().SetUsePoisson(true); // How to lump for Poisson calculator - 64 strips, 4 regions task->GetDensityCalculator().SetLumping(64,4); // Set whether or not to include sigma in cut task->GetDensityCalculator().SetCuts(cDensity); // Set whether or not to use the phi acceptance // AliFMDDensityCalculator::kPhiNoCorrect // AliFMDDensityCalculator::kPhiCorrectNch // AliFMDDensityCalculator::kPhiCorrectELoss task->GetDensityCalculator(). SetUsePhiAcceptance(AliFMDDensityCalculator::kPhiNoCorrect); // --- Set limits on fits the energy ------------------------------- // Maximum relative error on parameters AliFMDCorrELossFit::ELossFit::fgMaxRelError = .12; // Least weight to use AliFMDCorrELossFit::ELossFit::fgLeastWeight = 1e-5; // Maximum value of reduced chi^2 AliFMDCorrELossFit::ELossFit::fgMaxChi2nu = 20; // --- Debug ------------------------------------------------------- // Set the overall debug level (1: some output, 3: a lot of output) task->SetDebug(3); // --- Make the output container and connect it -------------------- // TString outputfile = ; // outputfile += ":PWGLFforwardDnDeta"; // Form(":%s",pars->GetDndetaAnalysisName()); AliAnalysisDataContainer* histOut = mgr->CreateContainer("Forward", TList::Class(), AliAnalysisManager::kOutputContainer, AliAnalysisManager::GetCommonFileName()); AliAnalysisDataContainer *output = mgr->CreateContainer("ForwardResults", TList::Class(), AliAnalysisManager::kParamContainer, "trending.root"); // AliAnalysisManager::GetCommonFileName()); mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer()); mgr->ConnectOutput(task, 1, histOut); mgr->ConnectOutput(task, 2, output); return task; } // // EOF //