update of readDigits

diff --git a/START/readDigits.C b/START/readDigits.C
index a37a6fb2e6b0fad7867cc200bb903e37ea9e3a57..4f43da924c6e16c0bfdb238527c9285e7f3c1414 100644 (file)
--- a/START/readDigits.C
+++ b/START/readDigits.C
@@ -2,128 +2,92 @@ void readDigits()
 {
  
   // Dynamically link some shared libs
-  if (gClassTable->GetID("AliRun") < 0) {
+  /*  if (gClassTable->GetID("AliRun") < 0) {
     gROOT->LoadMacro("loadlibs.C");
     loadlibs();
   }
-  char filename[100];
-  sprintf(filename,"galice.root");
-  AliRunLoader* rl = AliRunLoader::Open("galice.root",AliConfig::GetDefaultEventFolderName(),"read");
-  if (rl == 0x0)
+  */
+    Float_t c = 0.0299792; // cm/ps
+    Float_t channelWidth = 25;
+    Float_t timeDelay = 150;
+    char filename[100];
+    sprintf(filename,"galice.root");
+    AliRunLoader* rl = AliRunLoader::Open("galice.root",AliConfig::GetDefaultEventFolderName(),"read");
+    if (rl == 0x0)
    {
      cerr<<"Can not open session for file galice.root\n";
      return;
    }
-
-  rl->LoadgAlice();
-  gAlice = rl->GetAliRun();
-  
-  AliSTART* START  = (AliSTART *)gAlice->GetDetector("START");
-  
-  rl->LoadHeader();
-  rl->LoadKinematics("READ");
-  Int_t retval;
-  AliLoader* lstart = rl->GetLoader("STARTLoader");
+    
+    rl->LoadgAlice();
+    gAlice = rl->GetAliRun();
+    
+    AliSTART* START  = (AliSTART *)gAlice->GetDetector("START");
+    
+    rl->LoadHeader();
+    Int_t retval;
+    AliLoader* lstart = rl->GetLoader("STARTLoader");
 
   Int_t iNevents=rl->GetNumberOfEvents();
   cout<<"  nevents   "<<iNevents<<endl;
 
-   char nameTD[8],nameTR[8];
 
-  TH1F *hTimediff = new TH1F("hTimediff","Time difference",100,0,200);
-  TH1F *hTimePs = new TH1F("hTimePs","Time in Ps",100,-200,200);
-  TH1F *hMeanTime = new TH1F("hMeanTime","Time in Ps",100,2000,2500);
-  TH1F *hallADC = new TH1F("hallright","ADC summary right",100,0,200);
-  TH1F *hADCright = new TH1F("hADCright","ADC right",100,0,200);
-  TH1F *hADCleft = new TH1F("hADCleft","ADC left",100,0,200);
-  TH1F *hTimeright = new TH1F("hTimeright","Time right",100,2000.,3200.);
-  TH1F *hTimeleft = new TH1F("hTimeleft","Time left",100,2000.,3200.);
+  TH1F *hTimediff = new TH1F("hTimediff","Time difference",100,5950,6050);
   TH1F *hBestTimeright = new TH1F("hBestTimeright","First time right",
-                                 100,2000.,3200);
+                                 100,8450.,8650);
   TH1F *hBestTimeleft = new TH1F("hBestTimeleft","First time left",
-                                100,11000.,13000.);
-   
+                                100,8450.,8650.);
+  TH1F *hRealVertex = new TH1F("hRealVertex","Real Vertex",100,-15,15);
+  
+  TH1F *hVertex = new TH1F("hVertex","Z position of vertex",   100,-15,15);
   //  digits = new AliSTARTdigit();
-  AliSTARTdigit *digits   ; // digits
-  digits = new AliSTARTdigit();
- 
-  timeRight = new TArrayI(12);
-  timeLeft = new TArrayI(12);
-  ADCRight = new TArrayI(12);
-  ADCLeft = new TArrayI(12);
+  AliSTARTdigit *fDigits   ; // digits
+  fDigits = new AliSTARTdigit();
 
  // Event ------------------------- LOOP  
   for (Int_t j=0; j<iNevents; j++){
-    //  gAlice->GetEvent(j);
     rl->GetEvent(j);
-    //   lstart->Dump();
-    lstart->LoadHits("READ");
+
+    AliHeader *header = gAlice->GetHeader();
+    AliGenEventHeader* genHeader = header->GenEventHeader();
+    TArrayF *o = new TArrayF(3); 
+    genHeader->PrimaryVertex(*o);
+    Float_t zRealVertex=o->At(2);
+    hRealVertex->Fill( zRealVertex);
+
     lstart->LoadDigits("READ");
-    sprintf(nameTD,"START_D_%d",j);
-    printf("%s\n",nameTD);
-  TObject *td = (TObject*)gDirectory->Get(nameTD);
-  // td->Dump();
-    //   cout<<" td "<<td<<endl;
-  td->Read(nameTD);
-   digits->Read(nameTD);
-   //   digits->Read(nameTD);
-  // digits->Print();
-  printf("time %d\n",digits->GetTimeDiff());
-   
-    if(digits->GetTimeDiff()!=999999){
-      Int_t timediff = digits->GetTimeDiff();
-      //     Double_t timePs=(timediff-128)*10.; // time in Ps channel_width =10ps
-      Int_t timePs=(512-timediff)*2.5.; // time in Ps channel_width =10ps
-      cout<<"timediff "<<timediff<<" timePs "<<timePs<<endl;
-      hTimediff->Fill(timediff);
-      hTimePs->Fill(timePs);
-      Int_t mean=digits->GetMeanTime();
-      cout<<" mean "<<mean<<endl;
-      mean=mean*2.5;
-      hMeanTime->Fill(mean);
-      Int_t br=digits->GetBestTimeRight(); 
-      Int_t bl=digits->GetBestTimeLeft();
-      cout<<"BestTimeRight "<<br*2.5<<" BestTimeLeft "<<bl*2.5<<endl;
-      hBestTimeright->Fill(br*2.5); 
-      hBestTimeleft ->Fill(bl*2.5); 
-      digits->GetTimeRight(*timeRight );
-      digits->GetTimeLeft(*timeLeft );
-      digits->GetADCRight(*ADCRight );
-      digits->GetADCLeft(*ADCLeft );
-       for (Int_t i=0; i<12; i++) 
-       {
-         Int_t t=timeRight.At(i);
-         Int_t ADC=ADCRight.At(i);
-         hTimeright->Fill(t*2.5);
-         hADCright->Fill(ADC);
-       }
-       for (Int_t i=0; i<12; i++) 
-       {
-         Int_t ADC=ADCRight.At(i);
-         Int_t t=timeLeft.At(i);
-         hTimeleft->Fill(t*2.5);
-         hADCleft->Fill(ADC);
-       }
+    TTree *digitsTree =  lstart->TreeD();
+    TBranch *brDigits=digitsTree->GetBranch("START");
+    AliSTARTdigit *fDigits = new AliSTARTdigit();
+    if (brDigits) {
+      brDigits->SetAddress(&fDigits);
+    }else{
+      cerr<<"EXEC Branch START digits not found"<<endl;
+      return;
     }
+    brDigits->GetEntry(0);
+    Int_t   besttimeright = fDigits->BestTimeRight();
+    Int_t   besttimeleft = fDigits->BestTimeLeft();
+    Int_t   timeDiff = fDigits->TimeDiff();
+    Int_t    sumMult=   fDigits->SumMult();
+    hTimediff->Fill(timeDiff);
+    hBestTimeright->Fill(besttimeright);
+    hBestTimeleft->Fill(besttimeleft );
+    Float_t vertex= (timeDiff* channelWidth - timeDelay*1000.)*c;
+    cout<<j<<" "<<besttimeright<<" "<< besttimeleft<<" "<<timeDiff<<" "<<vertex<<endl;
+    hVertex->Fill(vertex);
   }
-  
   Hfile = new TFile("Figdigits.root","RECREATE","Histograms for START digits");
   printf("Writting histograms to root file \n");
   Hfile->cd();
   //Create a canvas, set the view range, show histograms
   gStyle->SetOptStat(111111);
   //  TCanvas *c1 = new TCanvas("c1","Alice START Time ",400,10,600,600);
-  hADCright->Write();
-  hADCleft->Write();
-  hTimeright->Write();
-  hTimeleft->Write();
-  hTimePs->SetXTitle("arriving time, ps");
-  hTimePs->SetYTitle("number of events");
-  hTimePs->Write();
-  hMeanTime->Write();
+   hTimediff->Write();
   hBestTimeright->Write(); 
   hBestTimeleft ->Write(); 
-  Hfile->Close();
+  hVertex->Write();
+  hRealVertex->Write();
     
  
 } // end of macro