From 2d065274fced2a4b03240b173616f8b34a92105d Mon Sep 17 00:00:00 2001
From: alla <alla@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 25 Oct 2007 12:25:53 +0000
Subject: [PATCH] remove warnings

---
 T0/AliT0CalibLaserData.cxx | 401 ++++++++++++++++++-------------------
 1 file changed, 200 insertions(+), 201 deletions(-)

diff --git a/T0/AliT0CalibLaserData.cxx b/T0/AliT0CalibLaserData.cxx
index d2c386ae19d..f08997f15b9 100644
--- a/T0/AliT0CalibLaserData.cxx
+++ b/T0/AliT0CalibLaserData.cxx
@@ -16,7 +16,6 @@
 #include "AliT0LookUpValue.h"
 #include "AliT0Parameters.h"
 #include "AliT0RawReader.h"
-#include "AliLog.h"		  
 
 ClassImp(AliT0CalibLaserData)
 
@@ -43,10 +42,10 @@ AliT0CalibLaserData &AliT0CalibLaserData::operator =(const AliT0CalibLaserData&
   return *this;
 }
 //________________________________________________________________
-//AliT0CalibLaserData::~AliT0CalibLaserData()
-//{
+AliT0CalibLaserData::~AliT0CalibLaserData()
+{
   //
-//}
+}
 */
 //________________________________________________________________
 
@@ -96,11 +95,11 @@ void AliT0CalibLaserData::ReadHistSize(Int_t rNumber)
     fEntries[8] = new TGNumberEntry(fMain, 700);
     fMain->AddFrame(fEntries[8]);
 // CFD axis Y 
-   fEntries[9] = new TGNumberEntry(fMain, 0);
+   fEntries[9] = new TGNumberEntry(fMain, 12000.5);
     fMain->AddFrame(fEntries[9]);
-    fEntries[10] = new TGNumberEntry(fMain, 10000);
+    fEntries[10] = new TGNumberEntry(fMain, 18000.5);
     fMain->AddFrame(fEntries[10]);
-    fEntries[11] = new TGNumberEntry(fMain, 10000);
+    fEntries[11] = new TGNumberEntry(fMain, 2000);
     fMain->AddFrame(fEntries[11]);
 //
     fMain->AddFrame( new TGLabel(fMain, "CFD vs LED-CFD " ) );
@@ -112,11 +111,11 @@ void AliT0CalibLaserData::ReadHistSize(Int_t rNumber)
     fEntries[14] = new TGNumberEntry(fMain, 500);
     fMain->AddFrame(fEntries[14]);
 // CFD axis Y
-    fEntries[15] = new TGNumberEntry(fMain, 3900);
+    fEntries[15] = new TGNumberEntry(fMain, 14600);
     fMain->AddFrame(fEntries[15]);
-    fEntries[16] = new TGNumberEntry(fMain, 4100);
+    fEntries[16] = new TGNumberEntry(fMain, 14700);
     fMain->AddFrame(fEntries[16]);
-    fEntries[17] = new TGNumberEntry(fMain, 200);
+    fEntries[17] = new TGNumberEntry(fMain, 100);
     fMain->AddFrame(fEntries[17]);
  
    fMain->AddFrame( new TGLabel(fMain, " Number of run") );
@@ -147,216 +146,216 @@ void AliT0CalibLaserData::DoOk()
     for( int i=0; i<18; i++ ) cout<<fHistLimits[i]<<" ";
     cout<<endl;
     ReadData();
-
 }
 
 void AliT0CalibLaserData::ReadData()
 {
 
-  TH1F *hChannel[105];  TH1F *hQTC[12];  
-  TH2F *hCFD_QTC[12]; TH2F *hCFD_LED[12]; TH1F *h1CFD_LED[12];
-  TH1F *hmpd[12];
-  Int_t allData[110][20];
-  Int_t numberOfHits[105];
-  
-  Char_t  buf1[20], buf2[20], buf3[20], buf4[20],  buf7[20];
-  
-  TTree* digitsTree = new TTree("testData","Tree of test data Digits");
-  TBranch *b[106];
-  
-  Int_t channels[106];
-  
-  TString names[106], type;
-  AliT0LookUpKey* lookkey= new AliT0LookUpKey();
-  AliT0LookUpValue*  lookvalue= new AliT0LookUpValue();
-  AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
-  AliCDBManager::Instance()->SetRun(0);
-  AliT0Parameters *fParam = AliT0Parameters::Instance();
-  fParam->Init();
-  TMap *lookup = fParam->GetMapLookup();
-  TMapIter *iter = new TMapIter(lookup);
-  
-  for( Int_t iline=0; iline<106; iline++)
-    {
-      lookvalue = ( AliT0LookUpValue*) iter->Next();
-      lookkey = (AliT0LookUpKey*) lookup->GetValue((TObject*)lookvalue);
-      if(lookkey){
-	Int_t key=lookkey->GetKey();
-	names[key]=lookkey->GetChannelName();
-	//	cout<<lookkey->GetKey()<<" "<<lookkey->GetChannelName()<<" trm "<<lookvalue->GetTRM()<<" tdc "<<lookvalue->GetTDC()<<" chain  "<<lookvalue->GetChain()<<" channel "<<lookvalue->GetChannel()<<endl;
-	hChannel[key] = new TH1F(names[key].Data(),names[key].Data(),30000,0,30000);
-	//	hitsname="xHits" + names[key];
-	//	hNumHits[key] = new TH1F(hitsname.Data(),hitsname.Data(),50,-0.25,24.25);
-	type =names[key] + "/I";
-	b[key]=digitsTree->Branch(names[key].Data(),&channels[key], type);
-      }
-      else
-	{cout<<iline<<" no such value "<<endl;}
-      
-    } 
-  for(Int_t ic=0; ic<12; ic++) {
-    {
-      sprintf(buf1,"QTC%i",ic+1);
-      sprintf(buf2,"CFDvsQTC%i",ic+1);
-      sprintf(buf3,"CFDvsLED%i",ic+1);
-      sprintf(buf4,"LED_CFD%i",ic+1);
-      sprintf(buf7,"mpd%i",ic+1);
-      
-      hQTC[ic] = new TH1F(buf1,"QTC",(Int_t)fHistLimits[2],fHistLimits[0],fHistLimits[1]);
-       h1CFD_LED[ic] = new TH1F(buf4,"LED - CFD",(Int_t)fHistLimits[5],fHistLimits[3],fHistLimits[4]);
-       hmpd[ic] = new TH1F(buf7,"mpd",20000,-10000.0,10000.0);
-      hCFD_QTC[ic] = new TH2F(buf2,"CFD vs	QTC",
-			      (Int_t)fHistLimits[8],fHistLimits[6],fHistLimits[7],
-			      (Int_t)fHistLimits[11],fHistLimits[9],fHistLimits[10]);
-      hCFD_LED[ic] = new TH2F(buf3,"CFD vs LED-CFD",
-			      (Int_t)fHistLimits[14],fHistLimits[12],fHistLimits[13],
-			      (Int_t)fHistLimits[17],fHistLimits[15],fHistLimits[16]);
-      
-      
-    }
+      TH1F *hChannel[105];  TH1F *hQTC[12];  
+      TH2F *hCFD_QTC[12]; TH2F *hCFD_LED[12]; TH1F *h1CFD_LED[12];
+      TH1F *hmpd[12];
+      Int_t allData[110][5];
+      Int_t numberOfHits[105];
 
+      Char_t  buf1[20], buf2[20], buf3[20], buf4[20], buf7[20];
 
-  }
-      //   cout<<" hist created "<<endl;
-  
-  TH1F*hEffCFD= new TH1F("hEffCFD","Effeciency",50,-0.25,24.25);
-  TH1F*hEffLED= new TH1F("hEffLED","Effeciency",50,-0.25,24.25);
-  TH1F*hEffQT0= new TH1F("hEffQT0","Effeciency",50,-0.25,24.25);
-  TH1F*hEffQT1= new TH1F("hEffQT1","Effeciency",50,-0.25,24.25);
-  
-  
-  Char_t filename[100];
-  sprintf(filename,"/home/t0/alice/testSep07/raw/t0%i.001.raw",fRunNumber);
-  AliRawReader *reader = new AliRawReaderDate(filename);
-  if(!reader) AliFatal(Form("Can not opne file ",filename));
-  // AliRawReader *reader = new AliRawReaderFile();
-  reader->LoadEquipmentIdsMap("T0map.txt");
-  //    reader->RequireHeader(kFALSE);
-  reader->RequireHeader(kTRUE);
-  AliT0RawReader *start = new AliT0RawReader(reader);
-  //  start->SetNumberOfTRM(1);
-  for (Int_t i0=0; i0<105; i0++)
-    {
-      for (Int_t j0=0; j0<5; j0++) allData[i0][j0]=0; 	
-      numberOfHits[i0]=0;
-    }
-  Int_t event=0;
-
-  while (reader->NextEvent()) {
-    start->Next();
-    for (Int_t i=0; i<105; i++) {
-      for (Int_t iHit=0; iHit<20; iHit++) 
+      TTree* digitsTree = new TTree("testData","Tree of test data Digits");
+      TBranch *b[106];
+
+       Int_t channels[106];
+
+      TString names[106], type;
+      AliT0LookUpKey* lookkey= new AliT0LookUpKey();
+      AliT0LookUpValue*  lookvalue= new AliT0LookUpValue();
+      AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
+      AliCDBManager::Instance()->SetRun(0);
+      AliT0Parameters *fParam = AliT0Parameters::Instance();
+      fParam->Init();
+      TMap *lookup = fParam->GetMapLookup();
+      TMapIter *iter = new TMapIter(lookup);
+
+      for( Int_t iline=0; iline<106; iline++)
 	{
-	  allData[i][iHit]= start->GetData(i,iHit);
+	  lookvalue = ( AliT0LookUpValue*) iter->Next();
+	   lookkey = (AliT0LookUpKey*) lookup->GetValue((TObject*)lookvalue);
+	  if(lookkey){
+	    Int_t key=lookkey->GetKey();
+	    names[key]=lookkey->GetChannelName();
+	    //	cout<<lookkey->GetKey()<<" "<<lookkey->GetChannelName()<<" trm "<<lookvalue->GetTRM()<<" tdc "<<lookvalue->GetTDC()<<" chain  "<<lookvalue->GetChain()<<" channel "<<lookvalue->GetChannel()<<endl;
+	    hChannel[key] = new TH1F(names[key].Data(),names[key].Data(),30000,0,30000);
+	    //	hitsname="xHits" + names[key];
+	    //	hNumHits[key] = new TH1F(hitsname.Data(),hitsname.Data(),50,-0.25,24.25);
+	    type =names[key] + "/I";
+	    b[key]=digitsTree->Branch(names[key].Data(),&channels[key], type);
+	   }
+	  else
+	    {cout<<iline<<" no such value "<<endl;}
+
+	} 
+      for(Int_t ic=0; ic<12; ic++) {
+	{
+	  sprintf(buf1,"QTC%i",ic+1);
+	  sprintf(buf2,"CFDvsQTC%i",ic+1);
+	  sprintf(buf3,"CFDvsLED%i",ic+1);
+	  sprintf(buf4,"LED_CFD%i",ic+1);
+	  sprintf(buf7,"mpd%i",ic+1);
+
+	  hQTC[ic] = new TH1F(buf1,"QTC",(Int_t)fHistLimits[2],fHistLimits[0],fHistLimits[1]);
+	  cout<<" hQTC "<<fHistLimits[2]<<" "<<fHistLimits[0]<<" "<<fHistLimits[1]<<endl;
+	  h1CFD_LED[ic] = new TH1F(buf4,"LED - CFD",(Int_t)fHistLimits[5],fHistLimits[3],fHistLimits[4]);
+	  cout<<" hLED-CFD  "<<fHistLimits[5]<<" "<<fHistLimits[3]<<" "<<fHistLimits[4]<<endl;
+	  hmpd[ic] = new TH1F(buf7,"mpd",20000,-10000.0,10000.0);
+	  hCFD_QTC[ic] = new TH2F(buf2,"CFD vs	QTC",
+	  	                 (Int_t)fHistLimits[8],fHistLimits[6],fHistLimits[7],
+		                 (Int_t)fHistLimits[11],fHistLimits[9],fHistLimits[10]);
+	  hCFD_LED[ic] = new TH2F(buf3,"CFD vs LED-CFD",
+	                         (Int_t)fHistLimits[14],fHistLimits[12],fHistLimits[13],
+	                         (Int_t)fHistLimits[17],fHistLimits[15],fHistLimits[16]);
+
+
 	}
-    }
-    
-    if(event%1000 == 0) printf("Event:%d\n",event);
-    
-    //        if(event > 100000) break;
-    
 
-    for (Int_t it = 0; it<24; it=it+2)
-      {
-	for (Int_t iHit=0; iHit<20; iHit++)
-	  {
-	    if(allData[it+25][iHit] != 0 && allData[it+26][iHit] !=0)
+
+      }
+      //   cout<<" hist created "<<endl;
+
+      TH1F*hEffCFD= new TH1F("hEffCFD","Effeciency",50,-0.25,24.25);
+      TH1F*hEffLED= new TH1F("hEffLED","Effeciency",50,-0.25,24.25);
+      TH1F*hEffQT0= new TH1F("hEffQT0","Effeciency",50,-0.25,24.25);
+      TH1F*hEffQT1= new TH1F("hEffQT1","Effeciency",50,-0.25,24.25);
+
+      
+      Char_t filename[13];
+      sprintf(filename,"t0%i.raw",fRunNumber);
+       AliRawReader *reader = new AliRawReaderDate(filename);
+     // AliRawReader *reader = new AliRawReaderFile();
+       reader->LoadEquipmentIdsMap("T0map.txt");
+      //    reader->RequireHeader(kFALSE);
+       reader->RequireHeader(kTRUE);
+       AliT0RawReader *start = new AliT0RawReader(reader);
+       //  start->SetNumberOfTRM(1);
+       for (Int_t i0=0; i0<105; i0++)
+	 {
+	   for (Int_t j0=0; j0<5; j0++) allData[i0][j0]=0; 	
+	   numberOfHits[i0]=0;
+	 }
+       Int_t event=0;
+
+       while (reader->NextEvent()) {
+	 start->Next();
+	 for (Int_t i=0; i<105; i++) {
+	   for (Int_t iHit=0; iHit<5; iHit++) 
+	    {
+	      allData[i][iHit]= start->GetData(i,iHit);
+	    }
+	}
+      
+       if(event%1000 == 0) printf("Event:%d\n",event);
+
+          if(event > 200000) break;
+	//  cout<<"!!!!  Event Number "<< event-2<<endl;
+
+	    for (Int_t it = 0; it<24; it=it+2)
 	      {
-		Int_t cc=it/2;
-		hQTC[cc]->Fill(allData[it+25][iHit]-allData[it+26][iHit]);
-		hmpd[cc]->Fill(allData[it+26][iHit]-allData[it+25][iHit]);
-		if(allData[cc+1][iHit] != 0 ) hCFD_QTC[cc]->Fill(allData[it+25][iHit]-allData[it+26][iHit],allData[cc+1][iHit]-allData[0][0]);
-		if(allData[cc+1][iHit] != 0 && allData[cc+13][iHit]!=0 ) 
+		for (Int_t iHit=0; iHit<5; iHit++)
+		  {
+		if(allData[it+25][iHit] != 0 && allData[it+26][iHit] !=0)
 		  {
-		    hCFD_LED[cc]->Fill(allData[cc+13][iHit]-allData[cc+1][iHit],allData[cc+1][iHit]-allData[0][0]);  
-		    h1CFD_LED[cc]->Fill(allData[cc+13][iHit]-allData[cc+1][iHit]);
+		    Int_t cc=it/2;
+		    // if( allData[56][0]-allData[0][0] > 0) 
+		    hQTC[cc]->Fill(allData[it+25][iHit]-allData[it+26][iHit]);
+//		    if( allData[55][0]-allData[0][0] > 0) hQTCc[cc]->Fill(allData[it+26][iHit]-allData[it+25][iHit]);
+		    hmpd[cc]->Fill(allData[it+26][iHit]-allData[it+25][iHit]);
+		    if(allData[cc+1][iHit] != 0 ) hCFD_QTC[cc]->Fill(allData[it+25][iHit]-allData[it+26][iHit],allData[cc+1][iHit]-allData[0][0]);
+		    if(allData[cc+1][iHit] != 0 && allData[cc+13][iHit]!=0 ) 
+		      {
+			hCFD_LED[cc]->Fill(allData[cc+13][iHit]-allData[cc+1][iHit],allData[cc+1][iHit]-allData[0][0]);  
+			h1CFD_LED[cc]->Fill(allData[cc+13][iHit]-allData[cc+1][iHit]);
+		      }
+		    //  cout<<allData[cc+1][iHit]-allData[0][0]<<" "<<cc<<endl;
 		  }
 	      }
 	  }
-      }
-    
-    for (Int_t iHit=0; iHit<20; iHit++) 
-      {
-	
-	for(Int_t ik=1; ik<105; ik++)
-	  { 
-	    channels[ik] = -100;
-	    if((allData[ik][iHit] - allData[0][0]) > 0 ) 
-	      {
-		numberOfHits[ik]++;
-		hChannel[ik] -> Fill(allData[ik][iHit] - allData[0][0]);
-		channels[ik] = allData[ik][iHit] - allData[0][0];
+
+	for (Int_t iHit=0; iHit<5; iHit++) 
+	  {
+
+	    for(Int_t ik=1; ik<105; ik++)
+	      { 
+		channels[ik] = -100;
+		if((allData[ik][iHit] - allData[0][0]) > 0 ) 
+		  {
+		    numberOfHits[ik]++;
+		    hChannel[ik] -> Fill(allData[ik][iHit] - allData[0][0]);
+		    channels[ik] = allData[ik][iHit] - allData[0][0];
+		  }
 	      }
+	      digitsTree->Fill();   
 	  }
-	//	digitsTree->Fill();   
-      }
-    
-    event++;
-    
-  } //event
-  
-  
-  if (event>1)
-    {
-      cout<<"efficiency for "<<event<<" events"<<endl;
-      for (Int_t i0=1; i0<13;  i0++)
+
+	event++;
+
+      } //event
+
+
+     if (event>1)
 	{
-	  
-	  cout<<names[i0].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<" ";
-	  cout<<names[i0+13].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<" ";
-	  cout<<names[i0+57].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<" ";
-	  cout<<names[i0+69].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<endl;
-	  
-	  hEffCFD->Fill(i0,Float_t(numberOfHits[i0])  / Float_t(event));
-	  hEffLED->Fill(i0,Float_t(numberOfHits[i0+13]) / Float_t(event));
-	  hEffCFD->Fill(i0+12,Float_t(numberOfHits[i0+57]) /Float_t(event));
-	  hEffLED->Fill(i0+12,Float_t(numberOfHits[i0+69]) /Float_t(event));
-	}
-      cout<<endl;
-      
-      for (Int_t i0=0; i0<24;  i0=i0+2)
+	  cout<<"efficiency for "<<event<<" events"<<endl;
+	  for (Int_t i0=1; i0<13;  i0++)
+	    {
+
+	      cout<<names[i0].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<" ";
+	      cout<<names[i0+13].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<" ";
+	      cout<<names[i0+57].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<" ";
+	      cout<<names[i0+69].Data()<<" "<<Float_t(numberOfHits[i0])/Float_t(event)<<endl;
+
+	      hEffCFD->Fill(i0,Float_t(numberOfHits[i0])  / Float_t(event));
+	      hEffLED->Fill(i0,Float_t(numberOfHits[i0+13]) / Float_t(event));
+	      hEffCFD->Fill(i0+12,Float_t(numberOfHits[i0+57]) /Float_t(event));
+	      hEffLED->Fill(i0+12,Float_t(numberOfHits[i0+69]) /Float_t(event));
+	    }
+	  cout<<endl;
+
+	  for (Int_t i0=0; i0<24;  i0=i0+2)
+	    {
+	      hEffQT1->Fill(i0, Float_t (numberOfHits[i0+25]) / Float_t(event));
+	      hEffQT0->Fill(i0, Float_t (numberOfHits[i0]+26) / Float_t(event));
+	      hEffQT1->Fill((i0+12), Float_t (numberOfHits[i0]+81) /  Float_t(event));
+	      hEffQT0->Fill((i0+12), Float_t (numberOfHits[i0]+82) /  Float_t(event));
+	      cout<<names[i0+25].Data()<<" "<<Float_t(numberOfHits[i0+25])/Float_t(event)<<" ";
+	      cout<<names[i0+26].Data()<<" "<<Float_t(numberOfHits[i0+26])/Float_t(event)<<" ";
+	      cout<<names[i0+81].Data()<<" "<<Float_t(numberOfHits[i0]+81)/Float_t(event)<<" ";
+	      cout<<names[i0+82].Data()<<" "<<Float_t(numberOfHits[i0]+82)/Float_t(event)<<endl;
+
+
+	    }
+	}	      
+
+
+      Char_t filehist[40]; 
+     sprintf(filehist,"t0tree%i.root",fRunNumber);
+     //  sprintf(filehist,"test.root",runNumber);
+     TFile *hist = new TFile(filehist,"RECREATE");
+      cout<<" writing hist in file "<<filehist<<endl;
+      hist->cd();
+      //  digitsTree->Write("",TObject::kOverwrite);
+
+      hEffCFD->Write();
+      hEffLED->Write();
+      hEffQT0->Write();
+      hEffQT1->Write();
+
+      for(Int_t ik=0; ik<105; ik++)	hChannel[ik] ->Write();
+
+      for (Int_t i=0; i<12; i++)
 	{
-	  hEffQT1->Fill(i0, Float_t (numberOfHits[i0+25]) / Float_t(event));
-	  hEffQT0->Fill(i0, Float_t (numberOfHits[i0]+26) / Float_t(event));
-	  hEffQT1->Fill((i0+12), Float_t (numberOfHits[i0]+81) /  Float_t(event));
-	  hEffQT0->Fill((i0+12), Float_t (numberOfHits[i0]+82) /  Float_t(event));
-	  cout<<names[i0+25].Data()<<" "<<Float_t(numberOfHits[i0+25])/Float_t(event)<<" ";
-	  cout<<names[i0+26].Data()<<" "<<Float_t(numberOfHits[i0+26])/Float_t(event)<<" ";
-	  cout<<names[i0+81].Data()<<" "<<Float_t(numberOfHits[i0]+81)/Float_t(event)<<" ";
-	  cout<<names[i0+82].Data()<<" "<<Float_t(numberOfHits[i0]+82)/Float_t(event)<<endl;
-
-	  
+	  hQTC[i]->Write();
+	  hmpd[i]->Write();
+	  hCFD_QTC[i]->Write();
+	  hCFD_LED[i]->Write();
+	  h1CFD_LED[i]->Write();
 	}
-    }	      
-  
-  
-  Char_t filehist[100]; 
-   sprintf(filehist,"/home/t0/alice/testSep07/tree/t0tree%i.root",fRunNumber);
-  //   sprintf(filehist,"test.root");
-  TFile *hist = new TFile(filehist,"RECREATE");
-  cout<<" writing hist in file "<<filehist<<endl;
-  
-  //  digitsTree->Write("",TObject::kOverwrite);
-  
-  hEffCFD->Write();
-  hEffLED->Write();
-  hEffQT0->Write();
-  hEffQT1->Write();
-  
-  for(Int_t ik=0; ik<105; ik++) {	
-    if (hChannel[ik]->GetEntries()>0 ) hChannel[ik] ->Write();
-  }
-  for (Int_t i=0; i<12; i++)
-    {
-      hQTC[i]->Write();
-      hmpd[i]->Write();
-      hCFD_QTC[i]->Write();
-      hCFD_LED[i]->Write();
-      h1CFD_LED[i]->Write();
-    }
-  
+     
+      cout<<" hist in file"<<endl;
 
-  hist->Close();   
-  cout<<" hist in file"<<endl;
- 
 }
-- 
2.43.0