calo/tower_slope_run23/do_eta_fit2_ohcal.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 do_eta_fit2_ohcal(const char * reffile, const char * infile, const char * modfile)
0008 {
0009   TStopwatch t;
0010   t.Start();
0011
0012   gSystem->Load("libLiteCaloEvalTowSlope.so");
0013   LiteCaloEval reflce, modlce;
0014
0015   reflce.CaloType(LiteCaloEval::HCALOUT);
0016   modlce.CaloType(LiteCaloEval::HCALOUT);
0017
0018   reflce.Get_Histos(reffile);
0019   modlce.Get_Histos(infile,modfile);
0020
0021   modlce.setFitMin(0.12);
0022   modlce.setFitMax(0.7);
0023
0024   modlce.FitRelativeShifts(&reflce,10);
0025
0026   t.Stop();
0027   t.Print();
0028
0029 }