class AliAnalysisTaskPartonDisc: public AliAnalysisTaskSE

 Constructor

 read the cross sections
 and number of trials from pyxsec.root
 from AliAnalysisTaskJetSpectrum2

 Create histograms
 Called once

 Main loop
 Called for each event

 Draw result to the screen
 Called once at the end of the query

 Get the event type from the pythia headers

 Get the event type from the phojet header

 Check if the jet is inside abs(eta)<=fJetAcceptance

 Get the jet flavour, using the definition:
 The flavour will be that of the parton with the highest energy
 within an angular distance <= 0.3
 This method also keeps track of the mother of the parton giving the flavor

 Return R between the two jets or particles

 Calculate the number of charged particles necessary to
 add the given percentage of the jet energy (transverse energy)
 for the MC case

 Calculate the number of charged particles necessary to
 add the given percentage of the jet energy (transverse energy)
 for the AOD track case

 Allocate the static container with the given dimensions

 Initialize the static container with the given dimensions

 Sort the contents of the array
 From lower to higher value

 Loop over the contents and count how many tracks are necessary to recover
 the given percentage of the energy.
 If all tracks did not sum the required fraction, it returns 500

 Return kTrue if the mc track is inside the area covered by the cone of the jet

 Return kTrue if the track is inside the area covered by the cone of the jet

 Return kTRUE if the cone covers the vector
 from the primary vertex to this vertex

 Get the eta value

 Get the theta angle related to these coordinates

 Get the phi angle related to these coordinates

 Get delta-phi MC jet-track

 Get delta-eta MC jet-track

 Get delta-phi^2 MC jet-track

 Get delta-eta^2 MC jet-track

 Get delta-phi Track jet-track

 Get delta-eta Track jet-track

 Get delta-phi^2 Track jet-track

 Get delta-eta^2 Track jet-track

 get the number of events read out to create this AOD
 from the njets distribution from UA1
 This is to called in Notify and should provide the path to the AOD/ESD file
 code from AliAnalysisTaskJetSpectrum2

 Check if there are jet cones that overlap on the current event
 for UA1 and SISCone, based on the cone axis and the radius.
 There can be maximum 7.85 jet in the phi acceptance, with no overlap in
 the phi direction (2pi/0.8).
 Plus there can be up to two jets in contained in the eta acceptance,
 if they are centered in -0.5 and 0.5 respectively, per phi interval
 In total there can be up to 15.7 jets
 Check up to 16 jets, inside the acceptance,
 set the flags for up to those 8 jets.

 Reset the flags used for tagging jets from the branches
 Use before calling HasOverlapedCones(jetbranch)

 Calculate the number of charged particles above threshold
 inside this jet for the MC case
 the threshold is fCurrentJetMinPtNT90

 Calculate the number of charged particles necessary to
 add the given percentage of the jet energy (transverse energy)
 after background substraction, for the MC case

 Loop over the contents and count how many tracks are necessary to recover
 the given percentage of the energy.
 If all tracks did not sum the required fraction, it returns 500
 this saves the minimum pT used, during the nt90 recalculation
 after jet energy correction

 Calculate the number of charged particles above threshold
 inside this jet for the MC case
 the threshold is fCurrentJetMinPtNT90Recalc

 Calculate the number of charged particles
 above the threshold set by the NTX calculation
 for the AOD track case, the threshold is fCurrentJetMinPtNT90
 the fCurrentJetMinPtNT90 was set when calling the NTX method before

 Calculate the number of charged particles necessary to
 add the given percentage of the jet energy (transverse energy)
 after pT recalculation, for the AOD track case

 Calculate the number of charged particles
 above the threshold set by the NTX calculation
 for the AOD track case, the threshold is fCurrentJetMinPtNT90Recalc
 the fCurrentJetMinPtNT90Recalc was set when calling the NTX method before

 Check if these track coordinates are in an excluded area

 Calculate the V0 MIP multiplicity that should be excluded due
 to possible jet activity

 Calculate the V0 like multiplicity that should be excluded due
 to possible jet activity in MC events

 Comparison of Double_t values