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" );
}