From 9a7352cc50d3fc91bacaeb97a8397be1c18ad25e Mon Sep 17 00:00:00 2001
From: cvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 12 Oct 2011 12:58:53 +0000
Subject: [PATCH] Adding methods to calculate the number of trigger flags on
 each side and fill the V0 trigger bits which should be contained in the
 raw-data header.

---
 VZERO/AliVZERORawStream.cxx | 88 +++++++++++++++++++++++++++++++++++++
 VZERO/AliVZERORawStream.h   |  4 ++
 2 files changed, 92 insertions(+)

diff --git a/VZERO/AliVZERORawStream.cxx b/VZERO/AliVZERORawStream.cxx
index c146c1cb7c4..7d0ed3c5108 100644
--- a/VZERO/AliVZERORawStream.cxx
+++ b/VZERO/AliVZERORawStream.cxx
@@ -263,3 +263,91 @@ void AliVZERORawStream::CalculateChargeForCentrTriggers(AliVZEROTriggerData *tri
     }
   }
 }
+
+//_____________________________________________________________________________
+void AliVZERORawStream::CalculateBBandBGFlags(AliVZEROTriggerData *triggerData,
+					      UChar_t &nBBA, UChar_t &nBBC,
+					      UChar_t &nBGA, UChar_t &nBGC) const
+{
+  // Use the raw-data payload
+  // in order to calculate the total
+  // number of beam-beam and beam-gas flags
+  // (which is used in the centrality and
+  // multiplicity triggers) on each side of V0
+  if (!triggerData) {
+    AliFatal("Trigger configuration data is not provided. Exiting...");
+    return;
+  }
+  nBBA = nBBC = nBGA = nBGC = 0;
+  for(Int_t iChannel=0; iChannel<64; iChannel++) {
+    Int_t offlineCh = GetOfflineChannel(iChannel);
+    Int_t board = AliVZEROCalibData::GetBoardNumber(offlineCh);
+    Int_t feeChannel = AliVZEROCalibData::GetFEEChannelNumber(offlineCh);
+    if (triggerData->GetEnableTiming(board,feeChannel)) {
+      if (offlineCh < 32) {
+	if (GetBBFlag(iChannel,10)) nBBC++;
+	if (GetBGFlag(iChannel,10)) nBGC++;
+      }
+      else {
+	if (GetBBFlag(iChannel,10)) nBBA++;
+	if (GetBGFlag(iChannel,10)) nBGA++;
+      }
+    }
+  }
+}
+
+//_____________________________________________________________________________
+void AliVZERORawStream::FillTriggerBits(AliVZEROTriggerData *triggerData)
+{
+  // Calculate the charge sums and
+  // number of trigger flags and then
+  // fill the V0 trigger bits word
+  // following the trigger logic implemented
+  // in the firmware
+  UShort_t chargeA,chargeC;
+  CalculateChargeForCentrTriggers(triggerData,chargeA,chargeC);
+  UChar_t nBBA,nBBC,nBGA,nBGC;
+  CalculateBBandBGFlags(triggerData,
+			nBBA,nBBC,
+			nBGA,nBGC);
+
+  fTrigger = 0;
+  //BBA and BBC
+  if((nBBC >= triggerData->GetBBCThreshold()) && (nBBA >= triggerData->GetBBAThreshold())) fTrigger |= 1;
+  //BBA or BBC
+  if((nBBC >= triggerData->GetBBCThreshold()) || (nBBA >= triggerData->GetBBAThreshold())) fTrigger |= (1<<1);
+  //BGA and BBC
+  if((nBBC >= triggerData->GetBBCForBGThreshold()) && (nBGA >= triggerData->GetBGAThreshold())) fTrigger |= (1<<2);
+  //BGA 
+  if((nBGA >= triggerData->GetBGAThreshold())) fTrigger |= (1<<3);
+  //BGC and BBA
+  if((nBGC >= triggerData->GetBGCThreshold()) && (nBBA >= triggerData->GetBBAForBGThreshold())) fTrigger |= (1<<4);
+  //BGC 
+  if((nBGC >= triggerData->GetBGCThreshold())) fTrigger |= (1<<5);
+  //CTA1 and CTC1
+  if((chargeC >= triggerData->GetCentralityV0CThrLow()) && (chargeA >= triggerData->GetCentralityV0AThrLow())) fTrigger |= (1<<6);
+  //CTA1 or CTC1
+  if((chargeC >= triggerData->GetCentralityV0CThrLow()) || (chargeA >= triggerData->GetCentralityV0AThrLow())) fTrigger |= (1<<7);
+  //CTA2 and CTC2
+  if((chargeC >= triggerData->GetCentralityV0CThrHigh()) && (chargeA >= triggerData->GetCentralityV0AThrHigh())) fTrigger |= (1<<8);
+  //CTA2 or CTC2
+  if((chargeC >= triggerData->GetCentralityV0CThrHigh()) || (chargeA >= triggerData->GetCentralityV0AThrHigh())) fTrigger |= (1<<9);
+  //MTA and MTC
+  if(((nBBC >= triggerData->GetMultV0CThrLow()) && (nBBC <= triggerData->GetMultV0CThrHigh())) &&
+     ((nBBA >= triggerData->GetMultV0AThrLow()) && (nBBA <= triggerData->GetMultV0AThrHigh()))) 
+    fTrigger |= (1<<10);
+  //MTA or MTC
+  if(((nBBC >= triggerData->GetMultV0CThrLow()) && (nBBC <= triggerData->GetMultV0CThrHigh())) ||
+     ((nBBA >= triggerData->GetMultV0AThrLow()) && (nBBA <= triggerData->GetMultV0AThrHigh()))) 
+    fTrigger |= (1<<11);
+  //BBA 
+  if((nBBA >= triggerData->GetBBAThreshold())) fTrigger |= (1<<12);
+  //BBC
+  if((nBBC >= triggerData->GetBBCThreshold())) fTrigger |= (1<<13);
+  //BGA or BGC 
+  if((nBGC >= triggerData->GetBGCThreshold()) || (nBGA >= triggerData->GetBGAThreshold())) fTrigger |= (1<<14);
+  //(BGA and BBC) or (BGC and BBA) 
+  if(((nBBC >= triggerData->GetBBCForBGThreshold()) && (nBGA >= triggerData->GetBGAThreshold())) ||
+     ((nBGC >= triggerData->GetBGCThreshold()) && (nBBA >= triggerData->GetBBAForBGThreshold()))) fTrigger |= (1<<15);
+
+}
diff --git a/VZERO/AliVZERORawStream.h b/VZERO/AliVZERORawStream.h
index 2cd6c537f4b..e8135e90efe 100644
--- a/VZERO/AliVZERORawStream.h
+++ b/VZERO/AliVZERORawStream.h
@@ -76,6 +76,10 @@ class AliVZERORawStream: public TObject {
 
     void              CalculateChargeForCentrTriggers(AliVZEROTriggerData *triggerData,
 						      UShort_t &chargeA, UShort_t &chargeC) const;
+    void              CalculateBBandBGFlags(AliVZEROTriggerData *triggerData,
+					    UChar_t &nBBA, UChar_t &nBBC,
+					    UChar_t &nBGA, UChar_t &nBGC) const;
+    void              FillTriggerBits(AliVZEROTriggerData *triggerData);
 
 // Getter of Offline Channel number as used in aliroot (defined by aliroot 
 // numbering convention) from FEE channel (electronic channel number given 
-- 
2.43.0