Update from Francesco Noferini: VZERO event plane calibration

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALTriggerRawDigitMaker.cxx

diff --git a/EMCAL/AliEMCALTriggerRawDigitMaker.cxx b/EMCAL/AliEMCALTriggerRawDigitMaker.cxx
index 41ac1e11fc09367d1e983af1854cd096a61f0641..0dee4f20fa348d62205b7555a7e263d27337b8dc 100644 (file)
--- a/EMCAL/AliEMCALTriggerRawDigitMaker.cxx
+++ b/EMCAL/AliEMCALTriggerRawDigitMaker.cxx
@@ -32,6 +32,9 @@ Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
 #include "AliCaloBunchInfo.h"
 #include "AliRawReader.h"
 #include "AliEMCALTriggerDCSConfigDB.h"
+#include "AliEMCALTriggerDCSConfig.h"
+#include "AliEMCALTriggerTRUDCSConfig.h"
+#include "AliEMCALTriggerSTUDCSConfig.h"
 #include "AliEMCALTriggerData.h"
 #include "AliEMCALTriggerPatch.h"
 #include "AliLog.h"
@@ -97,7 +100,8 @@ AliEMCALTriggerRawDigitMaker::~AliEMCALTriggerRawDigitMaker()
 //_______________
 void AliEMCALTriggerRawDigitMaker::SetIO(AliRawReader* reader, AliCaloRawStreamV3& in, AliEMCALTriggerSTURawStream& inSTU, TClonesArray* digits, AliEMCALTriggerData* data)
 {
-       //
+       // Connect I/O
+       
        fRawReader     = reader;
        fCaloRawStream = ∈
        fRawDigits     = digits;
@@ -108,7 +112,8 @@ void AliEMCALTriggerRawDigitMaker::SetIO(AliRawReader* reader, AliCaloRawStreamV
 //_______________
 void AliEMCALTriggerRawDigitMaker::Add(const std::vector<AliCaloBunchInfo> &bunchlist)
 {
-       //
+       // Add bunch list
+       
        Int_t    hwAdd   = fCaloRawStream->GetHWAddress();
        UShort_t iRCU    = fCaloRawStream->GetDDLNumber() % 2; // 0/1
        UShort_t iBranch = ( hwAdd >> 11 ) & 0x1;              // 0/1
@@ -277,40 +282,43 @@ void AliEMCALTriggerRawDigitMaker::Add(const std::vector<AliCaloBunchInfo> &bunc
 //_______________
 void AliEMCALTriggerRawDigitMaker::PostProcess()
 {      
-       //
+       // Post process digits
+       
         AliDebug(2,"Start post processing the raw digit maker");
        Int_t idx;
        
        AliEMCALTriggerRawDigit* dig = 0x0;
        
-       Int_t sizeL1gsubr[2], sizeL1gpatch[2], sizeL1jsubr[2], sizeL1jpatch[2];
+       TVector2 sizeL1gsubr, sizeL1gpatch, sizeL1jsubr, sizeL1jpatch;
        
-       fDCSConfig->GetSTUSegmentation(sizeL1gsubr, sizeL1gpatch, sizeL1jsubr, sizeL1jpatch);
+       fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->GetSegmentation(sizeL1gsubr, sizeL1gpatch, sizeL1jsubr, sizeL1jpatch);
        
        fRawReader->Reset();
        fRawReader->Select("EMCAL",44); 
 
-       Bool_t STUin = kFALSE;
+       Bool_t isSTUin = kFALSE;
        
        Int_t nSubEv = fRawReader->GetEvent()->GetNSubEvents();
        
        for ( Int_t iSubEv=0; iSubEv<nSubEv; iSubEv++)
        {
-               AliRawVEvent *SubEv = ((AliRawEvent*)fRawReader->GetEvent())->GetSubEvent(iSubEv);
-               if ( !SubEv ) continue;
+               AliRawVEvent *subEv = ((AliRawEvent*)fRawReader->GetEvent())->GetSubEvent(iSubEv);
+               if ( !subEv ) continue;
                
-               for (Int_t iEquip = 0; iEquip < SubEv->GetNEquipments(); iEquip++)
+               for (Int_t iEquip = 0; iEquip < subEv->GetNEquipments(); iEquip++)
                {
-                       Int_t eqId = SubEv->GetEquipment(iEquip)->GetEquipmentHeader()->GetId();
+                       Int_t eqId = subEv->GetEquipment(iEquip)->GetEquipmentHeader()->GetId();
                        
-                       if (eqId == kSTUEqId) STUin = kTRUE;
+                       if (eqId == kSTUEqId) isSTUin = kTRUE;
                }
        }
        
        fRawReader->Reset();
        
-       if (STUin && fSTURawStream && fSTURawStream->ReadPayLoad())
+       if (isSTUin && fSTURawStream && fSTURawStream->ReadPayLoad())
        {
+               fTriggerData->SetL1DataDecoded(1);
+               
                fTriggerData->SetL1GammaThreshold(fSTURawStream->GetL1GammaThreshold());
                fTriggerData->SetL1JetThreshold(  fSTURawStream->GetL1JetThreshold()  );
                
@@ -326,12 +334,24 @@ void AliEMCALTriggerRawDigitMaker::PostProcess()
                        fSTURawStream->GetJC(),
                        fSTURawStream->GetRegionEnable(), 
                        fSTURawStream->GetFwVersion()
-               };
+               };              
+
+               // Modify DCS config from STU payload content
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetGA(type[0]);
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetGB(type[1]);
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetGC(type[2]);
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetJA(type[3]);
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetJB(type[4]);
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetJC(type[5]);
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetFw(type[7]);
+               fDCSConfig->GetTriggerDCSConfig()->GetSTUDCSConfig()->SetRawData(fSTURawStream->GetRawData());
                
                fTriggerData->SetL1FrameMask(fSTURawStream->GetFrameReceived());
                fTriggerData->SetL1V0(v0);
                fTriggerData->SetL1TriggerType(type);
                
+               fTriggerData->SetL1RawData(fSTURawStream->GetRawData());
+               
                Int_t iTRU, x, y;
 
                if (fSTURawStream->GetRawData())
@@ -345,20 +365,9 @@ void AliEMCALTriggerRawDigitMaker::PostProcess()
                                UInt_t adc[96]; for (Int_t j = 0; j < 96; j++) adc[j] = 0;
                                
                                fSTURawStream->GetADC(i, adc);
-                               /*
-                               ofstream outfile(Form("data_TRU%d.txt",i),ios_base::trunc);
-                               
-                               for (Int_t j = 0; j < 96; j++) 
-                               {
-                                       outfile << adc[j] << endl;
-                               }
-                               
-                               outfile.close();
-                               */
+
                                for (Int_t j = 0; j < 96; j++)
                                {
-                                       //if (adc[j] < 5) continue;
-                                       
                                        if (AliDebugLevel()) printf("| STU => TRU# %2d raw data: ADC# %2d: %d\n", iTRU, j, adc[j]);
                                        
                                        fGeometry->GetAbsFastORIndexFromTRU(iTRU, j, idx);
@@ -389,13 +398,16 @@ void AliEMCALTriggerRawDigitMaker::PostProcess()
                for (Int_t i = 0; i < fSTURawStream->GetNL0GammaPatch(); i++)
                {
                        fSTURawStream->GetL0GammaPatch(i, iTRU, x);
-                       
+
                        iTRU = fGeometry->GetTRUIndexFromSTUIndex(iTRU);
                        
-                       if (AliDebugLevel()) printf("| STU => Found L0 patch id: %2d in TRU# %2d\n", x, iTRU);
-                       
-                       const Int_t sizePatchL0 = fDCSConfig->GetTRUSegmentation(iTRU) * fDCSConfig->GetTRUSegmentation(iTRU);
+                       const Int_t sizePatchL0 = 
+                       ((AliEMCALTriggerTRUDCSConfig*)fDCSConfig->GetTriggerDCSConfig()->GetTRUArr()->At(fGeometry->GetOnlineIndexFromTRUIndex(iTRU)))->GetSegmentation() 
+                       * 
+                       ((AliEMCALTriggerTRUDCSConfig*)fDCSConfig->GetTriggerDCSConfig()->GetTRUArr()->At(fGeometry->GetOnlineIndexFromTRUIndex(iTRU)))->GetSegmentation();
                        
+                       if (AliDebugLevel()) printf("| STU => Found L0 patch id: %2d in TRU# %2d\n", x, iTRU);
+                                               
                        Int_t idFastOR[4];
                        for (Int_t j = 0; j < 4; j++) idFastOR[j] = -1;
                        
@@ -437,9 +449,9 @@ void AliEMCALTriggerRawDigitMaker::PostProcess()
                                
                                if (iTRU % 2) vx += 24; // C side
                                
-                               vx = vx - sizeL1gsubr[0] * sizeL1gpatch[0] + 1;
-                               
-                               if (vx >= 0) 
+                               vx = vx - int(sizeL1gsubr.X()) * int(sizeL1gpatch.X()) + 1;
+
+                               if (vx >= 0 && vy < 63) 
                                {
                                        if (fGeometry->GetAbsFastORIndexFromPositionInEMCAL(vx, vy, idx))
                                        {
@@ -469,9 +481,9 @@ void AliEMCALTriggerRawDigitMaker::PostProcess()
                        {
                                if (AliDebugLevel()) printf("| STU => Found L1 jet patch at (%2d , %2d)\n", x, y);
                                
-                               Int_t ix = sizeL1jsubr[0] * (11 - y - sizeL1jpatch[0] + 1);
-
-                               Int_t iy = sizeL1jsubr[1] * (15 - x - sizeL1jpatch[1] + 1);
+                               Int_t ix = int(sizeL1jsubr.X()) * (11 - y - int(sizeL1jpatch.X()) + 1);
+                               
+                               Int_t iy = int(sizeL1jsubr.Y()) * (15 - x - int(sizeL1jpatch.Y()) + 1);
                                
                                // FIXME: x = 0 || y = 0 (Olivier's CS) patches a lost?
                                
@@ -504,7 +516,8 @@ void AliEMCALTriggerRawDigitMaker::PostProcess()
 //_______________
 void AliEMCALTriggerRawDigitMaker::Reset()
 {
-       //      
+       // Reset
+       
        for (Int_t i = 0; i < 3072; i++) fRawDigitIndex[i] = -1;
 }