Mods in QuadraticRoots

[u/mrichter/AliRoot.git] / FMD / AliFMDBaseDA.cxx

diff --git a/FMD/AliFMDBaseDA.cxx b/FMD/AliFMDBaseDA.cxx
index 35f68165b5e5762db760c73b372430f54454cd1a..2cb422c7a4d1859d17f961ad48fc9662d404a8f2 100644 (file)
--- a/FMD/AliFMDBaseDA.cxx
+++ b/FMD/AliFMDBaseDA.cxx
@@ -28,7 +28,6 @@
 
 #include "AliFMDBaseDA.h"
 #include "iostream"
-
 #include "AliFMDRawReader.h"
 #include "AliFMDCalibSampleRate.h"
 #include "AliFMDCalibStripRange.h"
@@ -92,9 +91,9 @@ AliFMDBaseDA::AliFMDBaseDA() :
   fConditionsFile(),
   fSaveHistograms(kFALSE),
   fDetectorArray(),
-  fPulseSize(16),
-  fPulseLength(16),
-  fRequiredEvents(5),
+  fPulseSize(10),
+  fPulseLength(10),
+  fRequiredEvents(0),
   fCurrentEvent(0)
  {
   fDetectorArray.SetOwner();
@@ -134,41 +133,43 @@ void AliFMDBaseDA::Run(AliRawReader* reader)
   
   
   
-  InitContainer(diagFile);
-  Init();
-
-
+  
   reader->Reset();
   
   AliFMDRawReader* fmdReader  = new AliFMDRawReader(reader,0);
   TClonesArray*    digitArray = new TClonesArray("AliFMDDigit",0);
-    
-  //reader->NextEvent(); // Read Start-of-Run event
-  // reader->NextEvent(); // Read Start-of-Files event
-  
-  // reader->NextEvent(); // Start-of-Data
   
+  Bool_t SOD_read = kFALSE;
   
-  int lastProgress = 0;
-  
-  for(Int_t n =1;n <= GetRequiredEvents(); n++) {
-    if(!reader->NextEvent()) continue;
+  for(Int_t i=0;i<3;i++) {
+    reader->NextEvent(); // Read Start-of-Run / Start-of-Files event
     
     UInt_t eventType = reader->GetType();
     if(eventType == AliRawEventHeaderBase::kStartOfData || 
        eventType == AliRawEventHeaderBase::kFormatError) { 
       
       WriteConditionsData(fmdReader);
-      
-    }
-    if(eventType != AliRawEventHeaderBase::kPhysicsEvent)  {
-      n--; 
-      continue; 
+      Init();
+      SOD_read = kTRUE;
+      break;
     }
-    
+  }
+  
+  InitContainer(diagFile);
+  
+  if(!SOD_read) 
+    AliWarning("No SOD event detected!");
+  
+  int lastProgress = 0;
+  
+  
+  
+  for(Int_t n =1;n <= GetRequiredEvents(); n++) {
+    if(!reader->NextEvent()) continue;
     SetCurrentEvent(n);
     digitArray->Clear();
     fmdReader->ReadAdcs(digitArray);
+    
     //std::cout<<"in event "<<*(reader->GetEventId())<<"   "<<n<<std::endl;
     //AliDebug(5, Form("In event # %d with %d entries", 
     //              *(reader->GetEventId()), digitArray->GetEntriesFast()));
@@ -178,13 +179,16 @@ void AliFMDBaseDA::Run(AliRawReader* reader)
       FillChannels(digit);
     }
     
+   
     FinishEvent();
+    
     int progress = int((n *100)/ GetRequiredEvents()) ;
     if (progress <= lastProgress) continue;
     lastProgress = progress;
     std::cout << "Progress: " << lastProgress << " / 100 " << std::endl;
+    
   }
-
+  
   AliInfo(Form("Looped over %d events",GetCurrentEvent()));
   WriteHeaderToFile();
   
@@ -203,24 +207,25 @@ void AliFMDBaseDA::Run(AliRawReader* reader)
        std::cout << '.' << std::flush;
       }
       if(fSaveHistograms)
-       diagFile->Flush();
+       diagFile->Flush();
       std::cout << "done" << std::endl;
     }
   }
-
+  
   if(fOutputFile.is_open()) {
     fOutputFile.write("# EOF\n",6);
     fOutputFile.close();
   }
   
+  Terminate(diagFile);
+    
   if(fSaveHistograms ) {
     
-    Terminate(diagFile);
-  
     AliInfo("Closing diagnostics file - please wait ...");
     // diagFile->Write();
     diagFile->Close();
     AliInfo("done");
+    
   }
 }
 //_____________________________________________________________________
@@ -336,10 +341,46 @@ void AliFMDBaseDA::WriteConditionsData(AliFMDRawReader* fmdReader)
   
   // Strip Range
   
- 
+  fConditionsFile.write("# Gain Events \n",15);
   
+  for(UShort_t det=1; det<=3;det++) {
+    UShort_t firstring = (det == 1 ? 1 : 0);
+    for(UShort_t iring = firstring; iring <=1;iring++) {
+      Char_t ring = (iring == 1 ? 'I' : 'O');
+      for(UShort_t board =0 ; board <=1; board++) {
+       
+       Int_t idx = GetHalfringIndex(det,ring,board);
+       
+       fConditionsFile << det                     << ','
+                       << ring                    << ','
+                       << board                   << ','
+                       << fPulseLength.At(idx)    << "\n";
+       
+      }
+    }
+  }
+  
+  fConditionsFile.write("# Gain Pulse \n",14);
+  
+  for(UShort_t det=1; det<=3;det++) {
+    UShort_t firstring = (det == 1 ? 1 : 0);
+    for(UShort_t iring = firstring; iring <=1;iring++) {
+      Char_t ring = (iring == 1 ? 'I' : 'O');
+      for(UShort_t board =0 ; board <=1; board++) {
+       
+       Int_t idx = GetHalfringIndex(det,ring,board);
+       
+       fConditionsFile << det                     << ','
+                       << ring                    << ','
+                       << board                   << ','
+                       << fPulseSize.At(idx)      << "\n";
+       
+      }
+    }
+  }
   
   
+
   // Gain Relevant stuff
   /*
   UShort_t defPulseSize = 32 ; 
@@ -379,7 +420,7 @@ Int_t AliFMDBaseDA::GetHalfringIndex(UShort_t det, Char_t ring, UShort_t board)
   
   Int_t index = (((det-1) << 2) | (iring << 1) | (board << 0));
   
-  return index;
+  return index-2;
   
 }