Double_t cr = ar->GetBinCenter(ir+1);
Double_t cphi = aphi->GetBinCenter(iphi+1);
Double_t cz = az->GetBinCenter(iz+1);
- hProj->SetNameTitle(Form("h_%02d_%02d_%d02",iz+1, iphi+1, ir+1),
- Form("z,phi,r: %02d,%02d,%d02 (%.2f, %.2f, %.2f)",iz+1,iphi+1,ir+1, cr, cphi, cz )
+ hProj->SetNameTitle(Form("h_%02d_%02d_%02d",iz+1, iphi+1, ir+1),
+ Form("z,phi,r: %02d,%02d,%02d (%.2f, %.2f, %.2f)",iz+1,iphi+1,ir+1, cz, cphi, cr )
);
hProj->Fit(&fg,"QR");
arrFits.Add(hProj);
fAverage.Close();
TObjArray *arrMaps = new TObjArray(2);
- arrMaps->Add(corrFluct); // distortion with the fluctuation Map
+ // !!!!! In AliTPCComposedCorrection::GetDistortion MakeInverseIterator is called !!!!
+ // for this reason we have to add the maps in the wrong order
+
arrMaps->Add(corrAverage); // correction with the average Map
+ arrMaps->Add(corrFluct); // distortion with the fluctuation Map
// create the composed correction
// if the weight are set to +1 and -1, the first map will be responsible for the distortions
// The second map for the corrections
AliTPCComposedCorrection *residualDistortion = new AliTPCComposedCorrection(arrMaps, AliTPCComposedCorrection::kQueueResidual);
TVectorD weights(2);
- weights(0)=+1.;
- weights(1)=-1.;
+ weights(1)=+1.;
+ weights(0)=-1.;
if (rescale) {
Float_t dummy=0;
- weights(1)=-AliToyMCEventGenerator::GetSCScalingFactor(corrFluct, corrAverage,dummy);
+ weights(0)=-AliToyMCEventGenerator::GetSCScalingFactor(corrFluct, corrAverage,dummy);
}
residualDistortion->SetWeights(&weights);