Mostly documentation changes.
[u/mrichter/AliRoot.git] / PWG2 / FORWARD / analysis2 / ForwardAODConfig.C
1 /**
2  * @file   ForwardAODConfig.C
3  * @author Christian Holm Christensen <cholm@dalsgaard.hehi.nbi.dk>
4  * @date   Wed Mar 23 13:56:02 2011
5  * 
6  * @brief  
7  * 
8  * @ingroup pwg2_forward_scripts_tasks
9  * 
10  */
11 /**
12  * Configuration script for forward multiplicity task.  
13  *
14  * You can copy this to your working directory or to some other
15  * directory up-front in your ROOT macro path, and edit it to suit your
16  * needs.
17  * 
18  * @ingroup pwg2_forward_aod
19  */
20 void
21 ForwardAODConfig(AliForwardMultiplicityBase* task)
22 {
23   if (!task) return;
24
25   Info("ForwardAODConfig", "Setting up task %s (%p)", task->GetName(), task);
26
27   // --- General parameters ------------------------------------------
28   // Whether to enable low flux specific code 
29   task->SetEnableLowFlux(kFALSE);
30
31   // --- Event inspector ---------------------------------------------
32   // Set the number of SPD tracklets for which we consider the event a
33   // low flux event
34   task->GetEventInspector().SetLowFluxCut(1000); 
35   // Set the maximum error on v_z [cm]
36   task->GetEventInspector().SetMaxVzErr(0.2);
37   
38   // --- Energy Loss Fitter ------------------------------------------
39   // Set the eta axis to use - note, this overrides whatever is used
40   // by the rest of the algorithms - but only for the energy fitter
41   // algorithm. 
42   task->GetEnergyFitter().SetEtaAxis(200, -4, 6);
43   // Set maximum energy loss to consider 
44   task->GetEnergyFitter().SetMaxE(10); 
45   // Set number of energy loss bins 
46   task->GetEnergyFitter().SetNEbins(300);
47   // Set whether to use increasing bin sizes 
48   task->GetEnergyFitter().SetUseIncreasingBins(true);
49   // Set whether to do fit the energy distributions 
50   task->GetEnergyFitter().SetDoFits(kFALSE);
51   // Set whether to make the correction object 
52   task->GetEnergyFitter().SetDoMakeObject(kFALSE);
53   // Set the low cut used for energy
54   task->GetEnergyFitter().SetLowCut(0.4);
55   // Set the number of bins to subtract from maximum of distributions
56   // to get the lower bound of the fit range
57   task->GetEnergyFitter().SetFitRangeBinWidth(4);
58   // Set the maximum number of landaus to try to fit (max 5)
59   task->GetEnergyFitter().SetNParticles(5);
60   // Set the minimum number of entries in the distribution before
61   // trying to fit to the data
62   task->GetEnergyFitter().SetMinEntries(1000);
63
64   // --- Sharing filter ----------------------------------------------
65   // Set the low cut used for sharing - overrides settings in eloss fits
66   task->GetSharingFilter().SetLowCut(0.3);
67   // Set the number of xi's (width of landau peak) to stop at 
68   task->GetSharingFilter().SetNXi(1);
69
70   // --- Density calculator 
71   // Set the maximum number of particle to try to reconstruct 
72   task->GetDensityCalculator().SetMaxParticles(3);
73   // Wet whether to use poisson statistics to estimate N_ch
74   task->GetDensityCalculator().SetUsePoisson(false);
75   // Set the lower multiplicity cut.  Overrides setting in energy loss fits.
76   task->GetDensityCalculator().SetMultCut(0.3); //was 0.3
77
78   // --- Corrector ---------------------------------------------------
79   // Whether to use the secondary map correction
80   task->GetCorrections().SetUseSecondaryMap(true);
81   // Whether to use the vertex bias correction
82   task->GetCorrections().SetUseVertexBias(false);
83   // Whether to use the vertex bias correction
84   task->GetCorrections().SetUseAcceptance(true);
85   // Whether to use the merging efficiency correction 
86   task->GetCorrections().SetUseMergingEfficiency(false);
87
88   // --- Histogram Collector -----------------------------------------
89   // Set the number of extra bins (beyond the secondary map border) 
90   task->GetHistCollector().SetNCutBins(2);
91   // Set the correction cut, that is, when bins in the secondary map 
92   // is smaller than this, they are considered empty 
93   task->GetHistCollector().SetCorrectionCut(0.5);
94   // How to calculate the value of overlapping bins. 
95   // Possible values are 
96   //    kStraightMean 
97   //    kStraightMeanNoZero 
98   //    kWeightedMean 
99   //    kLeastError 
100   task->GetHistCollector().SetMergeMethod(AliFMDHistCollector::kStraightMean);
101   // How to find the fiducial area of the secondary maps 
102   // Possible values are 
103   //   kByCut    Only bins larger that cut are trusted 
104   //   kDistance Only bins that are more than half the size of it neighbors
105   task->GetHistCollector().SetFiducialMethod(AliFMDHistCollector::kByCut);
106
107   // --- Debug -------------------------------------------------------
108   // Set the overall debug level (1: some output, 3: a lot of output)
109   task->SetDebug(0);
110   // Set the debug level of a single algorithm 
111   // task->GetEventInspector().SetDebug(4);
112
113   // --- Set limits on fits the energy -------------------------------
114   // Maximum relative error on parameters 
115   AliFMDCorrELossFit::ELossFit::fgMaxRelError = .12;
116   // Least weight to use 
117   AliFMDCorrELossFit::ELossFit::fgLeastWeight = 1e-5;
118   // Maximum value of reduced chi^2 
119   AliFMDCorrELossFit::ELossFit::fgMaxChi2nu   = 20;
120 }
121 //
122 // EOF
123 //