Fixed memory leaks for #86360: High memory consumption in 2.76TeV p+p RAW reco jobs
[u/mrichter/AliRoot.git] / PHOS / AliPHOSRawDigiProducer.cxx
index d21da2c..5e54017 100644 (file)
@@ -155,8 +155,18 @@ AliPHOSRawDigiProducer::~AliPHOSRawDigiProducer()
 //--------------------------------------------------------------------------------------
 void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawFitterv0* fitter) 
 {
+  // Create a temporary array of LG digits and then make digits from raw data
+
+  TClonesArray *tmpLG = new TClonesArray("AliPHOSDigit",10000) ;
+  MakeDigits(digits, tmpLG, fitter);
+  tmpLG->Delete();
+  delete tmpLG;
+}
+//--------------------------------------------------------------------------------------
+void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, TClonesArray *tmpDigLG, AliPHOSRawFitterv0* fitter) 
+{
   //Makes the job.
-  //TClonesArray *digits and raw data fitter should be provided by calling function.
+  //TClonesArray *digits, *tmpDigLG and raw data fitter should be provided by calling function.
 
   digits->Clear();
  
@@ -173,8 +183,8 @@ void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawFitterv0
   else
     fSampleToSec=fRawStream->GetTSample() ;
   
-  //Temporary array for LowGain digits
-  TClonesArray tmpLG("AliPHOSDigit",10000) ;
+  // Clear a temporary array for LowGain digits
+  tmpDigLG->Clear();
   Int_t ilgDigit=0 ;
 
   //Let fitter subtract pedestals in case of ZS
@@ -207,9 +217,9 @@ void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawFitterv0
       fGeom->RelToAbsNumbering(relId, absId);
       
       fitter->SetNBunches(0);
-      Int_t sigStart  ;
-      Int_t sigLength ;
-      while (fRawStream->NextBunch()) { //Take first in time banch
+      Int_t sigStart =0 ;
+      Int_t sigLength=0 ;
+      while (fRawStream->NextBunch()) { //Take the first in time bunch
        const UShort_t *sig = fRawStream->GetSignals();
        sigStart  = fRawStream->GetStartTimeBin();
        sigLength = fRawStream->GetBunchLength();
@@ -244,17 +254,23 @@ void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawFitterv0
        continue ;
       
       energy = CalibrateE(energy,relId,!caloFlag) ;
-//       time   = CalibrateT(time,relId,lowGainFlag) ;
 
       //convert time from sample bin units to s
       time*=fSampleToSec ;
+//CalibrateT moved to Clusterizer
+//      time = CalibrateT(time,relId,!caloFlag) ;
+      // subtract RCU L1 phase (L1Phase is in seconds) w.r.t. L0:
+      time -= fRawStream->GetL1Phase();
+
+
       
       if(energy <= 0.) 
        continue;
       
       if (caloFlag == AliCaloRawStreamV3::kLowGain) {
-       new(tmpLG[ilgDigit]) AliPHOSDigit(-1,absId,(Float_t)energy,(Float_t)time);
-       dynamic_cast<AliPHOSDigit*>(tmpLG.At(ilgDigit))->SetALTROSamplesHG(sigLength,fADCValuesLG);
+       new((*tmpDigLG)[ilgDigit]) AliPHOSDigit(-1,absId,(Float_t)energy,(Float_t)time);
+       if (sigLength>0 && fADCValuesLG!=0)
+         static_cast<AliPHOSDigit*>(tmpDigLG->At(ilgDigit))->SetALTROSamplesLG(sigLength,fADCValuesLG);
        ilgDigit++ ; 
       }
       else if (caloFlag == AliCaloRawStreamV3::kHighGain) {
@@ -263,7 +279,8 @@ void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawFitterv0
          new((*digits)[iDigit]) AliPHOSDigit(-1,absId,-1.f,(Float_t)time);
        else
          new((*digits)[iDigit]) AliPHOSDigit(-1,absId,(Float_t)energy,(Float_t)time);
-       dynamic_cast<AliPHOSDigit*>(digits->At(iDigit))->SetALTROSamplesHG(sigLength,fADCValuesLG);
+       if (sigLength>0 && fADCValuesHG!=0)
+         static_cast<AliPHOSDigit*>(digits->At(iDigit))->SetALTROSamplesHG(sigLength,fADCValuesHG);
        iDigit++;
       }
     } // End of NextChannel()
@@ -272,17 +289,17 @@ void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawFitterv0
     //replace energy of HighGain digits only if there is overflow
     //negative energy (overflow)
     digits->Sort() ;
-    tmpLG.Sort() ;
+    tmpDigLG->Sort() ;
     Int_t iLG = 0;
-    Int_t nLG1 = tmpLG.GetEntriesFast()-1 ;
+    Int_t nLG1 = tmpDigLG->GetEntriesFast()-1 ;
     
     for(Int_t iDig=0 ; iDig < digits->GetEntriesFast() ; iDig++) { 
       AliPHOSDigit * digHG = dynamic_cast<AliPHOSDigit*>(digits->At(iDig)) ;
       if (!digHG) continue;
-      AliPHOSDigit * digLG = dynamic_cast<AliPHOSDigit*>(tmpLG.At(iLG)) ;
+      AliPHOSDigit * digLG = dynamic_cast<AliPHOSDigit*>(tmpDigLG->At(iLG)) ;
       while(digLG && iLG<nLG1 && digHG->GetId()> digLG->GetId()){
        iLG++ ;
-       digLG = dynamic_cast<AliPHOSDigit*>(tmpLG.At(iLG)) ;
+       digLG = dynamic_cast<AliPHOSDigit*>(tmpDigLG->At(iLG)) ;
       }
       absId=digHG->GetId() ;
       fGeom->AbsToRelNumbering(absId,relId) ;
@@ -325,7 +342,6 @@ Double_t AliPHOSRawDigiProducer::CalibrateE(Double_t amp, Int_t* relId, Bool_t i
 Double_t AliPHOSRawDigiProducer::CalibrateT(Double_t time, Int_t * relId, Bool_t /* isLowGain */)
 {
   //Calibrate time
-  time*=fPulseGenerator->GetRawFormatTimeTrigger() ;
   if(fgCalibData){
     Int_t   module = relId[0];
     Int_t   column = relId[3];