o same here In AliTPCComposedCorrection::GetDistortion MakeInverseIterator is called
authorwiechula <wiechula@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 14 Nov 2013 17:27:31 +0000 (17:27 +0000)
committerwiechula <wiechula@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 14 Nov 2013 17:27:31 +0000 (17:27 +0000)
TPC/Upgrade/macros/AnaDelta.C

index d11e26b33d82ac4c88663c5effd709b636d76dfa..6f71a78e1fd46687fc11f9a8fa0a2bb114ec61ab 100644 (file)
@@ -228,17 +228,20 @@ void AnaDeltaResiduals(TString fluctuationMap, TString averageMap, TString outFi
   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(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
+  // !!!!! 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);
-  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.);