test/MuonSeedPTAnalysis/RunAllScale.C

0001 #include "TROOT.h"
0002 #include "TSystem.h"
0003 #include "TApplication.h"
0004
0005 void RunScale(int type, int st =0, double h1=0, double h2=0, int idx=0, int np=3, int mode=0 ) {
0006
0007   // type = 0 && mode = 0   : do all scaling

0008   // mode = 1               : scaling for a specific case

0009   // mode = 0 && type != 0  : scaling for whole CSC/DT

0010
0011   gROOT->ProcessLine(".L SeedPtFunction.C+");
0012   gROOT->ProcessLine(".L SeedPtScale.C+");
0013   SeedPtScale *scaler = new SeedPtScale() ;
0014
0015   if ( mode != 0 ) scaler->PtScale(type, st, h1, h2, idx, np );
0016
0017   if ( mode == 0 ) {
0018      np = 2;
0019      if ( type  == 1 || type == 0 ) {
0020         scaler->PtScale(1,01, 1.55, 1.65, 1, np);
0021         scaler->PtScale(1,12, 1.20, 1.30, 1, np);
0022         scaler->PtScale(1,12, 1.70, 1.90, 2, np);
0023         scaler->PtScale(1,12, 1.70, 1.90, 3, np);
0024         scaler->PtScale(1,13, 1.30, 1.40, 2, np);
0025         scaler->PtScale(1,13, 1.80, 2.00, 3, np);
0026         scaler->PtScale(1,14, 1.80, 2.00, 3, np);
0027         scaler->PtScale(1,23, 1.20, 1.25, 1, np);
0028         scaler->PtScale(1,23, 1.80, 2.00, 2, np);
0029         scaler->PtScale(1,24, 1.85, 2.05, 1, np);
0030         scaler->PtScale(1,34, 1.82, 2.02, 1, np);
0031         cout<<" CSC DONE!!"<<endl;
0032      }
0033      if ( type  == 2 || type == 0 ) {
0034         scaler->PtScale(2,1213, 0.90, 0.95, 0, np);
0035         scaler->PtScale(2,1222, 1.04, 1.09, 0, np);
0036         scaler->PtScale(2,1232, 1.1 , 1.2 , 0, np);
0037         scaler->PtScale(2,2213, 0.9 , 0.96, 0, np);
0038         scaler->PtScale(2,2222, 1.0 , 1.04, 0, np);
0039      }
0040      if ( type  == 3 || type == 0 ) {
0041         scaler->PtScale(3,12, 0.0 , 0.1 , 1, np);
0042         scaler->PtScale(3,12, 0.85, 0.9 , 2, np);
0043         scaler->PtScale(3,13, 0.0 , 0.1 , 1, np);
0044         scaler->PtScale(3,13, 0.65, 0.75, 2, np);
0045         scaler->PtScale(3,14, 0.0 , 0.1 , 1, np);
0046         scaler->PtScale(3,14, 0.55, 0.65, 2, np);
0047         scaler->PtScale(3,23, 0.0 , 0.1 , 1, np);
0048         scaler->PtScale(3,23, 0.75, 0.85, 2, np);
0049         scaler->PtScale(3,24, 0.0 , 0.1 , 1, np);
0050         scaler->PtScale(3,24, 0.58, 0.68, 2, np);
0051         scaler->PtScale(3,34, 0.0,  0.1,  1, np);
0052         scaler->PtScale(3,34, 0.65, 0.75, 2, np);
0053      }
0054      if ( type  == 4 || type == 0 ) {
0055         scaler->PtScale(4,10, 0.00, 0.20, 0, np);
0056         scaler->PtScale(4,11, 0.40, 0.60, 0, np);
0057         scaler->PtScale(4,12, 0.85, 0.95, 0, np);
0058         scaler->PtScale(4,20, 0.00, 0.20, 0, np);
0059         scaler->PtScale(4,21, 0.40, 0.60, 0, np);
0060         scaler->PtScale(4,22, 0.72, 0.82, 0, np);
0061         scaler->PtScale(4,30, 0.00, 0.18, 0, np);
0062         scaler->PtScale(4,31, 0.30, 0.50, 0, np);
0063         scaler->PtScale(4,32, 0.60, 0.80, 0, np);
0064      }
0065      if ( type  == 5 || type == 0 ) {
0066         scaler->PtScale(5,11, 1.78, 1.83, 0, np);
0067         scaler->PtScale(5,12, 1.36, 1.44, 0, np);
0068         scaler->PtScale(5,13, 1.03, 1.09, 0, np);
0069         scaler->PtScale(5,21, 1.82, 2.00, 0, np);
0070         scaler->PtScale(5,22, 1.28, 1.36, 0, np);
0071      }
0072
0073
0074   }
0075   gROOT->Reset();
0076
0077 }