do_eta_fit_emc.C

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file do_eta_fit_emc.C

#include <GlobalVariables.C>

#include "litecaloeval/LiteCaloEval.h"
R__LOAD_LIBRARY(libLiteCaloEvalTowSlope.so)

// File parameters - not modified histos, modified histos, output
void do_eta_fit_emc(const char * reffile = "condor/combine_out/out_sum.root", const char * infile="condor/combine_out/out_23746.root", const char * modfile="calib.root")
{

  TStopwatch t;
  t.Start();


  gSystem->Load("libLiteCaloEvalTowSlope.so");

  LiteCaloEval reflce, modlce;


  reflce.CaloType(LiteCaloEval::CEMC);
  modlce.CaloType(LiteCaloEval::CEMC);


  reflce.Get_Histos(reffile);

  modlce.Get_Histos(infile,modfile);

  modlce.setFitMin(0.12);
  modlce.setFitMax(0.7);


  /*Fit histos,obtain fit values. Second parameter determines smoothing/include phi loop (individual towers). 
    If set to e.g. 10 - do 1 smoothing and run over phi. 
    If set to 11, do one smoothing and dont do phi loop. Can do just 0 or 1. Which would be no smoothing, and 
    include phi if 0, and dont include if 1.
  */
  modlce.FitRelativeShifts(&reflce,10);
  
  t.Stop();
  t.Print();

}