d2eeaf60ab2c8036fd2dae0684f41f7d060bf034
[u/mrichter/AliRoot.git] / PWGDQ / dielectron / TakuAlberica / pair / run_all.C
1 void run_all(int id=0, int end=50, int pro=0, int icut=0, int itrig=0){
2
3   gStyle->SetPalette(1);
4   gStyle->SetOptStat(0);
5   gStyle->SetOptFit(111111);
6   gStyle->SetOptFit(0);
7   gStyle->SetOptTitle(0);
8   //gStyle->SetFillColor(10);
9   gStyle->SetCanvasColor(10);
10   gStyle->SetFrameBorderMode(0);
11   gStyle->SetFrameFillColor(0);
12   gStyle->SetCanvasColor(0);
13   gStyle->SetPadBorderSize(0);
14   gStyle->SetCanvasBorderSize(0);
15   //gStyle->SetPadLeftMargin(0.15);
16   gStyle->SetPadLeftMargin(0.125);
17   gStyle->SetPadBottomMargin(0.125);
18   gStyle->SetPadTopMargin(0.1);
19   gStyle->SetTitleYOffset(1.3);
20   //gStyle->SetPadLeftMargin(0.1);
21   cout<<"physics is always fun! "<<endl; 
22
23
24   gSystem->Load("libCore");// no
25   gSystem->Load("libTree");
26   gSystem->Load("libGeom");
27   gSystem->Load("libVMC");
28   gSystem->Load("libXMLIO");// no
29   gSystem->Load("libPhysics");
30   gSystem->Load("libXMLParser");
31   gSystem->Load("libProof");
32   gSystem->Load("libMinuit");
33
34
35   gSystem->Load("libana_sgl.so");
36   ana_sgl *ana=new ana_sgl();
37
38   char outputname[100];
39   sprintf(outputname,"ntpair_rev2_%d_seg%d_%d_cut%d_trig%d_v2.root", pro, id,end, icut, itrig);
40
41
42   ana->ana_init(outputname);
43
44  //ana->set_tof_cuts(-3,3);
45   //ana->enable_pair_emc_cut(0.7, 1.3);
46   //ana->enable_pait_pt_cut(0.4, 20);
47
48   if(pro==1){
49     ana->set_veto_for_proton(-2, 3.5);
50     ana->set_veto_for_kaon(-2, 3.5);
51   }
52
53   //0->MB, 1->SemiCent, 2-->Cent
54   if(itrig==0){
55     ana->select_trigger(2);
56   }else if(itrig==1){
57     ana->select_trigger(1);
58   }
59
60   if(icut==0){
61     ana->reject_conversion(true);
62     //ana->set_tpc_dedx_cuts(75, 90);
63     ana->set_tpc_dedx_cuts(65, 90);
64     ana->enable_pair_phiv_cut(0.6);
65   }else if(icut==1){
66     //ana->reject_conversion(true);
67     //ana->set_tpc_dedx_cuts(80, 85);
68     //ana->set_tpc_dedx_cuts(68, 75);
69     ana->set_tpc_dedx_cuts(65, 90);
70     ana->enable_pair_phiv_cut(0.6);
71   }else if(icut==2){
72     ana->reject_conversion(true);
73     ana->set_tpc_dedx_cuts(75, 90);
74     ana->enable_pair_phiv_cut(0.6);
75   }else if(icut==3){
76     ana->reject_conversion(true);
77     ana->set_tpc_dedx_cuts(80, 85);
78     ana->enable_pair_phiv_cut(0.6);
79   }
80
81   ana->print_cuts();
82
83   if(pro==0){
84     ifstream f("list_rev1_0.all");
85     int nline = 1418;
86   }else if(pro==1){
87     ifstream f("list_rev2_6.all");
88     int nline = 50;
89   }else if(pro==2){
90     ifstream f("list_rev2_7.all");
91     int nline = 50;
92   }
93
94
95   char inputfilename[100];
96   for(int i=0;i<nline;i++){
97     f >> inputfilename; 
98     if(i>=id && i<end){
99       sprintf(inputfilename,"%s",inputfilename);
100       ana->loop_a_file(inputfilename);
101       cout<<" end of run : "<<i<<endl;
102     }
103   }
104
105   ana->ana_end();
106
107 }
108