// Fills QA histograms
- Short_t vCharge;
Double_t vEta;
- Double_t vY;
Double_t vPhi;
Double_t vPt;
Double_t vDCAconstrainedxy;
}
// track parameters
- vCharge = aodTrack->Charge();
vEta = aodTrack->Eta();
- vY = aodTrack->Y();
vPhi = aodTrack->Phi();// * TMath::RadToDeg();
vPt = aodTrack->Pt();
vDCAconstrainedxy = aodTrack->DCA();
labelArray.AddAt(label,labelCounter);
labelCounter += 1;
- Bool_t iFound = kFALSE;
Int_t mcGoods = nMCLabelCounter;
for (Int_t k = 0; k < mcGoods; k++) {
Int_t mcLabel = labelMCArray.At(k);
- iFound = kFALSE;
if (mcLabel != TMath::Abs(label)) continue;
if(mcLabel != label) continue;
const Float_t kLimit = 0.02 * 3;
Float_t dphistarminabs = 1e5;
- Float_t dphistarmin = 1e5;
if (TMath::Abs(dphistar1) < kLimit || TMath::Abs(dphistar2) < kLimit || dphistar1 * dphistar2 < 0 )
{
if (dphistarabs < dphistarminabs)
{
- dphistarmin = dphistar;
dphistarminabs = dphistarabs;
}
}
const Float_t kLimit = 0.02 * 3;
Float_t dphistarminabs = 1e5;
- Float_t dphistarmin = 1e5;
if (TMath::Abs(dphistar1) < kLimit || TMath::Abs(dphistar2) < kLimit || dphistar1 * dphistar2 < 0 )
{
if (dphistarabs < dphistarminabs)
{
- dphistarmin = dphistar;
dphistarminabs = dphistarabs;
}
}
//____________________________________________________________________//
void AliBalanceEventMixing::PrintResults(Int_t iAnalysisType, TH1D *gHistBalance) {
//Prints the calculated width of the BF and its error
- Double_t x[MAXIMUM_NUMBER_OF_STEPS];
Double_t gSumXi = 0.0, gSumBi = 0.0, gSumBiXi = 0.0;
Double_t gSumBiXi2 = 0.0, gSumBi2Xi2 = 0.0;
Double_t gSumDeltaBi2 = 0.0, gSumXi2DeltaBi2 = 0.0;
cout<<"=================================================="<<endl;
for(Int_t i = 1; i <= fNumberOfBins[iAnalysisType]; i++) {
- x[i-1] = fP2Start[iAnalysisType] + fP2Step[iAnalysisType]*i + fP2Step[iAnalysisType]/2;
//cout<<"B: "<<gHistBalance->GetBinContent(i)<<"\t Error: "<<gHistBalance->GetBinError(i)<<"\t bin: "<<gHistBalance->GetBinCenter(i)<<endl;
}
//cout<<"=================================================="<<endl;
binVertexMin = fHistPN->GetGrid(0)->GetGrid()->GetAxis(5)->FindBin(vertexZMin);
binVertexMax = fHistPN->GetGrid(0)->GetGrid()->GetAxis(5)->FindBin(vertexZMax-0.00001);
}
- Double_t binPsiLowEdge = 0.;
- Double_t binPsiUpEdge = 0.;
- Double_t binVertexLowEdge = 0.;
- Double_t binVertexUpEdge = 0.;
TH2D* h1 = NULL;
TH2D* h2 = NULL;
cout<<"In the balance function (2D->1D) loop: "<<iBinPsi<<" (psiBin), "<<iBinVertex<<" (vertexBin) "<<endl;
- // determine the bin edges for this bin
- binPsiLowEdge = fHistPN->GetGrid(0)->GetGrid()->GetAxis(0)->GetBinLowEdge(iBinPsi);
- binPsiUpEdge = fHistPN->GetGrid(0)->GetGrid()->GetAxis(0)->GetBinUpEdge(iBinPsi);
- if(fVertexBinning){
- binVertexLowEdge = fHistPN->GetGrid(0)->GetGrid()->GetAxis(5)->GetBinLowEdge(iBinVertex);
- binVertexUpEdge = fHistPN->GetGrid(0)->GetGrid()->GetAxis(5)->GetBinUpEdge(iBinVertex);
- }
-
// Psi_2
fHistP->GetGrid(0)->GetGrid()->GetAxis(0)->SetRangeUser(psiMin,psiMax-0.00001);
fHistN->GetGrid(0)->GetGrid()->GetAxis(0)->SetRangeUser(psiMin,psiMax-0.00001);