X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FAliMUONTriggerCircuit.cxx;h=016f2067f4389a3d08581d7e4ef9b68991660124;hb=b8f92f9dec22bca57816e5765e3216b63199ae46;hp=0dc2a678b6d93d968a664cd49e994dd0b06458c0;hpb=d9a3473dbf35b79b7c40a1cd6c43db0a09803d9a;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/AliMUONTriggerCircuit.cxx b/MUON/AliMUONTriggerCircuit.cxx index 0dc2a678b6d..016f2067f43 100644 --- a/MUON/AliMUONTriggerCircuit.cxx +++ b/MUON/AliMUONTriggerCircuit.cxx @@ -1,494 +1,560 @@ /************************************************************************** - * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * - * * - * Author: The ALICE Off-line Project. * - * Contributors are mentioned in the code where appropriate. * - * * - * Permission to use, copy, modify and distribute this software and its * - * documentation strictly for non-commercial purposes is hereby granted * - * without fee, provided that the above copyright notice appears in all * - * copies and that both the copyright notice and this permission notice * - * appear in the supporting documentation. The authors make no claims * - * about the suitability of this software for any purpose. It is * - * provided "as is" without express or implied warranty. * - **************************************************************************/ +* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * +* * +* Author: The ALICE Off-line Project. * +* Contributors are mentioned in the code where appropriate. * +* * +* Permission to use, copy, modify and distribute this software and its * +* documentation strictly for non-commercial purposes is hereby granted * +* without fee, provided that the above copyright notice appears in all * +* copies and that both the copyright notice and this permission notice * +* appear in the supporting documentation. The authors make no claims * +* about the suitability of this software for any purpose. It is * +* provided "as is" without express or implied warranty. * +**************************************************************************/ /* $Id$ */ -#include "AliRun.h" -#include "AliMUON.h" +//----------------------------------------------------------------------------- +/// \class AliMUONTriggerCircuit +/// Contains as data members the Y positions of the X declusturized strips and +/// the X positions of the (doubled or not) Y strips. +/// This is used to associate the global positions to the fired strips of the +/// local trigger output (see AliMUONTrackReconstructor::MakeTriggerTrack) +/// +/// \author Philippe Crochet (LPCCFd) +//----------------------------------------------------------------------------- + #include "AliMUONTriggerCircuit.h" -#include "AliMUONTriggerConstants.h" -#include "AliSegmentation.h" -#include "AliMUONResponse.h" -#include "AliMUONChamber.h" -#include "TMath.h" -#include "Riostream.h" +#include "AliMUONConstants.h" +#include "AliMUONGeometryTransformer.h" + +#include "AliMpTrigger.h" +#include "AliMpSlat.h" +#include "AliMpPCB.h" +#include "AliMpSegmentation.h" +#include "AliMpVSegmentation.h" +#include "AliMpCathodType.h" +#include "AliMpDDLStore.h" +#include "AliMpLocalBoard.h" +#include "AliMpConstants.h" +#include "AliMpPad.h" +#include "AliMpEncodePair.h" + +#include "AliRun.h" +#include "AliLog.h" +#include +#include + +using std::cout; +using std::endl; +/// \cond CLASSIMP ClassImp(AliMUONTriggerCircuit) +/// \endcond //---------------------------------------------------------------------- -AliMUONTriggerCircuit::AliMUONTriggerCircuit() +AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONGeometryTransformer* transformer) + : TObject(), + fkTransformer(transformer), + fkCurrentSeg(0x0), + fCurrentDetElem(0x0), + fCurrentLocalBoard(0x0) { -// Constructor - fSegmentation=0; - fIdCircuit=0; - fX2m=0; - fX2ud=0; - fOrMud[0]=fOrMud[1]=0; - Int_t i; - for (i=0; i<4; i++) { - for (Int_t j=0; j<32; j++) { - fXcode[i][j]=0; - fYcode[i][j]=0; +/// Constructor + + for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) { + fXpos11[i].Set(16); + fYpos11[i].Set(32); + fYpos21[i].Set(64); + fZpos11[i].Set(32); + fZpos21[i].Set(64); + fXwidth11[i].Set(16); + fYwidth11[i].Set(32); + fYwidth21[i].Set(64); + } + + for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) { // board begins at 1 + + AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(i); + + if (!localBoard) + { + AliError(Form("Did not get localboard %d",i)); + continue; } + + LoadXPos(localBoard); + LoadYPos(localBoard); + /* + printf("LocalBoard %03d \n",i); + printf("fXpos11 \n"); + for (Int_t i1 = 0; i1 < 16; i1++) + printf("%02d %7.2f \n",i1,fXpos11[i][i1]); + printf("fYpos11 \n"); + for (Int_t i2 = 0; i2 < 32; i2++) + printf("%02d %7.2f \n",i2,fYpos11[i][i2]); + printf("fYpos21 \n"); + for (Int_t i3 = 0; i3 < 64; i3++) + printf("%02d %7.2f \n",i3,fYpos21[i][i3]); + printf("fZpos11 \n"); + for (Int_t i4 = 0; i4 < 32; i4++) + printf("%02d %8.2f \n",i4,fZpos11[i][i4]); + printf("fZpos21 \n"); + for (Int_t i5 = 0; i5 < 64; i5++) + printf("%02d %8.2f \n",i5,fZpos21[i][i5]); + */ } - for (i=0; i<16; i++) { fXpos11[i]=0.; } - for (i=0; i<31; i++) { fYpos11[i]=0.; } - for (i=0; i<63; i++) { fYpos21[i]=0.; } -} -//---------------------------------------------------------------------- -AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONTriggerCircuit& MUONTriggerCircuit):TObject(MUONTriggerCircuit) -{ -// Dummy copy constructor } //---------------------------------------------------------------------- -AliMUONTriggerCircuit & AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& /*MUONTriggerCircuit*/) +AliMUONTriggerCircuit::~AliMUONTriggerCircuit() { -// Dummy assignment operator - return *this; -} +/// Destructor + for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) { + fXpos11[i].Reset(); + fYpos11[i].Reset(); + fYpos21[i].Reset(); + fZpos11[i].Reset(); + fZpos21[i].Reset(); + fXwidth11[i].Reset(); + fYwidth11[i].Reset(); + fYwidth21[i].Reset(); + } -//---------------------------------------------------------------------- -void AliMUONTriggerCircuit::Init(Int_t iCircuit) { -// initialize circuit characteristics - fIdCircuit=AliMUONTriggerConstants::CircuitId(iCircuit); - LoadX2(); - LoadXCode(); - LoadYCode(); - LoadXPos(); - LoadYPos(); -} +} //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::CircuitNumber(Int_t idCircuit){ -// returns circuit number iCircuit (0-234) corresponding to circuit idCircuit - Int_t iCircuit=0; - for (Int_t i=0; i<234; i++) { - if (AliMUONTriggerConstants::CircuitId(i)==idCircuit) { - iCircuit=i; - break; +AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONTriggerCircuit& circuit) + : TObject(circuit), + fkTransformer(circuit.fkTransformer), // do not copy, just pointed to + fkCurrentSeg(circuit.fkCurrentSeg), + fCurrentDetElem(circuit.fCurrentDetElem), + fCurrentLocalBoard(circuit.fCurrentLocalBoard) +{ +/// Copy constructor + + for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) { + fXpos11[i] = circuit.fXpos11[i]; + fYpos11[i] = circuit.fYpos11[i]; + fYpos21[i] = circuit.fYpos21[i]; + fZpos11[i] = circuit.fZpos11[i]; + fZpos21[i] = circuit.fZpos21[i]; + fXwidth11[i] = circuit.fXwidth11[i]; + fYwidth11[i] = circuit.fYwidth11[i]; + fYwidth21[i] = circuit.fYwidth21[i]; } - } - return iCircuit; + } //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::ModuleNumber(Int_t idModule){ -// returns module number imod (from 0 to 63) corresponding to module idmodule - Int_t absidModule=TMath::Abs(idModule); - Int_t iModule=0; - for (Int_t i=0; i<63; i++) { - if (AliMUONTriggerConstants::ModuleId(i)==absidModule) { - iModule=i; - break; - } +AliMUONTriggerCircuit& AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& circuit) +{ +/// Assignment operator + + if (this == &circuit) return *this; + + fkTransformer = circuit.fkTransformer; + fkCurrentSeg = circuit.fkCurrentSeg; + fCurrentDetElem = circuit.fCurrentDetElem; + fCurrentLocalBoard = circuit.fCurrentLocalBoard; + + for (Int_t i = 1; i < AliMpConstants::NofLocalBoards()+1; ++i) { + fXpos11[i] = circuit.fXpos11[i]; + fYpos11[i] = circuit.fYpos11[i]; + fYpos21[i] = circuit.fYpos21[i]; + fZpos11[i] = circuit.fZpos11[i]; + fZpos21[i] = circuit.fZpos21[i]; + fXwidth11[i] = circuit.fXwidth11[i]; + fYwidth11[i] = circuit.fYwidth11[i]; + fYwidth21[i] = circuit.fYwidth21[i]; } - return iModule; -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::Module(Int_t idCircuit) { -// returns ModuleId where Circuit idCircuit is sitting - return Int_t(idCircuit/10); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::Position(Int_t idCircuit) { -// returns position of idCircuit in correcponding Module - return TMath::Abs(idCircuit)-TMath::Abs(Module(idCircuit))*10; + return *this; + } -//---------------------------------------------------------------------- -void AliMUONTriggerCircuit::LoadX2() { -// initialize fX2m, fX2ud and fOrMud - - Int_t idModule=Module(fIdCircuit); // corresponding module Id. -// and its number of X strips - Int_t nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule)); -// and its number of Y strips - Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule)); - Int_t iPosCircuit=Position(fIdCircuit); // position of circuit in module - -// first step : look at lower part - if (iPosCircuit==1) { // need to scan lower module - if(idModule<91&&TMath::Abs(idModule)!=41&&idModule>-91) { - fOrMud[0]=1; - Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule); - Int_t nStrD=AliMUONTriggerConstants::NstripY(ModuleNumber(idModuleD)); - - if (nStrY!=nStrD - &&TMath::Abs(idModule)!=42&&TMath::Abs(idModule)!=52) { - if (nStrY==8) fX2m=1; - if (nStrD==8) fX2ud=1; - } - } - - } else { // lower strips within same module - fOrMud[0]=0; - } +//--------------------------------------------------------------------- +void AliMUONTriggerCircuit::LoadYPos(AliMpLocalBoard* const localBoard) +{ +/// fill fYpos11 and fYpos21 -> y position of X declusterized strips -// second step : look at upper part - if ((iPosCircuit==1&&nStrX==16)||(iPosCircuit==2&&nStrX==32)|| - (iPosCircuit==3&&nStrX==48)||(iPosCircuit==4&&nStrX==64)) { - if ((idModule>17||idModule<-17)&&TMath::Abs(idModule)!=61) { - fOrMud[1]=1; - Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule); - Int_t nStrU=AliMUONTriggerConstants::NstripY(ModuleNumber(idModuleU)); - - if (nStrY!=nStrU - &&TMath::Abs(idModule)!=62&&TMath::Abs(idModule)!=52) { - if (nStrY==8) fX2m=1; - if (nStrU==8) fX2ud=1; - } - } + fCurrentLocalBoard = localBoard->GetId(); + Int_t ichamber = 0; + Int_t icathode = 0; - } else { // upper strips within same module - fOrMud[1]=0; - } -} + Int_t zeroDown = localBoard->GetSwitch(AliMpLocalBoard::kZeroDown); + Int_t zeroUp = localBoard->GetSwitch(AliMpLocalBoard::kZeroUp); + + Int_t iline = AliMp::PairFirst(localBoard->GetPosition()); + Int_t icol = AliMp::PairSecond(localBoard->GetPosition()); + if ( iline == 5 ) --icol; + + //--- first plane + ichamber = 10; + fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber); + + fkCurrentSeg = AliMpSegmentation::Instance() + ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode)); + + Int_t iFirstStrip = FirstStrip(localBoard); + Int_t iLastStrip = iFirstStrip + 16; + Int_t iStripCircuit = 0; + + FillXstrips(icol, iFirstStrip, iLastStrip, + iStripCircuit, kTRUE); + + //--- second plane + ichamber = 12; + fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber); -//---------------------------------------------------------------------- -void AliMUONTriggerCircuit::LoadXCode(){ -// assign a Id. number to each X strip of current circuit -// Id.=(corresponding module Id.)*100+(Id. strip of module) - -// first part : fill XMC11 XMC12 and strips 8 to 24 (middle) XMC21 XMC22 - Int_t iStripCircMT1=0, iStripCircMT2=8; - Int_t idModule=Module(fIdCircuit); // corresponding module Id. -// and its number of strips - Int_t nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule)); - Int_t iPosCircuit=Position(fIdCircuit); // position of circuit in module - Int_t sign=TMath::Abs(idModule)/idModule; // left or right - Int_t istrip; - - for (istrip=(iPosCircuit-1)*16; - istrip<(iPosCircuit-1)*16+16; istrip++) { - - fXcode[0][iStripCircMT1]=sign*(TMath::Abs(idModule)*100+istrip); - fXcode[1][iStripCircMT1]=sign*(TMath::Abs(idModule)*100+istrip); - fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip); - fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip); - iStripCircMT1++; - iStripCircMT2++; - } + fkCurrentSeg = AliMpSegmentation::Instance() + ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode)); -// second part -// XMC21 XMC22 strips 0 to 7 and 24 to 31 - Int_t idModuleD, idModuleU; - Int_t nStrD, nStrU; - - idModule=Module(fIdCircuit); // corresponding module Id. -// number of X strips - nStrX=AliMUONTriggerConstants::NstripX(ModuleNumber(idModule)); - sign=TMath::Abs(idModule)/idModule; - -// fill lower part (0 to 7) - if (iPosCircuit==1) { // need to scan lower module - if(idModule<91&&TMath::Abs(idModule)!=41&&idModule>-91) { // non-existing - idModuleD=sign*(TMath::Abs(idModule)+10); // lower module Id -// and its number of strips - nStrD=AliMUONTriggerConstants::NstripX(ModuleNumber(idModuleD)); - - iStripCircMT2=0; - for (istrip=nStrD-8; istripPadByIndices(icol-1,iLastStripMiddle+1,kFALSE); + + if (pad.IsValid()) { // upper strips within same detElemId + iFirstStripUp = iLastStripMiddle; + iLastStripUp = iFirstStripUp + 8; + + } else { // upper strips in another detElemId + fCurrentDetElem = AliMpDDLStore::Instance()-> + GetNextDEfromLocalBoard(fCurrentLocalBoard, ichamber); + + fkCurrentSeg = AliMpSegmentation::Instance() + ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode)); + + iFirstStripUp = 0; + iLastStripUp = iFirstStripUp + 8; + if (iline == 4) icolUp = icol - 1; // special case + if (iline == 5) icolUp = icol + 1; // special case } - - } else { // lower strips within same module - iStripCircMT2=0; - for (istrip=(iPosCircuit-1)*16-8; - istrip<(iPosCircuit-1)*16; istrip++) { - fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip); - fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip); - iStripCircMT2++; - } - } + iStripCircuit = 24; + FillXstrips(icolUp, iFirstStripUp, iLastStripUp, + iStripCircuit, kFALSE); + + // fill strip between middle and upper part + fYpos21[fCurrentLocalBoard][47] = (fYpos21[fCurrentLocalBoard][46] + + fYpos21[fCurrentLocalBoard][48])/2.; + fZpos21[fCurrentLocalBoard][47] = (fZpos21[fCurrentLocalBoard][46] + + fZpos21[fCurrentLocalBoard][48])/2.; + fYwidth21[fCurrentLocalBoard][47] = (fYwidth21[fCurrentLocalBoard][46] + + fYwidth21[fCurrentLocalBoard][48])/2.; + } // end of something up -// fill upper part (24 to 31) - if ((iPosCircuit==1&&nStrX==16)||(iPosCircuit==2&&nStrX==32)|| - (iPosCircuit==3&&nStrX==48)||(iPosCircuit==4&&nStrX==64)) { - if ((idModule>17||idModule<-17)&&TMath::Abs(idModule)!=61) { - idModuleU=sign*(TMath::Abs(idModule)-10); // upper module Id -// and its number of strips - nStrU=AliMUONTriggerConstants::NstripX(ModuleNumber(idModuleU)); - - iStripCircMT2=24; - for (istrip=0; istrip<8; istrip++) { - fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModuleU)*100+istrip); - fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModuleU)*100+istrip); - iStripCircMT2++; - } - } + // restore current detElemId & segmentation + fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber); + fkCurrentSeg = AliMpSegmentation::Instance() + ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode)); + + // second plane lower part + if (zeroDown == 0) { // something down + Int_t iFirstStripDo; + Int_t iLastStripDo; + Int_t icolDo = icol; - } else if ((iPosCircuit==1&&nStrX>16)||(iPosCircuit==2&&nStrX>32)|| - (iPosCircuit==3&&nStrX>48)) { // upper strips within same mod. + // check if we need to move to another detElemId + AliMpPad pad = fkCurrentSeg->PadByIndices(icol-1,iFirstStripMiddle-1,kFALSE); + if (pad.IsValid()) { // lower strips within same detElemId + iFirstStripDo = iFirstStripMiddle - 8; + iLastStripDo = iFirstStripDo + 8; + + } else { // lower strips in another detElemId + fCurrentDetElem = AliMpDDLStore::Instance() + ->GetPreviousDEfromLocalBoard(fCurrentLocalBoard, ichamber); + + fkCurrentSeg = AliMpSegmentation::Instance() + ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode)); + + // get iFirstStrip in this module + const AliMpTrigger* t = AliMpSegmentation::Instance()->GetTrigger(fkCurrentSeg); + const AliMpSlat* slat = t->GetLayer(0); + + if (iline == 5) icolDo = icol + 1; // special case + if (iline == 6) icolDo = icol - 1; // special case + + const AliMpPCB* pcb = slat->GetPCB(icolDo-1); + iFirstStripDo = (pcb->Iymax() + 1) - 8; + iLastStripDo = iFirstStripDo + 8; + } - iStripCircMT2=24; - for (istrip=(iPosCircuit-1)*16+16; - istrip<(iPosCircuit-1)*16+24; istrip++) { - fXcode[2][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip); - fXcode[3][iStripCircMT2]=sign*(TMath::Abs(idModule)*100+istrip); - iStripCircMT2++; - } - } + iStripCircuit = 0; + FillXstrips(icolDo, iFirstStripDo, iLastStripDo, + iStripCircuit, kFALSE); + + // fill strip between middle and upper part + fYpos21[fCurrentLocalBoard][15] = (fYpos21[fCurrentLocalBoard][14] + + fYpos21[fCurrentLocalBoard][16])/2.; + fZpos21[fCurrentLocalBoard][15] = (fZpos21[fCurrentLocalBoard][14] + + fZpos21[fCurrentLocalBoard][16])/2.; + fYwidth21[fCurrentLocalBoard][15] = (fYwidth21[fCurrentLocalBoard][14] + + fYwidth21[fCurrentLocalBoard][16])/2.; + } // end of something down + } //---------------------------------------------------------------------- -void AliMUONTriggerCircuit::LoadYCode(){ -// assign a Id. number to each Y strip of current circuit -// Id.=(corresponding module Id.)*100+(Id. strip of module) -// note : for Y plane fill only "central part" of circuit -// (upper and lower parts are filled in PreHandlingY of AliMUONTriggerDecision) +void AliMUONTriggerCircuit::FillXstrips(const Int_t icol, + const Int_t iFirstStrip, const Int_t iLastStrip, + Int_t liStripCircuit, const Bool_t is11) +{ +/// fill + TArrayF& ypos = (is11) ? fYpos11[fCurrentLocalBoard] : fYpos21[fCurrentLocalBoard]; + TArrayF& zpos = (is11) ? fZpos11[fCurrentLocalBoard] : fZpos21[fCurrentLocalBoard]; + TArrayF& ywidth = (is11) ? fYwidth11[fCurrentLocalBoard] : fYwidth21[fCurrentLocalBoard]; + + Double_t xyGlobal[3] = {0.}; + for (Int_t istrip = iFirstStrip; istrip < iLastStrip; ++istrip) { + + AliMpPad pad = fkCurrentSeg->PadByIndices(icol-1,istrip,kTRUE); + if ( !pad.IsValid() ) { + StdoutToAliError(cout << "Pad not found in seg " << endl; + fkCurrentSeg->Print(); + cout << " ix,iy=" << icol-1 << "," << istrip << endl; + ); + } + Float_t yDim = pad.GetDimensionY(); // half size! + + XYGlobal(pad,xyGlobal); - Int_t idModule=Module(fIdCircuit); // corresponding module Id. -// and its number of Y strips - Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule)); - Int_t sign=TMath::Abs(idModule)/idModule; // left or right - - for (Int_t istrip=0; istrip y position of X declusterized strips - - Int_t chamber, cathode; - Int_t code, idModule, idStrip, idSector; - Float_t x, y, z, width; - Int_t istrip; - - AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON"); - AliMUONChamber* iChamber; - AliSegmentation* segmentation; - -// first plane (11) - chamber=11; - cathode=1; - iChamber = &(pMUON->Chamber(chamber-1)); - segmentation=iChamber->SegmentationModel(cathode); +void AliMUONTriggerCircuit::LoadXPos(AliMpLocalBoard* const localBoard) +{ +/// fill fXpos11 -> x position of Y strips for the first plane only +/// fXpos11 contains the x position of Y strip for the current circuit +/// taking into account whether or nor not part(s) of the circuit +/// (middle, up or down) has(have) 16 strips (handdled by means of switchs) + + fCurrentLocalBoard = localBoard->GetId(); + + Int_t ichamber = 10; + Int_t icathode = 1; - for (istrip=0; istrip<16; istrip++) { - code=fXcode[0][istrip]; // decode current strip - idModule=Int_t(code/100); // corresponding module Id. - idStrip=TMath::Abs(code-idModule*100); // corresp. strip number in module - idSector=segmentation->Sector(idModule,idStrip); // corresponding sector - width=segmentation->Dpy(idSector); // corresponding strip width - segmentation->GetPadC(idModule,idStrip,x,y,z); // get strip real position - - fYpos11[2*istrip]=y; - if (istrip!=15) fYpos11[2*istrip+1]=y+width/2.; - } - -// second plane (21) - chamber=13; - cathode=1; - iChamber = &(pMUON->Chamber(chamber-1)); - segmentation=iChamber->SegmentationModel(cathode); + Int_t x2u = localBoard->GetSwitch(AliMpLocalBoard::kX2u); + Int_t x2m = localBoard->GetSwitch(AliMpLocalBoard::kX2m); + Int_t x2d = localBoard->GetSwitch(AliMpLocalBoard::kX2d); + Int_t zeroAllYLSB = localBoard->GetSwitch(AliMpLocalBoard::kZeroAllYLSB); + + Int_t iStripCircuit = 0; + Int_t iFirstStrip = 0; + Int_t iLastStrip = 0; + Bool_t doubling = kFALSE; - for (istrip=0; istrip<32; istrip++) { - code=fXcode[2][istrip]; // decode current strip - idModule=Int_t(code/100); // corresponding module Id. - idStrip=TMath::Abs(code-idModule*100); // corresp. strip number in module - idSector=segmentation->Sector(idModule,idStrip); // corresponding sector - width=segmentation->Dpy(idSector); // corresponding strip width - segmentation->GetPadC(idModule,idStrip,x,y,z); // get strip real position - -// using idModule!=0 prevents to fill garbage in case of circuits -// in the first and last rows - if (idModule!=0) { - fYpos21[2*istrip]=y; - if (istrip!=31) fYpos21[2*istrip+1]=y+width/2.; - } - } -} + Int_t iline = AliMp::PairFirst(localBoard->GetPosition()); + Int_t icol = AliMp::PairSecond(localBoard->GetPosition()); + if ( iline == 5 ) --icol; -//---------------------------------------------------------------------- -void AliMUONTriggerCircuit::LoadXPos(){ -// fill fXpos11 -> x position of Y strips for the first plane only -// fXpos11 contains the x position of Y strip for the current circuit -// taking into account whether or nor not part(s) of the circuit -// (middle, up or down) has(have) 16 strips + fCurrentDetElem = AliMpDDLStore::Instance()->GetDEfromLocalBoard(fCurrentLocalBoard, ichamber); + + fkCurrentSeg = AliMpSegmentation::Instance() + ->GetMpSegmentation(fCurrentDetElem, AliMp::GetCathodType(icathode)); + + // check if one needs a strip doubling or not + if ( (x2u || x2m || x2d ) && x2m ) doubling = kTRUE; - Float_t x, y, z; - Int_t istrip; - - Int_t chamber=11; - Int_t cathode=2; - AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON"); - AliMUONChamber* iChamber; - AliSegmentation* segmentation; - iChamber = &(pMUON->Chamber(chamber-1)); - segmentation=iChamber->SegmentationModel(cathode); + + // check if one starts at strip = 0 or 8 (boards 26-29 and 143-146) + if (zeroAllYLSB) iStripCircuit = 8; - Int_t idModule=Module(fIdCircuit); // corresponding module Id. -// number of Y strips - Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule)); - Int_t idSector=segmentation->Sector(idModule,0); // corresp. sector - Float_t width=segmentation->Dpx(idSector); // corresponding strip width + // get iFirstStrip in this module + const AliMpTrigger* t = AliMpSegmentation::Instance()->GetTrigger(fkCurrentSeg); + const AliMpSlat* slat = t->GetLayer(0); -// first case : up middle and down parts have all 8 or 16 strip - if ((nStrY==16)||(nStrY==8&&fX2m==0&&fX2ud==0)) { - for (istrip=0; istripGetPadC(idModule,istrip,x,y,z); - fXpos11[istrip]=x; - } -// second case : mixing 8 and 16 strips within same circuit - } else { - for (istrip=0; istripGetPadC(idModule,istrip,x,y,z); - fXpos11[2*istrip]=x-width/4.; - fXpos11[2*istrip+1]=fXpos11[2*istrip]+width/2.; - } - } + const AliMpPCB* pcb = slat->GetPCB(icol-1); + iFirstStrip = pcb->Ixmin(); + + + if (doubling || zeroAllYLSB == 1) iLastStrip = iFirstStrip + 8; + else iLastStrip = iFirstStrip + 16; + + FillYstrips(iFirstStrip, iLastStrip, iStripCircuit, doubling); } //---------------------------------------------------------------------- -Float_t AliMUONTriggerCircuit::PtCal(Int_t istripX, Int_t idev, Int_t istripY){ -// returns calculated pt for circuit/istripX/idev/istripY according -// to the formula of the TRD. Note : idev (input) is in [0+30] +void AliMUONTriggerCircuit::FillYstrips(const Int_t iFirstStrip, const Int_t iLastStrip, + Int_t liStripCircuit, + const Bool_t doubling) +{ +/// fill + Double_t xyGlobal[3] = {0.}; - // Int_t jdev = idev - 15; // jdev in [-15+15] - Int_t istripX2=istripX+idev+1; // find istripX2 using istripX and idev + for (Int_t istrip = iFirstStrip; istrip < iLastStrip; ++istrip) { - Float_t yPosX1=fYpos11[istripX]; - Float_t yPosX2=fYpos21[istripX2]; - Float_t xPosY1=fXpos11[istripY]; - - Float_t zf=975., z1=1603.5, z2=1703.5; - Float_t thetaDev=(1./zf)*(yPosX1*z2-yPosX2*z1)/(z2-z1); - Float_t xf=xPosY1*zf/z1; - Float_t yf=yPosX2-((yPosX2-yPosX1)*(z2-zf))/(z2-z1); - return (3.*0.3/TMath::Abs(thetaDev)) * TMath::Sqrt(xf*xf+yf*yf)/zf; -} + AliMpPad pad = fkCurrentSeg->PadByIndices(istrip,0,kTRUE); -//---------------------------------------------------------------------- -//--- methods which return member data related info -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdCircuit(){ -// returns circuit Id - return fIdCircuit; -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdModule(){ -// returns module Id - return Module(fIdCircuit); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetNstripX() { -// returns the number of X strips in the module where the circuit is sitting - return AliMUONTriggerConstants::NstripX(ModuleNumber(Module(fIdCircuit))); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetNstripY() { -// returns the number of Y strips in the module where the circuit is sitting - return AliMUONTriggerConstants::NstripY(ModuleNumber(Module(fIdCircuit))); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetPosCircuit() { -// returns the position of the circuit in its module - return Position(fIdCircuit); + if ( !pad.IsValid() ) + { + StdoutToAliError(cout << "Pad not found in seg " << endl; + fkCurrentSeg->Print(); + cout << " ix,iy=" << istrip << "," << 0 << endl; + ); + } + Float_t xDim = pad.GetDimensionX(); // half size! + + XYGlobal(pad,xyGlobal); + + if (!doubling) { + fXpos11[fCurrentLocalBoard].AddAt(xyGlobal[0], liStripCircuit); + fXwidth11[fCurrentLocalBoard].AddAt(2. * xDim, liStripCircuit); + } else if (doubling) { + + fXpos11[fCurrentLocalBoard].AddAt(TMath::Sign(1.,xyGlobal[0]) * + (TMath::Abs(xyGlobal[0]) - xDim/2.), 2*liStripCircuit); + fXwidth11[fCurrentLocalBoard].AddAt(2. * xDim, 2*liStripCircuit); + + fXpos11[fCurrentLocalBoard].AddAt(TMath::Sign(1.,xyGlobal[0]) * + (TMath::Abs(fXpos11[fCurrentLocalBoard][2*liStripCircuit]) + xDim), + 2*liStripCircuit + 1); + fXwidth11[fCurrentLocalBoard].AddAt(2. * xDim, 2*liStripCircuit + 1); + } + + liStripCircuit++; + } } + //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdCircuitD(){ -// returns the Id of the circuit down - Int_t idModule=Module(fIdCircuit); - Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule); - return (TMath::Abs(idModuleD)*10+1)*(TMath::Abs(idModule)/idModule); +void AliMUONTriggerCircuit::XYGlobal(const AliMpPad& pad, + Double_t* xyGlobal) +{ +/// returns pad x & y positions and x & y pad dimensions in global coordinates +/// note: no need for transformation for pad dimensions + + // get the pad position and dimensions + Double_t xl1 = pad.GetPositionX(); + Double_t yl1 = pad.GetPositionY(); + + // positions from local to global + fkTransformer->Local2Global(fCurrentDetElem, xl1, yl1, 0, + xyGlobal[0], xyGlobal[1], xyGlobal[2]); } + + //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetICircuitD(){ -// returns the number of the circuit down - Int_t idModule=Module(fIdCircuit); - Int_t idModuleD=(TMath::Abs(idModule)+10)*(TMath::Abs(idModule)/idModule); - Int_t idCircuitD= - (TMath::Abs(idModuleD)*10+1)*(TMath::Abs(idModule)/idModule); - return CircuitNumber(idCircuitD); -} +//--- methods which return member data related info //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdCircuitU(){ -// returns the Id of the circuit up - Int_t idModule=Module(fIdCircuit); - Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule); - return (TMath::Abs(idModuleU)*10+1)*(TMath::Abs(idModule)/idModule); -} //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetICircuitU(){ -// returns the number of the circuit up - Int_t idModule=Module(fIdCircuit); - Int_t idModuleU=(TMath::Abs(idModule)-10)*(TMath::Abs(idModule)/idModule); - Int_t idCircuitU= - (TMath::Abs(idModuleU)*10+1)*(TMath::Abs(idModule)/idModule); - return CircuitNumber(idCircuitU); +Float_t AliMUONTriggerCircuit::GetX11Pos(Int_t localBoardId, Int_t istrip) const +{ +/// returns X position of Y strip istrip in MC11 + return fXpos11[localBoardId][istrip]; } -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetX2m(){ -// returns fX2m - return fX2m; +Float_t AliMUONTriggerCircuit::GetY11Pos(Int_t localBoardId, Int_t istrip) const +{ +/// returns Y position of X strip istrip in MC11 + return fYpos11[localBoardId][istrip]; } //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetX2ud(){ -// returns fX2ud - return fX2ud; +Float_t AliMUONTriggerCircuit::GetY21Pos(Int_t localBoardId, Int_t istrip) const +{ +/// returns Y position of X strip istrip in MC21 + return fYpos21[localBoardId][istrip]; } //---------------------------------------------------------------------- -void AliMUONTriggerCircuit::GetOrMud(Int_t orMud[2]){ -// returns fOrMud - orMud[0]=fOrMud[0]; - orMud[1]=fOrMud[1]; +Float_t AliMUONTriggerCircuit::GetZ11Pos(Int_t localBoardId, Int_t istrip) const +{ +/// returns Z position of X strip istrip in MC11 + return fZpos11[localBoardId][istrip]; } //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetXcode(Int_t chamber, Int_t istrip){ -// returns X code of circuit/chamber/istrip (warning : chamber in [0,3]) - return fXcode[chamber][istrip]; +Float_t AliMUONTriggerCircuit::GetZ21Pos(Int_t localBoardId, Int_t istrip) const +{ +/// returns Z position of X strip istrip in MC21 + return fZpos21[localBoardId][istrip]; } -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetYcode(Int_t chamber, Int_t istrip){ -// returns Y code of circuit/chamber/istrip (warning : chamber in [0,3]) - return fYcode[chamber][istrip]; +Float_t AliMUONTriggerCircuit::GetX11Width(Int_t localBoardId, Int_t istrip) const +{ +/// returns width of Y strip istrip in MC11 + return fXwidth11[localBoardId][istrip]; } -//---------------------------------------------------------------------- -Float_t AliMUONTriggerCircuit::GetY11Pos(Int_t istrip){ -// returns Y position of X strip istrip in MC11 - return fYpos11[istrip]; +Float_t AliMUONTriggerCircuit::GetY11Width(Int_t localBoardId, Int_t istrip) const +{ +/// returns width of X strip istrip in MC11 + return fYwidth11[localBoardId][istrip]; } //---------------------------------------------------------------------- -Float_t AliMUONTriggerCircuit::GetY21Pos(Int_t istrip){ -// returns Y position of X strip istrip in MC21 - return fYpos21[istrip]; +Float_t AliMUONTriggerCircuit::GetY21Width(Int_t localBoardId, Int_t istrip) const +{ +/// returns width of X strip istrip in MC21 + return fYwidth21[localBoardId][istrip]; } + //---------------------------------------------------------------------- -Float_t AliMUONTriggerCircuit::GetX11Pos(Int_t istrip){ -// returns X position of Y strip istrip in MC11 - return fXpos11[istrip]; +Int_t AliMUONTriggerCircuit::FirstStrip(AliMpLocalBoard* localBoard) +{ +/// returns the first strip from mapping for board boardName +/// take care of special case for boards RC1L6B12 & LC1L6B12 + Int_t iFirstStrip = -1; + Int_t boardNumber = atoi(localBoard->GetName()+6); + + Int_t iline = AliMp::PairFirst(localBoard->GetPosition()); + Int_t icol = AliMp::PairSecond(localBoard->GetPosition()); + if ( iline == 5 ) --icol; + + switch (boardNumber) + { + case 12: + iFirstStrip = 0; + break; + case 34: + iFirstStrip = 16; + break; + case 56: + iFirstStrip = 32; + break; + case 78: + iFirstStrip = 48; + break; + } + if (icol == 1 && iline == 6) iFirstStrip += 16; // special case + return iFirstStrip; } -//---------------------------------------------------------------------- -//--- end of methods which return member data related info -//---------------------------------------------------------------------- - - - +//---------------------------------------------------------------------- +Float_t AliMUONTriggerCircuit::PtCal(Int_t localBoardId, Int_t istripX, Int_t idev, Int_t istripY) const{ +/// returns calculated pt for circuit/istripX/idev/istripY according +/// to the formula of the TRD. Note : idev (input) is in [0+30] + Int_t istripX2=istripX+idev+1; // find istripX2 using istripX and idev + Float_t yPosX1=fYpos11[localBoardId][istripX]; + Float_t yPosX2=fYpos21[localBoardId][istripX2]; + Float_t xPosY1=fXpos11[localBoardId][istripY]; +// Z distance between IP and center of dipole + Float_t zf= 0.5 *(AliMUONConstants::CoilZ() + AliMUONConstants::YokeZ()); + Float_t z1=fZpos11[localBoardId][istripX]; + Float_t z2=fZpos21[localBoardId][istripX2]; + Float_t thetaDev=(1./TMath::Abs(zf))*(yPosX1*z2-yPosX2*z1)/(z2-z1); + Float_t xf=xPosY1*zf/z1; + Float_t yf=yPosX2-((yPosX2-yPosX1)*(z2-zf))/(z2-z1); + return (3.*0.3/TMath::Abs(thetaDev)) * TMath::Sqrt(xf*xf+yf*yf)/TMath::Abs(zf); +}