o same here In AliTPCComposedCorrection::GetDistortion MakeInverseIterator is called

diff --git a/TPC/Upgrade/macros/AnaDelta.C b/TPC/Upgrade/macros/AnaDelta.C
index d11e26b33d82ac4c88663c5effd709b636d76dfa..6f71a78e1fd46687fc11f9a8fa0a2bb114ec61ab 100644 (file)
--- a/TPC/Upgrade/macros/AnaDelta.C
+++ b/TPC/Upgrade/macros/AnaDelta.C
@@ -228,17 +228,20 @@ void AnaDeltaResiduals(TString fluctuationMap, TString averageMap, TString outFi
   fAverage.Close();
   
   TObjArray *arrMaps = new TObjArray(2);
   fAverage.Close();
   
   TObjArray *arrMaps = new TObjArray(2);

-  arrMaps->Add(corrFluct);   // distortion with the fluctuation Map

   arrMaps->Add(corrAverage); // correction with the average Map
   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
   
   // 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

+  // !!!!! In AliTPCComposedCorrection::GetDistortion MakeInverseIterator is called !!!!
+  // for this reason we have to add the maps in the wrong order
+  

   AliTPCComposedCorrection *residualDistortion = new AliTPCComposedCorrection(arrMaps, AliTPCComposedCorrection::kQueueResidual);
   Float_t dummy=0;
   TVectorD weights(2);
   AliTPCComposedCorrection *residualDistortion = new AliTPCComposedCorrection(arrMaps, AliTPCComposedCorrection::kQueueResidual);
   Float_t dummy=0;
   TVectorD weights(2);

-  weights(0)=+1.;
-  weights(1)=-AliToyMCEventGenerator::GetSCScalingFactor(corrFluct, corrAverage,dummy);
+  weights(1)=+1.;
+  weights(0)=-AliToyMCEventGenerator::GetSCScalingFactor(corrFluct, corrAverage,dummy);

   residualDistortion->SetWeights(&weights);
 
   TVectorD *vR   = MakeLinBinning(10,86.,250.);
   residualDistortion->SetWeights(&weights);
 
   TVectorD *vR   = MakeLinBinning(10,86.,250.);