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   // --- Sharing filter ----------------------------------------------
  39   // Set the low cut used for sharing - overrides settings in eloss fits
  40   task->GetSharingFilter().SetLowCut(0.3);
  41   // Set the number of xi's (width of landau peak) to stop at
  42   task->GetSharingFilter().SetNXi(1);
  43
  44   // --- Density calculator
  45   // Set the maximum number of particle to try to reconstruct
  46   task->GetDensityCalculator().SetMaxParticles(3);
  47   // Wet whether to use poisson statistics to estimate N_ch
  48   task->GetDensityCalculator().SetUsePoisson(false);
  49   // Set the lower multiplicity cut.  Overrides setting in energy loss fits.
  50   task->GetDensityCalculator().SetMultCut(0.3); //was 0.3
  51
  52   // --- Corrector ---------------------------------------------------
  53   // Whether to use the secondary map correction
  54   task->GetCorrections().SetUseSecondaryMap(true);
  55   // Whether to use the vertex bias correction
  56   task->GetCorrections().SetUseVertexBias(false);
  57   // Whether to use the vertex bias correction
  58   task->GetCorrections().SetUseAcceptance(true);
  59   // Whether to use the merging efficiency correction
  60   task->GetCorrections().SetUseMergingEfficiency(false);
  61
  62   // --- Histogram Collector -----------------------------------------
  63   // Set the number of extra bins (beyond the secondary map border)
  64   task->GetHistCollector().SetNCutBins(2);
  65   // Set the correction cut, that is, when bins in the secondary map
  66   // is smaller than this, they are considered empty
  67   task->GetHistCollector().SetCorrectionCut(0.5);
  68   // How to calculate the value of overlapping bins.
  69   // Possible values are
  70   //    kStraightMean
  71   //    kStraightMeanNoZero
  72   //    kWeightedMean
  73   //    kLeastError
  74   task->GetHistCollector().SetMergeMethod(AliFMDHistCollector::kStraightMean);
  75   // How to find the fiducial area of the secondary maps
  76   // Possible values are
  77   //   kByCut    Only bins larger that cut are trusted
  78   //   kDistance Only bins that are more than half the size of it neighbors
  79   task->GetHistCollector().SetFiducialMethod(AliFMDHistCollector::kByCut);
  80
  81   // --- Debug -------------------------------------------------------
  82   // Set the overall debug level (1: some output, 3: a lot of output)
  83   task->SetDebug(0);
  84   // Set the debug level of a single algorithm
  85   // task->GetEventInspector().SetDebug(4);
  86
  87   // --- Set limits on fits the energy -------------------------------
  88   // Maximum relative error on parameters
  89   AliFMDCorrELossFit::ELossFit::fgMaxRelError = .12;
  90   // Least weight to use
  91   AliFMDCorrELossFit::ELossFit::fgLeastWeight = 1e-5;
  92   // Maximum value of reduced chi^2
  93   AliFMDCorrELossFit::ELossFit::fgMaxChi2nu   = 20;
  94 }
  95 //
  96 // EOF
  97 //