]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MUON/AliMUONTriggerElectronics.cxx
Fixed calculation of fluctuations of energy loss in absrber.
[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerElectronics.cxx
index 88dcaf2ad4e9b3bb0ac911cc29e77346d1b04d5a..9104af8f87ebe7d6d98132fcc22d3d137c7e2018 100644 (file)
@@ -28,7 +28,6 @@
 
 #include "AliLoader.h"
 #include "AliLog.h"
-#include "AliMUON.h" 
 #include "AliMUONCalibrationData.h"
 #include "AliMUONVDigit.h"
 #include "AliMUONVDigitStore.h"
@@ -47,6 +46,7 @@
 #include "AliMUONVTriggerStore.h"
 #include "AliMUONVCalibParam.h"
 #include "AliMpCathodType.h"
+#include "AliMpCDB.h"
 #include "AliMpDEManager.h"
 #include "AliMpSegmentation.h"
 #include "AliMpVSegmentation.h"
@@ -63,6 +63,8 @@
 #include <TBits.h>
 #include <TSystem.h>
 
+#include "AliCodeTimer.h"
+
 
 /// \cond CLASSIMP
 ClassImp(AliMUONTriggerElectronics)
@@ -77,14 +79,13 @@ AliMUONTriggerElectronics::AliMUONTriggerElectronics(AliMUONCalibrationData* cal
  /// CONSTRUCTOR
 ///
 
-  for (Int_t i = 0; i < 2; ++i) {
-    fCopyXInput[i] = new TList();
-    fCopyXInput[i]->SetOwner();
-    fCopyYInput[i] = new TList();
-    fCopyYInput[i]->SetOwner();
-  }
 
-  SetCopyInput();
+  // force loading of mapping if not already done
+  if ( !AliMpDDLStore::Instance(kFALSE) )
+  {
+    AliMpCDB::LoadDDLStore();
+  }
+  
   
   Factory(calibData);
   LoadMasks(calibData); 
@@ -97,61 +98,7 @@ AliMUONTriggerElectronics::~AliMUONTriggerElectronics()
 ///
   delete fGlobalTriggerBoard;
   delete fCrates;
-  for (Int_t i = 0; i < 2; ++i) {
-    delete fCopyXInput[i];
-    delete fCopyYInput[i];
-  }
-
-}
-
-//___________________________________________
-void AliMUONTriggerElectronics::SetCopyInput()
-{  
-  /// set list of copy input
-  
-    for (Int_t iDDL = 0; iDDL < 2; ++iDDL) { 
     
-      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));
-          
-         }
-       }
-      }
-    }
 }
 
 //___________________________________________
@@ -168,6 +115,8 @@ void AliMUONTriggerElectronics::Feed(const AliMUONVDigitStore& digitStore)
 {
   /// FILL INPUTS
   ///
+
+  AliCodeTimerAuto("",0);
   
   TIter next(digitStore.CreateTriggerIterator());
   AliMUONVDigit* mdig;
@@ -207,69 +156,40 @@ void AliMUONTriggerElectronics::Feed(const AliMUONVDigitStore& digitStore)
         
         if (b) 
         {
-          if (cathode && b->GetSwitch(6)) ibitxy += 8;
+          if (cathode && b->GetSwitch(AliMpLocalBoard::kZeroAllYLSB)) ibitxy += 8;
           
           b->SetbitM(ibitxy,cathode,ichamber-10);
+          
+          if ( cathode == 0 ) {
+            // Particular case of the columns with 22 local boards (2R(L) 3R(L))
+            // Fill copy boards
+            AliMpLocalBoard* mpLocalBoard = AliMpDDLStore::Instance()->GetLocalBoard(nboard);
+            Int_t nboardCopy = mpLocalBoard->GetInputXto();
+            if ( nboardCopy > 0 ) {
+              AliMUONLocalTriggerBoard* copyBoard = fCrates->LocalBoard(nboardCopy);
+              copyBoard->SetbitM(ibitxy,cathode,ichamber-10);
+            }
+          }
         }
         else
         {
-          AliError(Form("Could not get local board number %d",b->GetNumber()));
+          AliError(Form("Could not get local board number %d",nboard));
         }
       }
     }          
   }
 
-  // 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;
-
-  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();
-
-    pair = (AliMpIntPair*)fCopyXInput[1]->At(i);
-    TString crateTo  =  AliMpExMap::GetString(pair->GetFirst());
-    Int_t   slotTo   =  pair->GetSecond();
-
-    AliDebug(3, Form("copy xInputs from local  %s_%d to %s_%d\n", crateFrom.Data(), slotFrom, 
-                    crateTo.Data(), slotTo));
-
-    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);
-
+  FeedCopyNeighbours();
+}
 
-  }
 
-  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);
-  }
+//___________________________________________
+void AliMUONTriggerElectronics::FeedCopyNeighbours()
+{
+  //
+  /// Feed the local copies
+  /// and complete the feed with the information of neighbours
+  //
   
   // FILL UP/DOWN OF CURRENT BOARD (DONE VIA J3 BUS IN REAL LIFE)
   AliMUONTriggerCrate* cr;
@@ -348,7 +268,7 @@ void AliMUONTriggerElectronics::DumpOS()
 }
 
 //___________________________________________
-void AliMUONTriggerElectronics::Scan(Option_t *option)
+void AliMUONTriggerElectronics::Scan(const Option_t *option)
 {
   /// SCAN
   ///
@@ -411,9 +331,10 @@ void AliMUONTriggerElectronics::LoadMasks(AliMUONCalibrationData* calibData)
   // Set mask
   
   AliMUONRegionalTriggerConfig* regionalConfig = calibData->RegionalTriggerConfig();
-  if (!regionalConfig)
-     AliWarning("No valid regional trigger configuration in CDB");
-
+  if (!regionalConfig) {
+     AliError("No valid regional trigger configuration in CDB");
+     return;
+  }   
   
   AliMUONTriggerCrate* cr;
   TIter next(fCrates->CreateCrateIterator());
@@ -444,7 +365,7 @@ void AliMUONTriggerElectronics::LoadMasks(AliMUONCalibrationData* calibData)
       
       Int_t cardNumber = b->GetNumber();
       
-      if (cardNumber) // interface board are not interested
+      if (cardNumber) // skip empty slots
       {
         AliMUONVCalibParam* localBoardMasks = calibData->LocalTriggerBoardMasks(cardNumber);
         for ( Int_t i = 0; i < localBoardMasks->Size(); ++i )
@@ -458,8 +379,10 @@ void AliMUONTriggerElectronics::LoadMasks(AliMUONCalibrationData* calibData)
   }
   
    AliMUONGlobalCrateConfig * globalConfig = calibData->GlobalTriggerCrateConfig();
-  if (!globalConfig)
-     AliWarning("No valid trigger crate configuration in CDB");
+   if (!globalConfig) {
+     AliError("No valid trigger crate configuration in CDB");
+     return;
+   }   
 
     UInt_t gmask = 0;
     for (Int_t i = 0; i < 4; i++) {
@@ -472,9 +395,13 @@ void AliMUONTriggerElectronics::LoadMasks(AliMUONCalibrationData* calibData)
 void AliMUONTriggerElectronics::LocalResponse()
 {
 /// Compute the response for local cards
+
+  AliCodeTimerAuto("",0);
        
   AliMUONTriggerCrate* cr;
   TIter next(fCrates->CreateCrateIterator());
+
+  UShort_t thisl[16];
   
   while ( ( cr = static_cast<AliMUONTriggerCrate*>(next()) ) )
   {            
@@ -483,7 +410,7 @@ void AliMUONTriggerElectronics::LocalResponse()
     
     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=0; j<16; ++j) thisl[j] = 0;
   
     for (Int_t j = 1; j < boards->GetEntries(); j++)
     {     
@@ -492,26 +419,23 @@ void AliMUONTriggerElectronics::LocalResponse()
        if (!o) break;
       
        AliMUONLocalTriggerBoard *board = (AliMUONLocalTriggerBoard*)o;
-      
-       if (board) // check if empty slot
-       {
-         board->Response();
+
+       board->Response();
                                
-         UShort_t response = board->GetResponse();            
-        
-         // 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()));
-
-         }
+       UShort_t response = board->GetResponse();            
         
-         thisl[j-1] = response;
+       // CRATE CONTAINING INTERFACE BOARD
+  if (!board->IsNotified()) // 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] = response;
+    }
     
     regb->SetLocalResponse(thisl);
   }
@@ -521,6 +445,9 @@ void AliMUONTriggerElectronics::LocalResponse()
 void AliMUONTriggerElectronics::RegionalResponse()
 {
   /// Compute the response for all regional cards.
+
+  AliCodeTimerAuto("",0);
+
   AliMUONTriggerCrate* cr;
   TIter next(fCrates->CreateCrateIterator());
   
@@ -529,11 +456,9 @@ void AliMUONTriggerElectronics::RegionalResponse()
       TObjArray *boards = cr->Boards();
 
       AliMUONRegionalTriggerBoard *regb = (AliMUONRegionalTriggerBoard*)boards->At(0);
-      
-      if (regb) 
-      {
-         regb->Response();
-      }  
+
+      regb->Response();
+
    }
 }
 
@@ -542,12 +467,14 @@ void AliMUONTriggerElectronics::GlobalResponse()
 {
   /// Compute the global response
 
+  AliCodeTimerAuto("",0);
+
   UShort_t regional[16];
   
   AliMUONTriggerCrate* cr;
   Int_t irb(0);
   
-  if ( !fCrates->NumberOfCrates() >= 16 ) 
+  if ( fCrates->NumberOfCrates() > 16 ) 
   {
     AliFatal(Form("Something is wrong : too many crates %d",
                   fCrates->NumberOfCrates()));
@@ -577,10 +504,15 @@ void AliMUONTriggerElectronics::GlobalResponse()
 void AliMUONTriggerElectronics::Digits2Trigger(const AliMUONVDigitStore& digitStore,
                                                AliMUONVTriggerStore& triggerStore)
 {
+  AliCodeTimerAuto("",0);
+
   /// Main method to go from digits to trigger decision
   AliMUONRegionalTrigger pRegTrig;
   
   triggerStore.Clear();
+
+  // NOW RESET ELECTRONICS
+  Reset();
   
   // RUN THE FULL BEE CHAIN
   Feed(digitStore);
@@ -623,8 +555,6 @@ void AliMUONTriggerElectronics::Digits2Trigger(const AliMUONVDigitStore& digitSt
           //     {
           
           Int_t icirc = board->GetNumber();
-          if (icirc != 0) { // pcrochet 181206: MOOD needs ALL boards
-            
             localTrigger.SetLoCircuit(icirc);
             localTrigger.SetLoStripX(board->GetStripX11());
             localTrigger.SetLoDev(board->GetDev());
@@ -662,7 +592,6 @@ void AliMUONTriggerElectronics::Digits2Trigger(const AliMUONVDigitStore& digitSt
             triggerStore.Add(localTrigger);  
             
           }
-          }
         }
       pRegTrig.SetId(iReg + 8*iSide);
       pRegTrig.SetLocalOutput(regInpLpt, 0);
@@ -683,7 +612,124 @@ void AliMUONTriggerElectronics::Digits2Trigger(const AliMUONVDigitStore& digitSt
   globalTrigger.SetFromGlobalInput(globalInput);
   // ADD A LOCAL TRIGGER IN THE LIST 
   triggerStore.SetGlobal(globalTrigger);
-  
-  // NOW RESET ELECTRONICS
+
+}
+
+//___________________________________________
+void AliMUONTriggerElectronics::Feed(const AliMUONVTriggerStore& triggerStore)
+{
+  //
+  /// 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();
+}
+
+//_______________________________________________________________________
+Bool_t AliMUONTriggerElectronics::ModifiedLocalResponse(Int_t loCircuit,
+                                                       Bool_t& bendingPlaneResp,
+                                                       Bool_t& nonBendingPlaneResp,
+                                                       Bool_t isCoinc44,
+                                                       Int_t removeChamber)
+{
+  //
+  /// 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()));
+  }
 }