ROOT logo
AliRoot » PWGHF » VERTEXINGHF » AliVertexingHFUtils

class AliVertexingHFUtils: public TObject


Class with functions useful for different D2H analyses
- event plane resolution
- <pt> calculation with side band subtraction
- tracklet multiplicity calculation
Origin: F.Prino, Torino, prino@to.infn.it


Function Members (Methods)

public:
AliVertexingHFUtils()
AliVertexingHFUtils(Int_t k)
AliVertexingHFUtils(const AliVertexingHFUtils&)
virtual~AliVertexingHFUtils()
voidTObject::AbstractMethod(const char* method) const
virtual voidTObject::AppendPad(Option_t* option = "")
static voidAveragePt(Float_t& averagePt, Float_t& errorPt, Float_t ptmin, Float_t ptmax, TH2F* hMassD, Float_t massFromFit, Float_t sigmaFromFit, TF1* funcB2, Float_t sigmaRangeForSig = 2.5, Float_t sigmaRangeForBkg = 4.5, Float_t minMass = 0., Float_t maxMass = 3., Int_t rebin = 1)
virtual voidTObject::Browse(TBrowser* b)
static Int_tCheckD0Decay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckD0Decay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
static Int_tCheckDplusDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckDplusDecay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
static Int_tCheckDplusK0spiDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckDplusKKpiDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckDsDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckDsDecay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
static Int_tCheckDsK0sKDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckDstarDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckDstarDecay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
static Int_tCheckLcpKpiDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckLcV0bachelorDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
static Int_tCheckOrigin(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Bool_t searchUpToQuark = kTRUE)
static Int_tCheckOrigin(AliStack* stack, TParticle* mcPart, Bool_t searchUpToQuark = kTRUE)
static Bool_tCheckT0TriggerFired(AliAODEvent* aodEv)
static TClass*Class()
virtual const char*TObject::ClassName() const
virtual voidTObject::Clear(Option_t* = "")
virtual TObject*TObject::Clone(const char* newname = "") const
virtual Int_tTObject::Compare(const TObject* obj) const
static voidComputeSignificance(Double_t signal, Double_t errsignal, Double_t background, Double_t errbackground, Double_t& significance, Double_t& errsignificance)
virtual voidTObject::Copy(TObject& object) const
virtual voidTObject::Delete(Option_t* option = "")MENU
virtual Int_tTObject::DistancetoPrimitive(Int_t px, Int_t py)
virtual voidTObject::Draw(Option_t* option = "")
virtual voidTObject::DrawClass() constMENU
virtual TObject*TObject::DrawClone(Option_t* option = "") constMENU
virtual voidTObject::Dump() constMENU
virtual voidTObject::Error(const char* method, const char* msgfmt) const
virtual voidTObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual voidTObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual voidTObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
virtual voidTObject::Fatal(const char* method, const char* msgfmt) const
Double_tFindChi() const
static Double_tFindChi(Double_t res, Int_t k = 1)
virtual TObject*TObject::FindObject(const char* name) const
virtual TObject*TObject::FindObject(const TObject* obj) const
static Double_tGetCorrectedNtracklets(TProfile* estimatorAvg, Double_t uncorrectedNacc, Double_t vtxZ, Double_t refMult)
virtual Option_t*TObject::GetDrawOption() const
static Long_tTObject::GetDtorOnly()
Double_tGetFullEvResol() const
static Double_tGetFullEvResol(Double_t resSub, Int_t k = 1)
static Double_tGetFullEvResol(const TH1F* hSubEvCorr, Int_t k = 1)
static Double_tGetFullEvResolHighLim(const TH1F* hSubEvCorr, Int_t k = 1)
static Double_tGetFullEvResolLowLim(const TH1F* hSubEvCorr, Int_t k = 1)
static Int_tGetGeneratedMultiplicityInEtaRange(TClonesArray* arrayMC, Double_t mineta, Double_t maxeta)
static Int_tGetGeneratedPhysicalPrimariesInEtaRange(TClonesArray* arrayMC, Double_t mineta, Double_t maxeta)
static Int_tGetGeneratedPrimariesInEtaRange(TClonesArray* arrayMC, Double_t mineta, Double_t maxeta)
static TStringGetGenerator(Int_t label, AliAODMCHeader* header)
virtual const char*TObject::GetIconName() const
Int_tGetK() const
virtual const char*TObject::GetName() const
Int_tGetNumberOfTrackletsInEtaRange(AliAODEvent* ev) const
static Int_tGetNumberOfTrackletsInEtaRange(AliAODEvent* ev, Double_t mineta, Double_t maxeta)
virtual char*TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_tTObject::GetObjectStat()
virtual Option_t*TObject::GetOption() const
Double_tGetSubEvResol() const
static Double_tGetSubEvResol(const TH1F* hSubEvCorr)
static Double_tGetSubEvResolHighLim(const TH1F* hSubEvCorr)
static Double_tGetSubEvResolLowLim(const TH1F* hSubEvCorr)
virtual const char*TObject::GetTitle() const
voidGetTrackPrimaryGenerator(AliAODTrack* track, AliAODMCHeader* header, TClonesArray* arrayMC, TString& nameGen)
static Double_tGetTrueImpactParameterDplus(AliAODMCHeader* mcHeader, TClonesArray* arrayMC, AliAODMCParticle* partDp)
static Double_tGetTrueImpactParameterDzero(AliAODMCHeader* mcHeader, TClonesArray* arrayMC, AliAODMCParticle* partDp)
virtual UInt_tTObject::GetUniqueID() const
static Double_tGetVZEROAEqualizedMultiplicity(AliAODEvent* ev)
static Double_tGetVZEROCEqualizedMultiplicity(AliAODEvent* ev)
virtual Bool_tTObject::HandleTimer(TTimer* timer)
Bool_tHasCascadeCandidateAnyDaughInjected(AliAODRecoCascadeHF* cand, AliAODMCHeader* header, TClonesArray* arrayMC)
virtual ULong_tTObject::Hash() const
virtual voidTObject::Info(const char* method, const char* msgfmt) const
virtual Bool_tTObject::InheritsFrom(const char* classname) const
virtual Bool_tTObject::InheritsFrom(const TClass* cl) const
virtual voidTObject::Inspect() constMENU
voidTObject::InvertBit(UInt_t f)
virtual TClass*IsA() const
Bool_tIsCandidateInjected(AliAODRecoDecayHF* cand, AliAODMCHeader* header, TClonesArray* arrayMC)
virtual Bool_tTObject::IsEqual(const TObject* obj) const
virtual Bool_tTObject::IsFolder() const
Bool_tTObject::IsOnHeap() const
virtual Bool_tTObject::IsSortable() const
Bool_tIsTrackInjected(AliAODTrack* track, AliAODMCHeader* header, TClonesArray* arrayMC)
Bool_tTObject::IsZombie() const
virtual voidTObject::ls(Option_t* option = "") const
voidTObject::MayNotUse(const char* method) const
virtual Bool_tTObject::Notify()
voidTObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const
static voidTObject::operator delete(void* ptr)
static voidTObject::operator delete(void* ptr, void* vp)
static voidTObject::operator delete[](void* ptr)
static voidTObject::operator delete[](void* ptr, void* vp)
void*TObject::operator new(size_t sz)
void*TObject::operator new(size_t sz, void* vp)
void*TObject::operator new[](size_t sz)
void*TObject::operator new[](size_t sz, void* vp)
AliVertexingHFUtils&operator=(const AliVertexingHFUtils&)
virtual voidTObject::Paint(Option_t* option = "")
Double_tPol(Double_t x) const
static Double_tPol(Double_t x, Int_t k)
virtual voidTObject::Pop()
virtual voidTObject::Print(Option_t* option = "") const
virtual Int_tTObject::Read(const char* name)
virtual voidTObject::RecursiveRemove(TObject* obj)
voidTObject::ResetBit(UInt_t f)
static Double_tResolK1(Double_t x)
virtual voidTObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU
virtual voidTObject::SavePrimitive(basic_ostream<char,char_traits<char> >& out, Option_t* option = "")
voidTObject::SetBit(UInt_t f)
voidTObject::SetBit(UInt_t f, Bool_t set)
virtual voidTObject::SetDrawOption(Option_t* option = "")MENU
static voidTObject::SetDtorOnly(void* obj)
voidSetEtaRangeForTracklets(Double_t mineta, Double_t maxeta)
voidSetK(Int_t k)
static voidTObject::SetObjectStat(Bool_t stat)
voidSetSubEventHisto(const TH1F* hSub)
voidSetSubEvResol(Double_t res)
virtual voidTObject::SetUniqueID(UInt_t uid)
virtual voidShowMembers(TMemberInspector&)
virtual voidStreamer(TBuffer&)
voidStreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
virtual voidTObject::SysError(const char* method, const char* msgfmt) const
Bool_tTObject::TestBit(UInt_t f) const
Int_tTObject::TestBits(UInt_t f) const
virtual voidTObject::UseCurrentStyle()
virtual voidTObject::Warning(const char* method, const char* msgfmt) const
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const
protected:
virtual voidTObject::DoError(int level, const char* location, const char* fmt, va_list va) const
voidTObject::MakeZombie()

Data Members

private:
Int_tfKratio of measured harmonic to event plane harmonic
Double_tfMaxEtaForTrackletsmin eta for counting tracklets
Double_tfMinEtaForTrackletsmin eta for counting tracklets
Double_tfSubRessub-event resolution = sqrt(<cos[n(phiA-phiB)] >)

Class Charts

Inheritance Chart:
TObject
AliVertexingHFUtils

Function documentation

AliVertexingHFUtils()
 Default contructor
AliVertexingHFUtils(Int_t k)
 Standard constructor
void ComputeSignificance(Double_t signal, Double_t errsignal, Double_t background, Double_t errbackground, Double_t& significance, Double_t& errsignificance)
 calculate significance from S, B and errors
Double_t Pol(Double_t x, Int_t k)
 compute chi from polynomial approximation
Double_t FindChi(Double_t res, Int_t k = 1)
 compute chi variable (=v2*sqrt(N)) from external values
Double_t GetFullEvResol(Double_t resSub, Int_t k = 1)
 computes event plane resolution starting from sub event resolution
Double_t GetFullEvResol(const TH1F* hSubEvCorr, Int_t k = 1)
 computes event plane resolution starting from sub event correlation histogram
Double_t GetFullEvResolLowLim(const TH1F* hSubEvCorr, Int_t k = 1)
 computes low limit event plane resolution starting from sub event correlation histogram
Double_t GetFullEvResolHighLim(const TH1F* hSubEvCorr, Int_t k = 1)
 computes high limit event plane resolution starting from sub event correlation histogram
Int_t GetNumberOfTrackletsInEtaRange(AliAODEvent* ev, Double_t mineta, Double_t maxeta)
 counts tracklets in given eta range
Int_t GetGeneratedMultiplicityInEtaRange(TClonesArray* arrayMC, Double_t mineta, Double_t maxeta)
 counts generated particles in fgiven eta range
Int_t GetGeneratedPrimariesInEtaRange(TClonesArray* arrayMC, Double_t mineta, Double_t maxeta)
 counts generated primary particles in given eta range
Int_t GetGeneratedPhysicalPrimariesInEtaRange(TClonesArray* arrayMC, Double_t mineta, Double_t maxeta)
 counts generated primary particles in given eta range
Double_t GetVZEROAEqualizedMultiplicity(AliAODEvent* ev)
 Method to get VZERO-A equalized multiplicty as done in AliCentralitySelectionTask
  getting the equalized VZERO factors from the tender or AOD
 http://git.cern.ch/pubweb/AliRoot.git/blob/HEAD:/ANALYSIS/AliCentralitySelectionTask.cxx#l1345
Double_t GetVZEROCEqualizedMultiplicity(AliAODEvent* ev)
 Method to get VZERO-C equalized multiplicty as done in AliCentralitySelectionTask
  getting the equalized VZERO factors from the tender or AOD
 http://git.cern.ch/pubweb/AliRoot.git/blob/HEAD:/ANALYSIS/AliCentralitySelectionTask.cxx#l1345
void AveragePt(Float_t& averagePt, Float_t& errorPt, Float_t ptmin, Float_t ptmax, TH2F* hMassD, Float_t massFromFit, Float_t sigmaFromFit, TF1* funcB2, Float_t sigmaRangeForSig = 2.5, Float_t sigmaRangeForBkg = 4.5, Float_t minMass = 0., Float_t maxMass = 3., Int_t rebin = 1)
Bool_t CheckT0TriggerFired(AliAODEvent* aodEv)
 check if T0VTX trigger was fired, based on a workaround suggested by Alla
Double_t GetTrueImpactParameterDzero(AliAODMCHeader* mcHeader, TClonesArray* arrayMC, AliAODMCParticle* partDp)
 true impact parameter calculation for Dzero
Double_t GetTrueImpactParameterDplus(AliAODMCHeader* mcHeader, TClonesArray* arrayMC, AliAODMCParticle* partDp)
 true impact parameter calculation for Dplus
Double_t GetCorrectedNtracklets(TProfile* estimatorAvg, Double_t uncorrectedNacc, Double_t vtxZ, Double_t refMult)
 Correct the number of accepted tracklets based on a TProfile Hist


TString GetGenerator(Int_t label, AliAODMCHeader* header)
 get the name of the generator that produced a given particle
void GetTrackPrimaryGenerator(AliAODTrack* track, AliAODMCHeader* header, TClonesArray* arrayMC, TString& nameGen)
Bool_t IsTrackInjected(AliAODTrack* track, AliAODMCHeader* header, TClonesArray* arrayMC)
 method to check if a track comes from the signal event or from the underlying Hijing event
Bool_t IsCandidateInjected(AliAODRecoDecayHF* cand, AliAODMCHeader* header, TClonesArray* arrayMC)
 method to check if a D meson candidate comes from the signal event or from the underlying Hijing event
Bool_t HasCascadeCandidateAnyDaughInjected(AliAODRecoCascadeHF* cand, AliAODMCHeader* header, TClonesArray* arrayMC)
 method to check if a cascade candidate comes from the signal event or from the underlying Hijing event
Int_t CheckOrigin(AliStack* stack, TParticle* mcPart, Bool_t searchUpToQuark = kTRUE)
 checking whether the mother of the particles come from a charm or a bottom quark
Int_t CheckOrigin(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Bool_t searchUpToQuark = kTRUE)
 checking whether the mother of the particles come from a charm or a bottom quark
Int_t CheckD0Decay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the D0 decay channel. Returns 1 for the D0->Kpi case, 2 for the D0->Kpipipi case, -1 in other cases
Int_t CheckD0Decay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
 Checks the D0 decay channel. Returns 1 for the D0->Kpi case, 2 for the D0->Kpipipi case, -1 in other cases
Int_t CheckDplusDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Dplus decay channel. Returns 1 for the non-resonant case, 2 for the resonant case, -1 in other cases
Int_t CheckDplusDecay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
 Checks the Dplus decay channel. Returns 1 for the non-resonant case, 2 for the resonant case, -1 in other cases
Int_t CheckDplusKKpiDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Ds decay channel. Returns 1 for Ds->phipi->KKpi, 2 for Ds->K0*K->KKpi, 3 for the non-resonant case
Int_t CheckDplusK0spiDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Dplus->V0+pion decay channel. Returns 1 if success, -1 otherwise
Int_t CheckDsDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Ds decay channel. Returns 1 for Ds->phipi->KKpi, 2 for Ds->K0*K->KKpi, 3 for the non-resonant case
Int_t CheckDsK0sKDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Ds->K0s+S decay channel. Returns 1 in case of success, otherwise -1
Int_t CheckDsDecay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
 Checks the Ds decay channel. Returns 1 for Ds->phipi->KKpi, 2 for Ds->K0*K->KKpi, 3 for the non-resonant case
Int_t CheckDstarDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Dstar decay channel. Returns 1 for D*->D0pi->Kpipi, -1 in other cases
Int_t CheckDstarDecay(TClonesArray* arrayMC, AliAODMCParticle* mcPart, Int_t* arrayDauLab)
 Checks the Dstar decay channel. Returns 1 for D*->D0pi->Kpipi, -1 in other cases
Int_t CheckLcpKpiDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Lc->pKpi decay channel. Returns 1 for non-resonant decays and 2, 3 or 4 for resonant ones, -1 in other cases
Int_t CheckLcV0bachelorDecay(AliStack* stack, Int_t label, Int_t* arrayDauLab)
 Checks the Lc->V0+bachelor decay channel. Returns 1 for pK0s, 2 for piLambda, 3 for pK0l -1 in other cases
AliVertexingHFUtils()
virtual ~AliVertexingHFUtils()
{}
void SetK(Int_t k)
 Functions for Event plane resolution
{fK=k;}
void SetSubEvResol(Double_t res)
{fSubRes=res;}
void SetSubEventHisto(const TH1F* hSub)
Int_t GetK() const
{return fK;}
Double_t GetSubEvResol() const
{return fSubRes;}
Double_t Pol(Double_t x) const
{return Pol(x,fK);}
Double_t FindChi() const
{return FindChi(fSubRes,fK);}
Double_t GetFullEvResol() const
Double_t ResolK1(Double_t x)
Double_t GetSubEvResol(const TH1F* hSubEvCorr)
Double_t GetSubEvResolLowLim(const TH1F* hSubEvCorr)
Double_t GetSubEvResolHighLim(const TH1F* hSubEvCorr)
void SetEtaRangeForTracklets(Double_t mineta, Double_t maxeta)
 Functions for tracklet multiplcity calculation
Int_t GetNumberOfTrackletsInEtaRange(AliAODEvent* ev, Double_t mineta, Double_t maxeta)