2 // Class AliRsnCutPrimaryVertex
4 // This cut implementation checks the quality of event primary vertex.
5 // It currently works only with ESD events (not AOD).
7 // authors: Martin Vala (martin.vala@cern.ch)
8 // Alberto Pulvirenti (alberto.pulvirenti@ct.infn.it)
11 #include "AliRsnCutPrimaryVertex.h"
13 ClassImp(AliRsnCutPrimaryVertex)
15 //_________________________________________________________________________________________________
16 AliRsnCutPrimaryVertex::AliRsnCutPrimaryVertex
17 (const char *name, Double_t maxVz, Int_t nContributors, Bool_t acceptTPC) :
18 AliRsnCut(name, AliRsnCut::kEvent, 0, nContributors - 1, 0.0, maxVz),
19 fAcceptTPC(acceptTPC),
24 // Defines the cut range between 0 and
25 // the minimum required number of contributors.
26 // The cut will be passed when if the event has a
27 // primary vertex with number of contributors outside this interval.
29 // If the 'acceptTPC' argument is true, events with TPC
30 // primary vertex will be checked, otherwise they will be
31 // rejected by default.
33 // Since the range check uses the '>=' and '<=', the high edge
34 // must be decreased by 1 to get the right behaviour, since this is integer.
41 //_________________________________________________________________________________________________
42 Bool_t AliRsnCutPrimaryVertex::IsSelected(TObject *object)
49 AliRsnEvent *rsn = dynamic_cast<AliRsnEvent*>(object);
50 if (!rsn) return kFALSE;
51 AliESDEvent *esd = rsn->GetRefESD();
52 AliAODEvent *aod = rsn->GetRefAOD();
59 if (esd->IsPileupFromSPD()) return kFALSE;
62 // get the best primary vertex:
63 // first try the one with tracks
64 const AliESDVertex *vTrk = esd->GetPrimaryVertexTracks();
65 const AliESDVertex *vSPD = esd->GetPrimaryVertexSPD();
66 const AliESDVertex *vTPC = esd->GetPrimaryVertexTPC();
70 Double_t vzTrk = 1000000.0;
71 Double_t vzSPD = 1000000.0;
72 Double_t vzTPC = 1000000.0;
73 if (vTrk) vzTrk = TMath::Abs(vTrk->GetZv());
74 if (vSPD) vzSPD = TMath::Abs(vSPD->GetZv());
75 if (vTPC) vzTPC = TMath::Abs(vTPC->GetZv());
76 if (vTrk) ncTrk = (Int_t)vTrk->GetNContributors();
77 if (vSPD) ncSPD = (Int_t)vSPD->GetNContributors();
78 if (vTPC) ncTPC = (Int_t)vTPC->GetNContributors();
84 else if (vSPD && ncSPD > 0)
89 else if (vTPC && ncTPC > 0)
104 // pile-up check is not yet available for AODs
106 // lines suggested by Andrea to reject TPC-only events
109 if (!aod->GetPrimaryVertexSPD()) return kFALSE;
110 else if(aod->GetPrimaryVertexSPD()->GetNContributors() < 1) return kFALSE;
113 AliAODVertex *prim = (AliAODVertex*)aod->GetPrimaryVertex();
114 if (!prim) return kFALSE;
116 fCutValueI = prim->GetNContributors();
117 fCutValueD = prim->GetZ();
123 Bool_t result = ((!OkRangeI()) && OkRangeD());
127 //_________________________________________________________________________________________________
128 void AliRsnCutPrimaryVertex::Print(const Option_t *) const
131 // Print information on this cut
134 AliInfo(Form("Cut name : %s", GetName()));
135 AliInfo(Form("Accepting TPC primary vertex : %s", (fAcceptTPC ? "YES" : "NO")));
136 AliInfo(Form("Contributors range (outside) : %d - %d", fMinI, fMaxI));
137 AliInfo(Form("Z-vertex range (inside) : %f - %f", fMinD, fMaxD));