class AliAnalysisTaskLocalRho: public AliAnalysisTaskEmcalJet

 Default constructor

 Constructor

 destructor

 Init the analysis

 Initialize the anaysis

 create output objects

 Book a TH1F and connect it to the output container

 Book a TH2F and connect it to the output container

 Execute once for each event

 Get the vzero event plane

 Grab the TPC event plane. if parameter fExcludeLeadingJetsFromFit is larger than 0,
 strip in eta of width fExcludeLeadingJetsFromFit * GetJetContainer()->GetJetRadius() around the leading jet (before
 subtraction of rho) will be exluded from the event plane estimate

 Grab the combined vzero event plane

 Get the second order q-cumulant, a -999 return will be caught in the qa routine of CorrectRho

 Get the fourth order q-cumulant, a -999 return will be caught in the qa routine of CorrectRho

 Get the weighted n-th order q-vector, pass real and imaginary part as reference

 Get the weighted ij-th order autocorrelation correction

 Get multiplicity for unweighted q-cumulants. function QCnQnk should be called first

 Get multiplicity weights for the weighted two particle cumulant

 Get multiplicity weights for the weighted four particle cumulant

 Decides how to deal with the situation where c2 or c3 is negative
 Returns kTRUE depending on whether or not a modulated rho is used for the jet background

 Get rho' -> rho(phi)
 three routines are available, 1 and 2 can be used with or without pt weights
  [1] get vn from q-cumulants
      in case of cumulants, both cumulants and vn values are stored. in both cases, v2 and v3
      are expected. a check is performed to see if rho has no negative local minimum
      for full description, see Phys. Rev. C 83, 044913
      since the cn distribution has negative values, vn = sqrt(cn) can be imaginary sometimes
      in this case one can either roll back to the 'original' fixed rho, do a fit for vn or take use
      vn = - sqrt(|cn|) note that because of this, use of q-cumulants is not safe !
  [2] fitting a fourier expansion to the de/dphi distribution
      the fit can be done with either v2, v3 or a combination.
      in all cases, a cut can be made on the p-value of the chi-squared value of the fit
      and a check can be performed to see if rho has no negative local minimum
  [3] get v2 and v3 from user supplied histograms
      in this way, a fixed value of v2 and v3 is subtracted w.r.t. whichever event plane is requested

 Fill the analysis summary histrogram, saves all relevant analysis settigns

 Fill event plane histograms

 Terminate

 Set function to fit modulation

 constructors, destructor

 the cdf of the chisquare distribution is the normalized lower incomplete gamma function
inline
 setters - setup how to run

 setters - analysis details

 getters

 numerical evaluations