2 * @file AddTaskForwarddNdeta.C
3 * @author Christian Holm Christensen <cholm@nbi.dk>
4 * @date Fri Jan 28 10:22:26 2011
6 * @brief Script to add a multiplicity task for the central
10 * @ingroup pwg2_forward_scripts_tasks
13 * Create the Forward @f$ dN/d\eta@f$ analysis task
15 * @param trig Trigger to use
16 * @param vzMin Smallest @f$ v_z@f$
17 * @param vzMax Biggest @f$ v_z@f$
18 * @param useCent Whether to use the centrality or not
19 * @param scheme Normalisation scheme
20 * @param cutEdges Whether to cut edges when rebinning
21 * @param mcanalysisfilename Take final MC corrections from this - if present
23 * @return Newly created and configured task
25 * @ingroup pwg2_forward_dndeta
28 AddTaskForwarddNdeta(const char* trig = "INEL",
31 Bool_t useCent = false,
32 const char* scheme = 0,
33 Bool_t cutEdges = false,
34 const char* mcanalysisfilename = "none")
36 // --- Load libraries ----------------------------------------------
37 gROOT->LoadClass("AliAODForwardMult", "libPWG2forward2");
39 // --- Analysis manager --------------------------------------------
40 AliAnalysisManager* mgr = AliAnalysisManager::GetAnalysisManager();
42 // --- Make our object ---------------------------------------------
43 AliForwarddNdetaTask* task = new AliForwarddNdetaTask("Forward");
44 //Set the filename of the corresponding MC analysis
45 task->SetMCFinalCorrFilename(mcanalysisfilename);
47 // Set the vertex range to use
48 task->SetVertexRange(vzMin, vzMax);
49 // Set the trigger mask to use (INEL,INEL>0,NSD)
50 task->SetTriggerMask(trig);
51 // Set the trigger efficiency
52 // task->SetTriggerEff(0.997535);
53 // Whether to cut edges when re-binning
54 task->SetCutEdges(cutEdges);
55 // Whether to correct for empty bins when projecting
56 task->SetCorrEmpty(true);
57 // Whether to use TH2::ProjectionX
58 task->SetUseROOTProjectX(false);
61 // kNone Normalise to accepted events
62 // kEventLevel Normalise to all events in selected range
63 // kAltEventLevel Normalise to all events in selected range
64 // kBackground Also correct for background triggers
65 // kShape Correct shape
67 // kNone, kEventLevel, and kAltEventLevel are mutually exclusive.
68 // If neither kEventLevel, nor kAltEventLevel is specified, then
69 // kNone is assumed. kBackground (when implemented) only makes
70 // sense with kEventLevel and kAltEventLevel. Furthermore, there
71 // are some constants that encode the common cases
73 // kFull = kEventLevel | kBackground | kShape
74 // kAltFull = kAltEventLevel | kBackground | kShape
77 task->SetNormalizationScheme(AliBasedNdetaTask::kFull);
78 if (scheme) task->SetNormalizationScheme(scheme);
79 // Set the centrality bins to use. These are mutually exclusive.
80 // Note, that a bin specified as a-b, covers the interval from a,
81 // inclusive to b exclusive. An upper bound of 100 is treated
82 // especially, and the upper bound is inclusive in that case .
84 Short_t bins[] = { 0, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100 };
85 task->SetCentralityAxis(11, bins);
89 // --- create containers for input/output --------------------------
90 AliAnalysisDataContainer *sums =
91 mgr->CreateContainer("ForwardSums", TList::Class(),
92 AliAnalysisManager::kOutputContainer,
93 AliAnalysisManager::GetCommonFileName());
94 AliAnalysisDataContainer *output =
95 mgr->CreateContainer("ForwardResults", TList::Class(),
96 AliAnalysisManager::kParamContainer,
97 AliAnalysisManager::GetCommonFileName());
99 // --- connect input/output ----------------------------------------
100 mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
101 mgr->ConnectOutput(task, 1, sums);
102 mgr->ConnectOutput(task, 2, output);
108 //________________________________________________________________________