draw_txt.cxx File Reference

#include <fstream>
#include <iostream>
#include <cmath>

Go to the source code of this file.

Functions
int	draw_txt ()

Function Documentation

draw_txt()

int draw_txt

(

)

Definition at line 9 of file draw_txt.cxx.

{
    const int nfile = 5;
    double momentum[nfile]  = {150, 200, 300, 500, 1000};
    double sigma[nfile][25] = {0};
    double error[nfile][25] = {0};
 
    int l(0),i,file,last(-1);
    double s(0),e(0);
    string strcom;
    ifstream fin("log",ios::in);
    //std::getline(fin, strcom);
    //fin >> l  >>  s >>  e;
    //cout<<l<<"  "<<s<<"  "<<e<<endl;
    while(std::getline(fin, strcom)){
        fin >> l  >>  s >>  e;
        if(l == last)continue;
        if(l==0)i=0;
        sigma[file][i] = s;
        error[file][i] = e;
        //cout<<l<<"  "<<s<<"  "<<e<<endl;
        //cout<<i<<"  "<<sigma[file][i]<<"  "<<error[file][i]<<endl;
        i++;
        if(l==24)file++;
        last = l;
        //cout<<l<<"  "<<s<<"  "<<e<<endl;
    }
///*
    stringstream ssname;
    string sname;
    const char * name;
    TAxis *x_axis,*y_axis;
    TGraphErrors *gr_sigma_pt;
    TCanvas *c_pt = new TCanvas("sigma_pt","sigma_pt",2500,2500 );
    c_pt->Divide(5,5);
    for(int i=0;i<25;i++){
        TGraphErrors *gr_sigma_pt = new TGraphErrors();
        for(int file=0;file<nfile;file++){
            //cout<<i<<"  "<<sigma[file][i]<<"  "<<error[file][i]<<endl;
            gr_sigma_pt->SetPoint(file,momentum[file],sigma[file][i]);
            gr_sigma_pt->SetPointError(file,0,error[file][i]);
        }
        c_pt->cd(i+1);
        gr_sigma_pt->Draw("APLsame");
        gr_sigma_pt->SetMarkerStyle(20);
        gr_sigma_pt->SetMarkerSize(1.0);
        gr_sigma_pt->SetMarkerColor(3);
        gr_sigma_pt->SetLineColor(3);
        x_axis = gr_sigma_pt->GetXaxis();
        y_axis = gr_sigma_pt->GetYaxis();
        x_axis->SetTitle("transverse momentum (MeV)");
        x_axis->CenterTitle();
        x_axis->SetTitleSize(0.05);
        x_axis->SetLabelSize(0.04);
        x_axis->SetTitleOffset(1.05);
        y_axis->SetTitle("#sigma (cm)");
        y_axis->CenterTitle();
        y_axis->SetTitleSize(0.05);
        y_axis->SetLabelSize(0.04);
        //y_axis->SetTitleOffset(1.05);
        ssname.str("");
        if(i<3){ssname <<"layer"<<i+1;}
        else if(i<15){ssname <<"layer"<<i+6;}
        else if(i<22){ssname <<"layer"<<i+22;}
        else if(i==22){ssname <<"1~3layer";}
        else if(i==23){ssname <<"9~20layer";}
        else if(i==24){ssname <<"37~43layer";}
        //if(l<=3){ssname <<"layer "<<l;cut=1;}
        //else if(l<=15){ssname <<"layer "<<l+5;cut=1;}
        //else if(l<=22){ssname <<"layer "<<l+21;cut=1;}
        //else if(l==23){ssname <<"1~3 layer";cut=1;}
        //else if(l==24){ssname <<"9~20 layer";cut=1;}
        //else if(l==25){ssname <<"37~43 layer";cut=1;}
        sname = ssname.str();
        name = sname.c_str();
        gr_sigma_pt->SetTitle(name);
    }
    c_pt->SaveAs("sigma_pt.png");
    c_pt->SaveAs("sigma_pt.pdf");
    c_pt->SaveAs("sigma_pt.eps");
    TCanvas *c_pt_layer = new TCanvas("sigma_pt_layer","sigma_pt_layer",1000,1000 );
    TLegend* leg = new TLegend(0.30,0.60,0.60,0.88);
    for(int file=0;file<nfile;file++){
        TGraphErrors *gr_sigma_pt_layer = new TGraphErrors();
        for(int i=0;i<22;i++){
            int layer;
            if(i<3){layer = i+1;}
            else if(i<15){layer = i+6;}
            else if(i<22){layer = i+22;}
            //cout<<i<<"  "<<layer<<"  "<<sigma[file][i]<<endl;
            gr_sigma_pt_layer->SetPoint(i,layer,sigma[file][i]);
            gr_sigma_pt_layer->SetPointError(i,0,error[file][i]);
        }
        if(file==0)gr_sigma_pt_layer->Draw("APLsame");
        else gr_sigma_pt_layer->Draw("PLsame");
        gr_sigma_pt_layer->SetMarkerStyle(20);
        gr_sigma_pt_layer->SetMarkerSize(1.0);
        gr_sigma_pt_layer->SetMarkerColor(file+1);
        gr_sigma_pt_layer->SetLineColor(file+1);
        x_axis = gr_sigma_pt_layer->GetXaxis();
        y_axis = gr_sigma_pt_layer->GetYaxis();
        x_axis->SetTitle("layer");
        x_axis->CenterTitle();
        x_axis->SetTitleSize(0.05);
        x_axis->SetLabelSize(0.04);
        x_axis->SetTitleOffset(1.05);
        y_axis->SetTitle("#sigma (cm)");
        y_axis->CenterTitle();
        y_axis->SetTitleSize(0.05);
        y_axis->SetLabelSize(0.04);
        //y_axis->SetTitleOffset(1.05);
        //ssname.str("");
        //if(i<3){ssname <<"layer"<<i+1;}
        //else if(i<15){ssname <<"layer"<<i+6;}
        //else if(i<22){ssname <<"layer"<<i+22;}
        //else if(i==22){ssname <<"1~3layer";}
        //else if(i==23){ssname <<"9~20layer";}
        //else if(i==24){ssname <<"37~43layer";}
        //if(l<=3){ssname <<"layer "<<l;cut=1;}
        //else if(l<=15){ssname <<"layer "<<l+5;cut=1;}
        //else if(l<=22){ssname <<"layer "<<l+21;cut=1;}
        //else if(l==23){ssname <<"1~3 layer";cut=1;}
        //else if(l==24){ssname <<"9~20 layer";cut=1;}
        //else if(l==25){ssname <<"37~43 layer";cut=1;}
        //sname = ssname.str();
        //name = sname.c_str();
        ssname.str("");
        ssname <<momentum[file]<<"MeV";
        sname = ssname.str();
        name = sname.c_str();
        gr_sigma_pt_layer->SetTitle("sigma_layer");
        leg->AddEntry(gr_sigma_pt_layer,name,"pl");
 
    }
    leg->Draw();
    leg->SetFillColor(kWhite);
    leg->SetLineColor(kWhite);
    c_pt_layer->SaveAs("sigma_pt_layer.png");
    c_pt_layer->SaveAs("sigma_pt_layer.pdf");
    c_pt_layer->SaveAs("sigma_pt_layer.eps");
//*/
}