- containerInput[0] = vectorMC[0];
- containerInput[1] = vectorMC[1] ;
- containerInput[2] = vectorMC[2] ;
- containerInput[3] = vectorMC[3] ;
- containerInput[4] = vectorMC[4] ;
- containerInput[5] = vectorMC[5] ; // in micron
- containerInput[6] = 0.; // dummy value, meaningless in MC, in micron
- containerInput[7] = 0.; // dummy value, meaningless in MC, in micron
- containerInput[8] = 0.; // dummy value, meaningless in MC, in micron
- containerInput[9] = -100000.; // dummy value, meaningless in MC, in micron^2
- containerInput[10] = 1.01; // dummy value, meaningless in MC
- containerInput[11] = vectorMC[6]; // dummy value, meaningless in MC
- fCFManager->GetParticleContainer()->Fill(containerInput,kStepGenerated);
+ containerInputMC[0] = vectorMC[0];
+ containerInputMC[1] = vectorMC[1] ;
+ containerInputMC[2] = vectorMC[2] ;
+ containerInputMC[3] = vectorMC[3] ;
+ containerInputMC[4] = vectorMC[4] ;
+ containerInputMC[5] = vectorMC[5] ; // in micron
+ containerInputMC[6] = 0.; // dummy value, meaningless in MC, in micron
+ containerInputMC[7] = 0.; // dummy value, meaningless in MC, in micron
+ containerInputMC[8] = 0.; // dummy value, meaningless in MC, in micron
+ containerInputMC[9] = -100000.; // dummy value, meaningless in MC, in micron^2
+ containerInputMC[10] = 1.01; // dummy value, meaningless in MC
+ containerInputMC[11] = vectorMC[6]; // dummy value, meaningless in MC
+ containerInputMC[12] = zMCVertex; // z of reconstructed of primary vertex
+ if (fUseWeight) fWeight = GetWeight(vectorMC[0]); // setting the weight according to the function defined in AliCFHeavyFlavourTaskMultiVarMultiStep::GetWeight(Float_t pt)
+ AliDebug(3,Form("weight = %f",fWeight));
+ if (!fCuts->IsInFiducialAcceptance(vectorMC[0],vectorMC[1])) continue;
+ if (TMath::Abs(vectorMC[1]) < 0.5) {
+ fCFManager->GetParticleContainer()->Fill(containerInputMC,kStepGeneratedLimAcc,fWeight);
+ }
+ fCFManager->GetParticleContainer()->Fill(containerInputMC,kStepGenerated,fWeight);