test/MuonSeedPTAnalysis/RunAllPara.C

0001 #include "TROOT.h"
0002 #include "TSystem.h"
0003 #include "TApplication.h"
0004
0005 void testPara(int type, int st =0, double h1=0, double h2=0, int np=0, int mode=0 ) {
0006
0007   gROOT->ProcessLine(".L SeedPtFunction.C+");
0008   gROOT->ProcessLine(".L SeedParaFit.C+");
0009
0010   SeedParaFit *fitter = new SeedParaFit() ;
0011
0012   if ( mode == 1 ) fitter->ParaFit(type,st,h1,h2,np);
0013   // CSC Pair

0014   if ( mode != 1 ) {
0015      if ( type == 0 ) fitter->PrintTitle();
0016      if ( type == 1 || type == 0 ) {
0017         fitter->ParaFit(1, 11, 1.57, 1.67, 1);
0018     fitter->ParaFit(1, 12, 1.00, 1.55, 3);
0019     fitter->ParaFit(1, 12, 1.72, 2.40, 3);
0020     fitter->ParaFit(1, 13, 1.20, 1.55, 3);
0021     fitter->ParaFit(1, 13, 1.76, 2.40, 3);
0022     fitter->ParaFit(1, 14, 1.80, 2.40, 3);
0023     fitter->ParaFit(1, 23, 1.15, 2.38, 3);
0024     fitter->ParaFit(1, 24, 1.85, 2.38, 3);
0025     fitter->ParaFit(1, 34, 1.85, 2.38, 3);
0026     }
0027     // DT Pair

0028     if ( type == 3  || type == 0 ) {
0029        fitter->ParaFit(3, 12, 0., 1.02, 3);
0030        fitter->ParaFit(3, 13, 0., 0.88, 3);
0031        fitter->ParaFit(3, 14, 0., 0.78, 3);
0032        fitter->ParaFit(3, 23, 0., 0.88, 3);
0033        fitter->ParaFit(3, 24, 0., 0.78, 3);
0034        fitter->ParaFit(3, 34, 0., 0.78, 3);
0035     }
0036
0037     // OL Pair

0038     if ( type == 2  || type == 0 ) {
0039        fitter->ParaFit(2,1213, 0.95, 1.10, 2);
0040        fitter->ParaFit(2,1222, 1.02, 1.16, 2);
0041        fitter->ParaFit(2,1232, 1.08, 1.18, 2);
0042        fitter->ParaFit(2,2213, 0.92, 1.02, 2);
0043        fitter->ParaFit(2,2222, 1.00, 1.06, 1);
0044     }
0045     // CSC Single

0046     if ( type == 5 || type == 0 ) {
0047        fitter->ParaFit(5, 11, 1.60, 2.38, 3);
0048        fitter->ParaFit(5, 12, 1.22, 1.62, 2);
0049        fitter->ParaFit(5, 13, 0.94, 1.10, 2);
0050        fitter->ParaFit(5, 21, 1.62, 2.40, 3);
0051        fitter->ParaFit(5, 22, 1.02, 1.56, 3);
0052        fitter->ParaFit(5, 31, 1.75, 2.40, 3);
0053        fitter->ParaFit(5, 32, 1.12, 1.66, 3);
0054        fitter->ParaFit(5, 41, 1.82, 2.38, 3);
0055     }
0056     // DT Single

0057     if ( type == 4 || type == 0 ) {
0058        fitter->ParaFit(4, 10, 0.0 , 0.26, 2);
0059        fitter->ParaFit(4, 11, 0.32, 0.78, 3);
0060        fitter->ParaFit(4, 12, 0.85, 1.16, 2);
0061        fitter->ParaFit(4, 20, 0.0 , 0.22, 2);
0062        fitter->ParaFit(4, 21, 0.30, 0.66, 2);
0063        fitter->ParaFit(4, 22, 0.78, 1.02, 2);
0064        fitter->ParaFit(4, 30, 0.0 , 0.18, 2);
0065        fitter->ParaFit(4, 31, 0.24, 0.58, 2);
0066        fitter->ParaFit(4, 32, 0.62, 0.88, 2);
0067     }
0068     if ( type == 0 ) fitter->PrintEnd();
0069   }
0070   //gROOT->Reset();

0071
0072 }