TPC module

[u/mrichter/AliRoot.git] / TPC / Upgrade / macros / GetComposedResidualDistortion.C

diff --git a/TPC/Upgrade/macros/GetComposedResidualDistortion.C b/TPC/Upgrade/macros/GetComposedResidualDistortion.C
deleted file mode 100644 (file)
index df0e234..0000000
--- a/TPC/Upgrade/macros/GetComposedResidualDistortion.C
+++ /dev/null
@@ -1,67 +0,0 @@
-#include <TString.h>
-#include <TVectorD.h>
-#include <TObjArray.h>
-#include <TFile.h>
-
-#include <AliTPCComposedCorrection.h>
-#include <AliTPCCorrectionLookupTable.h>
-
-#include <AliToyMCEventGenerator.h>
-
-
-AliTPCComposedCorrection* GetComposedResidualDistortion(TString fluctuationMap, TString averageMap, Bool_t rescale=kTRUE)
-{
-  //
-  //
-  //
-
-
-  TFile fFluct(fluctuationMap);
-  AliTPCCorrectionLookupTable *corrFluct = (AliTPCCorrectionLookupTable*)fFluct.Get("map");
-  fFluct.Close();
-
-  TFile fAverage(averageMap);
-  AliTPCCorrectionLookupTable *corrAverage = (AliTPCCorrectionLookupTable*)fAverage.Get("map");
-  fAverage.Close();
-
-  TObjArray *arrMaps = new TObjArray(2);
-  // !!!!! 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.;
-  if (rescale) {
-    Float_t dummy=0;
-    weights(1)=-AliToyMCEventGenerator::GetSCScalingFactor(corrFluct, corrAverage,dummy);
-  }
-  residualDistortion->SetWeights(&weights);
-
-  return residualDistortion;
-}
-
-AliTPCCorrectionLookupTable* GetResidualTable(TString fluctuationMap, TString averageMap, Bool_t rescale=kTRUE)
-{
-  TFile fFluct(fluctuationMap);
-  AliTPCCorrectionLookupTable *corrFluct = (AliTPCCorrectionLookupTable*)fFluct.Get("map");
-  fFluct.Close();
-  
-  TFile fAverage(averageMap);
-  AliTPCCorrectionLookupTable *corrAverage = (AliTPCCorrectionLookupTable*)fAverage.Get("map");
-  fAverage.Close();
-
-  Double_t scale=AliToyMCEventGenerator::GetSCScalingFactor(corrFluct, corrAverage,dummy);
-
-  corrAverage->SetCorrScaleFactor(scale);
-
-  AliTPCCorrectionLookupTable *tab=new AliTPCCorrectionLookupTable;
-  tab->CreateResidual(corrFluct, corrAverage);
-  
-}
\ No newline at end of file