MUON/AliMUONTriggerDecision.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                         *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15 /*
  16 $Log$
  17 Revision 1.5  2000/10/02 16:58:29  egangler
  18 Cleaning of the code :
  19 -> coding conventions
  20 -> void Streamers
  21 -> some useless includes removed or replaced by "class" statement
  22
  23 Revision 1.4  2000/07/03 11:54:57  morsch
  24 AliMUONSegmentation and AliMUONHitMap have been replaced by AliSegmentation and AliHitMap in STEER
  25 The methods GetPadIxy and GetPadXxy of AliMUONSegmentation have changed name to GetPadI and GetPadC.
  26
  27 Revision 1.3  2000/06/25 17:02:19  pcrochet
  28 scope problem on HP, i declared once, pow replaced by TMath::Power (PH)
  29
  30 Revision 1.2  2000/06/15 07:58:49  morsch
  31 Code from MUON-dev joined
  32
  33 Revision 1.1.2.8  2000/06/14 14:54:34  morsch
  34 Complete redesign, make use of TriggerCircuit and TriggerLut (PC)
  35
  36 Revision 1.1.2.5  2000/04/26 19:59:57  morsch
  37 Constructor added.
  38
  39 Revision 1.1.2.4  2000/04/26 12:31:30  morsch
  40 Modifications by P. Crochet:
  41 - adapted to the new Trigger chamber geometry
  42 - condition on soft background added
  43 - contructor added in AliMUONTriggerDecision.h
  44 - single-undefined taken into account in the output of GlobalTrigger()
  45 - some bugs fixed
  46
  47 Revision 1.1.2.3  2000/03/21 09:29:58  morsch
  48 Put back comments
  49
  50 Revision 1.1.2.2  2000/03/21 09:24:34  morsch
  51 Author and responsible for the code: Philippe Crochet
  52 */
  53
  54 #include "AliMUONTriggerCircuit.h"
  55 #include "AliMUONTriggerDecision.h"
  56 #include "AliMUONTriggerLut.h"
  57 #include "AliMUONHitMapA1.h"
  58 #include "AliRun.h"
  59 #include "AliMUON.h"
  60 #include "AliMUONPoints.h"
  61 #include "AliSegmentation.h"
  62 #include "AliMUONResponse.h"
  63 #include "AliMUONChamber.h"
  64 #include "AliMUONDigit.h"
  65
  66
  67 #include <TF1.h>
  68 #include <TTree.h>
  69 #include <TCanvas.h>
  70 #include <TH1.h>
  71 #include <TPad.h>
  72 #include <TGraph.h>
  73 #include <TPostScript.h>
  74 #include <TMinuit.h>
  75 #include <iostream.h>
  76
  77 //----------------------------------------------------------------------
  78 ClassImp(AliMUONTriggerDecision)
  79
  80 //----------------------------------------------------------------------
  81 AliMUONTriggerDecision::AliMUONTriggerDecision(Int_t iprint)
  82 {
  83 // Constructor
  84   fDebug = iprint;            // print option
  85 // iprint = 0 : don't print anything
  86 // iprint = 1 : print Global Trigger Output
  87 // iprint = 2 : print Local and Global Trigger Outputs
  88 // iprint = 3 : iprint = 2 + detailed info on X strips
  89 // iprint = 4 : iprint = 2 + detailed info on Y strip
  90 // iprint = 5 : iprint = 2 + detailed info on X and Y strips
  91 // Note : with iprint>2, the strips detailed info is given for all circuits
  92
  93 // Global Trigger information
  94   Int_t i;
  95   Int_t icirc;
  96   Int_t istrip;
  97
  98   for (i=0; i<3; i++) {   // [0] : Low pt, [1] : High pt, [2] : All pt
  99     fGlobalSinglePlus[i]=0;     // tot num of single plus
 100     fGlobalSingleMinus[i]=0;    // tot num of single minus
 101     fGlobalSingleUndef[i]=0;    // tot num of single undefined
 102     fGlobalPairUnlike[i]=0;     // tot num of unlike-sign pairs
 103     fGlobalPairLike[i]=0;       // tot num of like-sign pairs
 104   }
 105   // Local Trigger information
 106   for (icirc=0; icirc<234; icirc++){
 107     fTrigger[icirc]=0;                   // trigger or not
 108     fStripX11[icirc]=0;                   // X strip in MC11 which triggers
 109     fDev[icirc]=0;                        // deviation which triggers
 110     fStripY11[icirc]=0;                   // Y strip in MC11 which triggers
 111     for (i=0; i<2; i++) {           // pt information via LuT
 112       fLutLpt[icirc][i]=fLutHpt[icirc][i]=fLutApt[icirc][i]=0;
 113     }
 114   }
 115   // bit pattern
 116   for (icirc=0; icirc<234; icirc++) {
 117     for (istrip=0; istrip<16; istrip++) {
 118       fXbit11[icirc][istrip]=fXbit12[icirc][istrip]=0;
 119       fYbit11[icirc][istrip]=fYbit12[icirc][istrip]=0;
 120       fYbit21[icirc][istrip]=fYbit22[icirc][istrip]=0;
 121       fYbit21U[icirc][istrip]=fYbit22U[icirc][istrip]=0;
 122       fYbit21D[icirc][istrip]=fYbit22D[icirc][istrip]=0;
 123     }
 124     for (istrip=0; istrip<32; istrip++) {
 125       fXbit21[icirc][istrip]=fXbit22[icirc][istrip]=0;
 126     }
 127   }
 128 }
 129
 130 //----------------------------------------------------------------------
 131 AliMUONTriggerDecision::~AliMUONTriggerDecision()
 132 {
 133 // Destructor
 134 }
 135
 136 //----------------------------------------------------------------------
 137 void AliMUONTriggerDecision::Trigger(){
 138 // main method of the class which calls the overall Trigger procedure
 139 //  cout << " In AliMUONTriggerDecision::Trigger " << "\n";
 140
 141   ResetBit();
 142   SetBit();
 143   SetBitUpDownY();
 144
 145   Int_t coinc44=0, resetMid=0; // initialize coincidence
 146
 147   AliMUON *pMUON  = (AliMUON*)gAlice->GetModule("MUON");
 148   AliMUONTriggerCircuit* triggerCircuit;
 149
 150   for (Int_t icirc=0; icirc<234; icirc++) {  // loop on circuits
 151     triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 152     //    Int_t idCircuit=triggerCircuit->GetIdCircuit();
 153
 154     Int_t minDevStrip[5], minDev[5], coordY[5];
 155     for (Int_t i=0; i<5; i++) {
 156       minDevStrip[i]=minDev[i]=coordY[i]=0;
 157     }
 158     Int_t x2m=triggerCircuit->GetX2m();
 159     Int_t x2ud=triggerCircuit->GetX2ud();
 160     Int_t orMud[2]={0,0};
 161     triggerCircuit->GetOrMud(orMud);
 162
 163 // call triggerX
 164     TrigX(fXbit11[icirc],fXbit12[icirc],fXbit21[icirc],fXbit22[icirc],
 165           coinc44, minDevStrip, minDev);
 166 // call triggerY
 167     TrigY(fYbit11[icirc],fYbit12[icirc],fYbit21[icirc],fYbit22[icirc],
 168           fYbit21U[icirc],fYbit21D[icirc],fYbit22U[icirc],fYbit22D[icirc],
 169           x2m,x2ud,orMud,resetMid,coinc44,coordY);
 170 // call LocalTrigger
 171     Int_t iTrigger=0;
 172     LocalTrigger(icirc, minDevStrip, minDev, coordY, iTrigger);
 173
 174     if (iTrigger==1&&fDebug>1) {
 175       PrintBitPatXInput(icirc);
 176       PrintBitPatYInput(icirc);
 177       PrintLocalOutput(minDevStrip, minDev, coordY);
 178     }
 179   }  //  end loop on circuits
 180
 181 // call Global Trigger
 182   GlobalTrigger();
 183   //  cout << " Leaving AliMUONTriggerDecision::Trigger " << "\n";
 184 }
 185
 186 //----------------------------------------------------------------------
 187 void AliMUONTriggerDecision::ResetBit(){
 188 // reset bit pattern, global and local trigger output tables to 0
 189
 190     Int_t i;
 191     Int_t icirc;
 192     Int_t istrip;
 193
 194   for (icirc=0; icirc<234; icirc++) {
 195     for (istrip=0; istrip<16; istrip++) {
 196       fXbit11[icirc][istrip]=fXbit12[icirc][istrip]=0;
 197       fYbit11[icirc][istrip]=fYbit12[icirc][istrip]=0;
 198       fYbit21[icirc][istrip]=fYbit22[icirc][istrip]=0;
 199       fYbit21U[icirc][istrip]=fYbit22U[icirc][istrip]=0;
 200       fYbit21D[icirc][istrip]=fYbit22D[icirc][istrip]=0;
 201     }
 202     for (istrip=0; istrip<32; istrip++) {
 203       fXbit21[icirc][istrip]=fXbit22[icirc][istrip]=0;
 204     }
 205   }
 206   for (i=0; i<3; i++) {
 207     fGlobalSinglePlus[i]=0;
 208     fGlobalSingleMinus[i]=0;
 209     fGlobalSingleUndef[i]=0;
 210     fGlobalPairLike[i]=0;
 211     fGlobalPairLike[i]=0;
 212   }
 213   for (icirc=0; icirc<234; icirc++){
 214     fTrigger[icirc]=0;
 215     fStripX11[icirc]=0;
 216     fDev[icirc]=0;
 217     fStripY11[icirc]=0;
 218     for (i=0; i<2; i++) {
 219       fLutLpt[icirc][i]=fLutHpt[icirc][i]=fLutApt[icirc][i]=0;
 220     }
 221   }
 222 }
 223
 224 //----------------------------------------------------------------------
 225 void AliMUONTriggerDecision::SetBit(){
 226 // 1) loop over chambers and cathodes
 227 // 2) load digits
 228 // 3) remove soft background
 229 // 4) set the bit patterns
 230
 231   AliMUON *pMUON  = (AliMUON*)gAlice->GetModule("MUON");
 232   AliMUONTriggerCircuit* triggerCircuit;
 233
 234   for (Int_t chamber=11; chamber<15; chamber++){
 235     for (Int_t cathode=1; cathode<3; cathode++){
 236
 237       AliMUONChamber*   iChamber;
 238       AliSegmentation*  segmentation;
 239
 240       TClonesArray *muonDigits  = pMUON->DigitsAddress(chamber-1);
 241       if (muonDigits == 0) return;
 242
 243       gAlice->ResetDigits();
 244
 245       Int_t nent=(Int_t)gAlice->TreeD()->GetEntries();
 246       gAlice->TreeD()->GetEvent(nent-2+cathode-1);
 247       Int_t ndigits = muonDigits->GetEntriesFast();
 248       if (ndigits == 0) return;
 249
 250       iChamber = &(pMUON->Chamber(chamber-1));
 251       segmentation=iChamber->SegmentationModel(cathode);
 252       AliMUONDigit  *mdig;
 253
 254       for (Int_t digit=0; digit<ndigits; digit++) {
 255         mdig    = (AliMUONDigit*)muonDigits->UncheckedAt(digit);
 256 // get the center of the pad Id
 257         Int_t ix=mdig->PadX();
 258         Int_t iy=mdig->PadY();
 259 // get the sum of the coded charge
 260 // see coding convention in AliMUONChamberTrigger::DisIntegration
 261         Int_t sumCharge=0;
 262         for (Int_t icharge=0; icharge<10; icharge++) {
 263           sumCharge=sumCharge+mdig->TrackCharge(icharge);
 264         }
 265 // apply condition on soft background
 266         Int_t testCharge=sumCharge-(Int_t(sumCharge/10))*10;
 267         testCharge=sumCharge-testCharge*10;
 268         if(sumCharge<=10||testCharge>0) {
 269 // code pad
 270           Int_t code=TMath::Abs(ix)*100+iy;
 271           if (ix<0) { code=-code; }
 272
 273           Int_t icirc;
 274           Int_t istrip;
 275           Int_t nStrip;
 276
 277           if (cathode==1) {
 278             switch (chamber)
 279               {
 280               case 11:
 281                 for (icirc=0; icirc<234; icirc++) {
 282                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 283                   for (istrip=0; istrip<16; istrip++) {
 284                     if (triggerCircuit->GetXcode(0,istrip)==code)
 285                       fXbit11[icirc][istrip]=1;
 286                   }
 287                 }
 288                 break;
 289               case 12:
 290                 for (icirc=0; icirc<234; icirc++) {
 291                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 292                   for (istrip=0; istrip<16; istrip++) {
 293                     if (triggerCircuit->GetXcode(1,istrip)==code)
 294                       fXbit12[icirc][istrip]=1;
 295                   }
 296                 }
 297                 break;
 298               case 13:
 299                 for (icirc=0; icirc<234; icirc++) {
 300                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 301                   for (istrip=0; istrip<32; istrip++) {
 302                     if (triggerCircuit->GetXcode(2,istrip)==code)
 303                       fXbit21[icirc][istrip]=1;
 304                   }
 305                 }
 306                 break;
 307               case 14:
 308                 for (icirc=0; icirc<234; icirc++) {
 309                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 310                   for (istrip=0; istrip<32; istrip++) {
 311                     if (triggerCircuit->GetXcode(3,istrip)==code)
 312                       fXbit22[icirc][istrip]=1;
 313                   }
 314                 }
 315                 break;
 316               }
 317
 318           } else {                // Y plane
 319             switch (chamber)
 320               {
 321               case 11:
 322                 for (icirc=0; icirc<234; icirc++) {
 323                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 324                   nStrip=triggerCircuit->GetNstripY();
 325                   for (istrip=0; istrip<nStrip; istrip++) {
 326                     if (triggerCircuit->GetYcode(0,istrip)==code)
 327                       fYbit11[icirc][istrip]=1;
 328                   }
 329                 }
 330                 break;
 331               case 12:
 332                 for (icirc=0; icirc<234; icirc++) {
 333                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 334                   nStrip=triggerCircuit->GetNstripY();
 335                   for (istrip=0; istrip<nStrip; istrip++) {
 336                     if (triggerCircuit->GetYcode(1,istrip)==code)
 337                       fYbit12[icirc][istrip]=1;
 338                   }
 339                 }
 340                 break;
 341               case 13:
 342                 for (icirc=0; icirc<234; icirc++) {
 343                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 344                   nStrip=triggerCircuit->GetNstripY();
 345                   for (istrip=0; istrip<nStrip; istrip++) {
 346                     if (triggerCircuit->GetYcode(2,istrip)==code)
 347                       fYbit21[icirc][istrip]=1;
 348                   }
 349                 }
 350                 break;
 351               case 14:
 352                 for (icirc=0; icirc<234; icirc++) {
 353                   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
 354                   nStrip=triggerCircuit->GetNstripY();
 355                   for (istrip=0; istrip<nStrip; istrip++) {
 356                     if (triggerCircuit->GetYcode(3,istrip)==code)
 357                       fYbit22[icirc][istrip]=1;
 358                   }
 359                 }
 360                 break;
 361               }
 362           } // if cathode
 363         }  // remove soft background
 364       }   // end loop on digit
 365     }    // end loop on cathode
 366   }     // end loop on chamber
 367 }
 368
 369 //----------------------------------------------------------------------
 370 void AliMUONTriggerDecision::SetBitUpDownY(){
 371 // Set Y bit for up and down parts of circuits
 372   Int_t idModule, nStripX, nStripY, iPosCircuit;
 373
 374   AliMUON *pMUON  = (AliMUON*)gAlice->GetModule("MUON");
 375
 376   for (Int_t icirc=0; icirc<234; icirc++) {
 377
 378     AliMUONTriggerCircuit* circuit;   // current circuit
 379     AliMUONTriggerCircuit* circuitD;  // circuit Down
 380     AliMUONTriggerCircuit* circuitU;  // circuit Up
 381
 382     circuit = &(pMUON->TriggerCircuit(icirc));
 383     idModule=circuit->GetIdModule();      // corresponding module Id.
 384     nStripX=circuit->GetNstripX();        // number of X strips
 385     nStripY=circuit->GetNstripY();        // number of Y strips
 386     iPosCircuit=circuit->GetPosCircuit(); // position of circuit in module
 387
 388 // fill lower part
 389     if (iPosCircuit==1) {               // need to scan lower module
 390       if(idModule<91&&TMath::Abs(idModule)!=41&&idModule>-91) {
 391         Int_t icircD=circuit->GetICircuitD();
 392         circuitD = &(pMUON->TriggerCircuit(icircD));
 393         Int_t nStripD=circuitD->GetNstripY();
 394
 395         if (TMath::Abs(idModule)==42) { // shift of +8 bits
 396           for (Int_t istrip=0; istrip<nStripD; istrip++) {
 397             fYbit21D[icirc][istrip+8]=fYbit21[icircD][istrip];
 398             fYbit22D[icirc][istrip+8]=fYbit22[icircD][istrip];
 399           }
 400         } else if (TMath::Abs(idModule)==52) { // shift of -8 bits
 401           for (Int_t istrip=0; istrip<nStripD; istrip++) {
 402             fYbit21D[icirc][istrip]=fYbit21[icircD][istrip+8];
 403             fYbit22D[icirc][istrip]=fYbit22[icircD][istrip+8];
 404           }
 405         } else {
 406           for (Int_t istrip=0; istrip<nStripD; istrip++) {
 407             fYbit21D[icirc][istrip]=fYbit21[icircD][istrip];
 408             fYbit22D[icirc][istrip]=fYbit22[icircD][istrip];
 409           }
 410         }
 411       }
 412     } else {                         // lower strips within same module
 413       for (Int_t istrip=0; istrip<nStripY; istrip++) {
 414         fYbit21D[icirc][istrip]=fYbit21[icirc][istrip];
 415         fYbit22D[icirc][istrip]=fYbit22[icirc][istrip];
 416       }
 417     }
 418
 419 // fill upper part
 420     if ((iPosCircuit==1&&nStripX==16)||(iPosCircuit==2&&nStripX==32)||
 421         (iPosCircuit==3&&nStripX==48)||(iPosCircuit==4&&nStripX==64)) {
 422       if ((idModule>17||idModule<-17)&&TMath::Abs(idModule)!=61) {
 423         Int_t icircU=circuit->GetICircuitU();
 424         circuitU = &(pMUON->TriggerCircuit(icircU));
 425         Int_t nStripU=circuitU->GetNstripY();
 426
 427         if (TMath::Abs(idModule)==62) { // shift of +8 bits
 428           for (Int_t istrip=0; istrip<nStripU; istrip++) {
 429             fYbit21U[icirc][istrip+8]=fYbit21[icircU][istrip];
 430             fYbit22U[icirc][istrip+8]=fYbit22[icircU][istrip];
 431           }
 432         } else if (TMath::Abs(idModule)==52) { // shift of -8 bits
 433           for (Int_t istrip=0; istrip<nStripU; istrip++) {
 434             fYbit21U[icirc][istrip]=fYbit21[icircU][istrip+8];
 435             fYbit22U[icirc][istrip]=fYbit22[icircU][istrip+8];
 436           }
 437         } else {
 438           for (Int_t istrip=0; istrip<nStripU; istrip++) {
 439             fYbit21U[icirc][istrip]=fYbit21[icircU][istrip];
 440             fYbit22U[icirc][istrip]=fYbit22[icircU][istrip];
 441           }
 442         }
 443       }
 444     } else {                       // upper strips within same module
 445       for (Int_t istrip=0; istrip<nStripY; istrip++) {
 446         fYbit21U[icirc][istrip]=fYbit21[icirc][istrip];
 447         fYbit22U[icirc][istrip]=fYbit22[icirc][istrip];
 448       }
 449     }
 450   } // loop on circuit
 451 }
 452
 453 //----------------------------------------------------------------------
 454 // x part of trigger Algo
 455 //----------------------------------------------------------------------
 456 //----------------------------------------------------------------------
 457 void AliMUONTriggerDecision::TrigX(Int_t ch1q[16], Int_t ch2q[16],
 458                                    Int_t ch3q[32], Int_t ch4q[32],
 459                                    Int_t coinc44, Int_t minDevStrip[5],
 460                                    Int_t minDev[5]){
 461 // note : coinc44 = flag 0 or 1 (0 coincidence -> 3/4, 1 coincidence -> 4/4)
 462 //---------------------------------------------------------
 463 // step # 1 : declustering, reduction DS, calculate sgle & dble
 464 //---------------------------------------------------------
 465   Int_t ch1e[19], ch2e[20], ch3e[35], ch4e[36];
 466   Int_t sgleHit1[31], sgleHit2[63];
 467   Int_t dbleHit1[31], dbleHit2[63];
 468
 469   Int_t i;
 470   Int_t j;
 471   Int_t istrip;
 472
 473   for (i=0; i<31; i++) {
 474     sgleHit1[i]=0;
 475     dbleHit1[i]=0;
 476   }
 477   for (i=0; i<63; i++) {
 478     sgleHit2[i]=0;
 479     dbleHit2[i]=0;
 480   }
 481
 482 //--- inititialize che using chq
 483   for (i=0; i<19; i++) {
 484     if (i<1||i>16)  ch1e[i]=0;
 485     else            ch1e[i]=ch1q[i-1];
 486   }
 487   for (i=0; i<20; i++) {
 488     if (i<2||i>17) ch2e[i]=0;
 489     else           ch2e[i]=ch2q[i-2];
 490   }
 491   for (i=0; i<35; i++) {
 492     if (i<1||i>32) ch3e[i]=0;
 493     else           ch3e[i]=ch3q[i-1];
 494   }
 495   for (i=0; i<36; i++) {
 496     if (i<2||i>33) ch4e[i]=0;
 497     else           ch4e[i]=ch4q[i-2];
 498   }
 499
 500
 501 //--- calculate dble & sgle first station
 502   for (i=0; i<=15; i++) {
 503     sgleHit1[2*i] = (!ch1e[i+1]|(ch1e[i]^ch1e[i+2])) &
 504       (!ch2e[i+2] | (ch2e[i+1]^ch2e[i+3]));
 505
 506     dbleHit1[2*i] = ch1e[i+1]&!(ch1e[i+2]^ch1e[i]) &
 507       (ch2e[i+2] | (!ch2e[i]&ch2e[i+1]) | (ch2e[i+3]&!ch2e[i+4]));
 508   }
 509
 510   for (i=0; i<=14; i++) {
 511     sgleHit1[2*i+1] = (!ch1e[i+1]|!ch1e[i+2]|(ch1e[i]^ch1e[i+3])) &
 512       (!ch2e[i+2] | !ch2e[i+3] | (ch2e[i+1]^ch2e[i+4]));
 513     dbleHit1[2*i+1] = ch1e[i+1]&ch1e[i+2]&!(ch1e[i]^ch1e[i+3]) &
 514       (ch2e[i+2]&(!ch2e[i+1]|!ch2e[i]) |
 515               ch2e[i+3]&(ch2e[i+2]|!ch2e[i+4]|!ch2e[i+5]));
 516   }
 517
 518 //--- calculate dble & sgle second station
 519   for (i=0; i<=31; i++) {
 520     sgleHit2[2*i] = (!ch3e[i+1]|(ch3e[i]^ch3e[i+2])) &
 521       (!ch4e[i+2] | (ch4e[i+1]^ch4e[i+3]));
 522     dbleHit2[2*i] = ch3e[i+1]&!(ch3e[i+2]^ch3e[i]) &
 523       (ch4e[i+2] | (!ch4e[i]&ch4e[i+1]) | (ch4e[i+3]&!ch4e[i+4]));
 524   }
 525
 526   for (i=0; i<=30; i++) {
 527     sgleHit2[2*i+1] = (!ch3e[i+1]|!ch3e[i+2]|(ch3e[i]^ch3e[i+3])) &
 528       (!ch4e[i+2] | !ch4e[i+3] | (ch4e[i+1]^ch4e[i+4]));
 529     dbleHit2[2*i+1] = ch3e[i+1]&ch3e[i+2]&!(ch3e[i]^ch3e[i+3]) &
 530       (ch4e[i+2]&(!ch4e[i+1]|!ch4e[i]) |
 531        ch4e[i+3]&(ch4e[i+2]|!ch4e[i+4]|!ch4e[i+5]));
 532   }
 533
 534 //---
 535   if(fDebug==3||fDebug==5) {
 536     cout << "===============================================================" << "\n";
 537     cout << " X plane after sgle and dble " << " \n";
 538     cout << "                       0987654321098765432109876543210";
 539     cout << "\n SGLE1                 ";
 540     for (istrip=30; istrip>=0; istrip--) { cout << (!sgleHit1[istrip]); }
 541     cout << "\n DBLE1                 ";
 542     for (istrip=30; istrip>=0; istrip--) { cout << dbleHit1[istrip]; }
 543     cout << "\n SGLE2 ";
 544     for (istrip=62; istrip>=0; istrip--) { cout << (!sgleHit2[istrip]); }
 545     cout << "\n DBLE2 ";
 546     for (istrip=62; istrip>=0; istrip--) { cout << dbleHit2[istrip]; }
 547     cout << "\n       210987654321098765432109876543210987654321098765432109876543210" << "\n";
 548   }
 549
 550 //---------------------------------------------------------
 551 // step # 2 : coincidence 3/4
 552 //---------------------------------------------------------
 553   Int_t rearImage[31][31];
 554   for (i=0; i<31; i++) {
 555     for (j=0; j<31; j++) {
 556       rearImage[i][j]=0;
 557     }
 558   }
 559
 560  Int_t notOr1=!dbleHit1[30] & !dbleHit1[29] & !dbleHit1[28] & !dbleHit1[27] &
 561  !dbleHit1[26] & !dbleHit1[25] & !dbleHit1[24] & !dbleHit1[23] &
 562  !dbleHit1[22] & !dbleHit1[21] & !dbleHit1[20] & !dbleHit1[19] &
 563  !dbleHit1[18] & !dbleHit1[17] & !dbleHit1[16] & !dbleHit1[15] &
 564  !dbleHit1[14] & !dbleHit1[13] & !dbleHit1[12] & !dbleHit1[11] &
 565  !dbleHit1[10] & !dbleHit1[9]  & !dbleHit1[8]  & !dbleHit1[7]  &
 566  !dbleHit1[6]  & !dbleHit1[5]  & !dbleHit1[4]  & !dbleHit1[3]  &
 567  !dbleHit1[2]  & !dbleHit1[1]  & !dbleHit1[0]  & !coinc44;
 568
 569  Int_t notOr2= !dbleHit2[62] & !dbleHit2[61] & !dbleHit2[60] & !dbleHit2[59] &
 570  !dbleHit2[58] & !dbleHit2[57] & !dbleHit2[56] & !dbleHit2[55] &
 571  !dbleHit2[54] & !dbleHit2[53] & !dbleHit2[52] & !dbleHit2[51] &
 572  !dbleHit2[50] & !dbleHit2[49] & !dbleHit2[48] & !dbleHit2[47] &
 573  !dbleHit2[46] & !dbleHit2[45] & !dbleHit2[44] & !dbleHit2[43] &
 574  !dbleHit2[42] & !dbleHit2[41] & !dbleHit2[40] & !dbleHit2[39] &
 575  !dbleHit2[38] & !dbleHit2[37] & !dbleHit2[36] & !dbleHit2[35] &
 576  !dbleHit2[34] & !dbleHit2[33] & !dbleHit2[32] & !dbleHit2[31] &
 577  !dbleHit2[30] & !dbleHit2[29] & !dbleHit2[28] & !dbleHit2[27] &
 578  !dbleHit2[26] & !dbleHit2[25] & !dbleHit2[24] & !dbleHit2[23] &
 579  !dbleHit2[22] & !dbleHit2[21] & !dbleHit2[20] & !dbleHit2[19] &
 580  !dbleHit2[18] & !dbleHit2[17] & !dbleHit2[16] & !dbleHit2[15] &
 581  !dbleHit2[14] & !dbleHit2[13] & !dbleHit2[12] & !dbleHit2[11] &
 582  !dbleHit2[10] & !dbleHit2[9]  & !dbleHit2[8]  & !dbleHit2[7]  &
 583  !dbleHit2[6]  & !dbleHit2[5]  & !dbleHit2[4]  & !dbleHit2[3]  &
 584  !dbleHit2[2]  & !dbleHit2[1]  & !dbleHit2[0]  & !coinc44;
 585
 586 // DS reduction
 587  for (i=0; i<31; i++) {
 588    sgleHit1[i] = !sgleHit1[i]&notOr1;
 589  }
 590  for (i=0; i<63; i++) {
 591    sgleHit2[i] = !sgleHit2[i]&notOr2;
 592  }
 593
 594 // extract rearImage
 595  for (i=0; i<31; i++){
 596    Int_t tmpSgleHit2[31];
 597    Int_t tmpDbleHit2[31];
 598    for (j=0; j<31; j++){
 599      tmpSgleHit2[j] = sgleHit2[i+j+1];
 600      tmpDbleHit2[j] = dbleHit2[i+j+1];
 601    }
 602
 603    for (Int_t k=0; k<31; k++) {
 604      rearImage[i][k]=(sgleHit1[i]&tmpDbleHit2[k])|
 605        (dbleHit1[i]&(tmpSgleHit2[k]|tmpDbleHit2[k]));
 606    }
 607  }
 608
 609   //-----------
 610  if(fDebug==3||fDebug==5) {
 611     cout << "===============================================================" << "\n";
 612    for (i=30; i>=0; i--) {
 613      cout << i << "\t ";
 614      for (istrip=31; istrip>=0; istrip--) {
 615        cout << rearImage[i][istrip];
 616      }
 617      cout << " " << "\n";
 618    }
 619  }
 620
 621
 622 //---------------------------------------------------------
 623 // step # 3 : calculate deviation
 624 //---------------------------------------------------------
 625  Int_t dev[31][6];
 626  for (i=0; i<31; i++) {
 627    for (j=0; j<6; j++) {
 628      dev[i][j]=0;
 629    }
 630  }
 631
 632  for (i=0; i<31; i++){
 633    Int_t leftDev[5], rightDev[5];
 634    Int_t orL1, andL1, andL2, orR1, orR2, andR1, andR2, andR3;
 635
 636 // calculate Left deviation
 637  orL1=rearImage[i][16]|rearImage[i][18]|rearImage[i][20]|rearImage[i][22];
 638  andL1=!rearImage[i][17]&!rearImage[i][19]&!rearImage[i][21] & !orL1;
 639  andL2=!rearImage[i][23]&!rearImage[i][24]&!rearImage[i][25]&!rearImage[i][26];
 640
 641  leftDev[0] = (rearImage[i][16]|!rearImage[i][17]) &
 642  (rearImage[i][16]|rearImage[i][18]|!rearImage[i][19]&
 643  (rearImage[i][20]|!rearImage[i][21])) &
 644  (orL1|!rearImage[i][23]&(rearImage[i][24]|!rearImage[i][25])) &
 645  (orL1|rearImage[i][24]|rearImage[i][26]|!rearImage[i][27]&
 646  (rearImage[i][28]|!rearImage[i][29]));
 647
 648  leftDev[1] = !rearImage[i][16] &
 649  !(!rearImage[i][17]&!rearImage[i][18]&!rearImage[i][21]&!rearImage[i][22] &
 650  (!rearImage[i][25]&!rearImage[i][26]&(rearImage[i][27]|rearImage[i][28]))) &
 651  (rearImage[i][17]|rearImage[i][18] | !rearImage[i][19]&!rearImage[i][20]) &
 652  (rearImage[i][17]|rearImage[i][18]|rearImage[i][21]|rearImage[i][22] |
 653  !rearImage[i][23]&!rearImage[i][24]);
 654
 655  leftDev[2] = (!rearImage[i][16]&!rearImage[i][17]&!rearImage[i][18]) &
 656  (rearImage[i][19]|rearImage[i][20]|rearImage[i][21]|rearImage[i][22] | andL2);
 657
 658  leftDev[3] = andL1;
 659
 660  leftDev[4] =
 661  !rearImage[i][27]&!rearImage[i][28]&!rearImage[i][29]&!rearImage[i][30] &
 662  andL1 & andL2;
 663
 664  // calculate Right deviation
 665  orR1=rearImage[i][8]|rearImage[i][10]|rearImage[i][12]|rearImage[i][14];
 666  orR2=rearImage[i][8]|rearImage[i][9]|rearImage[i][10]|rearImage[i][11];
 667  andR1=!rearImage[i][12]&!rearImage[i][13]&!rearImage[i][14]&!rearImage[i][15];
 668  andR2=
 669  !rearImage[i][8]&!rearImage[i][9]&!rearImage[i][10]&!rearImage[i][11] & andR1;
 670  andR3=!rearImage[i][4]&!rearImage[i][5]&!rearImage[i][6]&!rearImage[i][7];
 671
 672  rightDev[0] = !rearImage[i][15]&(rearImage[i][14]|!rearImage[i][13]) &
 673  ((rearImage[i][12]|rearImage[i][14]|!rearImage[i][11]&
 674  (rearImage[i][10]|!rearImage[i][9])) &
 675  ((orR1|!rearImage[i][7]&(rearImage[i][6]|!rearImage[i][5])) &
 676  (orR1|rearImage[i][4]|rearImage[i][6]|!rearImage[i][3]&(rearImage[i][2]|
 677  !rearImage[i][1]))));
 678
 679  rightDev[1] = !rearImage[i][15]&!rearImage[i][14] &
 680  !(!rearImage[i][4]&!rearImage[i][5]&!rearImage[i][8]&!rearImage[i][9] &
 681  (!rearImage[i][12]&!rearImage[i][13]&(rearImage[i][2]|rearImage[i][3]))) &
 682  (rearImage[i][12]|rearImage[i][13] | !rearImage[i][10]&!rearImage[i][11]) &
 683  (rearImage[i][8]|rearImage[i][9]|rearImage[i][12]|rearImage[i][13] |
 684  !rearImage[i][6]&!rearImage[i][7]);
 685
 686  rightDev[2] = andR1 & (orR2 | andR3);
 687  rightDev[3] = andR2;
 688  rightDev[4] =
 689  !rearImage[i][0]&!rearImage[i][1]&!rearImage[i][2]&!rearImage[i][3] &
 690  andR2 & andR3 ;
 691
 692  // compare Left & Right deviations
 693  Int_t tmpLeftDev=0, tmpRightDev=0;
 694  for (j=0; j<5; j++){
 695    tmpLeftDev  = tmpLeftDev + Int_t(leftDev[j]*TMath::Power(2,j));
 696    tmpRightDev = tmpRightDev + Int_t(rightDev[j]*TMath::Power(2,j));
 697  }
 698
 699  // assign mimimum deviation do dev[][]
 700  if (tmpLeftDev < tmpRightDev ){
 701    for (j=0; j<5; j++){ dev[i][j]=leftDev[j];}
 702    dev[i][5]=1;
 703  } else {
 704    for (j=0; j<5; j++){ dev[i][j]=rightDev[j];}
 705    dev[i][5]=0;
 706  }
 707   }
 708
 709 //---
 710   if(fDebug==3||fDebug==5) {
 711     cout << "===============================================================" << "\n";
 712     for (i=30; i>=0; i--) {
 713       cout << i << "\t ";
 714       for (istrip=5; istrip>=0; istrip--) { cout << dev[i][istrip]; }
 715       cout << " " << "\n";
 716     }
 717   }
 718
 719 //---------------------------------------------------------
 720 // step # 4 : sort deviation
 721 //---------------------------------------------------------
 722  Int_t bga1[16], bga2[8], bga3[4], bga4[2], bga5;
 723  Int_t tmpbga1[16][6], tmpbga2[8][6], tmpbga3[4][6], tmpbga4[2][6], tmpbga5[6];
 724  Int_t tmpMax[6]={1,1,1,1,1,0};
 725
 726   for (i=0; i<15; i++) {
 727     Sort2x5(dev[2*i],dev[2*i+1],tmpbga1[i],bga1[i]);
 728   }
 729     Sort2x5(dev[30],tmpMax,tmpbga1[15],bga1[15]);
 730
 731 //--
 732   if(fDebug==3||fDebug==5) {
 733     cout << "===============================================================" << "\n";
 734     cout << " sorting : 1st level " << "\n";
 735     for (i=15; i>=0; i--) {
 736       cout << i << "\t " << bga1[i] << "\t";
 737       for (j=5; j>=0; j--) {
 738         cout << tmpbga1[i][j];
 739       }
 740       cout << " " << "\n";
 741     }
 742   }
 743
 744   for (i=0; i<8; i++) {
 745     Sort2x5(tmpbga1[2*i],tmpbga1[2*i+1],tmpbga2[i],bga2[i]);
 746   }
 747
 748 //--
 749   if(fDebug==3||fDebug==5) {
 750     cout << "===============================================================" << "\n";
 751     cout << " sorting : 2nd level " << "\n";
 752     for (i=7; i>=0; i--) {
 753       cout << i << "\t " << bga2[i] << "\t";
 754       for (j=5; j>=0; j--) {
 755         cout << tmpbga2[i][j];
 756       }
 757       cout << " " << "\n";
 758     }
 759   }
 760
 761   for (i=0; i<4; i++) {
 762     Sort2x5(tmpbga2[2*i],tmpbga2[2*i+1],tmpbga3[i],bga3[i]);
 763   }
 764
 765 //--
 766   if(fDebug==3||fDebug==5) {
 767     cout << "===============================================================" << "\n";
 768     cout << " sorting : 3rd level " << "\n";
 769     for (i=3; i>=0; i--) {
 770       cout << i << "\t " << bga3[i] << "\t";
 771       for (j=5; j>=0; j--) {
 772         cout << tmpbga3[i][j];
 773       }
 774       cout << " " << "\n";
 775     }
 776   }
 777
 778   for (i=0; i<2; i++) {
 779     Sort2x5(tmpbga3[2*i],tmpbga3[2*i+1],tmpbga4[i],bga4[i]);
 780   }
 781
 782 //--
 783   if(fDebug==3||fDebug==5) {
 784     cout << "===============================================================" << "\n";
 785     cout << " sorting : 4th level " << "\n";
 786     for (i=1; i>=0; i--) {
 787       cout << i << "\t " << bga4[i] << "\t";
 788       for (j=5; j>=0; j--) {
 789         cout << tmpbga4[i][j];
 790       }
 791       cout << " " << "\n";
 792     }
 793   }
 794
 795     Sort2x5(tmpbga4[0],tmpbga4[1],tmpbga5,bga5);
 796
 797  // coding from 6 to 5 bits
 798     minDev[4] = tmpbga5[5] | tmpbga5[4];
 799     for (i=0; i<4; i++) {
 800       minDev[i]=tmpbga5[i] & !tmpbga5[4];
 801     }
 802
 803  // find address of strip with minimum deviation
 804     minDevStrip[4]=bga5;
 805     if (bga5<=1) minDevStrip[3]=bga4[bga5];
 806
 807     Int_t tmpAd=minDevStrip[3]+minDevStrip[4]*2;
 808     if (tmpAd<=3) minDevStrip[2]=bga3[tmpAd];
 809
 810     tmpAd=minDevStrip[2]+minDevStrip[3]*2+minDevStrip[4]*4;
 811     if (tmpAd<=7) minDevStrip[1]=bga2[tmpAd];
 812
 813     tmpAd=minDevStrip[1]+minDevStrip[2]*2+minDevStrip[3]*4+minDevStrip[4]*8;
 814     if (tmpAd<=15) minDevStrip[0]=bga1[tmpAd];
 815
 816   if(fDebug==3||fDebug==5) {
 817     cout << "===============================================================" << "\n";
 818     cout << "minDevStrip = ";
 819     for  (i=4; i>=0; i--) {cout << minDevStrip[i];}
 820     cout << " minDev = ";
 821     for  (i=4; i>=0; i--) {cout << minDev[i];}
 822     cout << " " << "\n";
 823     cout << "===============================================================" << "\n";
 824   }
 825
 826 }
 827
 828 //---------------------------------------------
 829 void AliMUONTriggerDecision::Sort2x5(Int_t dev1[6], Int_t dev2[6],
 830                                      Int_t minDev[6], Int_t &dev1GTdev2){
 831 // returns minimun between dev1 and dev2
 832  Int_t tmpDev1=0, tmpDev2=0;
 833  for (Int_t j=0; j<5; j++){
 834    tmpDev1 = tmpDev1 + Int_t(dev1[j]*TMath::Power(2,j));
 835    tmpDev2 = tmpDev2 + Int_t(dev2[j]*TMath::Power(2,j));
 836  }
 837  if (tmpDev1 <= tmpDev2 ){
 838    for (Int_t j=0; j<=5; j++) { minDev[j]=dev1[j];}
 839    dev1GTdev2=0;
 840  } else {
 841    for (Int_t j=0; j<=5; j++) { minDev[j]=dev2[j];}
 842    dev1GTdev2=1;
 843  }
 844 }
 845
 846 //----------------------------------------------------------------------
 847 // y part of trigger Algo
 848 //----------------------------------------------------------------------
 849 //----------------------------------------------------------------------
 850 void AliMUONTriggerDecision::TrigY(Int_t y1[16], Int_t y2[16],
 851                                    Int_t y3[16], Int_t y4[16],
 852                                    Int_t y3u[16], Int_t y3d[16],
 853                                    Int_t y4u[16], Int_t y4d[16],
 854                                    Int_t x2m, Int_t x2ud, Int_t orMud[2],
 855                                    Int_t resetMid, Int_t coinc44,
 856                                    Int_t coordY[5]){
 857 // note : resMid = 1 -> cancel
 858 //---------------------------------------------------------
 859 // step # 1 : prehandling Y
 860 //---------------------------------------------------------
 861     Int_t i;
 862     Int_t istrip;
 863
 864   for (i=0; i<16; i++){
 865     y3[i]=y3[i]&!resetMid;
 866     y4[i]=y4[i]&!resetMid;
 867   }
 868
 869   Int_t ch1[16], ch2[16], ch3[16], ch4[16];
 870
 871   Int_t tmpy3to16[16], tmpy4to16[16];
 872   Int_t tmpy3uto16[16], tmpy3dto16[16], tmpy4uto16[16], tmpy4dto16[16];
 873   for (i=0; i<8; i++){
 874     ch1[2*i]   = y1[i]&x2m | y1[2*i]&!x2m;
 875     ch1[2*i+1] = y1[i]&x2m | y1[2*i+1]&!x2m;
 876
 877     ch2[2*i]   = y2[i]&x2m | y2[2*i]&!x2m;
 878     ch2[2*i+1] = y2[i]&x2m | y2[2*i+1]&!x2m;
 879
 880     tmpy3to16[2*i]   = y3[i]&x2m | y3[2*i]&!x2m;
 881     tmpy3to16[2*i+1] = y3[i]&x2m | y3[2*i+1]&!x2m;
 882
 883     tmpy4to16[2*i]   = y4[i]&x2m | y4[2*i]&!x2m;
 884     tmpy4to16[2*i+1] = y4[i]&x2m | y4[2*i+1]&!x2m;
 885
 886     tmpy3uto16[2*i]   = y3u[i]&x2ud | y3u[2*i]&!x2ud;
 887     tmpy3uto16[2*i+1] = y3u[i]&x2ud | y3u[2*i+1]&!x2ud;
 888
 889     tmpy4uto16[2*i]   = y4u[i]&x2ud | y4u[2*i]&!x2ud;
 890     tmpy4uto16[2*i+1] = y4u[i]&x2ud | y4u[2*i+1]&!x2ud;
 891
 892     tmpy3dto16[2*i]   = y3d[i]&x2ud | y3d[2*i]&!x2ud;
 893     tmpy3dto16[2*i+1] = y3d[i]&x2ud | y3d[2*i+1]&!x2ud;
 894
 895     tmpy4dto16[2*i]   = y4d[i]&x2ud | y4d[2*i]&!x2ud;
 896     tmpy4dto16[2*i+1] = y4d[i]&x2ud | y4d[2*i+1]&!x2ud;
 897   }
 898
 899   if (orMud[0]==0&&orMud[1]==0){
 900     for (i=0; i<16; i++){
 901       ch3[i] = tmpy3to16[i];
 902       ch4[i] = tmpy4to16[i];
 903     }
 904   }
 905   if (orMud[0]==0&&orMud[1]==1){
 906       for (i=0; i<16; i++){
 907         ch3[i] = tmpy3uto16[i]|tmpy3to16[i];
 908         ch4[i] = tmpy4uto16[i]|tmpy4to16[i];
 909       }
 910   }
 911   if (orMud[0]==1&&orMud[1]==0){
 912       for (i=0; i<16; i++){
 913         ch3[i] = tmpy3dto16[i]|tmpy3to16[i];
 914         ch4[i] = tmpy4dto16[i]|tmpy4to16[i];
 915       }
 916   }
 917   if (orMud[0]==1&&orMud[1]==1){
 918       for (i=0; i<16; i++){
 919         ch3[i] = tmpy3dto16[i]|tmpy3to16[i]|tmpy3uto16[i];
 920         ch4[i] = tmpy4dto16[i]|tmpy4to16[i]|tmpy4uto16[i];
 921       }
 922   }
 923
 924 // debug
 925   if(fDebug==4||fDebug==5) {
 926     cout << "===============================================================" << "\n";
 927     cout << " Y plane after PreHandling x2m x2ud orMud "
 928          << x2m << " , " << x2ud << " , " << orMud[0] << orMud[1] << "\n";
 929     cout << "                            ";
 930     for (istrip=15; istrip>=0; istrip--) {
 931       if (istrip>9) cout << istrip-10*Int_t(istrip/10);
 932       if (istrip<10) cout << istrip;
 933     }
 934     cout << "\n YMC11                      ";
 935     for (istrip=15; istrip>=0; istrip--) {
 936       cout << ch1[istrip];
 937     }
 938     cout << "\n YMC12                      ";
 939     for (istrip=15; istrip>=0; istrip--) {
 940       cout << ch2[istrip];
 941     }
 942     cout << "\n YMC21                      ";
 943     for (istrip=15; istrip>=0; istrip--) {
 944       cout << ch3[istrip];
 945     }
 946     cout << "\n YMC22                      ";
 947     for (istrip=15; istrip>=0; istrip--) {
 948       cout << ch4[istrip];
 949     }
 950     cout << " \n";
 951   }
 952 //debug
 953
 954 //---------------------------------------------------------
 955 // step # 2 : calculate sgle and dble, apply DS reduction
 956 //---------------------------------------------------------
 957   Int_t sgle1[16], dble1[16];
 958   Int_t sgle2[16], dble2[16];
 959
 960   // Calculate simple and double hits
 961   for (i=0; i<16; i++) {
 962     dble1[i] = ch1[i] & ch2[i];
 963     dble2[i] = ch3[i] & ch4[i];
 964
 965     sgle1[i] = (ch1[i]|ch2[i]);
 966     sgle2[i] = (ch3[i]|ch4[i]);
 967   }
 968
 969   //debug
 970   if(fDebug==4||fDebug==5) {
 971     cout << "===============================================================" << "\n";
 972     cout << " Y plane after sgle dble " << "\n";
 973     cout << "                            ";
 974     for (istrip=15; istrip>=0; istrip--) {
 975       if (istrip>9) { cout << istrip-10*Int_t(istrip/10);}
 976       if (istrip<10) { cout << istrip;}
 977     }
 978     cout << "\n SGLE1                      ";
 979     for (istrip=15; istrip>=0; istrip--) {
 980       cout << sgle1[istrip];
 981     }
 982     cout << "\n DBLE1                      ";
 983     for (istrip=15; istrip>=0; istrip--) {
 984       cout << dble1[istrip];
 985     }
 986     cout << "\n SGLE2                      ";
 987     for (istrip=15; istrip>=0; istrip--) {
 988       cout << sgle2[istrip];
 989     }
 990     cout << "\n DBLE2                      ";
 991     for (istrip=15; istrip>=0; istrip--) {
 992       cout << dble2[istrip];
 993     }
 994     cout << " \n";
 995   }
 996   //debug
 997
 998   // DS Reduction
 999   Int_t notOr1, notOr2;
1000
1001   notOr1=!dble1[15] & !dble1[14] & !dble1[13] & !dble1[12] &
1002          !dble1[11] & !dble1[10] & !dble1[9]  & !dble1[8]  &
1003          !dble1[7]  & !dble1[6]  & !dble1[5]  & !dble1[4]  &
1004          !dble1[3]  & !dble1[2]  & !dble1[1]  & !dble1[0];
1005
1006   notOr2=!dble2[15] & !dble2[14] & !dble2[13] & !dble2[12] &
1007          !dble2[11] & !dble2[10] & !dble2[9]  & !dble2[8]  &
1008          !dble2[7]  & !dble2[6]  & !dble2[5]  & !dble2[4]  &
1009          !dble2[3]  & !dble2[2]  & !dble2[1]  & !dble2[0];
1010
1011   for (i=0; i<16; i++) {
1012     sgle1[i] = sgle1[i] & notOr1 & !coinc44;
1013     sgle2[i] = sgle2[i] & notOr2 & !coinc44;
1014   }
1015
1016 //---------------------------------------------------------
1017 // step # 3 : 3/4 coincidence
1018 //---------------------------------------------------------
1019   Int_t frontImage[16];
1020
1021   for (i=1; i<15; i++) {
1022   frontImage[i] = (dble1[i] | sgle1[i]) &
1023     (dble2[i+1] | dble2[i] | dble2[i-1]) |
1024      dble1[i] & (sgle2[i+1] | sgle2[i] | sgle2[i-1]);
1025   }
1026   frontImage[0] = (dble1[0] | sgle1[0]) &
1027     (dble2[1] | dble2[0]) | dble1[0] & (sgle2[1] | sgle2[0]);
1028
1029   frontImage[15] = (dble1[15] | sgle1[15]) &
1030     (dble2[15] | dble2[14]) | dble1[15] & (sgle2[15] | sgle2[14]);
1031
1032
1033 //debug
1034   if(fDebug==4||fDebug==5) {
1035     cout << "===============================================================" << "\n";
1036     cout << " Y plane frontImage\n";
1037     cout << "                            ";
1038   for (istrip=15; istrip>=0; istrip--) {
1039     if (istrip>9) cout << istrip-10*Int_t(istrip/10);
1040     if (istrip<10) cout << istrip;
1041   }
1042     cout << "\n                            ";
1043     for (istrip=15; istrip>=0; istrip--) {
1044       cout << frontImage[istrip];
1045     }
1046     cout << "\n";
1047   }
1048 //debug
1049
1050 //---------------------------------------------------------
1051 // step # 4 : Y position
1052 //---------------------------------------------------------
1053   Int_t or1, or2, and1, and2, and3;
1054
1055  or1  = frontImage[7]|frontImage[5]|frontImage[3]|frontImage[1];
1056  or2  = frontImage[7]|frontImage[6]|frontImage[5]|frontImage[4];
1057  and1 = !frontImage[3]&!frontImage[2]&!frontImage[1]&!frontImage[0];
1058  and2 = !frontImage[7]&!frontImage[6]&!frontImage[5]&!frontImage[4] & and1;
1059  and3 = !frontImage[11]&!frontImage[10]&!frontImage[9]&!frontImage[8];
1060
1061  coordY[0] = !frontImage[0]&(frontImage[1]|!frontImage[2]) &
1062 (frontImage[3]|frontImage[1]|!frontImage[4]&(frontImage[5]|!frontImage[6])) &
1063 (or1|!frontImage[8]&(frontImage[9]|!frontImage[10])) &
1064 (or1|frontImage[11]|frontImage[9]|!frontImage[12]&(frontImage[13]|!frontImage[14]));
1065
1066  coordY[1] = !frontImage[0]&!frontImage[1] &
1067 !(!frontImage[11]&!frontImage[10]&!frontImage[7]&!frontImage[6] &
1068   !frontImage[3]&!frontImage[2]&(frontImage[13]|frontImage[12])) &
1069   (frontImage[3]|frontImage[2] | !frontImage[5]&!frontImage[4]) &
1070   (frontImage[7]|frontImage[6]|frontImage[3]|frontImage[2] |
1071 !frontImage[9]&!frontImage[8]);
1072
1073  coordY[2] = and1 & (or2 | and3);
1074
1075  coordY[3] = and2;
1076
1077  coordY[4] = !frontImage[15]&!frontImage[14]&!frontImage[13]&!frontImage[12] &
1078  and2 & and3 ;
1079
1080 }
1081 //----------------------------------------------------------------------
1082 // end of trigger Algo
1083 //----------------------------------------------------------------------
1084
1085 //----------------------------------------------------------------------
1086 void AliMUONTriggerDecision::LocalTrigger(Int_t icirc,
1087                                           Int_t minDevStrip[5],
1088                                           Int_t minDev[5], Int_t coordY[5],
1089                                           Int_t &iTrigger){
1090 // returns local trigger answer for circuit icirc
1091   Int_t i;
1092
1093   AliMUON *pMUON  = (AliMUON*)gAlice->GetModule("MUON");
1094   AliMUONTriggerCircuit* triggerCircuit;
1095   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
1096   Int_t idCircuit=triggerCircuit->GetIdCircuit();
1097
1098   Int_t signDev=minDev[4];
1099   Int_t deviation=0;
1100   for (i=0; i<4; i++) {          // extract deviation
1101     deviation = deviation+Int_t(minDev[i]*TMath::Power(2,i));
1102   }
1103
1104   Int_t istripX1Circ=0;
1105   for (i=0; i<5; i++) {          // extract X1 strip fired
1106     istripX1Circ = istripX1Circ+Int_t(minDevStrip[i]*TMath::Power(2,i));
1107   }
1108
1109   Int_t iStripY=0;
1110   for (i=0; i<4; i++) {          // extract Y strip fired
1111     iStripY = iStripY+Int_t(coordY[i]*TMath::Power(2,i));
1112   }
1113
1114 // trigger or not
1115   if (signDev==1&&deviation==0) {      // something in X ?
1116     iTrigger=0;
1117   } else {
1118     if (coordY[4]==1&&iStripY==15) {   // something in Y ?
1119       iTrigger=0;
1120     } else {
1121       iTrigger=1;
1122     }
1123   }
1124
1125   if (iTrigger==1) {
1126 // fill fTrigger fStripX11 fStripY11
1127     fTrigger[icirc] = 1;
1128     fStripX11[icirc] = istripX1Circ;
1129     fStripY11[icirc] = iStripY;
1130
1131 // calculate deviation in [0+30]
1132     Int_t sign=0;
1133     if (signDev==0&&deviation!=0) sign=-1;
1134     if (signDev==0&&deviation==0) sign=0;
1135     if (signDev==1)               sign=1;
1136     fDev[icirc] = sign * deviation + 15; // fill fDev
1137
1138 // get Lut output for circuit/istripX/idev/istripY
1139     AliMUONTriggerLut* lut = new AliMUONTriggerLut;
1140     //    lut->StartEvent();
1141     lut->GetLutOutput(icirc,fStripX11[icirc],fDev[icirc],fStripY11[icirc],
1142                       fLutLpt[icirc],fLutHpt[icirc],fLutApt[icirc]);
1143     //    lut->FinishEvent();
1144     delete lut;
1145
1146     if (fDebug>1) {
1147       Float_t pt= // get ptCal corresponding to istripX1Circ/idev/iStripY
1148       triggerCircuit->PtCal(fStripX11[icirc],fDev[icirc],fStripY11[icirc]);
1149       cout << "-------------------------------------------" << "\n";
1150       cout << " Local Trigger info for circuit Id " << idCircuit
1151            << " (number " << icirc << ")" << "\n";
1152       cout << " istripX1 signDev deviation istripY = "
1153            << istripX1Circ << " , " << signDev
1154            << " , " << deviation << " , " << iStripY << "\n";
1155       cout << " pt = " << pt << " (GeV/c) " << "\n";
1156       cout << "-------------------------------------------" << "\n";
1157       cout << " Local Trigger Lut Output = Lpt : " ;
1158       for (i=1; i>=0; i--) { cout << fLutLpt[icirc][i] ; }
1159       cout << " Hpt : ";
1160       for (i=1; i>=0; i--) { cout << fLutHpt[icirc][i] ; }
1161       cout << " Apt : ";
1162       for (i=1; i>=0; i--) { cout << fLutApt[icirc][i] ; }
1163       cout << "\n";
1164       cout << "-------------------------------------------" << "\n";
1165     } // fDebug > 1
1166   }  // local trigger = 1
1167 }
1168
1169 //----------------------------------------------------------------------
1170 void AliMUONTriggerDecision::GlobalTrigger(){
1171 // loop on Lut[icirc] and give Global Trigger output
1172     Int_t i;
1173
1174   for (Int_t icirc=0; icirc<234; icirc++){
1175     if (fLutLpt[icirc][0]==1&&fLutLpt[icirc][1]==1)
1176       fGlobalSingleUndef[0] = fGlobalSingleUndef[0] + 1;
1177     if (fLutHpt[icirc][0]==1&&fLutHpt[icirc][1]==1)
1178       fGlobalSingleUndef[1] = fGlobalSingleUndef[1] + 1;
1179     if (fLutApt[icirc][0]==1&&fLutApt[icirc][1]==1)
1180       fGlobalSingleUndef[2] = fGlobalSingleUndef[2] + 1;
1181
1182     if (fLutLpt[icirc][0]==0&&fLutLpt[icirc][1]==1)
1183       fGlobalSinglePlus[0] = fGlobalSinglePlus[0] + 1;
1184     if (fLutHpt[icirc][0]==0&&fLutHpt[icirc][1]==1)
1185       fGlobalSinglePlus[1] = fGlobalSinglePlus[1] + 1;
1186     if (fLutApt[icirc][0]==0&&fLutApt[icirc][1]==1)
1187       fGlobalSinglePlus[2] = fGlobalSinglePlus[2] + 1;
1188
1189     if (fLutLpt[icirc][0]==1&&fLutLpt[icirc][1]==0)
1190       fGlobalSingleMinus[0] = fGlobalSingleMinus[0] + 1;
1191     if (fLutHpt[icirc][0]==1&&fLutHpt[icirc][1]==0)
1192       fGlobalSingleMinus[1] = fGlobalSingleMinus[1] + 1;
1193     if (fLutApt[icirc][0]==1&&fLutApt[icirc][1]==0)
1194       fGlobalSingleMinus[2] = fGlobalSingleMinus[2] + 1;
1195   }
1196
1197   // like sign low, high and all pt
1198   for (i=0; i<3; i++) {
1199     fGlobalPairLike[i]=fGlobalSingleMinus[i]*(fGlobalSingleMinus[i]-1)/2 +
1200       fGlobalSinglePlus[i]*(fGlobalSinglePlus[i]-1)/2 +
1201       fGlobalSingleUndef[i]*(fGlobalSingleUndef[i]-1)/2 +
1202       fGlobalSingleUndef[i]*fGlobalSinglePlus[i] +
1203       fGlobalSingleUndef[i]*fGlobalSingleMinus[i];
1204   }
1205
1206   // unlike sign low, high and all pt
1207   for (i=0; i<3; i++) {
1208     fGlobalPairUnlike[i]=fGlobalSingleMinus[i]*fGlobalSinglePlus[i] +
1209       fGlobalSingleUndef[i]*(fGlobalSingleUndef[i]-1)/2 +
1210       fGlobalSingleUndef[i]*fGlobalSinglePlus[i] +
1211       fGlobalSingleUndef[i]*fGlobalSingleMinus[i];
1212   }
1213
1214   if (fDebug>=1) {
1215     cout << "\n";
1216     cout << "===================================================" << "\n";
1217     cout << " Global Trigger output       " << "Low pt  High pt   All"  << "\n";
1218     cout << " number of Single Plus      :\t";
1219     for (i=0; i<3; i++) { cout << fGlobalSinglePlus[i] <<"\t";}
1220     cout << "\n";
1221     cout << " number of Single Minus     :\t";
1222     for (i=0; i<3; i++) { cout << fGlobalSingleMinus[i] <<"\t";}
1223     cout << "\n";
1224     cout << " number of Single Undefined :\t";
1225     for (i=0; i<3; i++) { cout << fGlobalSingleUndef[i] <<"\t";}
1226     cout << "\n";
1227     cout << " number of UnlikeSign pair  :\t";
1228     for (i=0; i<3; i++) { cout << fGlobalPairUnlike[i] <<"\t";}
1229     cout << "\n";
1230     cout << " number of LikeSign pair    :\t";
1231     for (i=0; i<3; i++) { cout << fGlobalPairLike[i] <<"\t";}
1232     cout << "\n";
1233     cout << "===================================================" << "\n";
1234   }
1235 }
1236
1237 //----------------------------------------------------------------------
1238 void AliMUONTriggerDecision::PrintBitPatXInput(Int_t icirc){
1239 // print bit pattern for X strips
1240
1241     Int_t istrip;
1242
1243   cout << "-------- TRIGGER INPUT ---------" << "\n";
1244   cout << "===============================================================" << "\n";
1245   cout << "                            5432109876543210";
1246   cout << "\n XMC11                      ";
1247   for (istrip=15; istrip>=0; istrip--) {
1248     cout << fXbit11[icirc][istrip];
1249   }
1250   cout << "\n XMC12                      ";
1251   for (istrip=15; istrip>=0; istrip--) {
1252     cout << fXbit12[icirc][istrip];
1253   }
1254   cout << "\n XMC21              ";
1255   for (istrip=31; istrip>=0; istrip--) {
1256     cout << fXbit21[icirc][istrip];
1257   }
1258   cout << "\n XMC22              ";
1259   for (istrip=31; istrip>=0; istrip--) {
1260     cout << fXbit22[icirc][istrip];
1261   }
1262   cout << "\n                    ";
1263   cout << "10987654321098765432109876543210" << "\n";
1264 }
1265
1266 //----------------------------------------------------------------------
1267 void AliMUONTriggerDecision::PrintBitPatYInput(Int_t icirc){
1268 // print bit pattern for Y strips
1269
1270     Int_t istrip;
1271
1272   AliMUON *pMUON  = (AliMUON*)gAlice->GetModule("MUON");
1273   AliMUONTriggerCircuit* triggerCircuit;
1274   triggerCircuit = &(pMUON->TriggerCircuit(icirc));
1275   Int_t idCircuit=triggerCircuit->GetIdCircuit();
1276   Int_t nStrip=triggerCircuit->GetNstripY();
1277
1278   cout << "---------------------------------------------------------------" << "\n";
1279   cout << "                            ";
1280   for (istrip=nStrip-1; istrip>=0; istrip--) {
1281     if (istrip>9) { cout << istrip-10*Int_t(istrip/10);}
1282     if (istrip<10) { cout << istrip;}
1283   }
1284   cout << "\n YMC11                      ";
1285   for (istrip=nStrip-1; istrip>=0; istrip--) {
1286     cout << fYbit11[icirc][istrip];
1287   }
1288   cout << "\n YMC12                      ";
1289   for (istrip=nStrip-1; istrip>=0; istrip--) {
1290     cout << fYbit12[icirc][istrip];
1291   }
1292   cout << "\n YMC21                      ";
1293   for (istrip=nStrip-1; istrip>=0; istrip--) {
1294     cout << fYbit21[icirc][istrip];
1295   }
1296   cout << "\n YMC22                      ";
1297   for (istrip=nStrip-1; istrip>=0; istrip--) {
1298     cout << fYbit22[icirc][istrip];
1299   }
1300   cout << "\n";
1301 // tmp
1302   cout << "---------------------------------------------------------------";
1303   cout << "\n upper part of circuit " << idCircuit ;
1304   cout << "\n UMC21                      ";
1305   for (istrip=15; istrip>=0; istrip--) {
1306     cout << fYbit21U[icirc][istrip];
1307   }
1308   cout << "\n UMC22                      ";
1309   for (istrip=15; istrip>=0; istrip--) {
1310     cout << fYbit22U[icirc][istrip];
1311   }
1312
1313   cout << "\n lower part of circuit " << idCircuit ;
1314   cout << "\n LMC21                      ";
1315   for (istrip=15; istrip>=0; istrip--) {
1316     cout << fYbit21D[icirc][istrip];
1317   }
1318   cout << "\n LMC22                      ";
1319   for (istrip=15; istrip>=0; istrip--) {
1320     cout << fYbit22D[icirc][istrip];
1321   }
1322   cout << "\n";
1323   cout << "===============================================================" << "\n";
1324 }
1325 //----------------------------------------------------------------------
1326 void AliMUONTriggerDecision::PrintLocalOutput(Int_t minDevStrip[5],
1327                                               Int_t minDev[5],
1328                                               Int_t coordY[5]){
1329 // print Local trigger output before the LuT step
1330
1331     Int_t i;
1332
1333   cout << "===============================================================" << "\n";
1334   cout << "-------- TRIGGER OUTPUT --------" << "\n";
1335   cout << "minDevStrip = ";
1336   for  (i=4; i>=0; i--) {cout << minDevStrip[i];}
1337   cout << " minDev = ";
1338   for  (i=4; i>=0; i--) {cout << minDev[i];}
1339   cout << " coordY = ";
1340   for  (i=4; i>=0; i--) {cout << coordY[i];}
1341   cout << " " << "\n";
1342 }
1343
1344 //----------------------------------------------------------------------
1345 //--- methods which return member data related info
1346 //----------------------------------------------------------------------
1347 Int_t AliMUONTriggerDecision::GetITrigger(Int_t icirc){
1348 // returns Local Trigger Status
1349   return fTrigger[icirc];
1350 }
1351 //----------------------------------------------------------------------
1352 Int_t AliMUONTriggerDecision::GetStripX11(Int_t icirc){
1353 // returns fStripX11
1354   return fStripX11[icirc];
1355 }
1356 //----------------------------------------------------------------------
1357 Int_t AliMUONTriggerDecision::GetDev(Int_t icirc){
1358 // returns idev
1359   return fDev[icirc];
1360 }
1361 //----------------------------------------------------------------------
1362 Int_t AliMUONTriggerDecision::GetStripY11(Int_t icirc){
1363 // returns fStripY11;
1364    return fStripY11[icirc];
1365 }
1366 //----------------------------------------------------------------------
1367 void AliMUONTriggerDecision::GetLutOutput(Int_t icirc, Int_t lpt[2],
1368                                           Int_t hpt[2], Int_t apt[2]){
1369 // returns Look up Table output
1370   for (Int_t i=0; i<2; i++) {
1371     lpt[i]=fLutLpt[icirc][i];
1372     hpt[i]=fLutHpt[icirc][i];
1373     apt[i]=fLutApt[icirc][i];
1374   }
1375 }
1376 //----------------------------------------------------------------------
1377 void AliMUONTriggerDecision::GetGlobalTrigger(Int_t singlePlus[3],
1378                                               Int_t singleMinus[3],
1379                                               Int_t singleUndef[3],
1380                                               Int_t pairUnlike[3],
1381                                               Int_t pairLike[3]){
1382 // returns Global Trigger information (0,1,2 : Lpt,Hpt,Apt)
1383   for (Int_t i=0; i<3; i++) {
1384     singlePlus[i]  = fGlobalSinglePlus[i];
1385     singleMinus[i] = fGlobalSingleMinus[i];
1386     singleUndef[i] = fGlobalSingleUndef[i];
1387     pairUnlike[i]  = fGlobalPairUnlike[i];
1388     pairLike[i]    = fGlobalPairLike[i];
1389   }
1390 }
1391 //----------------------------------------------------------------------
1392 //--- end of methods which return member data related info
1393 //----------------------------------------------------------------------
1394 //----------------------------------------------------------------------
1395 /*
1396 void AliMUONTriggerDecision::AddLocalTrigger(const AliMUONLocalTrigger c){
1397 // Add a Local Trigger copy to the list
1398   AliMUON *MUON=(AliMUON*)gAlice->GetModule("MUON");
1399   MUON->AddLocalTrigger(c);
1400   fNLocalTriggers++;
1401 }
1402 */