DataBase/tau_mode.c

Go to the documentation of this file.

{
    gROOT->Reset();
 
    TFile *f1 = new TFile("YYYY/m_root_dir/LumTau_XXXX.root");
    TTree *data =event;
 
    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);
    }
 
    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("YYYY/m_txt_dir/LumTau_XXXX.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();
 
    gROOT->ProcessLine(".q");
}