calib_nov23/towerslope/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.m_myminbin =  8;       // this really is min bin (no -1/+1 needed)
0024                                 //  but not sure if for summer 23 this should
0025                                 // be 8 or 9 I think 8 according to Joey Clement
0026                                 // on ~Nov 1, mattermost post
0027   modlce.m_mymaxbin =  95 + 1 ;  // this actually means maxbin is 1 less
0028   //  modlce.m_mymaxbin =  95+1;  // so that means maxbin is 95 which is
0029   //  really the max eta bin for emc, indexed 0 through 95
0030
0031
0032   modlce.setFitMin(0.12);
0033   modlce.setFitMax(0.7);
0034
0035   modlce.FitRelativeShifts(&modlce,110);
0036
0037   t.Stop();
0038   t.Print();
0039
0040 }