calo/tower_slope_run23/do_eta_fit2_ihcal.C

0001 #include <GlobalVariables.C>
0002
0003 #include "LiteCaloEval.h"
0004 R__LOAD_LIBRARY(libLiteCaloEvalTowSlope.so)
0005
0006 /// typical parameters, despite the naming here, is modified histo, non modified, output file
0007 void do_eta_fit2_ihcal(const char * reffile, const char * infile, const char * modfile)
0008 {
0009
0010   TStopwatch t;
0011   t.Start();
0012
0013   gSystem->Load("libLiteCaloEvalTowSlope.so");
0014
0015
0016   LiteCaloEval reflce, modlce;
0017
0018   reflce.CaloType(LiteCaloEval::HCALIN);
0019   modlce.CaloType(LiteCaloEval::HCALIN);
0020
0021   reflce.Get_Histos(reffile); // get histos for modified files
0022   modlce.Get_Histos(infile,modfile); // get histos for non mod files
0023
0024   /// we want to fit our nonmod/reference histogram and then
0025   /// adjust the modified histo according to such fit
0026   modlce.setFitMin(0.1);
0027   modlce.setFitMax(0.4);
0028
0029   /// Modlce holds histos with no energy modifications. What we then want to
0030   /// do is give the modified files, fit them, and compare with the nonmod
0031   modlce.FitRelativeShifts(&reflce,10);
0032
0033   t.Stop();
0034   t.Print();
0035
0036 }