class AliAnalysisTaskPID: public AliAnalysisTaskPIDV0base
Function Members (Methods)
public:
protected:
| void | CheckDoAnyStematicStudiesOnTheExpectedSignal() |
| void | AliAnalysisTaskPIDV0base::ClearV0PIDlist() |
| void | AliAnalysisTaskSE::ConnectMultiHandler() |
| Double_t | ConvolutedGaus(const Double_t* xx, const Double_t* par) const |
| void | AliAnalysisTask::DefineInput(Int_t islot, TClass* type) |
| void | AliAnalysisTask::DefineOutput(Int_t islot, TClass* type) |
| void | AliAnalysisTaskSE::DisconnectMultiHandler() |
| virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| void | AliAnalysisTask::EnableBranch(Int_t islot, const char* bname) const |
| Double_t | FastGaus(Double_t x, Double_t mean, Double_t sigma) const |
| Double_t | FastNormalisedGaus(Double_t x, Double_t mean, Double_t sigma) const |
| void | AliAnalysisTaskPIDV0base::FillV0PIDlist(AliESDEvent* esdEvent = 0x0) |
| Int_t | FindBinWithinRange(TAxis* axis, Double_t value) const |
| Int_t | FindFirstBinAboveIn3dSubset(const TH3* hist, Double_t threshold, Int_t yValue, Int_t zValue) const |
| Int_t | FindLastBinAboveIn3dSubset(const TH3* hist, Double_t threshold, Int_t yValue, Int_t zValue) const |
| char* | AliAnalysisTask::GetBranchAddress(Int_t islot, const char* branch) const |
| void | TObject::MakeZombie() |
| TFile* | AliAnalysisTask::OpenFile(Int_t iout, Option_t* option = "RECREATE") const |
| Bool_t | AliAnalysisTask::PostData(Int_t iout, TObject* data, Option_t* option = "") |
| Bool_t | AliAnalysisTask::SetBranchAddress(Int_t islot, const char* branch, void* address) const |
| virtual void | SetUpBinZeroStudyHist(THnSparse* hist, const Double_t* binsCent, const Double_t* binsPt) const |
| virtual void | SetUpDeDxCheckHist(THnSparse* hist, const Double_t* binsPt, const Double_t* binsJetPt, const Double_t* binsEtaAbs) const |
| virtual void | SetUpGenHist(THnSparse* hist, Double_t* binsPt, Double_t* binsDeltaPrime, Double_t* binsCent, Double_t* binsJetPt) const |
| virtual void | SetUpGenYieldHist(THnSparse* hist, Double_t* binsPt, Double_t* binsCent, Double_t* binsJetPt) const |
| virtual void | SetUpHist(THnSparse* hist, Double_t* binsPt, Double_t* binsDeltaPrime, Double_t* binsCent, Double_t* binsJetPt) const |
| virtual void | SetUpPIDcombined() |
| virtual void | SetUpPtResHist(THnSparse* hist, Double_t* binsPt, Double_t* binsJetPt, Double_t* binsCent) const |
| virtual void | SetUpSharedClsHist(THnSparse* hist, Double_t* binsPt, Double_t* binsJetPt) const |
private:
| AliAnalysisTaskPID(const AliAnalysisTaskPID&) | |
| AliAnalysisTaskPID& | operator=(const AliAnalysisTaskPID&) |
Data Members
public:
| enum ErrorCode { | kNoErrors | |
| kWarning | ||
| kError | ||
| }; | ||
| enum dataAxes { | kDataMCID | |
| kDataSelectSpecies | ||
| kDataPt | ||
| kDataDeltaPrimeSpecies | ||
| kDataCentrality | ||
| kDataJetPt | ||
| kDataZ | ||
| kDataXi | ||
| kDataCharge | ||
| kDataTOFpidInfo | ||
| kDataNumAxes | ||
| }; | ||
| enum genAxes { | kGenMCID | |
| kGenSelectSpecies | ||
| kGenPt | ||
| kGenDeltaPrimeSpecies | ||
| kGenCentrality | ||
| kGenJetPt | ||
| kGenZ | ||
| kGenXi | ||
| kGenCharge | ||
| kGenTOFpidInfo | ||
| kGenNumAxes | ||
| }; | ||
| enum genYieldAxes { | kGenYieldMCID | |
| kGenYieldPt | ||
| kGenYieldCentrality | ||
| kGenYieldJetPt | ||
| kGenYieldZ | ||
| kGenYieldXi | ||
| kGenYieldCharge | ||
| kGenYieldNumAxes | ||
| }; | ||
| enum ptResolutionAxes { | kPtResJetPt | |
| kPtResGenPt | ||
| kPtResRecPt | ||
| kPtResCharge | ||
| kPtResCentrality | ||
| kPtResNumAxes | ||
| }; | ||
| enum qaSharedClsAxes { | kQASharedClsJetPt | |
| kQASharedClsPt | ||
| kQASharedClsNumSharedCls | ||
| kQASharedClsPadRow | ||
| kQASharedClsNumAxes | ||
| }; | ||
| enum dEdxCheckAxes { | kDeDxCheckPID | |
| kDeDxCheckP | ||
| kDeDxCheckJetPt | ||
| kDeDxCheckEtaAbs | ||
| kDeDxCheckDeDx | ||
| kDeDxCheckNumAxes | ||
| }; | ||
| enum binZeroStudyAxes { | kBinZeroStudyCentrality | |
| kBinZeroStudyGenPt | ||
| kBinZeroStudyGenEta | ||
| kBinZeroStudyNumAxes | ||
| }; | ||
| enum efficiencyAxes { | kEffMCID | |
| kEffTrackPt | ||
| kEffTrackEta | ||
| kEffTrackCharge | ||
| kEffCentrality | ||
| kEffJetPt | ||
| kEffZ | ||
| kEffXi | ||
| kEffNumAxes | ||
| }; | ||
| enum EffSteps { | kStepGenWithGenCuts | |
| kStepRecWithGenCuts | ||
| kStepRecWithGenCutsMeasuredObs | ||
| kStepRecWithRecCutsMeasuredObs | ||
| kStepRecWithRecCutsMeasuredObsPrimaries | ||
| kStepRecWithRecCutsMeasuredObsStrangenessScaled | ||
| kStepRecWithRecCutsPrimaries | ||
| kNumSteps | ||
| }; | ||
| enum TOFpidInfo { | kNoTOFinfo | |
| kNoTOFpid | ||
| kTOFpion | ||
| kTOFkaon | ||
| kTOFproton | ||
| kNumTOFspecies | ||
| kNumTOFpidInfoBins | ||
| }; | ||
| enum EventCounterType { | kTriggerSel | |
| kTriggerSelAndVtxCut | ||
| kTriggerSelAndVtxCutAndZvtxCutNoPileUpRejection | ||
| kTriggerSelAndVtxCutAndZvtxCut | ||
| }; | ||
| enum CutHistoType { | kMCPtHardCut | |
| }; | ||
| enum AliAnalysisTaskPIDV0base::PileUpRejectionType { | kPileUpRejectionOff | |
| kPileUpRejectionSPD | ||
| kPileUpRejectionMV | ||
| }; | ||
| enum AliAnalysisTaskPIDV0base::TPCcutType { | kNoCut | |
| kTPCCutMIGeo | ||
| kTPCnclCut | ||
| }; | ||
| enum AliAnalysisTask::EAnalysisTaskFlags { | kTaskUsed | |
| kTaskZombie | ||
| kTaskChecked | ||
| kTaskPostEventLoop | ||
| }; | ||
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
protected:
| Bool_t | TTask::fActive | true if task is active |
| AliAnalysisUtils* | AliAnalysisTaskPIDV0base::fAnaUtils | ! Object to use analysis utils like pile-up rejection |
| TString | AliAnalysisTask::fBranchNames | List of input branches that need to be loaded for this task |
| Int_t | TTask::fBreakin | =1 if a break point set at task extry |
| Int_t | TTask::fBreakout | =1 if a break point set at task exit |
| Int_t | AliAnalysisTaskSE::fCurrentRunNumber | ! Current run number |
| Int_t | AliAnalysisTaskSE::fDebug | Debug flag |
| AliESDEvent* | AliAnalysisTaskPIDV0base::fESD | ! ESDEvent object, if ESD |
| AliESDfriend* | AliAnalysisTaskSE::fESDfriend | ! ESD friend |
| Int_t | AliAnalysisTaskSE::fEntry | Current entry in the chain |
| Double_t | AliAnalysisTaskPIDV0base::fEtaCut | Eta cut |
| AliVEvent* | AliAnalysisTaskPIDV0base::fEvent | ! VEvent object |
| Bool_t | TTask::fHasExecuted | True if task has executed |
| TList* | AliAnalysisTaskSE::fHistosQA | ! Output histos for QA |
| Bool_t | AliAnalysisTask::fInitialized | True if Init() was called |
| AliVEvent* | AliAnalysisTaskSE::fInputEvent | ! VEvent Input |
| AliInputEventHandler* | AliAnalysisTaskSE::fInputHandler | ! Input Handler |
| TObjArray* | AliAnalysisTask::fInputs | Array of input slots |
| Bool_t | AliAnalysisTaskPIDV0base::fIsPbpOrpPb | Pbp/pPb collision or something else? |
| AliMCEvent* | AliAnalysisTaskPIDV0base::fMC | ! MC object |
| AliMCEvent* | AliAnalysisTaskSE::fMCEvent | ! MC |
| AliInputEventHandler* | AliAnalysisTaskSE::fMCEventHandler | ! pointer to MCEventHandler |
| AliMultiInputEventHandler* | AliAnalysisTaskSE::fMultiInputHandler | ! pointer to multihandler |
| TString | TNamed::fName | object identifier |
| Int_t | AliAnalysisTask::fNinputs | Number of inputs |
| Int_t | AliAnalysisTask::fNoutputs | Number of outputs |
| Int_t | AliAnalysisTaskPIDV0base::fNumTagsStored | Number of entries of fV0tags |
| UInt_t | AliAnalysisTaskSE::fOfflineTriggerMask | Task processes collision candidates only |
| TString | TTask::fOption | Option specified in ExecuteTask |
| AliAODEvent* | AliAnalysisTaskSE::fOutputAOD | ! AOD out |
| Bool_t* | AliAnalysisTask::fOutputReady | [fNoutputs] Flags for output readyness |
| TObjArray* | AliAnalysisTask::fOutputs | Array of output slots |
| AliPIDResponse* | AliAnalysisTaskPIDV0base::fPIDResponse | ! PID response Handler |
| TF1* | AliAnalysisTaskPIDV0base::fPhiCutHigh | phi prime cut, high |
| TF1* | AliAnalysisTaskPIDV0base::fPhiCutLow | phi prime cut, low |
| TObject* | AliAnalysisTask::fPublishedData | ! published data |
| TRandom3* | AliAnalysisTaskPIDV0base::fRandom | ! Can be used to statistically determine the shape in the pt bins e.g. |
| Bool_t | AliAnalysisTask::fReady | Flag if the task is ready |
| Bool_t | AliAnalysisTaskPIDV0base::fStoreMotherIndex | Switch on/off storing the mother indices of V0 daughters |
| AliAnalysisTaskPIDV0base::TPCcutType | AliAnalysisTaskPIDV0base::fTPCcutType | Type of TPC cut to be used |
| TList* | TTask::fTasks | List of Tasks |
| TString | TNamed::fTitle | object title |
| AliAnalysisFilter* | AliAnalysisTaskPIDV0base::fTrackFilter | Track Filter |
| TTree* | AliAnalysisTaskSE::fTreeA | AOD output Tree |
| Bool_t | AliAnalysisTaskPIDV0base::fUsePhiCut | Use cut on phi (useful for TPC) |
| AliESDv0KineCuts* | AliAnalysisTaskPIDV0base::fV0KineCuts | ! ESD V0 kine cuts |
| Int_t* | AliAnalysisTaskPIDV0base::fV0motherIndex | ! Pointer to array with index of the mother V0 |
| Char_t* | AliAnalysisTaskPIDV0base::fV0tags | ! Pointer to array with tags for identified particles from V0 decays |
| Double_t | AliAnalysisTaskPIDV0base::fZvtxCutEvent | Vertex z cut for the event (cm) |
| static TClonesArray* | AliAnalysisTaskSE::fgAODCaloClusters | ! CaloClusters for replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODDimuons | ! Dimuons replication |
| static AliAODCaloTrigger* | AliAnalysisTaskSE::fgAODEMCALTrigger | ! Emcal Trigger for replication |
| static AliAODCaloCells* | AliAnalysisTaskSE::fgAODEmcalCells | ! Emcal Cell replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODFMDClusters | ! FMDClusters for replication |
| static AliVHeader* | AliAnalysisTaskSE::fgAODHeader | ! Header for replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODHmpidRings | ! HMPID replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODJets | ! Jets for replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODMCParticles | ! MC Particles for replicatio |
| static AliAODCaloTrigger* | AliAnalysisTaskSE::fgAODPHOSTrigger | ! Phos Trigger for replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODPMDClusters | ! PMDClusters for replication |
| static AliAODCaloCells* | AliAnalysisTaskSE::fgAODPhosCells | ! Phos Cell replication |
| static AliAODTracklets* | AliAnalysisTaskSE::fgAODTracklets | ! Tracklets for replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODTracks | ! Tracks for replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODV0s | ! V0s for replication |
| static AliAODVZERO* | AliAnalysisTaskSE::fgAODVZERO | ! VZERO for replication |
| static TClonesArray* | AliAnalysisTaskSE::fgAODVertices | ! Vertices for replication |
| static TTask* | TTask::fgBeginTask | pointer to task initiator |
| static TTask* | TTask::fgBreakPoint | pointer to current break point |
| static Double_t | AliAnalysisTaskPIDV0base::fgCutGeo | Cut variable for TPCCutMIGeo concerning geometry |
| static Double_t | AliAnalysisTaskPIDV0base::fgCutNcl | Cut variable for TPCCutMIGeo concerning num clusters |
| static Double_t | AliAnalysisTaskPIDV0base::fgCutNcr | Cut variable for TPCCutMIGeo concerning num crossed rows |
| static UShort_t | AliAnalysisTaskPIDV0base::fgCutPureNcl | Cut variable for TPCnclCut |
| static AliTOFHeader* | AliAnalysisTaskSE::fgTOFHeader | ! TOFHeader for replication |
| static const Double_t | fgkEpsilon | Double_t threshold above zero |
| static const Int_t | fgkMaxNumGenEntries | Maximum number of generated detector responses per track and delta(Prime) and associated species |
| static const Int_t | fgkNumJetAxes | Number of additional axes for jets |
| static const Double_t | fgkOneOverSqrt2 | = 1. / TMath::Sqrt2(); |
private:
| Double_t | fAccuracyNonGaussianTail | Accuracy of the non-gaussian tail (fraction of the maximum) |
| TString | fCentralityEstimator | Estimator for the centrality (e.g. V0A, V0M) |
| THnSparseD* | fChargedGenPrimariesTriggerSel | ! Histo holding the generated charged primary yields for triggered events |
| THnSparseD* | fChargedGenPrimariesTriggerSelVtxCut | ! Histo holding the generated charged primary yields for triggered events passing vertex cuts |
| THnSparseD* | fChargedGenPrimariesTriggerSelVtxCutZ | ! Histo holding the generated charged primary yields for triggered events passing vertex cuts (including cut on z) |
| THnSparseD* | fChargedGenPrimariesTriggerSelVtxCutZPileUpRej | ! Histo holding the generated charged primary yields for triggered events passing vertex cuts (including cut on z) and pile-up rejection |
| AliCFContainer* | fContainerEff | ! Container for efficiency determination |
| TF1* | fConvolutedGausDeltaPrime | Gaus convoluted with exponential tail to generate detector response (deltaPrime) |
| Double_t | fConvolutedGaussTransitionPars[3] | Parameter for transition from gaussian parameters to asymmetric shape |
| THnSparseD* | fDeDxCheck | ! dEdx check |
| TAxis* | fDeltaPrimeAxis | ! Axis holding the deltaPrime binning |
| Bool_t | fDoAnySystematicStudiesOnTheExpectedSignal | Internal flag indicating whether any systematic studies are going to be performed |
| Bool_t | fDoBinZeroStudy | Do bin zero study, if flag is set to kTRUE |
| Bool_t | fDoDeDxCheck | Check dEdx, if flag set to kTRUE |
| Bool_t | fDoEfficiency | Do efficiency processing (and post the output), if flag is set to kTRUE |
| Bool_t | fDoPID | Do PID processing (and post the output), if flag is set to kTRUE |
| Bool_t | fDoPtResolution | Do pT resolution processing (and post the output), if flag is set to kTRUE |
| Double_t | fEtaAbsCutLow | Lower cut value on |eta| |
| Double_t | fEtaAbsCutUp | Upper cut value on |eta| |
| TH3D* | fFractionHists[5] | 3D histos of particle fraction as function with trackPt, jetPt (-1 for inclusive spectra), centralityPercentile (-1 for pp) |
| TH3D* | fFractionSysErrorHists[5] | 3D histos of sys. error of particle fraction as function with trackPt, jetPt (-1 for inclusive spectra), centralityPercentile (-1 for pp) |
| Double_t* | fGenRespElDeltaPrimeEl | ! Generated responses for a single track |
| Double_t* | fGenRespElDeltaPrimeKa | ! Generated responses for a single track |
| Double_t* | fGenRespElDeltaPrimePi | ! Generated responses for a single track |
| Double_t* | fGenRespElDeltaPrimePr | ! Generated responses for a single track |
| Double_t* | fGenRespKaDeltaPrimeEl | ! Generated responses for a single track |
| Double_t* | fGenRespKaDeltaPrimeKa | ! Generated responses for a single track |
| Double_t* | fGenRespKaDeltaPrimePi | ! Generated responses for a single track |
| Double_t* | fGenRespKaDeltaPrimePr | ! Generated responses for a single track |
| Double_t* | fGenRespMuDeltaPrimeEl | ! Generated responses for a single track |
| Double_t* | fGenRespMuDeltaPrimeKa | ! Generated responses for a single track |
| Double_t* | fGenRespMuDeltaPrimePi | ! Generated responses for a single track |
| Double_t* | fGenRespMuDeltaPrimePr | ! Generated responses for a single track |
| Double_t* | fGenRespPiDeltaPrimeEl | ! Generated responses for a single track |
| Double_t* | fGenRespPiDeltaPrimeKa | ! Generated responses for a single track |
| Double_t* | fGenRespPiDeltaPrimePi | ! Generated responses for a single track |
| Double_t* | fGenRespPiDeltaPrimePr | ! Generated responses for a single track |
| Double_t* | fGenRespPrDeltaPrimeEl | ! Generated responses for a single track |
| Double_t* | fGenRespPrDeltaPrimeKa | ! Generated responses for a single track |
| Double_t* | fGenRespPrDeltaPrimePi | ! Generated responses for a single track |
| Double_t* | fGenRespPrDeltaPrimePr | ! Generated responses for a single track |
| Bool_t | fInputFromOtherTask | If set to kTRUE, no events are processed and the input must be fed in from another task. If set to kFALSE, normal event processing |
| TObjArray* | fOutputContainer | ! output data container |
| AliPIDCombined* | fPIDcombined | ! PID combined object |
| AliAnalysisTaskPIDV0base::PileUpRejectionType | fPileUpRejectionType | Which pile-up rejection is used (if any) |
| THnSparseD* | fPtResolution[5] | ! Pt Resolution for the individual species |
| TObjArray* | fQAContainer | ! output data container for QA |
| THnSparseD* | fQASharedCls | ! QA for shared clusters |
| Int_t | fRun | Current run number |
| Bool_t | fStoreAdditionalJetInformation | If set to kTRUE, additional jet information like jetPt, z, xi will be stored in the THnSparses |
| Bool_t | fStoreCentralityPercentile | If set to kTRUE, store centrality percentile for each event. In case of kFALSE (appropriate for pp), centrality percentile will be set to -1 for every event |
| Double_t | fSystematicScalingEtaCorrectionHighMomenta | Systematic scale factor for the eta correction (1. = no systematics) at high momenta |
| Double_t | fSystematicScalingEtaCorrectionLowMomenta | Systematic scale factor for the eta correction (1. = no systematics) at low momenta |
| Double_t | fSystematicScalingEtaCorrectionMomentumThr | Momentum threshold for the systematic scale factor for the eta correction (separates low-p from high-p |
| Double_t | fSystematicScalingEtaSigmaParaAboveThreshold | Systematic scale factor for the parametrisation of the relative signal width (1. = no systematics) above threshold |
| Double_t | fSystematicScalingEtaSigmaParaBelowThreshold | Systematic scale factor for the parametrisation of the relative signal width (1. = no systematics) below threshold |
| Double_t | fSystematicScalingEtaSigmaParaThreshold | dEdx threshold for the systematic scale factor for the parametrisation of the relative signal width |
| Double_t | fSystematicScalingMultCorrection | Systematic scale factor for the multiplicity correction (1. = no systematics) |
| Double_t | fSystematicScalingSplinesAboveThreshold | Systematic scale factor for the splines (1. = no systematics) above threshold |
| Double_t | fSystematicScalingSplinesBelowThreshold | Systematic scale factor for the splines (1. = no systematics) below threshold |
| Double_t | fSystematicScalingSplinesThreshold | beta-gamma threshold for the systematic spline scale factor |
| Int_t | fTOFmode | TOF mode used for TOF PID info (affects num sigma inclusion/exclusion) |
| Bool_t | fTPCDefaultPriors | Use TPC default priors for PID combined probabilities, if priors are enabled |
| Bool_t | fTakeIntoAccountMuons | Also take into account muons for the generation of the expected response and the most probable PID |
| Bool_t | fUseConvolutedGaus | Use convoluted gaus to generate detector response instead of pure gaus |
| Bool_t | fUseITS | Use ITS for PID combined probabilities |
| Bool_t | fUseMCidForGeneration | If MC, use MCid instead of PIDcombined to generate the signals |
| Bool_t | fUsePriors | Use priors for PID combined probabilities |
| Bool_t | fUseTOF | Use TOF for PID combined probabilities |
| static const Double_t | fgkSigmaReferenceForTransitionPars | Reference sigma chosen to calculate transition parameters |
| TH1F* | fh1EvtsPtHardCut | ! Number events before and after the cut on MC pT hard |
| TH1D* | fh1Trials | ! sum of trials |
| TProfile* | fh1Xsec | ! pythia cross section and trials |
| TH2D* | fh2FFJetPtGen | ! Number of generated jets vs. jetPt and centrality |
| TH2D* | fh2FFJetPtRec | ! Number of reconstructed jets vs. jetPt and centrality |
| TH1D* | fhEventsProcessed | ! Histo holding the number of processed events (i.e. passing trigger selection, vtx and zvtx cuts and (if enabled) pile-up rejection) |
| TH1D* | fhEventsProcessedNoPileUpRejection | ! Histo holding the number of processed events before pile-up rejection |
| TH1D* | fhEventsTriggerSel | ! Histo holding the number of events passing trigger selection |
| TH1D* | fhEventsTriggerSelVtxCut | ! Histo holding the number of events passing trigger selection and vtx cut |
| THnSparseD* | fhGenEl | ! Generated response for el |
| THnSparseD* | fhGenKa | ! Generated response for ka |
| THnSparseD* | fhGenMu | ! Generated response for mu |
| THnSparseD* | fhGenPi | ! Generated response for pi |
| THnSparseD* | fhGenPr | ! Generated response for pr |
| THnSparseD* | fhMCgeneratedYieldsPrimaries | ! Histo holding the generated (no reco, no cuts) primary particle yields in considered eta range |
| TH1D* | fhMaxEtaVariation | ! Histo holding the maximum deviation of the eta correction factor from unity vs. 1/dEdx(splines) |
| THnSparseD* | fhPIDdataAll | ! Data histo |
| const Int_t | fkConvolutedGausNPar | Number of parameters for convoluted gaus |
| const Double_t | fkDeltaPrimeLowLimit | Lower deltaPrime limit |
| const Double_t | fkDeltaPrimeUpLimit | Upper deltaPrime limit |
Class Charts
| Inheritance Chart: | ||||||||||||||||||||
|
Function documentation
AliAnalysisTaskPID()
, fGenRespElDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespElDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespElDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespElDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaPr(new Double_t[fgkMaxNumGenEntries]) default Constructor
AliAnalysisTaskPID(const char* name)
, fGenRespElDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespElDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespElDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespElDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespKaDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespPiDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespMuDeltaPr(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaEl(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaKa(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaPi(new Double_t[fgkMaxNumGenEntries]) , fGenRespPrDeltaPr(new Double_t[fgkMaxNumGenEntries]) Constructor
Bool_t CalculateMaxEtaVariationMapFromPIDResponse()
Calculate the maximum deviation from unity of the eta correction factors for each row in 1/dEdx(splines) from the eta correction map of the TPCPIDResponse. The result is stored in fhMaxEtaVariation.
void CheckDoAnyStematicStudiesOnTheExpectedSignal()
Check whether at least one scale factor indicates the ussage of systematic studies and set internal flag accordingly.
void ConfigureTaskForCurrentEvent(AliVEvent* event)
Configure the task for the current event. In particular, this is needed if the run number changes
void GetJetTrackObservables(Double_t trackPt, Double_t jetPt, Double_t& z, Double_t& xi)
Uses trackPt and jetPt to obtain z and xi.
Double_t ConvolutedGaus(const Double_t* xx, const Double_t* par) const
Convolutes gauss with an exponential tail which describes dEdx-response better than pure gaussian
Double_t FastGaus(Double_t x, Double_t mean, Double_t sigma) const
Calculate an unnormalised gaussian function with mean and sigma.
Double_t FastNormalisedGaus(Double_t x, Double_t mean, Double_t sigma) const
Calculate a normalised (divided by sqrt(2*Pi)*sigma) gaussian function with mean and sigma.
Int_t FindBinWithinRange(TAxis* axis, Double_t value) const
Find the corresponding bin of the axis. Values outside the range (also under and overflow) will be set to the first/last available bin
Int_t FindFirstBinAboveIn3dSubset(const TH3* hist, Double_t threshold, Int_t yValue, Int_t zValue) const
Kind of projects a TH3 to 1 bin combination in y and z and looks for the first x bin above a threshold for this projection. If no such bin is found, -1 is returned.
Int_t FindLastBinAboveIn3dSubset(const TH3* hist, Double_t threshold, Int_t yValue, Int_t zValue) const
Kind of projects a TH3 to 1 bin combination in y and z and looks for the last x bin above a threshold for this projection. If no such bin is found, -1 is returned.
Bool_t GetParticleFraction(Double_t trackPt, Double_t jetPt, Double_t multiplicity, AliPID::EParticleType species, Double_t& fraction, Double_t& fractionErrorStat, Double_t& fractionErrorSys) const
Computes the particle fraction for the corresponding species for the given trackPt, jetPt and centrality. Use jetPt = -1 for inclusive spectra and centralityPercentile = -1 for pp. On success (return value kTRUE), fraction contains the particle fraction, fractionErrorStat(Sys) the sigma of its statistical (systematic) error
Bool_t GetParticleFractions(Double_t trackPt, Double_t jetPt, Double_t centralityPercentile, Double_t* prob, Int_t smearSpeciesByError, Int_t takeIntoAccountSpeciesSysError, Bool_t uniformSystematicError = kFALSE) const
Fills the particle fractions for the given trackPt, jetPt and centrality into "prob". Use jetPt = -1 for inclusive spectra and centralityPercentile = -1 for pp. If smearSpeciesByError is >= 0 && < AliPID::kSPECIES, the returned fractions will be a random number distributed with a gauss with mean being the corresponding particle fraction and sigma it's error for the considered species "smearSpeciesByError". Note that in this case the fractions for all species will NOT sum up to 1! Thus, all other species fractions will be re-scaled weighted with their corresponding statistical error. A similar procedure is used for "takeIntoAccountSpeciesSysError": The systematic error of the corresponding species is used to generate a random number with uniform distribution in [mean - sysError, mean + sysError] for the new mean (in cace of uniformSystematicError = kTRUE, otherwise it will be a gaus(mean, sysError)), then the other species will be re-scaled according to their systematic errors. First, the systematic error uncertainty procedure will be performed (that is including re-scaling), then the statistical uncertainty procedure. On success, kTRUE is returned.
const TH3D* GetParticleFractionHisto(Int_t species, Bool_t sysError = kFALSE) const
Double_t GetMCStrangenessFactorCMS(Int_t motherPDG, Double_t motherGenPt)
Strangeness ratio MC/data as function of mother pt from CMS data in |eta|<2.0 -> Based on function in PWGJE/AliAnalysisTaskFragmentationFunction, which uses the following data from CMS pp @ 7 TeV inclusive (JHEP 05 (2011) 064)
Double_t GetMCStrangenessFactorCMS(AliMCEvent* mcEvent, AliMCParticle* daughter)
Strangeness ratio MC/data as function of mother pt from CMS data in |eta|<2.0 -> Based on function in PWGJE/AliAnalysisTaskFragmentationFunction
Bool_t IsSecondaryWithStrangeMotherMC(AliMCEvent* mcEvent, Int_t partLabel)
Bool_t SetParticleFractionHisto(const TH3D* hist, Int_t species, Bool_t sysError = kFALSE)
Bool_t SetParticleFractionHistosFromFile(const TString filePathName, Bool_t sysError = kFALSE)
Loads particle fractions for all species from the desired file and returns kTRUE on success. The maps are assumed to be of Type TH3D, to sit in the main directory and to have names Form("hFraction_%e", AliPID::ParticleName(i)) for sysError = kFALSE and Form("hFractionSysError_%e", AliPID::ParticleName(i)) for sysError = kTRUE.
Double_t GetMaxEtaVariation(Double_t dEdxSplines)
Returns the maximum eta variation (i.e. deviation of eta correction factor from unity) for the given (spline) dEdx
Int_t GetRandomParticleTypeAccordingToParticleFractions(Double_t trackPt, Double_t jetPt, Double_t centralityPercentile, Bool_t smearByError, Bool_t takeIntoAccountSysError = kFALSE) const
Uses the stored histograms with the particle fractions to generate a random particle type according to these fractions. In case of problems (e.g. histo missing), AliPID::kUnknown is returned. If smearByError is kTRUE, the used fractions will be random numbers distributed with a gauss with mean being the corresponding particle fraction and sigma it's error. Note that in this case only the fraction of a random species is varied in this way. The other fractions will be re-normalised according their statistical errors. The same holds for the systematic error of species "takeIntoAccountSpeciesSysError", but the random number will be uniformly distributed within [mean - sys, mean + sys] and the re-normalisation will be weighted with the systematic errors. Note that the fractions will be calculated first with only the systematic error taken into account (if desired), including re-normalisation. Then, the resulting fractions will be used to calculate the final fractions - either with statistical error or without. The species, for which the error will be used for smearing, is the same for sys and stat error.
Bool_t ProcessTrack(const AliVTrack* track, Int_t particlePDGcode, Double_t centralityPercentile, Double_t jetPt)
Process the track (generate expected response, fill histos, etc.). particlePDGcode == 0 means data. Otherwise, the corresponding MC ID will be assumed.
Bool_t SetConvolutedGaussLambdaParameter(Double_t lambda)
Set the lambda parameter of the convolution to the desired value. Automatically calculates the parameters for the transition (restricted) gauss -> convoluted gauss.
void SetUpGenHist(THnSparse* hist, Double_t* binsPt, Double_t* binsDeltaPrime, Double_t* binsCent, Double_t* binsJetPt) const
Sets bin limits for axes which are not standard binned and the axes titles.
void SetUpGenYieldHist(THnSparse* hist, Double_t* binsPt, Double_t* binsCent, Double_t* binsJetPt) const
Sets bin limits for axes which are not standard binned and the axes titles.
void SetUpHist(THnSparse* hist, Double_t* binsPt, Double_t* binsDeltaPrime, Double_t* binsCent, Double_t* binsJetPt) const
Sets bin limits for axes which are not standard binned and the axes titles.
void SetUpPtResHist(THnSparse* hist, Double_t* binsPt, Double_t* binsJetPt, Double_t* binsCent) const
Sets bin limits for axes which are not standard binned and the axes titles.
void SetUpSharedClsHist(THnSparse* hist, Double_t* binsPt, Double_t* binsJetPt) const
Sets bin limits for axes which are not standard binned and the axes titles.
void SetUpDeDxCheckHist(THnSparse* hist, const Double_t* binsPt, const Double_t* binsJetPt, const Double_t* binsEtaAbs) const
Sets bin limits for axes which are not standard binned and the axes titles.
void SetUpBinZeroStudyHist(THnSparse* hist, const Double_t* binsCent, const Double_t* binsPt) const
Sets bin limits for axes which are not standard binned and the axes titles.
Bool_t FillEfficiencyContainer(const Double_t* values, AliAnalysisTaskPID::EffSteps step, Double_t weight = 1.0)
Bool_t IncrementEventCounter(Double_t centralityPercentile, AliAnalysisTaskPID::EventCounterType type)
Int_t GetIndexOfChargeAxisData() const
{ return fStoreAdditionalJetInformation ? kDataCharge : kDataCharge - fgkNumJetAxes; }Int_t GetIndexOfChargeAxisGen() const
{ return fStoreAdditionalJetInformation ? kGenCharge : kGenCharge - fgkNumJetAxes; }Int_t GetIndexOfChargeAxisGenYield() const
{ return fStoreAdditionalJetInformation ? kGenYieldCharge : kGenYieldCharge - fgkNumJetAxes; }Int_t GetIndexOfTOFpidInfoAxisData() const
{ return fStoreAdditionalJetInformation ? kDataTOFpidInfo : kDataTOFpidInfo - fgkNumJetAxes; }Int_t GetIndexOfTOFpidInfoAxisGen() const
{ return fStoreAdditionalJetInformation ? kGenTOFpidInfo : kGenTOFpidInfo - fgkNumJetAxes; }Bool_t FillPythiaTrials(Double_t avgTrials)
ErrorCode GenerateDetectorResponse(AliAnalysisTaskPID::ErrorCode errCode, Double_t mean, Double_t sigma, Double_t* responses, Int_t nResponses, Bool_t usePureGaus = kFALSE)
const TString GetPPCentralityEstimator() const
Bool_t IsInAcceptedEtaRange(Double_t etaAbs) const
{ return (etaAbs >= fEtaAbsCutLow && etaAbs <= fEtaAbsCutUp); }AliAnalysisTaskPIDV0base::PileUpRejectionType GetPileUpRejectionType() const
{ return fPileUpRejectionType; }void SetPileUpRejectionType(AliAnalysisTaskPIDV0base::PileUpRejectionType newType)
{ fPileUpRejectionType = newType; }Double_t GetSystematicScalingSplinesThreshold() const
{ return fSystematicScalingSplinesThreshold; }void SetSystematicScalingSplinesThreshold(Double_t threshold)
{ fSystematicScalingSplinesThreshold = threshold; }Double_t GetSystematicScalingSplinesBelowThreshold() const
{ return fSystematicScalingSplinesBelowThreshold; }void SetSystematicScalingSplinesBelowThreshold(Double_t scaleFactor)
{ fSystematicScalingSplinesBelowThreshold = scaleFactor; CheckDoAnyStematicStudiesOnTheExpectedSignal(); }Double_t GetSystematicScalingSplinesAboveThreshold() const
{ return fSystematicScalingSplinesAboveThreshold; }void SetSystematicScalingSplinesAboveThreshold(Double_t scaleFactor)
{ fSystematicScalingSplinesAboveThreshold = scaleFactor; CheckDoAnyStematicStudiesOnTheExpectedSignal(); }Double_t GetSystematicScalingEtaCorrectionMomentumThr() const
{ return fSystematicScalingEtaCorrectionMomentumThr; }void SetSystematicScalingEtaCorrectionMomentumThr(Double_t threshold)
{ fSystematicScalingEtaCorrectionMomentumThr = threshold; }Double_t GetSystematicScalingEtaCorrectionLowMomenta() const
{ return fSystematicScalingEtaCorrectionLowMomenta; }void SetSystematicScalingEtaCorrectionLowMomenta(Double_t scaleFactor)
{ fSystematicScalingEtaCorrectionLowMomenta = scaleFactor; CheckDoAnyStematicStudiesOnTheExpectedSignal(); }Double_t GetSystematicScalingEtaCorrectionHighMomenta() const
{ return fSystematicScalingEtaCorrectionHighMomenta; }void SetSystematicScalingEtaCorrectionHighMomenta(Double_t scaleFactor)
{ fSystematicScalingEtaCorrectionHighMomenta = scaleFactor; CheckDoAnyStematicStudiesOnTheExpectedSignal(); }Double_t GetSystematicScalingEtaSigmaParaThreshold() const
{ return fSystematicScalingEtaSigmaParaThreshold; }void SetSystematicScalingEtaSigmaParaThreshold(Double_t threshold)
{ fSystematicScalingEtaSigmaParaThreshold = threshold; }Double_t GetSystematicScalingEtaSigmaParaBelowThreshold() const
{ return fSystematicScalingEtaSigmaParaBelowThreshold; }void SetSystematicScalingEtaSigmaParaBelowThreshold(Double_t scaleFactor)
{ fSystematicScalingEtaSigmaParaBelowThreshold = scaleFactor; CheckDoAnyStematicStudiesOnTheExpectedSignal(); }Double_t GetSystematicScalingEtaSigmaParaAboveThreshold() const
{ return fSystematicScalingEtaSigmaParaAboveThreshold; }void SetSystematicScalingEtaSigmaParaAboveThreshold(Double_t scaleFactor)
{ fSystematicScalingEtaSigmaParaAboveThreshold = scaleFactor; CheckDoAnyStematicStudiesOnTheExpectedSignal(); }void SetSystematicScalingMultCorrection(Double_t scaleFactor)
{ fSystematicScalingMultCorrection = scaleFactor; CheckDoAnyStematicStudiesOnTheExpectedSignal(); }TOFpidInfo GetTOFType(const AliVTrack* track, Int_t tofMode) const
AliAnalysisTaskPID& operator=(const AliAnalysisTaskPID& )