change axis title style to italics

[u/mrichter/AliRoot.git] / EVGEN / AliGenHBTosl.cxx

diff --git a/EVGEN/AliGenHBTosl.cxx b/EVGEN/AliGenHBTosl.cxx
index cf0610377b22ead3cfe335842151a30d5c8632c8..7d763d215903175dbb3a9c74580e58f3778ce40b 100644 (file)
--- a/EVGEN/AliGenHBTosl.cxx
+++ b/EVGEN/AliGenHBTosl.cxx
@@ -87,6 +87,10 @@
 
 
 /***********************************************************/
+using std::cout;
+using std::endl;
+using std::ofstream;
+using std::ios;
 ClassImp(AliGenHBTosl)
 
 AliGenHBTosl::AliGenHBTosl():
@@ -373,7 +377,7 @@ void AliGenHBTosl::Generate()
   }
  
  fGenerator->Generate();
- Int_t j = 0, ntr = 0;
+ Int_t ntr = 0;
  if ( genstack->GetNtrack() < fNpart/2)
   {
     Warning("Generate","************************************************************");
@@ -426,7 +430,7 @@ void AliGenHBTosl::Generate()
  
  char msg[1000];
  logfile<<endl;
- sprintf(msg,"\n");
+ snprintf(msg,1000, "\n");
  Int_t middlebin = fQNBins/2;
  
  for (Int_t k = middlebin-5; k < middlebin+5; k++)
@@ -434,21 +438,21 @@ void AliGenHBTosl::Generate()
      Double_t tx = work->GetXaxis()->GetBinCenter(30);
      Double_t ty = work->GetYaxis()->GetBinCenter(30);
      Double_t tz = work->GetZaxis()->GetBinCenter(k);
-     sprintf(msg,"% 6.5f ",GetQOutQSideQLongCorrTheorValue(tx,ty,tz));
+     snprintf(msg,1000, "% 6.5f ",GetQOutQSideQLongCorrTheorValue(tx,ty,tz));
      logfile<<msg;
    }
  logfile<<endl;
   
  for (Int_t k = middlebin-5; k < middlebin+5; k++)
   {
-    sprintf(msg,"% 6.5f ",work->GetBinContent(30,30,k));
+    snprintf(msg,1000, "% 6.5f ",work->GetBinContent(30,30,k));
     logfile<<msg;
   }
  logfile<<endl;
 
  for (Int_t k = middlebin-5; k < middlebin+5; k++)
   {
-    sprintf(msg,"% 6.5f ",chiarray[30][30][k]);
+    snprintf(msg, 1000, "% 6.5f ",chiarray[30][30][k]);
     logfile<<msg;
   }
  logfile<<endl;
@@ -520,7 +524,7 @@ void AliGenHBTosl::Generate()
      Double_t qoutc  = work->GetXaxis()->GetBinCenter(xmax);
      
 
-     sprintf(msg,"Generate Fill bin chi2(%d,%d,%d)=%f",xmax,ymax,zmax,chiarray[xmax][ymax][zmax]);
+     snprintf(msg,1000, "Generate Fill bin chi2(%d,%d,%d)=%f",xmax,ymax,zmax,chiarray[xmax][ymax][zmax]);
      logfile<<msg;
      
      qout  = gRandom->Uniform(qoutc -binwdh, qoutc +binwdh);
@@ -528,9 +532,11 @@ void AliGenHBTosl::Generate()
      qlong = gRandom->Uniform(qlongc-binwdh, qlongc+binwdh);
 
      TParticle* first = 0;
-     while (j < genstack->GetNtrack())
+     Int_t jj = 0;
+     
+     while (jj < genstack->GetNtrack())
       {
-        TParticle* tmpp = genstack->Particle(j++);
+        TParticle* tmpp = genstack->Particle(jj++);
         if (tmpp->GetPdgCode() == fPID)
          {
            if (CheckParticle(tmpp,0x0,stack) == kFALSE)
@@ -896,21 +902,21 @@ void AliGenHBTosl::StartSignal()
     {
       c1->cd();
       char buff[50];
-      sprintf(buff,"QTWorkPass2.%3d.root",counter);
+      snprintf(buff,50, "QTWorkPass2.%3d.root",counter);
       TFile* file = TFile::Open(buff,"update");
       work->Write();
       work->SetDirectory(0x0);
       file->Close();
       delete file;
 
-      sprintf(buff,"QTBackgroundPass2.%3d.root",counter);
+      snprintf(buff,50, "QTBackgroundPass2.%3d.root",counter);
       file = TFile::Open(buff,"update");
       fQBackground->Write();
       fQBackground->SetDirectory(0x0);
       file->Close();
       delete file;
 
-      sprintf(buff,"QTSignalPass2.%3d.root",counter);
+      snprintf(buff,50, "QTSignalPass2.%3d.root",counter);
       file = TFile::Open(buff,"update");
       fQSignal->Write();
       fQSignal->SetDirectory(0x0);
@@ -958,21 +964,21 @@ void AliGenHBTosl::StartSignal()
       Double_t tx = work->GetXaxis()->GetBinCenter(30);
       Double_t ty = work->GetYaxis()->GetBinCenter(30);
       Double_t tz = work->GetZaxis()->GetBinCenter(k);
-      sprintf(msg,"% 6.5f ",GetQOutQSideQLongCorrTheorValue(tx,ty,tz));
+      snprintf(msg,1000, "% 6.5f ",GetQOutQSideQLongCorrTheorValue(tx,ty,tz));
       logfile<<msg;
     }
    logfile<<endl;
 
    for (Int_t k = 25; k < 36; k++)
     {
-      sprintf(msg,"%6.5f ",work->GetBinContent(30,30,k));
+      snprintf(msg, 1000, "%6.5f ",work->GetBinContent(30,30,k));
       logfile<<msg;
     }
    logfile<<endl;
 
    for (Int_t k = 25; k < 36; k++)
     {
-      sprintf(msg,"% 6.5f ",chiarray[30][30][k]);
+      snprintf(msg,1000, "% 6.5f ",chiarray[30][30][k]);
       logfile<<msg;
     }
    logfile<<endl;
@@ -1290,7 +1296,7 @@ void AliGenHBTosl::StartSignalPass1()
     
     for (Int_t k = 1; k<=fQNBins; k++)
       {
-       for (Int_t j = 1; j<=fQNBins; j++)
+       for (j = 1; j<=fQNBins; j++)
          {  
            for (Int_t i = 1; i<=fQNBins; i++)
              {
@@ -1563,7 +1569,7 @@ void  AliGenHBTosl::GetQOutQSideQLong(TParticle* f, TParticle* s,Double_t& out,
   }
 
  Double_t beta = pzsum/esum;
- Double_t gamma = 1.0/TMath::Sqrt(1.0 - beta*beta);
+ Double_t gamma = 1.0/TMath::Sqrt((1.-beta)*(1.+beta));
 
  lon = gamma * ( pzdiff - beta*ediff );