calo/tower_slope_run23/rings_do_eta_fit2_emc.C

0001 #include <GlobalVariables.C>
0002
0003 #include "litecaloeval/LiteCaloEval.h"
0004 R__LOAD_LIBRARY(libLiteCaloEvalTowSlope.so)
0005
0006 // regular file, modified file, output file
0007 void rings_do_eta_fit2_emc(const char * infile, const char * modfile)
0008 {
0009   TStopwatch t;
0010   t.Start();
0011
0012   gSystem->Load("libLiteCaloEvalTowSlope.so");
0013   LiteCaloEval modlce;
0014
0015   //reflce.CaloType(LiteCaloEval::CEMC);
0016   modlce.CaloType(LiteCaloEval::CEMC);
0017
0018   /// we really dont need to retrieve histos here bc we arent using ref object
0019   /// reflce.Get_Histos(reffile);
0020
0021   modlce.Get_Histos(infile,modfile);
0022
0023   modlce.setFitMin(0.12);
0024   modlce.setFitMax(0.7);
0025
0026   modlce.FitRelativeShifts(&modlce,110);
0027
0028   t.Stop();
0029   t.Print();
0030
0031 }