tau_mode.c

Go to the documentation of this file.

{
    gROOT->Reset();
 
    TFile *f1 = new TFile("/bes3fs/groups/psip/psipgroup/temp/user/liuh/BbLum/psipp/BbEmc_00XXXX.root");
    TTree *data =event;
 
    TCanvas *myCanvas = new TCanvas();   
    myCanvas->Divide(1,1);
    TPad *c1_1 = new TPad("c1_1", "c1_1",0.01,0.01,0.99,0.99);
    c1_1->Draw();
    c1_1->cd();
    c1_1->Range(0.425458,-114.842,0.674993,802.951);
    c1_1->SetBorderSize(2);
    c1_1->SetBottomMargin(0.125129);
    c1_1->SetFrameFillColor(0);
 
    Double_t time = 0.0,dltphi = 0.0,costht1 = 0.0,costht2 = 0.0,e1 = 0.0,e2 = 0.0,etot = 0.0,phi1 = 0.0,phi2 = 0.0;
    Int_t nentries = 0;
 
    data->SetBranchAddress("time",&time);
    data->SetBranchAddress("dltphi",&dltphi);
    data->SetBranchAddress("costht1",&costht1);
    data->SetBranchAddress("costht2",&costht2);
    data->SetBranchAddress("e1",&e1);
    data->SetBranchAddress("e2",&e2);
    data->SetBranchAddress("etot",&etot);
    data->SetBranchAddress("phi1",&phi1);
    data->SetBranchAddress("phi2",&phi2);
    TH1D  *hdltphi[10];
 
    for(Int_t i=0;i<10;i++)
    {
        hdltphi[i] = new TH1D(" ","dltphi distribution",50,-50.,50.);
    }
 
    TH1 *htime = new TH1D("htime","time", 80, -10., 10.);
    Double_t timemin = 0.0,timemax = 0.0;
    nentries = (Int_t)data->GetEntries();
 
    data->GetEntry(0);
    timemin = time;
    data->GetEntry((nentries - 1));
    timemax = time;
 
    Double_t difft = timemax-timemin;
    Double_t xtime[11],lum[10];
    for(Int_t i=0;i<11;i++)
    {
        xtime[i]=timemin+(difft/10.)*i;
    }
    for(Int_t i=0;i<10;i++)
    {
        for(Int_t j=0;j<nentries;j++)
        {
            data->GetEntry(j);
            bool lumcut = e1>1.2&&e2>0.8&&fabs(costht1)<0.8&&
                fabs(costht2)<0.8&&fabs(dltphi+1.055)>2.&&
                fabs(dltphi+1.055)<25.&&etot<4.0;
            if(lumcut&&time>xtime[i]&&time<xtime[i+1])
            {
                hdltphi[i]->Fill(dltphi);
            }
        }
    }
 
    Int_t nevt[10];
    Double_t x[10];
    for(Int_t i=0;i<10;i++)
    {
        nevt[i] = hdltphi[i]->GetEntries();
        x[i]=xtime[i+1]-xtime[0];
        lum[i]=nevt[i]/(0.2*2514);
    }
 
    TCanvas *c1 = new TCanvas("c1"," ",700,500);
    c1->SetFillColor(kWhite);
    c1->SetGrid();
 
    const Int_t n = 10;
    TF1 *g1    = new TF1("g1","[0]*exp(-x/[1])",0,10000.);
    g1->SetParameters(lum[0],1e+4.);
    g1->SetLineColor(2);
    TGraph* gr = new TGraph(n,x,lum);
    gr->SetLineColor(2);
    gr->SetLineWidth(2);
    gr->SetMarkerColor(4);
    gr->SetMarkerStyle(21);
    gr->SetTitle("BbLum_000XXXX");
    gr->GetXaxis()->SetTitle("Time");
    gr->GetYaxis()->SetTitle("Luminosity");
    gr->Fit("g1","R");
    gr->Draw("ap");
 
    std::ofstream m_outputFile;
    m_outputFile.open("Lum_Tau.txt", ios_base::app);
    m_outputFile<<XXXX<<"       "<<difft<<"     "<<g1->GetParameter(0)<<"       "<<(g1->GetParameter(0) * exp(- 1.0 * difft / (g1->GetParameter(1))))<<"        "<<g1->GetParameter(1)<<endl;
    m_outputFile.close();
 
    TString result  = Form("tau = %6.2f",g1->GetParameter(1));
    TText* text01 = new TText(0.7,0.8,result);
    text01->SetNDC();
    text01->Draw("same");
    c1->Update();
    c1->GetFrame()->SetFillColor(kWhite);
    c1->GetFrame()->SetBorderSize(1);
    c1->Modified();
    c1->Print("lum_0000XXXX.ps");
}