DrawCorrELoss.C

Go to the documentation of this file.

void Setup(Bool_t compile)
{
  const char* post = (compile ? "+g" : "");
  const char* fwd  = "$ALICE_PHYSICS/PWGLF/FORWARD/analysis2";
  if (!gROOT->GetClass("AliFMDCorrELossFit"))
    gROOT->Macro(Form("%s/scripts/LoadLibs.C", fwd));
  gSystem->AddIncludePath(Form("-I%s/scripts -I%s/corrs -I%s "
             "-I$ALICE_PHYSICS/include "
             "-I$ALICE_ROOT/include",
             fwd, fwd, fwd));
  gROOT->LoadMacro(Form("%s/scripts/SummaryDrawer.C%s", fwd, post));
  gROOT->LoadMacro(Form("%s/corrs/CorrDrawer.C%s", fwd, post));

}
void
DrawCorrELoss(ULong_t runNo, UShort_t sys, UShort_t sNN, Short_t field,
        Bool_t mc=false, Bool_t sat=false, 
        const char* fname=0, Bool_t details=true)
{
  // Load libraries and object 
  // const char* fwd = "$ALICE_PHYSICS/PWGLF/FORWARD/analysis2";
  Setup(false);

  CorrDrawer d;
  d.Summarize(AliForwardCorrectionManager::kELossFits, runNo, sys, sNN, field, 
        mc, sat, "", fname);
}
void
DrawCorrELoss(Bool_t      mc, 
        const char* file="forward_eloss.root", 
        const char* local="fmd_corrections.root")
{
  Setup(true);

  CorrDrawer::Summarize(AliForwardCorrectionManager::kELossFits, 
      mc, file, local, "LANDSCAPE FEW");
}
void 
DrawCorrELoss(Bool_t mc, Bool_t dummy, 
        const char* file="forward_eloss_rerun.root")
{
  Setup(true);

  Printf("Drawing fit results from %s", file);
  TFile* hist = TFile::Open(file, "READ");
  if (!hist) { 
    Error("DrawCorrELoss", "Failed to open %s", file);
    return;
  }
  TList* res = static_cast<TList*>(hist->Get("ForwardELossResults"));
  if (!res) { 
    Error("DrawCorrEloss", "Failed to get ForwardELossResults from %s", file);
    return;
  }
  TList* ef = static_cast<TList*>(res->FindObject("fmdEnergyFitter"));
  if (!ef) { 
    Error("DrawCorrEloss", "Failed to get fmdEnergyFitter from %s:/%s", 
    file, res->GetName());
    return;
  }

  AliFMDCorrELossFit* fits = 
    static_cast<AliFMDCorrELossFit*>(ef->FindObject("AliFMDCorrELossFit"));
  if (!fits) {
    Error("DrawCorrEloss", "Failed to get AliFMDCorrELossFit from %s:/%s/%s", 
    file, res->GetName(), ef->GetName());
    return;
  }
  
  CorrDrawer* cd = new CorrDrawer;
  cd->fELossExtra = file;
  cd->fMinQuality = AliFMDCorrELossFit::kDefaultQuality;
  // cd->fLandscape = true;
  cd->Summarize(const_cast<const AliFMDCorrELossFit*>(fits), true);
}

//
// EOF
//