StRoot: StEEmcPool/macros/LCP/plCut.C Source File

00001 TFile *fd=0;
00002 
00003 plCut(TString run="R305", TString path="./wrkLcpX/") {
00004   gStyle->SetPalette(1,0);
00005   // gStyle->SetOptStat(0);
00006 
00007     //  run="R616";
00008 
00009   TString fname=path+"/"+run;  
00010   //   fname+="-D";  fname+=dca;
00011   fname+=".tree.root";
00012 
00013   // TString run1="bXoff/R3001002.hist.root";
00014 
00015   fd=new TFile(fname);
00016 
00017   if(!fd->IsOpen()) {
00018     printf("plCut: %s  not exits \n",fname.Data());
00019     return;
00020   }
00021 
00022   printf("plCut: %s  opened \n",fname.Data());
00023   fd->ls();
00024 
00025   ln = new TLine(10,260,40,260);
00026   ln = new TLine(0,340,4,340);
00027 
00028   ln->SetLineStyle(3);
00029 
00030   //  TString cname="maxDCA=";  cname+=dca;
00031   //TString cname="maxEta=";  cname+=eta;
00032   TString cname=run;
00033 
00034   TCanvas *c1=new TCanvas(cname,cname,600,350);
00035   TH1F* h=0;
00036 
00037   //  h=(TH1F*)fd->Get("pT"); h->Draw();   return;
00038      
00039 
00040   c1->Divide(2,1);
00041   c1->cd(1);
00042   mL=new TLegend(0.15,.5,.7,.7);
00043 
00044   // fd->ls();
00045   h=(TH1F*)fd->Get("CnFP-A"); 
00046   h->Draw();  mL->AddEntry(h,"#Delta#phi <1/8 #pi (data X 1/10) ");
00047   h->Scale(0.1);
00048   
00049   //  h->SetMaximum(1401);
00050   //  h->SetMaximum(1800);
00051  
00052   h=(TH1F*)fd->Get("CnFP-C"); h->SetLineColor(kBlue);
00053   h->Draw("same");  mL->AddEntry(h," #Delta#phi #in [ 3/8 #pi, 5/8 #pi ] ");
00054 
00055   h=(TH1F*)fd->Get("CnFP-E"); h->SetLineColor(kRed);
00056   h->Draw("same");  mL->AddEntry(h,"#Delta#phi > 7/8 #pi  ");
00057 
00058   mL->Draw(); 
00059   ln->Draw();
00060 
00061   // right plot
00062 
00063   c1->cd(2);
00064  
00065   mL=new TLegend(0.15,.45,.7,.70);
00066 
00067   h=(TH1F*)fd->Get("CnFP-A"); 
00068   h->Draw();  mL->AddEntry(h,"#Delta#phi <1/8 #pi  (data x 1/10)");
00069   
00070   h=(TH1F*)fd->Get("CnFP-B"); h->SetLineColor(kBlue);
00071   h->SetLineStyle(2);
00072   h->Draw("same");  mL->AddEntry(h," #Delta#phi #in [ 1/8 #pi, 3/8 #pi ] ");
00073 
00074   h=(TH1F*)fd->Get("CnFP-D"); h->SetLineColor(kRed);
00075   h->Draw("same");  mL->AddEntry(h," #Delta#phi #in [ 5/8 #pi, 7/8 #pi ] ");
00076 
00077   h=(TH1F*)fd->Get("CnFP-L"); h->SetLineColor(kGreen);
00078   h->Draw("same");  mL->AddEntry(h,"lost LCP");
00079 
00080   h=(TH1F*)fd->Get("CnFP-W"); h->SetLineColor(kRed);
00081   h->SetLineStyle(2);
00082   h->Draw("same");  mL->AddEntry(h,"won LCP");
00083 
00084   mL->Draw(); 
00085   ln->Draw();
00086 
00087   
00088 
00089 #if 0
00090   c1->Divide(2,4);
00091   c1->cd(1); (fd->Get("CnFP-A"))->Draw();
00092   c1->cd(2); (fd->Get("CnFP-B"))->Draw();
00093   c1->cd(3); (fd->Get("CnFP-C"))->Draw();
00094   c1->cd(4); (fd->Get("CnFP-D"))->Draw();
00095   c1->cd(5); (fd->Get("CnFP-E"))->Draw();
00096   c1->cd(7); (fd->Get("CnFP-L"))->Draw();
00097    c1->cd(8);  (fd->Get("CnFP-W"))->Draw();
00098 #endif
00099   
00100 }
00101