Fixing FORWARD_NULL defects reported by Coverity
[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerElectronics.cxx
index aaf996e..388c344 100644 (file)
 
 /* $Id$ */
 
-//*-- Author: Rachid Guernane (LPCCFd)
-//*   Manager class for muon trigger electronics
-//*   Client of trigger board classes
-//*
-//*
+//-----------------------------------------------------------------------------
+// Class AliMUONTriggerElectronics
+//--------------------------------
+// Manager class for muon trigger electronics
+// Client of trigger board classes
+// Debugged by Ph. Crochet & Ch. Finck
+// Interfaced with new mapping Ch. Finck
+//
+// Author: Rachid Guernane (LPCCFd)
+//-----------------------------------------------------------------------------
 
-#include "AliMUONTriggerElectronics.h"
-#include "AliMUONTriggerCrate.h"
-#include "AliMUONTriggerCrateStore.h"
-#include "AliMUONConstants.h"
-#include "AliMUONLocalTriggerBoard.h"
-#include "AliMUONRegionalTriggerBoard.h"
+#include "AliLoader.h"
+#include "AliLog.h"
+#include "AliMUONCalibrationData.h"
+#include "AliMUONVDigit.h"
+#include "AliMUONVDigitStore.h"
+#include "AliMUONGlobalTrigger.h"
 #include "AliMUONGlobalTriggerBoard.h"
 #include "AliMUONLocalTrigger.h"
+#include "AliMUONLocalTriggerBoard.h"
 #include "AliMUONRegionalTrigger.h"
-#include "AliMUONGlobalTrigger.h"
-#include "AliMUON.h" 
-#include "AliMUONData.h" 
-#include "AliMUONDigit.h"
-#include "AliMUONSegmentation.h"
-#include "AliMUONCalibrationData.h"
+#include "AliMUONRegionalTriggerBoard.h"
+#include "AliMUONTriggerCrate.h"
+#include "AliMUONTriggerCrateStore.h"
+#include "AliMUONTriggerElectronics.h"
+#include "AliMUONTriggerCrateConfig.h"
+#include "AliMUONRegionalTriggerConfig.h"
+#include "AliMUONGlobalCrateConfig.h"
+#include "AliMUONVTriggerStore.h"
 #include "AliMUONVCalibParam.h"
-
+#include "AliMpCathodType.h"
+#include "AliMpCDB.h"
+#include "AliMpDEManager.h"
+#include "AliMpSegmentation.h"
 #include "AliMpVSegmentation.h"
+#include "AliMpCathodType.h"
+#include "AliMpTriggerCrate.h"
+#include "AliMpLocalBoard.h"
+#include "AliMpDDLStore.h"
+#include "AliMpExMap.h"
+#include "AliMpIntPair.h"
 
 #include "AliLog.h"
 #include "AliLoader.h"
 #include "AliRun.h"
+#include <TBits.h>
+#include <TSystem.h>
+
+#include "AliCodeTimer.h"
 
-//#include "Riostream.h"
-#include "TBits.h"
-#include "TSystem.h"
 
+/// \cond CLASSIMP
 ClassImp(AliMUONTriggerElectronics)
+/// \endcond
 
 //___________________________________________
-AliMUONTriggerElectronics::AliMUONTriggerElectronics(AliMUONData *Data, AliMUONCalibrationData* calibData) 
-: TTask("AliMUONTriggerElectronics",
-        "From trigger digits to Local and Global Trigger objects"),
+AliMUONTriggerElectronics::AliMUONTriggerElectronics(AliMUONCalibrationData* calibData) 
+: TObject(),
   fCrates(new AliMUONTriggerCrateStore),
-  fGlobalTriggerBoard(new AliMUONGlobalTriggerBoard),
-  fMUONData(Data)
+  fGlobalTriggerBoard(new AliMUONGlobalTriggerBoard)
 {
-//* CONSTRUCTOR
-//*
-  if (!fMUONData)
-  {  
-    AliFatal("NO MUON TRIGGER DATA");
+ /// CONSTRUCTOR
+///
+
+  for (Int_t i = 0; i < 2; ++i) {
+    fCopyXInput[i] = new TList();
+    fCopyXInput[i]->SetOwner();
+    fCopyYInput[i] = new TList();
+    fCopyYInput[i]->SetOwner();
   }
-    
-  SetDataSource();
-  Factory(calibData);
-  LoadMasks(calibData);
-}
-
-//______________________________________________________________________________
-AliMUONTriggerElectronics::AliMUONTriggerElectronics(const AliMUONTriggerElectronics& right) 
-  : TTask(right) 
-{  
-/// Protected copy constructor (not implemented)
 
-  AliFatal("Copy constructor not provided.");
+  // force loading of mapping if not already done
+  if ( !AliMpDDLStore::Instance(kFALSE) )
+  {
+    AliMpCDB::LoadDDLStore();
+  }
+  
+  SetCopyInput();
+  
+  Factory(calibData);
+  LoadMasks(calibData); 
 }
 
 //___________________________________________
 AliMUONTriggerElectronics::~AliMUONTriggerElectronics()
 {
-//* DESTRUCTOR
-//*
+/// DESTRUCTOR
+///
   delete fGlobalTriggerBoard;
   delete fCrates;
-}
-
-//______________________________________________________________________________
-AliMUONTriggerElectronics& 
-AliMUONTriggerElectronics::operator=(const AliMUONTriggerElectronics& right)
-{
-/// Protected assignement operator (not implemented)
+  for (Int_t i = 0; i < 2; ++i) {
+    delete fCopyXInput[i];
+    delete fCopyYInput[i];
+  }
 
-  // check assignement to self
-  if (this == &right) return *this;
+}
 
-  AliFatal("Assignement operator not provided.");
+//___________________________________________
+void AliMUONTriggerElectronics::SetCopyInput()
+{  
+  /// set list of copy input
+  
+    for (Int_t iDDL = 0; iDDL < 2; ++iDDL) { 
     
-  return *this;  
-}    
+      for(Int_t iReg = 0; iReg < 8; ++iReg){   //reg loop
+      
+       AliMpTriggerCrate* crateMapping = AliMpDDLStore::Instance()->GetTriggerCrate(iDDL, iReg);
+      
+       for(Int_t iLocal = 0; iLocal < crateMapping->GetNofLocalBoards(); ++iLocal) { 
+        
+         Int_t localBoardFromId = crateMapping->GetLocalBoardId(iLocal);
+         if (!localBoardFromId) continue; //empty slot, should not happen
+        
+          AliMpLocalBoard* localBoardFrom = AliMpDDLStore::Instance()->GetLocalBoard(localBoardFromId);
+         Int_t localBoardToId;
+         if ((localBoardToId = localBoardFrom->GetInputXto())) {
+             AliMpLocalBoard* localBoardTo = AliMpDDLStore::Instance()->GetLocalBoard(localBoardToId);
+             TString crateFrom = localBoardFrom->GetCrate();
+             Int_t   slotFrom  = localBoardFrom->GetSlot();
+             TString crateTo   = localBoardTo->GetCrate();
+             Int_t   slotTo    = localBoardTo->GetSlot();
+          
+             fCopyXInput[0]->Add(new AliMpIntPair(AliMpExMap::GetIndex(crateFrom), slotFrom));
+             fCopyXInput[1]->Add(new AliMpIntPair(AliMpExMap::GetIndex(crateTo), slotTo));
+             AliDebug(3, Form("copy xInputs from local  %s_%d to %s_%d\n", crateFrom.Data(), slotFrom, 
+                              crateTo.Data(), slotTo));
+         }
+        
+         if ((localBoardToId = localBoardFrom->GetInputYto())) {
+             AliMpLocalBoard* localBoardTo = AliMpDDLStore::Instance()->GetLocalBoard(localBoardToId);
+             TString crateFrom = localBoardFrom->GetCrate();
+             Int_t   slotFrom  = localBoardFrom->GetSlot();
+             TString crateTo   = localBoardTo->GetCrate();
+             Int_t   slotTo    = localBoardTo->GetSlot();
+          
+             fCopyYInput[0]->Add(new AliMpIntPair(AliMpExMap::GetIndex(crateFrom), slotFrom));
+             fCopyYInput[1]->Add(new AliMpIntPair(AliMpExMap::GetIndex(crateTo), slotTo));
+             AliDebug(3, Form("copy yInputs from local  %s_%d to %s_%d\n", crateFrom.Data(), slotFrom, 
+                              crateTo.Data(), slotTo));
+          
+         }
+       }
+      }
+    }
+}
 
 //___________________________________________
 void AliMUONTriggerElectronics::Factory(AliMUONCalibrationData* calibData)
 {  
- //* BUILD ALL ELECTRONICS
- //*
-
-// get coinc44 from AliMUON (added 12/09/06)
-  AliMUON *pMUON  = (AliMUON*)gAlice->GetModule("MUON");
-  Int_t coinc44 = pMUON->GetTriggerCoinc44();
-  if (coinc44 != 0 && coinc44 != 1) {
-      AliFatal("Coinc 44 should be equal to 0 or 1");
-      return;
-  }
+ /// BUILD ALL ELECTRONICS
+ ///
 
-  fCrates->ReadFromFile(gSystem->ExpandPathName(fSourceFileName.Data()));
-  
-  if ( !calibData ) return;
-  
-  AliMUONTriggerLut* lut = calibData->TriggerLut();
-  
-  if (!lut) return;
-  
-  AliMUONLocalTriggerBoard* localBoard;
-  
-  fCrates->FirstLocalBoard();
-  
-  while ( (localBoard=fCrates->NextLocalBoard()) )
-  {
-    localBoard->SetLUT(lut);
-    localBoard->SetCoinc44(coinc44);
-  }
+    fCrates->ReadFromFile(calibData);
 }
 
 //___________________________________________
-void AliMUONTriggerElectronics::FeedM()
+void AliMUONTriggerElectronics::Feed(const AliMUONVDigitStore& digitStore)
 {
-//* FILL INPUTS
-//*
-       for (Int_t ichamber=10; ichamber<14; ichamber++) 
-       {
-      TClonesArray *muonDigits = fMUONData->Digits(ichamber);
-      Int_t ndigits = muonDigits->GetEntriesFast();
-
-      for (Int_t digit=0; digit<ndigits; digit++)
-               {
-                       AliMUONDigit *mdig = static_cast<AliMUONDigit*>(muonDigits->UncheckedAt(digit));
-
-//       CHECKME ! The TrackCharge is not ok with new digitizerV3 !
-//                     for (Int_t ichg=0; ichg<10; ichg++) schg += mdig->TrackCharge(ichg);
-//       assert(schg==mdig->Signal());
-                       Int_t schg = mdig->Signal();
-         
-//       APPLY CONDITION ON SOFT BACKGROUND    
-                       Int_t tchg = schg - (Int_t(schg/10))*10;        
-
-                       if (schg<=10 || tchg>0) 
-                       {
-//                             mdig->Print();
+  /// FILL INPUTS
+  ///
 
-                               Int_t digitindex = digit;
-                               Int_t detElemId  = mdig->DetElemId();
-                               Int_t cathode    = mdig->Cathode();
-
-                               const AliMpVSegmentation *seg = ((AliMUON*)gAlice->GetDetector("MUON"))->GetSegmentation()->GetMpSegmentation(detElemId,cathode);
-
-                               Int_t ix = mdig->PadX(), iy = mdig->PadY();
-                               
-                               AliDebug(3,Form("cathode %d ix %d iy %d ",cathode,ix,iy));
+  AliCodeTimerAuto("",0);
+  
+  TIter next(digitStore.CreateTriggerIterator());
+  AliMUONVDigit* mdig;
+  
+  while ( ( mdig = static_cast<AliMUONVDigit*>(next()) ) )
+  {      
+    //       CHECKME ! The TrackCharge is not ok with new digitizerV3 !
+    //                 for (Int_t ichg=0; ichg<10; ichg++) schg += mdig->TrackCharge(ichg);
+    Int_t ichamber = AliMpDEManager::GetChamberId(mdig->DetElemId());
+    Int_t schg = (Int_t)(mdig->Charge() + 0.5);
+    
+    //       APPLY CONDITION ON SOFT BACKGROUND        
+    Int_t tchg = schg - (Int_t(schg/10))*10;   
+    
+    if (schg<=10 || tchg>0) 
+    {
+      Int_t detElemId  = mdig->DetElemId();
+      Int_t cathode    = mdig->Cathode();
+    
+      const AliMpVSegmentation* seg = 
+         AliMpSegmentation::Instance()
+         ->GetMpSegmentation(detElemId,AliMp::GetCathodType(cathode));
+  
+      Int_t ix = mdig->PadX(), iy = mdig->PadY();
+      
+      AliDebug(3,Form("cathode %d ix %d iy %d ",cathode,ix,iy));
 
-                               AliMpPad pad = seg->PadByIndices(AliMpIntPair(ix,iy),kTRUE);
-                               
-                               for (Int_t i=0; i<pad.GetNofLocations(); i++) 
-                               {
-                                       AliMpIntPair location = pad.GetLocation(i);
-                                       
-                                       Int_t nboard = location.GetFirst();
-
-                                       Int_t ibitxy = location.GetSecond();
-
-                                       AliMUONLocalTriggerBoard *b = fCrates->LocalBoard(nboard);
-
-                                       if (b) 
-                                       {
-                                               if (cathode && b->GetSwitch(6)) ibitxy += 8;
-                                               
-                                               b->SetbitM(ibitxy,cathode,ichamber-10);
-                                               
-                                               DigitFiredCircuit(b->GetI(), cathode, ichamber, digitindex);
-                                       }
-          else
-          {
-            AliError(Form("Could not get local board number %d",b->GetNumber()));
-          }
-                               }
-                       }                       
-               }
-       }
+      AliMpPad pad = seg->PadByIndices(ix,iy,kTRUE);
+      
+      for (Int_t i=0; i<pad.GetNofLocations(); i++) 
+      {
+        Int_t nboard = pad.GetLocalBoardId(i);
+        
+        Int_t ibitxy = pad.GetLocalBoardChannel(i);
+        
+        AliMUONLocalTriggerBoard *b = fCrates->LocalBoard(nboard);
+        
+        if (b) 
+        {
+          if (cathode && b->GetSwitch(AliMpLocalBoard::kZeroAllYLSB)) ibitxy += 8;
+          
+          b->SetbitM(ibitxy,cathode,ichamber-10);
+        }
+        else
+        {
+          AliError(Form("Could not get local board number %d",nboard));
+        }
+      }
+    }          
+  }
 
-// Particular case of the columns with 22 local boards (2R(L) 3R(L))   
-       AliMUONTriggerCrate *crate = 0x0; TObjArray *bs = 0x0;
+  FeedCopyNeighbours();
+}
 
-       char *scratess[4] = {  "2R",   "2L",   "3L",   "3R"}; 
-       char *scratesd[4] = {"2-3R", "2-3L", "2-3L", "2-3R"}; 
-       Int_t    slotf[4] = {     2,      2,     10,     10}; 
-       Int_t    slotd[4] = {     1,      1,      9,      9}; 
 
-       for (Int_t i=0; i<4; i++)
-       {
-      crate = fCrates->Crate(scratess[i]); 
-      bs = crate->Boards();
-      AliMUONLocalTriggerBoard *desybb = (AliMUONLocalTriggerBoard*)bs->At(14);
-      AliMUONLocalTriggerBoard *fromcb = (AliMUONLocalTriggerBoard*)bs->At(15);
-      AliMUONLocalTriggerBoard *desxbb = (AliMUONLocalTriggerBoard*)bs->At(16);
+//___________________________________________
+void AliMUONTriggerElectronics::FeedCopyNeighbours()
+{
+  //
+  /// Feed the local copies
+  /// and complete the feed with the information of neighbours
+  //
 
-      crate = fCrates->Crate(scratesd[i]); 
-      bs = crate->Boards();
-      AliMUONLocalTriggerBoard *frombb = (AliMUONLocalTriggerBoard*)bs->At(slotf[i]);
-      AliMUONLocalTriggerBoard *desycb = (AliMUONLocalTriggerBoard*)bs->At(slotd[i]);
+  // Particular case of the columns with 22 local boards (2R(L) 3R(L))   
+  // fill copy input from mapping instead of hardcoded valued (Ch.F)
+  AliMUONTriggerCrate *crate = 0x0; TObjArray *bs = 0x0;
 
-      UShort_t cX[2];
+  for (Int_t i = 0; i < fCopyXInput[0]->GetEntries(); ++i) 
+  {
+    AliMpIntPair* pair = (AliMpIntPair*)fCopyXInput[0]->At(i);
+    TString crateFrom  =  AliMpExMap::GetString(pair->GetFirst());
+    Int_t   slotFrom   =  pair->GetSecond();
 
-//    COPY X3-4 FROM BOARD  2 OF CRATE 2-3 TO BOARD 16 OF CRATE 2
-//    COPY X3-4 FROM BOARD 10 OF CRATE 2-3 TO BOARD 16 OF CRATE 3
-      frombb->GetX34(cX); desxbb->SetX34(cX);
+    pair = (AliMpIntPair*)fCopyXInput[1]->At(i);
+    TString crateTo  =  AliMpExMap::GetString(pair->GetFirst());
+    Int_t   slotTo   =  pair->GetSecond();
 
-//    COPY X3-4 FROM BOARD 15 OF CRATE 2 TO BOARD 1 OF CRATE 2-3
-//    COPY X3-4 FROM BOARD 15 OF CRATE 3 TO BOARD 9 OF CRATE 2-3
-      fromcb->GetX34(cX); desycb->SetX34(cX);
+    AliDebug(3, Form("copy xInputs from local  %s_%d to %s_%d\n", crateFrom.Data(), slotFrom, 
+                    crateTo.Data(), slotTo));
 
-      UShort_t cY[4];
+    UShort_t cX[2];
+    crate = fCrates->Crate(crateFrom); 
+    bs = crate->Boards();
+    AliMUONLocalTriggerBoard *fromxb = (AliMUONLocalTriggerBoard*)bs->At(slotFrom);
+    crate = fCrates->Crate(crateTo); 
+    bs = crate->Boards();
+    AliMUONLocalTriggerBoard *desxb = (AliMUONLocalTriggerBoard*)bs->At(slotTo);
+    fromxb->GetX34(cX); desxb->SetX34(cX);
 
-      desybb->GetY(cY); frombb->SetY(cY);
 
-      frombb->GetY(cY); desxbb->SetY(cY);
-      fromcb->GetY(cY); desycb->SetY(cY);
-       }
+  }
 
-// FILL UP/DOWN OF CURRENT BOARD (DONE VIA J3 BUS IN REAL LIFE)
- AliMUONTriggerCrate* cr;
- fCrates->FirstCrate();
- while ( ( cr = fCrates->NextCrate() ) )
-       {            
-               TObjArray *boards = cr->Boards();
-               
-               for (Int_t j=1; j<boards->GetEntries()-1; j++)
-               {
-                       TObject *o = boards->At(j);
+  for (Int_t i = 0; i < fCopyYInput[0]->GetEntries(); ++i) 
+  {
+    AliMpIntPair* pair = (AliMpIntPair*)fCopyYInput[0]->At(i);
+    TString crateFrom  =  AliMpExMap::GetString(pair->GetFirst());
+    Int_t   slotFrom   =  pair->GetSecond();
+
+    pair = (AliMpIntPair*)fCopyYInput[1]->At(i);
+    TString crateTo  =  AliMpExMap::GetString(pair->GetFirst());
+    Int_t   slotTo   =  pair->GetSecond();
+
+    AliDebug(3, Form("copy yInputs from local  %s_%d to %s_%d\n", crateFrom.Data(), slotFrom, 
+                    crateTo.Data(), slotTo));
+
+    UShort_t cY[4];
+    crate = fCrates->Crate(crateFrom); 
+    bs = crate->Boards();
+    AliMUONLocalTriggerBoard *fromyb = (AliMUONLocalTriggerBoard*)bs->At(slotFrom);
+    crate = fCrates->Crate(crateTo); 
+    bs = crate->Boards();
+    AliMUONLocalTriggerBoard *desyb = (AliMUONLocalTriggerBoard*)bs->At(slotTo);
+    fromyb->GetY(cY); desyb->SetY(cY);
+  }
+  
+  // FILL UP/DOWN OF CURRENT BOARD (DONE VIA J3 BUS IN REAL LIFE)
+  AliMUONTriggerCrate* cr;
+  TIter next2(fCrates->CreateCrateIterator());
+  
+  while ( ( cr = static_cast<AliMUONTriggerCrate*>(next2()) ) )
+  {            
+    TObjArray *boards = cr->Boards();
+    
+    for (Int_t j = 1; j < boards->GetEntries()-1; j++)
+    {
+      TObject *o = boards->At(j);
                        
-                       if (!o) break;
+      if (!o) break;
                        
-                       AliMUONLocalTriggerBoard *currboard = (AliMUONLocalTriggerBoard*)o;
+      AliMUONLocalTriggerBoard *currboard = (AliMUONLocalTriggerBoard*)o;
                        
-                       AliMUONLocalTriggerBoard *neighbour = (AliMUONLocalTriggerBoard*)boards->At(j+1);
+      AliMUONLocalTriggerBoard *neighbour = (AliMUONLocalTriggerBoard*)boards->At(j+1);
                        
-                       UShort_t cXY[2][4];
+      UShort_t cXY[2][4];
                        
-                       if (j==1) {neighbour->GetXY(cXY); currboard->SetXYU(cXY);}
+      if (j==1) {neighbour->GetXY(cXY); currboard->SetXYU(cXY);}
                        
-//       LAST BOARD IN THE CRATE HAS NO UP EXCEPT FOR CRATES 2 & 3
-                       if (j<boards->GetEntries()-2)  
-                       {
-                               AliMUONLocalTriggerBoard *nextboard = (AliMUONLocalTriggerBoard*)boards->At(j+2);
+      //       LAST BOARD IN THE CRATE HAS NO UP EXCEPT FOR CRATES 2 & 3
+      if (j < boards->GetEntries()-2)  
+      {
+             AliMUONLocalTriggerBoard *nextboard = (AliMUONLocalTriggerBoard*)boards->At(j+2);
                                
-                               currboard->GetXY(cXY); neighbour->SetXYD(cXY);
-                               nextboard->GetXY(cXY); neighbour->SetXYU(cXY);
+             currboard->GetXY(cXY); neighbour->SetXYD(cXY);
+             nextboard->GetXY(cXY); neighbour->SetXYU(cXY);
                                
-                               if (j==boards->GetEntries()-3) {neighbour->GetXY(cXY); nextboard->SetXYD(cXY);}
-                       }
-               }
-       }
+             if (j==boards->GetEntries()-3) {neighbour->GetXY(cXY); nextboard->SetXYD(cXY);}
+      }
+    }
+  }
 }
 
+
 //___________________________________________
 void AliMUONTriggerElectronics::Feed(UShort_t pattern[2][4])
 {
-  //* FILL INPUTS
-  //*
+  /// FILL INPUTS
+  ///
   AliMUONTriggerCrate* cr;
+  TIter next(fCrates->CreateCrateIterator());
    
-   fCrates->FirstCrate();
-   
-   while ( ( cr = fCrates->NextCrate() ) )
+   while ( ( cr = static_cast<AliMUONTriggerCrate*>(next()) ) )
    {                 
      TObjArray *boards = cr->Boards();
      
-     for (Int_t j=1; j<boards->GetEntries(); j++)
+     for (Int_t j = 1; j < boards->GetEntries(); j++)
      {
        TObject *o = boards->At(j);
        
@@ -305,31 +359,30 @@ void AliMUONTriggerElectronics::Feed(UShort_t pattern[2][4])
 //___________________________________________
 void AliMUONTriggerElectronics::DumpOS()
 {
-//* DUMP IN THE OLD WAY
-//*
-   for (Int_t i=0;i<234;i++)
+/// DUMP IN THE OLD WAY
+///
+   for (Int_t i= 0; i < 234;i++)
    {
       AliMUONLocalTriggerBoard *board = fCrates->LocalBoard(i);
 
-                       if (board) board->Scan("ALL");
+      if (board) board->Scan("ALL");
    }
 }
 
 //___________________________________________
-void AliMUONTriggerElectronics::Scan(Option_t *option)
+void AliMUONTriggerElectronics::Scan(const Option_t *option)
 {
-  //* SCAN
-  //*
+  /// SCAN
+  ///
 
   AliMUONTriggerCrate* cr;
+  TIter next(fCrates->CreateCrateIterator());  
   
-  fCrates->FirstCrate();
-  
-  while ( ( cr = fCrates->NextCrate() ) )
-       {                
+  while ( ( cr = static_cast<AliMUONTriggerCrate*>(next()) ) )
+  {                
     TObjArray *boards = cr->Boards();
     
-    for (Int_t j=0; j<boards->GetEntries(); j++)
+    for (Int_t j = 0; j < boards->GetEntries(); j++)
     {
       TObject *o = boards->At(j);
       
@@ -353,14 +406,12 @@ void AliMUONTriggerElectronics::Scan(Option_t *option)
 //___________________________________________
 void AliMUONTriggerElectronics::Reset()
 {
-  //* RESET
-  //*
+  /// RESET
+  ///
   
    AliMUONTriggerCrate* cr;
-   
-   fCrates->FirstCrate();
-   
-   while ( ( cr = fCrates->NextCrate() ) )
+  TIter next(fCrates->CreateCrateIterator());
+   while ( ( cr = static_cast<AliMUONTriggerCrate*>(next()) ) )
    {            
       TObjArray *boards = cr->Boards();
             
@@ -373,36 +424,44 @@ void AliMUONTriggerElectronics::Reset()
    }
 }
 
+
 //_______________________________________________________________________
 void AliMUONTriggerElectronics::LoadMasks(AliMUONCalibrationData* calibData)
 {
-  // LOAD MASKS FROM CDB
+  /// Load mask from config in CDB 
   
-
-  // SET MASKS
+  // Set mask
   
-  AliMUONTriggerCrate* cr;
+  AliMUONRegionalTriggerConfig* regionalConfig = calibData->RegionalTriggerConfig();
+  if (!regionalConfig) {
+     AliError("No valid regional trigger configuration in CDB");
+     return;
+  }   
   
-  fCrates->FirstCrate();
+  AliMUONTriggerCrate* cr;
+  TIter next(fCrates->CreateCrateIterator());
   
   Int_t irb(0);
   
-  while ( ( cr = fCrates->NextCrate() ) )
-       {            
+  while ( ( cr = static_cast<AliMUONTriggerCrate*>(next()) ) )
+  {            
     TObjArray *boards = cr->Boards();
     
-    AliMUONRegionalTriggerBoard *regb =
-      (AliMUONRegionalTriggerBoard*)boards->At(0);
+    AliMUONRegionalTriggerBoard *regb = (AliMUONRegionalTriggerBoard*)boards->At(0);
 
-    AliMUONVCalibParam* regionalBoardMasks = calibData->RegionalTriggerBoardMasks(irb);
+    AliMUONTriggerCrateConfig* crateConfig = regionalConfig->FindTriggerCrate(cr->GetName());
     
-    for ( Int_t i = 0; i < regionalBoardMasks->Size(); ++i )
+    if (!crateConfig)
     {
-      UShort_t rmask = static_cast<UShort_t>(regionalBoardMasks->ValueAsInt(i) & 0x3F);
-      regb->Mask(i,rmask);
+      AliError(Form("Crate %s not present in configuration !!!", cr->GetName()));
+      return;
     }
     
-    for (Int_t j=1; j<boards->GetEntries(); j++)
+    UShort_t rmask= crateConfig->GetMask();
+
+    regb->Mask(rmask);
+    
+    for (Int_t j = 1; j < boards->GetEntries(); j++)
     {
       AliMUONLocalTriggerBoard *b = (AliMUONLocalTriggerBoard*)boards->At(j);
       
@@ -421,87 +480,63 @@ void AliMUONTriggerElectronics::LoadMasks(AliMUONCalibrationData* calibData)
     ++irb;
   }
   
-  AliMUONVCalibParam* globalBoardMasks = calibData->GlobalTriggerBoardMasks();
-  for ( Int_t i = 0; i < globalBoardMasks->Size(); ++i )
-  {
-    UShort_t gmask = static_cast<UShort_t>(globalBoardMasks->ValueAsInt(i) & 0xFFF);
-    fGlobalTriggerBoard->Mask(i,gmask);
-  }
-}
+   AliMUONGlobalCrateConfig * globalConfig = calibData->GlobalTriggerCrateConfig();
+   if (!globalConfig) {
+     AliError("No valid trigger crate configuration in CDB");
+     return;
+   }   
 
+    UInt_t gmask = 0;
+    for (Int_t i = 0; i < 4; i++) {
+      gmask = globalConfig->GetGlobalMask(i);
+      fGlobalTriggerBoard->Mask(i,gmask);
+    }
+}
 
 //___________________________________________
 void AliMUONTriggerElectronics::LocalResponse()
 {
-  // INTERFACE BOARDS
-       struct crates_t 
-  {
-    TString name;
-    Int_t slots[5];
-    Int_t ns;
-  } crate[6];
-
-       crate[0].name = "2R";   crate[0].ns = 1; crate[0].slots[0] = 16;
-       crate[1].name = "2L";   crate[1].ns = 1; crate[1].slots[0] = 16;
-       crate[2].name = "3L";   crate[2].ns = 1; crate[2].slots[0] = 16;
-       crate[3].name = "3R";   crate[3].ns = 1; crate[3].slots[0] = 16;
-       crate[4].name = "2-3R"; crate[4].ns = 2; crate[4].slots[0] = 1;  crate[4].slots[1] = 9;
-       crate[5].name = "2-3L"; crate[5].ns = 2; crate[5].slots[0] = 1;  crate[5].slots[1] = 9; 
+/// Compute the response for local cards
+
+  AliCodeTimerAuto("",0);
        
   AliMUONTriggerCrate* cr;
+  TIter next(fCrates->CreateCrateIterator());
+
+  UShort_t thisl[16];
   
-  fCrates->FirstCrate();
-  
-  while ( ( cr = fCrates->NextCrate() ) )
-       {            
-    Int_t iib = -1;
-    
-    for (Int_t icr=0; icr<6; icr++) 
-    {
-                       const char *n = (crate[icr].name).Data();
-                       
-      AliMUONTriggerCrate *dcr = fCrates->Crate(n);
-      
-      //       THIS CRATE CONTAINS AN INTERFACE BOARD
-      if ( dcr && !strcmp(cr->GetName(),dcr->GetName()) ) iib = icr;
-    }
+  while ( ( cr = static_cast<AliMUONTriggerCrate*>(next()) ) )
+  {            
     
     TObjArray *boards = cr->Boards();
     
     AliMUONRegionalTriggerBoard *regb = (AliMUONRegionalTriggerBoard*)boards->At(0);
     
-    UShort_t thisl[16]; for (Int_t j=0; j<16; j++) thisl[j] = 0;
-        
-    for (Int_t j=1; j<boards->GetEntries(); j++)
+    for (Int_t j=0; j<16; ++j) thisl[j] = 0;
+  
+    for (Int_t j = 1; j < boards->GetEntries(); j++)
     {     
-      TObject *o = boards->At(j);
-      
-      if (!o) break;
+       TObject *o = boards->At(j);
       
-      AliMUONLocalTriggerBoard *board = (AliMUONLocalTriggerBoard*)o;
+       if (!o) break;
       
-      if (board) 
-      {
-        board->Response();
+       AliMUONLocalTriggerBoard *board = (AliMUONLocalTriggerBoard*)o;
+
+       board->Response();
                                
-        UShort_t tmp = board->GetResponse();            
+       UShort_t response = board->GetResponse();            
         
-        //          CRATE CONTAINING INTERFACE BOARD
-        if ( iib>-1 ) 
-        {
-          for (Int_t iid = 0; iid<crate[iib].ns; iid++) 
-                                       {
-            if ( j == crate[iib].slots[iid] )
-                                               {
-              if ( tmp != 0 ) 
-                AliWarning(Form("Interface board %s in slot %d of crate %s has a non zero response",
-                                board->GetName(),j,cr->GetName()));
-                                               }
-                                       }                                       
-        }
+       // CRATE CONTAINING INTERFACE BOARD
+       if (board->GetNumber() == 0) // copy boards
+       {
+         if ( response != 0 ) 
+           AliWarning(Form("Interface board %s in slot %d of crate %s has a non zero response",
+                           board->GetName(),j,cr->GetName()));
+         AliDebug(1, Form("local slot %d, number %d in crate %s\n", j, board->GetNumber(), cr->GetName()));
+         
+       }
         
-        thisl[j-1] = tmp;
-      }
+       thisl[j-1] = response;
     }
     
     regb->SetLocalResponse(thisl);
@@ -512,20 +547,20 @@ void AliMUONTriggerElectronics::LocalResponse()
 void AliMUONTriggerElectronics::RegionalResponse()
 {
   /// Compute the response for all regional cards.
+
+  AliCodeTimerAuto("",0);
+
   AliMUONTriggerCrate* cr;
+  TIter next(fCrates->CreateCrateIterator());
   
-  fCrates->FirstCrate();
-  
-  while ( ( cr = fCrates->NextCrate() ) )
-       {            
+  while ( ( cr = static_cast<AliMUONTriggerCrate*>(next()) ) )
+  {            
       TObjArray *boards = cr->Boards();
 
       AliMUONRegionalTriggerBoard *regb = (AliMUONRegionalTriggerBoard*)boards->At(0);
-      
-      if (regb) 
-      {
-         regb->Response();
-      }  
+
+      regb->Response();
+
    }
 }
 
@@ -534,11 +569,11 @@ void AliMUONTriggerElectronics::GlobalResponse()
 {
   /// Compute the global response
 
+  AliCodeTimerAuto("",0);
+
   UShort_t regional[16];
   
   AliMUONTriggerCrate* cr;
-  
-  fCrates->FirstCrate();
   Int_t irb(0);
   
   if ( !fCrates->NumberOfCrates() >= 16 ) 
@@ -546,260 +581,260 @@ void AliMUONTriggerElectronics::GlobalResponse()
     AliFatal(Form("Something is wrong : too many crates %d",
                   fCrates->NumberOfCrates()));
   }
+
+  // send regional responses to the global trigger in right order
+  // do not used iterator order
   
-  while ( ( cr = fCrates->NextCrate() ) )
-       {            
-    AliMUONTriggerBoard* rb = 
-      static_cast<AliMUONTriggerBoard*>(cr->Boards()->At(0));
-    regional[irb] = rb->GetResponse();
-    ++irb;
+  for (Int_t iSide = 0; iSide < 2; iSide++) // right & left side
+  {            
+    for (Int_t iReg = 0; iReg < 8; iReg++) // 8 crates/regional boards for each side.
+    {
+      cr = fCrates->Crate(iSide, iReg);     
+
+      AliMUONTriggerBoard* rb = 
+       static_cast<AliMUONTriggerBoard*>(cr->Boards()->At(0));
+      regional[irb] = rb->GetResponse();
+      ++irb;
+    }
   }
-  
+
   fGlobalTriggerBoard->SetRegionalResponse(regional);
   fGlobalTriggerBoard->Response();
 }
 
-//___________________________________________
-void AliMUONTriggerElectronics::BoardName(Int_t ix, Int_t iy, char *name)
-{
-//* BOARD NAME FROM PAD INFO (OLD MAPPING)
-//*
-   TString s = (ix>0) ? "R" : "L"; 
-
-   Int_t board = iy / 16, bid[4] = {12,34,56,78}; 
-
-   ix = abs(ix);
-
-   Int_t line = ix / 10, column = ix - 10 * line;
-
-// old scheme: line==1 is line==9
-   line -= 9; line = TMath::Abs(line); line++;
-
-   sprintf(name,"%sC%dL%dB%d", s.Data(), column, line, bid[board]);
-   
-   AliDebug(3, Form("Strip ( %d , %d ) connected to board %s ", ix, iy, name));
-}
-
-//___________________________________________
-void AliMUONTriggerElectronics::BuildName(Int_t icirc, char name[20])
-{
-//* GET BOARD NAME FROM OLD NUMBERING
-//*
-   const Int_t kCircuitId[234] = 
-      {
-          111,  121,  131,  141,  151,  161,  171,
-          211,  212,  221,  222,  231,  232,  241,  242,  251,  252,  261,  262,  271,
-          311,  312,  321,  322,  331,  332,  341,  342,  351,  352,  361,  362,  371,
-          411,  412,  413,  421,  422,  423,  424,  431,  432,  433,  434,  441,  442,  451,  452,  461,  462,  471,
-          521,  522,  523,  524,  531,  532,  533,  534,  541,  542,  551,  552,  561,  562,  571, 
-          611,  612,  613,  621,  622,  623,  624,  631,  632,  633,  634,  641,  642,  651,  652,  661,  662,  671,
-          711,  712,  721,  722,  731,  732,  741,  742,  751,  752,  761,  762,  771,
-          811,  812,  821,  822,  831,  832,  841,  842,  851,  852,  861,  862,  871,
-          911,  921,  931,  941,  951,  961,  971,
-         -111, -121, -131, -141, -151, -161, -171,
-         -211, -212, -221, -222, -231, -232, -241, -242, -251, -252, -261, -262, -271,
-         -311, -312, -321, -322, -331, -332, -341, -342, -351, -352, -361, -362, -371,
-         -411, -412, -413, -421, -422, -423, -424, -431, -432, -433, -434, -441, -442, -451, -452, -461, -462, -471,
-         -521, -522, -523, -524, -531, -532, -533, -534, -541, -542, -551, -552, -561, -562, -571, 
-         -611, -612, -613, -621, -622, -623, -624, -631, -632, -633, -634, -641, -642, -651, -652, -661, -662, -671,
-         -711, -712, -721, -722, -731, -732, -741, -742, -751, -752, -761, -762, -771,
-         -811, -812, -821, -822, -831, -832, -841, -842, -851, -852, -861, -862, -871,
-         -911, -921, -931, -941, -951, -961, -971 
-      };
-
-   Int_t b[4] = {12, 34, 56, 78};
-
-   Int_t code = TMath::Abs(kCircuitId[icirc]);
-
-   Int_t lL = code / 100;
-
-   Int_t cC = ( code - 100 * lL ) / 10;
-   
-   Int_t bB = code - 100 * lL - 10 * cC;
-   
-   const char *side = (kCircuitId[icirc]>0) ? "R" : "L";
-
-// lL=1 AT TOP
-   lL -= 9; lL = abs(lL); lL++;
-
-   sprintf(name,"%sC%dL%dB%d",side,cC,lL,b[bB-1]);
-}
-
-//_______________________________________________________________________
-void 
-AliMUONTriggerElectronics::Exec(Option_t*)
-{
-//*
-//*
-  Digits2Trigger();
-}
-
 //_______________________________________________________________________
-void AliMUONTriggerElectronics::Trigger()
+void AliMUONTriggerElectronics::Digits2Trigger(const AliMUONVDigitStore& digitStore,
+                                               AliMUONVTriggerStore& triggerStore)
 {
-//*
-//*
-   FeedM();
-   LocalResponse();
-   RegionalResponse();      
-   GlobalResponse();
-}
+  AliCodeTimerAuto("",0);
 
-//_______________________________________________________________________
-void AliMUONTriggerElectronics::Digits2Trigger()
-{
   /// Main method to go from digits to trigger decision
-
-  AliMUONRegionalTrigger *pRegTrig = new AliMUONRegionalTrigger();
-
-  ClearDigitNumbers();
+  AliMUONRegionalTrigger pRegTrig;
   
-  fMUONData->ResetTrigger(); 
+  triggerStore.Clear();
+
+  // NOW RESET ELECTRONICS
+  Reset();
   
   // RUN THE FULL BEE CHAIN
-  Trigger();
-  //   DumpOS();
+  Feed(digitStore);
+  LocalResponse();
+  RegionalResponse();      
+  GlobalResponse();
+  //    DumpOS();
        
   AliMUONTriggerCrate* cr;
+  AliMUONLocalTrigger localTrigger;
   
-  fCrates->FirstCrate();
+  // stored in right order
+  // do not used iterator order
   
-  while ( ( cr = fCrates->NextCrate() ) )
+  for (Int_t iSide = 0; iSide < 2; iSide++) // right & left side
   {            
-    TObjArray *boards = cr->Boards();
-
-    UInt_t regInpLpt = 0;
-    UInt_t regInpHpt = 0;
-    UShort_t localMask = 0x0;
-
-    AliMUONRegionalTriggerBoard *regBoard = (AliMUONRegionalTriggerBoard*)boards->At(0);
-
-    for (Int_t j = 1; j < boards->GetEntries(); j++)
-    {     
-      TObject *o = boards->At(j);
+    for (Int_t iReg = 0; iReg < 8; iReg++) // 8 crates/regional boards for each side.
+    {
+      cr = fCrates->Crate(iSide, iReg);     
+      TObjArray *boards = cr->Boards();
       
-      if (!o) break;
+      UInt_t regInpLpt = 0;
+      UInt_t regInpHpt = 0;
       
-      AliMUONLocalTriggerBoard *board = (AliMUONLocalTriggerBoard*)o;
+      AliMUONRegionalTriggerBoard *regBoard = (AliMUONRegionalTriggerBoard*)boards->At(0);
       
-      if (board) 
-      {
-        //          L0 TRIGGER
-        if (board->Triggered())
+      for (Int_t j = 1; j < boards->GetEntries(); j++)
+      {     
+        TObject *o = boards->At(j);
+        
+        if (!o) break;
+        
+        AliMUONLocalTriggerBoard *board = (AliMUONLocalTriggerBoard*)o;
+        
+        if (board) 
         {
-          Int_t localtr[15] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ,0};
+          //          L0 TRIGGER
+          // pcrochet 181206: MOOD needs ALL boards
+          //     if (board->Triggered())
+          //     {
           
           Int_t icirc = board->GetNumber();
-          
-          localtr[0] = icirc;
-          localtr[1] = board->GetStripX11();
-          localtr[2] = board->GetDev();
-          localtr[3] = board->GetStripY11();
-          
-          //             SAVE LUT OUTPUT 
-          UShort_t response = board->GetResponse();
-          localtr[4] = (response & 12) >> 2;
-          localtr[5] = (response & 48) >> 4;
-          localtr[6] = (response &  3);
-
-         // calculates regional inputs from local for the moment
-         UInt_t hPt = (response >> 4) & 0x3;
-         UInt_t lPt = (response >> 2) & 0x3;
-           
-         regInpHpt |= hPt << (30 - (j-1)*2);
-         regInpLpt |= lPt << (30 - (j-1)*2);
-         localMask |= (0x1 << (j-1)); // local mask
-
-
-          TBits rrr;
-          rrr.Set(6,&response);
-          
-          //             SAVE BIT PATTERN
-          localtr[7]  = board->GetXY(0,0);
-          localtr[8]  = board->GetXY(0,1);
-          localtr[9]  = board->GetXY(0,2);
-          localtr[10] = board->GetXY(0,3);
-          
-          localtr[11] = board->GetXY(1,0);
-          localtr[12] = board->GetXY(1,1);
-          localtr[13] = board->GetXY(1,2);
-          localtr[14] = board->GetXY(1,3);
-          
-          //             ADD A NEW LOCAL TRIGGER
-          AliMUONLocalTrigger *pLocTrig = new AliMUONLocalTrigger(localtr, fDigitNumbers[icirc]);
-          
-          fMUONData->AddLocalTrigger(*pLocTrig);  
-         delete pLocTrig;
+          if (icirc != 0) { // pcrochet 181206: MOOD needs ALL boards
+            
+            localTrigger.SetLoCircuit(icirc);
+            localTrigger.SetLoStripX(board->GetStripX11());
+            localTrigger.SetLoDev(board->GetDev());
+            localTrigger.SetLoSdev(board->GetSdev());
+            localTrigger.SetLoTrigY(board->GetTrigY());
+            localTrigger.SetLoStripY(board->GetStripY11());
+            
+            //             SAVE LUT OUTPUT 
+            UShort_t response = board->GetResponse();
+            localTrigger.SetLoHpt((response & 12) >> 2);
+            localTrigger.SetLoLpt(response &  3);
+            
+            // calculates regional inputs from local for the moment
+            UInt_t hPt = (response >> 2) & 0x3;
+            UInt_t lPt =  response       & 0x3;
+            
+            regInpHpt |= hPt << (30 - (j-1)*2);
+            regInpLpt |= lPt << (30 - (j-1)*2);
+            
+            TBits rrr;
+            rrr.Set(6,&response);        
+            
+            //             SAVE BIT PATTERN
+            localTrigger.SetX1Pattern(board->GetXY(0,0));
+            localTrigger.SetX2Pattern(board->GetXY(0,1));
+            localTrigger.SetX3Pattern(board->GetXY(0,2));
+            localTrigger.SetX4Pattern(board->GetXY(0,3));
+            
+            localTrigger.SetY1Pattern(board->GetXY(1,0));
+            localTrigger.SetY2Pattern(board->GetXY(1,1));
+            localTrigger.SetY3Pattern(board->GetXY(1,2));
+            localTrigger.SetY4Pattern(board->GetXY(1,3));
+            
+            //             ADD A NEW LOCAL TRIGGER          
+            triggerStore.Add(localTrigger);  
+            
+          }
+          }
         }
+      pRegTrig.SetId(iReg + 8*iSide);
+      pRegTrig.SetLocalOutput(regInpLpt, 0);
+      pRegTrig.SetLocalOutput(regInpHpt, 1);
+      pRegTrig.SetOutput(regBoard->GetResponse());
+      
+      triggerStore.Add(pRegTrig);  
       }
     }
-    pRegTrig->SetLocalOutput(regInpLpt, 0);
-    pRegTrig->SetLocalOutput(regInpHpt, 1);
-    pRegTrig->SetLocalMask(localMask);
-    pRegTrig->SetOutput((regBoard->GetResponse() >> 4) & 0xF); // to be uniformized (oct06 ?)
-
-    fMUONData->AddRegionalTrigger(*pRegTrig);  
-
-  }
-  delete pRegTrig;
-  
-  // GLOBAL TRIGGER INFORMATION: [0] -> LOW PT 
-  //                             [1] -> HIGH PT
-  //                             [2] -> ALL PT 
-  Int_t globalSinglePlus[3], globalSingleMinus[3], globalSingleUndef[3]; 
-  Int_t globalPairUnlike[3], globalPairLike[3];   
   
+  // GLOBAL TRIGGER INFORMATION
   UShort_t global = fGlobalTriggerBoard->GetResponse();
+  UInt_t *globalInput = fGlobalTriggerBoard->GetGlobalInput();  
+
+  AliMUONGlobalTrigger globalTrigger;
   
-  globalPairUnlike[0] = (global &  16) >> 4;
-  globalPairUnlike[1] = (global & 256) >> 8;
-  globalPairUnlike[2] = (global &   1);
-  
-  globalPairLike[0] = (global &  32) >> 5;
-  globalPairLike[1] = (global & 512) >> 9;
-  globalPairLike[2] = (global &   2) >> 1;
-  
-  globalSinglePlus[0] = ((global &  192) >>  6) == 2;
-  globalSinglePlus[1] = ((global & 3072) >> 10) == 2;
-  globalSinglePlus[2] = ((global &   12) >>  2) == 2;
-  
-  globalSingleMinus[0] = ((global &  192) >>  6) == 1;
-  globalSingleMinus[1] = ((global & 3072) >> 10) == 1;
-  globalSingleMinus[2] = ((global &   12) >>  2) == 1;
-  
-  globalSingleUndef[0] = ((global &  192) >>  6) == 3;
-  globalSingleUndef[1] = ((global & 3072) >> 10) == 3;
-  globalSingleUndef[2] = ((global &   12) >>  2) == 3;
-  
-  AliMUONGlobalTrigger *pGloTrig = new AliMUONGlobalTrigger(globalSinglePlus, globalSingleMinus,
-                                                            globalSingleUndef, globalPairUnlike, 
-                                                            globalPairLike);
-  
+  globalTrigger.SetFromGlobalResponse(global);
+  globalTrigger.SetFromGlobalInput(globalInput);
   // ADD A LOCAL TRIGGER IN THE LIST 
-  fMUONData->AddGlobalTrigger(*pGloTrig);
-  delete pGloTrig;
+  triggerStore.SetGlobal(globalTrigger);
 
-  // NOW RESET ELECTRONICS
-  Reset();
 }
 
-//_______________________________________________________________________
-void AliMUONTriggerElectronics::ClearDigitNumbers()
+//___________________________________________
+void AliMUONTriggerElectronics::Feed(const AliMUONVTriggerStore& triggerStore)
 {
-// RESET fDigitNumbers
-       for (Int_t i=0; i<AliMUONConstants::NTriggerCircuit(); i++) fDigitNumbers[i].Set(0);
+  //
+  /// Fill inputs from reconstructed local trigger store
+  //
+  AliMUONLocalTrigger* locTrg;
+  TIter next(triggerStore.CreateLocalIterator());
+  TArrayS xyPattern[2];
+  UShort_t xy[2][4];
+  Int_t loCircuit;
+  while ( ( locTrg = static_cast<AliMUONLocalTrigger*>( next() )) != NULL ){
+    locTrg->GetXPattern(xyPattern[0]);
+    locTrg->GetYPattern(xyPattern[1]);
+    loCircuit = locTrg->LoCircuit();
+    AliMUONLocalTriggerBoard* localBoard = fCrates->LocalBoard(loCircuit);
+    for (Int_t icath = 0; icath<2; ++icath){
+      for (Int_t ich = 0; ich < 4; ++ich){
+       xy[icath][ich] = xyPattern[icath][ich];
+      }
+    }
+    localBoard->SetXY(xy);
+  }
+
+  FeedCopyNeighbours();
 }
 
 //_______________________________________________________________________
-void AliMUONTriggerElectronics::DigitFiredCircuit(Int_t circuit, Int_t cathode,
-                                                  Int_t chamber, Int_t digit)
+Bool_t AliMUONTriggerElectronics::ModifiedLocalResponse(Int_t loCircuit,
+                                                       Bool_t& bendingPlaneResp,
+                                                       Bool_t& nonBendingPlaneResp,
+                                                       Bool_t isCoinc44,
+                                                       Int_t removeChamber)
 {
-// REGISTERS THAT THE SPECIFIED DIGIT FIRED THE SPECIFIED CIRCUIT
-// THIS DIGIT GETS ADDED TO AN ARRAY WHICH WILL BE COPIED TO
-// AliMUONLocalTrigger WHEN SUCH AN OBJECT IS CREATED FOR EACH CIRCUIT
-       Int_t digitnumber = AliMUONLocalTrigger::EncodeDigitNumber(chamber, cathode, digit);
-       Int_t last = fDigitNumbers[circuit].GetSize();
-       fDigitNumbers[circuit].Set(last + 1);
-       fDigitNumbers[circuit][last] = digitnumber;
+  //
+  /// Re-compute the local trigger response
+  /// with some modifications (i.e. setting coinc44 or after removing one chamber)
+  //
+
+  bendingPlaneResp = kFALSE;
+  nonBendingPlaneResp = kFALSE;
+
+  Bool_t isTriggered = kFALSE;
+
+  AliMUONLocalTriggerBoard* currBoard = fCrates->LocalBoard(loCircuit);
+
+  if ( ! currBoard ) return isTriggered;
+
+  AliMUONLocalTriggerBoard localBoard (*currBoard);
+
+  if (removeChamber>=0 && removeChamber<=3){
+
+    // Set the bit pattern of selected chamber to 0
+    UShort_t xy[2][4];
+    UShort_t xyu[2][4];
+    UShort_t xyd[2][4];
+
+    localBoard.GetXY(xy);
+    localBoard.GetXYU(xyu);
+    localBoard.GetXYD(xyd);
+
+    for(Int_t icath=0; icath<2; icath++){
+      xy[icath][removeChamber] = 0;
+      xyu[icath][removeChamber] = 0;
+      xyd[icath][removeChamber] = 0;
+    }
+
+    localBoard.SetXY(xy);
+    localBoard.SetXYU(xyu);
+    localBoard.SetXYD(xyd);
+  }
+
+  localBoard.ResetResponse();
+
+  localBoard.SetCoinc44((Int_t)isCoinc44);
+  localBoard.Response();
+
+  bendingPlaneResp = localBoard.IsTrigX();
+  nonBendingPlaneResp = localBoard.IsTrigY();
+  isTriggered = localBoard.Triggered();
+
+  return isTriggered;
 }
 
+
+//_______________________________________________________________________
+void AliMUONTriggerElectronics::ResponseRemovingChambers(AliMUONVTriggerStore& triggerStore)
+{
+  /// Update local board information with the trigger response after removing each chamber
+
+  AliCodeTimerAuto("", 0);
+
+  Reset();
+  Feed(triggerStore);
+
+  AliMUONLocalTrigger* locTrg;
+  TIter next(triggerStore.CreateLocalIterator());
+  Int_t loCircuit;
+  Bool_t planeResp[2], isTrig44;
+  Bool_t bendPlaneRespNoCh, nonBendPlaneRespNoCh, isTrigWithoutCh;
+  while ( ( locTrg = static_cast<AliMUONLocalTrigger*>( next() )) != NULL ){
+    if ( ! ( locTrg->IsTrigX() && locTrg->IsTrigY() ) ) continue;
+    loCircuit = locTrg->LoCircuit();
+    isTrig44 = ModifiedLocalResponse(loCircuit, planeResp[0], planeResp[1], kTRUE);
+    for (Int_t ich=0; ich<4; ++ich){
+      if ( ! isTrig44 ){
+       isTrigWithoutCh = ModifiedLocalResponse(loCircuit, bendPlaneRespNoCh, nonBendPlaneRespNoCh, kFALSE, ich);
+       if ( ! isTrigWithoutCh ) continue;
+       for (Int_t icath=0; icath<2; icath++){
+         if ( ! planeResp[icath] )
+           locTrg->SetNoHitInPlane(icath, ich);
+       } // loop on cathodes
+      }
+      locTrg->SetTriggerWithoutChamber(ich);
+    } // loop on chambers
+    AliDebug(1, Form("Is44 %i  triggers %i  pattern %i", isTrig44, locTrg->GetTriggerWithoutChamber(), locTrg->GetHitPatternFromResponse()));
+  }
+}