fKeepingOnlyPYTHIABkg(kFALSE),
fHistoMCLcK0SpGen(0x0),
fHistoMCLcK0SpGenAcc(0x0),
- fHistoMCLcK0SpGenLimAcc(0x0)
+ fHistoMCLcK0SpGenLimAcc(0x0),
+ fTriggerMask(0)
{
//
// Default ctor
fKeepingOnlyPYTHIABkg(kFALSE),
fHistoMCLcK0SpGen(0x0),
fHistoMCLcK0SpGenAcc(0x0),
- fHistoMCLcK0SpGenLimAcc(0x0)
-
+ fHistoMCLcK0SpGenLimAcc(0x0),
+ fTriggerMask(0)
{
//
// Constructor. Initialization of Inputs and Outputs
if ( !fUseMCInfo && fIspA) {
fAnalCuts->SetTriggerClass("");
- fAnalCuts->SetTriggerMask(AliVEvent::kINT7);
+ fAnalCuts->SetTriggerMask(fTriggerMask);
}
Int_t runnumber = aodEvent->GetRunNumber();
AliError("AliAnalysisTaskSELc2V0bachelorTMVA::UserExec: MC header branch not found!\n");
return;
}
+
+ Double_t zMCVertex = mcHeader->GetVtxZ();
+ if (TMath::Abs(zMCVertex) > fAnalCuts->GetMaxVtxZ()){
+ AliDebug(3,Form("z coordinate of MC vertex = %f, it was required to be within [-%f, +%f], skipping event", zMCVertex, fAnalCuts->GetMaxVtxZ(), fAnalCuts->GetMaxVtxZ()));
+ AliInfo(Form("z coordinate of MC vertex = %f, it was required to be within [-%f, +%f], skipping event", zMCVertex, fAnalCuts->GetMaxVtxZ(), fAnalCuts->GetMaxVtxZ()));
+ return;
+ }
+
//Printf("Filling MC histo");
FillMCHisto(mcArray);
}
-
+
// AOD primary vertex
fVtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
if (!fVtx1) return;