// get reconstructed vertex
const AliESDVertex* vtxESD = 0;
- Bool_t isRecVertex = kTRUE;
if(evtCuts->IsRecVertexRequired())
{
-// Bool_t bRedoTPCVertex = evtCuts->IsRedoTPCVertex();
- //Bool_t bUseConstraints = evtCuts->IsUseBeamSpotConstraint();
- //vtxESD = AlidNdPtHelper::GetVertex(esdEvent,evtCuts,accCuts,esdTrackCuts,GetAnalysisMode(),kFALSE,bRedoTPCVertex,bUseConstraints);
- //isRecVertex = AlidNdPtHelper::TestRecVertex(vtxESD, esdEvent->GetPrimaryVertexSPD(), GetAnalysisMode(), kFALSE);
if(GetAnalysisMode() == AlidNdPtHelper::kTPC) {
vtxESD = esdEvent->GetPrimaryVertexTPC();
- //isRecVertex = AlidNdPtHelper::TestRecVertex(vtxESD, esdEvent->GetPrimaryVertexSPD(), GetAnalysisMode(), kFALSE);
+ }
+ else if(GetAnalysisMode() == AlidNdPtHelper::kTPCITS) {
+ vtxESD = esdEvent->GetPrimaryVertexTracks();
+ }
+ else {
+ return;
}
}
- Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD) && isRecVertex;
+ Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD);
//printf("isEventOK %d, isEventTriggered %d \n",isEventOK, isEventTriggered);
//printf("GetAnalysisMode() %d \n",GetAnalysisMode());
Int_t multMBTracks = 0;
if(GetAnalysisMode() == AlidNdPtHelper::kTPC)
{
- if(vtxESD->GetStatus() && isRecVertex) {
- //multMBTracks = AlidNdPtHelper::GetTPCMBTrackMult(esdEvent,evtCuts,accCuts,esdTrackCuts);
+ if(vtxESD->GetStatus()) {
multMBTracks = vtxESD->GetNContributors();
}
}
+ else if(GetAnalysisMode() == AlidNdPtHelper::kTPCITS) {
+ if(vtxESD->GetStatus()) {
+ multMBTracks = vtxESD->GetNContributors();
+ }
+ }
else {
AliDebug(AliLog::kError, Form("Found analysis type %d", GetAnalysisMode()));
return;
//\r
//nClust:chi2PerClust:nClust/nFindableClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromMaterial:isPrim:charge\r
Int_t binsRecMCTrackHist[12]={160,20,40,100,100,30,90,ptNbins, 2, 2, 2, 3};\r
- Double_t minRecMCTrackHist[12]={0., 0., 0., -10.,-10.,-1.5, 0., ptMin, 0., 0., 0., -1.};\r
- Double_t maxRecMCTrackHist[12]={160.,10.,1.2, 10.,10.,1.5, 2.*TMath::Pi(), ptMax, 2., 2., 2., 2.};\r
+ //Double_t minRecMCTrackHist[12]={0., 0., 0., -10.,-10.,-1.5, 0., ptMin, 0., 0., 0., -1.};\r
+ //Double_t maxRecMCTrackHist[12]={160.,10.,1.2, 10.,10.,1.5, 2.*TMath::Pi(), ptMax, 2., 2., 2., 2.};\r
+ Double_t minRecMCTrackHist[12]={0., 0., 0., -1.,-1.,-1.5, 0., ptMin, 0., 0., 0., -1.};\r
+ Double_t maxRecMCTrackHist[12]={160.,10.,1.2, 1.,1.,1.5, 2.*TMath::Pi(), ptMax, 2., 2., 2., 2.};\r
\r
fRecMCTrackHist = new THnSparseF("fRecMCTrackHist","nClust:chi2PerClust:nClust/nFindableClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromMaterial:isPrim:charge",12,binsRecMCTrackHist,minRecMCTrackHist,maxRecMCTrackHist);\r
fRecMCTrackHist->SetBinEdges(7,binsPt);\r
\r
// get reconstructed vertex \r
const AliESDVertex* vtxESD = 0; \r
- Bool_t isRecVertex = kTRUE;\r
if(evtCuts->IsRecVertexRequired()) \r
{\r
if(GetAnalysisMode() == AlidNdPtHelper::kTPC) {\r
vtxESD = esdEvent->GetPrimaryVertexTPC();\r
}\r
+ else if(GetAnalysisMode() == AlidNdPtHelper::kTPCITS) {\r
+ vtxESD = esdEvent->GetPrimaryVertexTracks();\r
+ }\r
+ else {\r
+ return;\r
+ }\r
}\r
- Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD) && isRecVertex; \r
+ Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD); \r
//printf("isEventOK %d, isEventTriggered %d \n",isEventOK, isEventTriggered);\r
//printf("GetAnalysisMode() %d \n",GetAnalysisMode());\r
\r
AliESDtrack *track = (AliESDtrack*)allChargedTracks->At(i);\r
if(!track) continue;\r
\r
- //if(!esdTrackCuts->AcceptTrack(track)) continue;\r
+ if(!esdTrackCuts->AcceptTrack(track)) continue;\r
FillHistograms(track, stack);\r
multAll++;\r
}\r
TString tag = "TPC+ITS combine tracking + DCAr(pt) + DCAz(pt)";
}
+ // TPC+ITS combine tracking + DCAr(pt) (2010)
+ if (cutMode == 72)
+ {
+ Int_t minclsTPC=70;
+ Double_t maxchi2perTPCcl=4.;
+ Double_t maxdcazITSTPC=2.0;
+
+ //
+ // TPC
+ //
+ esdTrackCuts->SetRequireTPCRefit(kTRUE);
+ esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
+ esdTrackCuts->SetMinNClustersTPC(minclsTPC);
+ esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
+ //
+ // ITS
+ //
+ esdTrackCuts->SetRequireITSRefit(kTRUE);
+ esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
+ //
+ // primary selection
+ //
+ esdTrackCuts->SetDCAToVertex2D(kFALSE);
+ esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
+ esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
+
+ // DCArphi parametrization (LHC10c pass2)
+ // 7*(0.0026+0.0050/pt^1.01)
+ esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01");
+
+ TString tag = "TPC+ITS combine tracking + DCAr(pt) (2010)";
+ }
+
// TPC-tracks + SPD point + ITS refit + DCAr(pt) 4-sigma
if (cutMode == 75)
{
TString tag = "Adam Kisiel track selection";
}
+ // TPC+ITS refit
+ // for cut studies
+ if (cutMode == 151)
+ {
+ //
+ // TPC
+ //
+ esdTrackCuts->SetRequireTPCRefit(kTRUE);
+ esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
+ //
+ // ITS
+ //
+ esdTrackCuts->SetRequireITSRefit(kTRUE);
+ esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
+ //
+
+ TString tag = "TPC+ITS refit required - for cut studies";
+ }
+
+ // TPC+ITS
+ // for cut studies
+ if (cutMode == 152)
+ {
+ //
+ // TPC
+ //
+ esdTrackCuts->SetRequireTPCRefit(kTRUE);
+ esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
+ //
+ // ITS
+ //
+ //esdTrackCuts->SetRequireITSRefit(kTRUE);
+ //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
+ //
+
+ TString tag = "TPC refit required - for cut studies";
+ }
+
+