class AliAnalysisVertexingHF: public TNamed
Implementation of the heavy-flavour vertexing analysis class
Candidates are stored in the AOD as objects deriving from AliAODRecoDecay.
To be used as a task of AliAnalysisManager by means of the interface
class AliAnalysisTaskSEVertexingHF.
An example of usage in the macro AliAnalysisTaskSEVertexingHFTest.C.
Contact: andrea.dainese@pd.infn.it
Contributors: E.Bruna, G.E.Bruno, A.Dainese, C.Di Gliglio,
F.Prino, R.Romita, X.M.Zhang
Function Members (Methods)
public:
protected:
| virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| void | TObject::MakeZombie() |
private:
Data Members
public:
| enum { | kBitDispl | |
| kBitSoftPi | ||
| kBit3Prong | ||
| kBitPionCompat | ||
| kBitKaonCompat | ||
| kBitProtonCompat | ||
| }; | ||
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
protected:
| TString | TNamed::fName | object identifier |
| TString | TNamed::fTitle | object title |
private:
| Bool_t | f3Prong | D+,Ds,Lc |
| Bool_t | f4Prong | D0->Kpipipi |
| Int_t* | fAODMap | [fAODMapSize] map between index and ID for AOD tracks |
| Int_t | fAODMapSize | size of fAODMap |
| Double_t | fBzkG | z componenent of field in kG |
| Bool_t | fCascades | cascades, Lc --> v0+track |
| AliRDHFCutsD0toKpi* | fCutsD0toKpi | D0->Kpi cuts |
| AliRDHFCutsD0toKpipipi* | fCutsD0toKpipipi | D0->Kpipipi cuts |
| AliRDHFCutsDStartoKpipi* | fCutsDStartoKpipi | Dstar->D0pi cuts |
| AliRDHFCutsDplustoKpipi* | fCutsDplustoKpipi | D+->Kpipi cuts |
| AliRDHFCutsDstoKKpi* | fCutsDstoKKpi | Ds->KKpi cuts |
| AliRDHFCutsJpsitoee* | fCutsJpsitoee | J/psi->ee cuts |
| AliRDHFCutsLctoV0* | fCutsLctoV0 | Lc --> v0 + bachelor cuts |
| AliRDHFCutsLctopKpi* | fCutsLctopKpi | Lc->pKpi cuts |
| Bool_t | fD0toKpi | D0->Kpi |
| Bool_t | fDstar | D*->D0pi |
| Bool_t | fFindVertexForCascades | reconstruct a secondary vertex or assume it's from the primary vertex |
| Bool_t | fFindVertexForDstar | reconstruct a secondary vertex or assume it's from the primary vertex |
| Bool_t | fInputAOD | input from AOD (kTRUE) or ESD (kFALSE) |
| Bool_t | fJPSItoEle | Jpsi->ee |
| Bool_t | fLikeSign | Like-sign pairs |
| Bool_t | fLikeSign3prong | Like-sign triplets |
| TList* | fListOfCuts | pointer to list of cuts for output file |
| AliAODRecoDecay* | fMassCalc2 | for 2 prong |
| AliAODRecoDecay* | fMassCalc3 | for 3 prong |
| AliAODRecoDecay* | fMassCalc4 | for 4 prong |
| Bool_t | fMassCutBeforeVertexing | to go faster in PbPb |
| Double_t | fMassDplus | |
| Double_t | fMassDs | |
| Double_t | fMassDstar | |
| Double_t | fMassDzero | |
| Double_t | fMassJpsi | |
| Double_t | fMassLambdaC | |
| Float_t | fMaxCentPercentileForTightCuts | max. centrality percentile for using tight cuts |
| Double_t | fMaxMomForTPCPid | upper momentum limit to apply TPC PID |
| Bool_t | fMixEvent | event mixing |
| Bool_t | fOKInvMassD0 | pair fullfilling D0 inv mass selection |
| Bool_t | fOKInvMassD0to4p | 4tracks fullfilling D0 inv mass selection |
| Bool_t | fOKInvMassDplus | triplet fullfilling D+ inv mass selection |
| Bool_t | fOKInvMassDs | triplet fullfilling Ds inv mass selection |
| Bool_t | fOKInvMassDstar | combination fullfilling D* inv mass selection |
| Bool_t | fOKInvMassJpsi | pair fullfilling Jpsi inv mass selection |
| Bool_t | fOKInvMassLc | triplet fullfilling Lc inv mass selection |
| Bool_t | fOKInvMassLctoV0 | triplet fullfilling Lc inv mass selection |
| AliPIDResponse* | fPidResponse | PID response |
| Bool_t | fRecoPrimVtxSkippingTrks | flag for primary vertex reco on the fly |
| Bool_t | fRmTrksFromPrimVtx | flag for fast removal of daughters from |
| Bool_t | fSecVtxWithKF | if kTRUE use KF vertexer, else AliVertexerTracks |
| AliAnalysisFilter* | fTrackFilter | Track Filter for displaced vertices |
| AliAnalysisFilter* | fTrackFilter2prongCentral | Track Filter for displaced vertices in PbPb central events (tighter cuts) for 2 prong (D0->Kpi) |
| AliAnalysisFilter* | fTrackFilter3prongCentral | Track Filter for displaced vertices in PbPb central events (tighter cuts) for 3 prong (D+, Ds, Lc) |
| AliAnalysisFilter* | fTrackFilterSoftPi | Track Filter for D* soft pion |
| Bool_t | fUseKaonPIDfor3Prong | Kaon PID usage for 3 prongs |
| Bool_t | fUseKaonPIDforDs | Kaon PID usage for Ds |
| Int_t | fUsePIDforLc | PID for Lambdac: 0=no, 1=proton, 2=p and pi |
| Bool_t | fUsePIDforLc2V0 | PID for Lambdac 2 V0: 0=no, 1=proton, |
| Bool_t | fUseTOFPID | switch use/not use TOF PID |
| Bool_t | fUseTPCPID | switch use/not use TPC PID |
| Bool_t | fUseTPCPIDOnlyIfNoTOF | use TPC PID only for tracks that without TOF |
| Int_t | fV0TypeForCascadeVertex | Select which V0 type we want to use for the cascas |
| AliESDVertex* | fV1 | primary vertex |
| AliVertexerTracks* | fVertexerTracks | vertexer, to compute secondary vertices |
| Int_t | fnSeleTrksTotal | |
| Double_t | fnSigmaTOFKaonHi | High cut value on n. of sigmas for K TOF PID |
| Double_t | fnSigmaTOFKaonLow | Low cut value on n. of sigmas for K TOF PID |
| Double_t | fnSigmaTOFPionHi | High cut value on n. of sigmas for pi TOF PID |
| Double_t | fnSigmaTOFPionLow | Low cut value on n. of sigmas for pi TOF PID |
| Double_t | fnSigmaTOFProtonHi | High cut value on n. of sigmas for p TOF PID |
| Double_t | fnSigmaTOFProtonLow | Low cut value on n. of sigmas for p TOF PID |
| Double_t | fnSigmaTPCKaonHi | High cut value on n. of sigmas for K TPC PID |
| Double_t | fnSigmaTPCKaonLow | Low cut value on n. of sigmas for K TPC PID |
| Double_t | fnSigmaTPCPionHi | High cut value on n. of sigmas for pi TPC PID |
| Double_t | fnSigmaTPCPionLow | Low cut value on n. of sigmas for pi TPC PID |
| Double_t | fnSigmaTPCProtonHi | High cut value on n. of sigmas for p TPC PID |
| Double_t | fnSigmaTPCProtonLow | Low cut value on n. of sigmas for p TPC PID |
| Int_t | fnTrksTotal |
Class Charts
| Inheritance Chart: | ||||||||
|
Function documentation
void FindCandidates(AliVEvent* event, TClonesArray* aodVerticesHFTClArr, TClonesArray* aodD0toKpiTClArr, TClonesArray* aodJPSItoEleTClArr, TClonesArray* aodCharm3ProngTClArr, TClonesArray* aodCharm4ProngTClArr, TClonesArray* aodDstarTClArr, TClonesArray* aodCascadesTClArr, TClonesArray* aodLikeSign2ProngTClArr, TClonesArray* aodLikeSign3ProngTClArr)
Find heavy-flavour vertex candidates
Input: ESD or AOD
Output: AOD (additional branches added)
AliCodeTimerAuto("",0);
void AddRefs(AliAODVertex* v, AliAODRecoDecayHF* rd, const AliVEvent* event, const TObjArray* trkArray) const
Add the AOD tracks as daughters of the vertex (TRef)
Also add the references to the primary vertex and to the cuts
AliCodeTimerAuto("",0);
void AddDaughterRefs(AliAODVertex* v, const AliVEvent* event, const TObjArray* trkArray) const
Add the AOD tracks as daughters of the vertex (TRef)
AliCodeTimerAuto("",0);
void FixReferences(AliAODEvent* aod)
Checks that the references to the daughter tracks are properly
assigned and reassigns them if needed
AliCodeTimerAuto("",0);
AliAODRecoCascadeHF* MakeCascade(TObjArray* twoTrackArray, AliVEvent* event, AliAODVertex* secVert, AliAODRecoDecayHF2Prong* rd2Prong, Double_t dca, Bool_t& okDstar)
Make the cascade as a 2Prong decay and check if it passes Dstar
reconstruction cuts
AliCodeTimerAuto("",0);
AliAODRecoCascadeHF* MakeCascade(TObjArray* twoTrackArray, AliVEvent* event, AliAODVertex* secVert, AliAODv0* v0, Double_t dca, Bool_t& okCascades)
Make the cascade as a 2Prong decay and check if it passes
cascades reconstruction cuts
AliCodeTimerAuto("",0);
AliAODRecoDecayHF2Prong * Make2Prong(TObjArray* twoTrackArray1, AliVEvent* event, AliAODVertex* secVert, Double_t dcap1n1, Bool_t& okD0, Bool_t& okJPSI, Bool_t& okD0fromDstar)
Make 2Prong candidates and check if they pass D0toKpi or BtoJPSI
reconstruction cuts
G.E.Bruno (J/psi), A.Dainese (D0->Kpi)
AliCodeTimerAuto("",0);
AliAODRecoDecayHF3Prong* Make3Prong(TObjArray* threeTrackArray, AliVEvent* event, AliAODVertex* secVert, Double_t dispersion, const AliAODVertex* vertexp1n1, const AliAODVertex* vertexp2n1, Double_t dcap1n1, Double_t dcap2n1, Double_t dcap1p2, Bool_t useForLc, Bool_t useForDs, Bool_t& ok3Prong)
Make 3Prong candidates and check if they pass Dplus or Ds or Lambdac reconstruction cuts E.Bruna, F.Prino
AliAODRecoDecayHF4Prong* Make4Prong(TObjArray* fourTrackArray, AliVEvent* event, AliAODVertex* secVert, const AliAODVertex* vertexp1n1, const AliAODVertex* vertexp1n1p2, Double_t dcap1n1, Double_t dcap1n2, Double_t dcap2n1, Double_t dcap2n2, Bool_t& ok4Prong)
Make 4Prong candidates and check if they pass D0toKpipipi
reconstruction cuts
G.E.Bruno, R.Romita
AliCodeTimerAuto("",0);
AliAODVertex* PrimaryVertex(const TObjArray* trkArray = 0x0, AliVEvent* event = 0x0) const
Returns primary vertex to be used for this candidate
AliCodeTimerAuto("",0);
AliAODVertex* ReconstructSecondaryVertex(TObjArray* trkArray, Double_t& dispersion, Bool_t useTRefArray = kTRUE) const
Secondary vertex reconstruction with AliVertexerTracks or AliKFParticle AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPt3prong(TObjArray* trkArray)
Invariant mass cut on tracks
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPt4prong(TObjArray* trkArray)
Invariant mass cut on tracks
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPtDstarD0pi(TObjArray* trkArray)
Invariant mass cut on tracks
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPtD0Kpi(Double_t* px, Double_t* py, Double_t* pz)
Check invariant mass cut and pt candidate cut
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPtJpsiee(Double_t* px, Double_t* py, Double_t* pz)
Check invariant mass cut and pt candidate cut
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPt3prong(Double_t* px, Double_t* py, Double_t* pz, Int_t pidLcStatus = 3)
Check invariant mass cut and pt candidate cut
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPtDstarD0pi(Double_t* px, Double_t* py, Double_t* pz)
Check invariant mass cut and pt candidate cut
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPt4prong(Double_t* px, Double_t* py, Double_t* pz)
Check invariant mass cut and pt candidate cut
AliCodeTimerAuto("",0);
Bool_t SelectInvMassAndPtCascade(Double_t* px, Double_t* py, Double_t* pz)
Check invariant mass cut and pt candidate cut
AliCodeTimerAuto("",0);
void SelectTracksAndCopyVertex(const AliVEvent* event, Int_t trkEntries, TObjArray& seleTrksArray, TObjArray& tracksAtVertex, Int_t& nSeleTrks, UChar_t* seleFlags, Int_t* evtNumber)
Apply single-track preselection. Fill a TObjArray with selected tracks (for displaced vertices or soft pion from D*). Selection flag stored in seleFlags. Create the AliESDVertex object (convert from AliAODVertex if necessary) In case of AOD input, also fill fAODMap for track index<->ID AliCodeTimerAuto("",0);
void SetSelectionBitForPID(AliRDHFCuts* cuts, AliAODRecoDecayHF* rd, Int_t bit)
Set the selection bit for PID
AliCodeTimerAuto("",0);
Bool_t SingleTrkCuts(AliESDtrack* trk, Float_t centralityperc, Bool_t& okDisplaced, Bool_t& okSoftPi, Bool_t& ok3prong) const
Check if track passes some kinematical cuts
AliAODv0* TransformESDv0toAODv0(AliESDv0* esdv0, TObjArray* twoTrackArrayV0)
Transform ESDv0 to AODv0
this function takes the ESDv0 vertex, computes the DCA variables from the ESDv0
and creates an AODv0 out of them
AliCodeTimerAuto("",0);
void SetParametersAtVertex(AliESDtrack* esdt, const AliExternalTrackParam* extpar) const
Set the stored track parameters at primary vertex into AliESDtrack AliCodeTimerAuto("",0);
void SetTrackFilter(AliAnalysisFilter* trackF)
switch off the TOF selection that cannot be applied with AODTracks
void SetTrackFilter2prongPbCentral(Float_t maxPercentile, AliAnalysisFilter* trackF)
switch off the TOF selection that cannot be applied with AODTracks
void SetTrackFilter3prongPbCentral(Float_t maxPercentile, AliAnalysisFilter* trackF)
switch off the TOF selection that cannot be applied with AODTracks
void SetTrackFilterSoftPi(AliAnalysisFilter* trackF)
switch off the TOF selection that cannot be applied with AODTracks
void SetnSigmaTPCforPionSel(Double_t nsl, Double_t nsh)
void SetnSigmaTOFforPionSel(Double_t nsl, Double_t nsh)
void SetnSigmaTPCforKaonSel(Double_t nsl, Double_t nsh)
void SetnSigmaTOFforKaonSel(Double_t nsl, Double_t nsh)
void SetnSigmaTPCforProtonSel(Double_t nsl, Double_t nsh)
void SetnSigmaTOFforProtonSel(Double_t nsl, Double_t nsh)
void GetnSigmaTOFforPionSel(Double_t& minnsigma, Double_t& maxnsigma) const
void GetnSigmaTPCforPionSel(Double_t& minnsigma, Double_t& maxnsigma) const
void GetnSigmaTOFforKaonSel(Double_t& minnsigma, Double_t& maxnsigma) const
void GetnSigmaTPCforKaonSel(Double_t& minnsigma, Double_t& maxnsigma) const
void GetnSigmaTOFforProtonSel(Double_t& minnsigma, Double_t& maxnsigma) const
void GetnSigmaTPCforProtonSel(Double_t& minnsigma, Double_t& maxnsigma) const