X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FAliMUONTriggerCircuit.cxx;h=1b72adb1ba19204c47c6bda1fb0805dd12629efb;hb=c49cb542f8fdd7654586f9c9d52080d90ce37992;hp=269e9567e1a5541df1e42309e36e856b4f52072d;hpb=e516b01dc9b050949d85a6e082e9213ee783c82c;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/AliMUONTriggerCircuit.cxx b/MUON/AliMUONTriggerCircuit.cxx index 269e9567e1a..1b72adb1ba1 100644 --- a/MUON/AliMUONTriggerCircuit.cxx +++ b/MUON/AliMUONTriggerCircuit.cxx @@ -1,506 +1,588 @@ /************************************************************************** - * 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 +/// \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) #include "AliMUONTriggerCircuit.h" -#include "AliRun.h" #include "AliMUON.h" -#include "AliMUONTriggerConstants.h" -#include "AliSegmentation.h" -#include "AliMUONChamber.h" #include "AliMUONConstants.h" +#include "AliMUONLocalTriggerBoard.h" +#include "AliMUONTriggerCrateStore.h" +#include "AliMUONTriggerCrate.h" +#include "AliMUONGeometryTransformer.h" + +#include "AliMpTriggerSegmentation.h" +#include "AliMpTrigger.h" +#include "AliMpSlat.h" +#include "AliMpPCB.h" +#include "AliMpSegmentation.h" +#include "AliMpVSegmentation.h" +#include "AliMpCathodType.h" + +#include "AliRun.h" +#include "AliLog.h" + +#include +/// \cond CLASSIMP ClassImp(AliMUONTriggerCircuit) +/// \endcond //---------------------------------------------------------------------- AliMUONTriggerCircuit::AliMUONTriggerCircuit() - : TObject() +: TObject(), + fILocalBoard(0), + fTransformer(0x0) { -// Constructor - fSegmentation=0; - fIdCircuit=0; - fX2m=0; - fX2ud=0; - fOrMud[0]=fOrMud[1]=0; +/// Constructor + 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; - } - } 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& theMUONTriggerCircuit) - : TObject(theMUONTriggerCircuit) +AliMUONTriggerCircuit::~AliMUONTriggerCircuit() { -// Protected copy constructor - - Fatal("AliMUONTriggerCircuit", "Not implemented."); -} +/// Destructor +} //---------------------------------------------------------------------- -AliMUONTriggerCircuit & -AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& rhs) +AliMUONTriggerCircuit::AliMUONTriggerCircuit(const AliMUONTriggerCircuit& circuit) + : TObject(circuit), + fILocalBoard(circuit.fILocalBoard), + fTransformer(0x0) { -// Protected assignement operator +/// Copy constructor - if (this == &rhs) return *this; + for (Int_t i = 0; i < 16; ++i) + fXpos11[i] = circuit.fXpos11[i]; - Fatal("operator=", "Not implemented."); - - return *this; -} + for (Int_t i = 0; i < 31; ++i) + fYpos11[i] = circuit.fYpos11[i]; -//---------------------------------------------------------------------- -void AliMUONTriggerCircuit::Init(Int_t iCircuit) { -// initialize circuit characteristics - fIdCircuit=AliMUONTriggerConstants::CircuitId(iCircuit); - LoadX2(); - LoadXCode(); - LoadYCode(); - LoadXPos(); - LoadYPos(); -} + for (Int_t i = 0; i < 63; ++i) + fYpos21[i] = circuit.fYpos21[i]; -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::CircuitNumber(Int_t idCircuit) const { -// 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; - } - } - return iCircuit; } //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::ModuleNumber(Int_t idModule) const { -// 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; - } - } - return iModule; -} +AliMUONTriggerCircuit& AliMUONTriggerCircuit::operator=(const AliMUONTriggerCircuit& circuit) +{ +/// Assignment operator + + if (this == &circuit) return *this; + + fILocalBoard = circuit.fILocalBoard; + + for (Int_t i = 0; i < 16; ++i) + fXpos11[i] = circuit.fXpos11[i]; + + for (Int_t i = 0; i < 31; ++i) + fYpos11[i] = circuit.fYpos11[i]; + + for (Int_t i = 0; i < 63; ++i) + fYpos21[i] = circuit.fYpos21[i]; + + return *this; -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::Module(Int_t idCircuit) const { -// returns ModuleId where Circuit idCircuit is sitting - return Int_t(idCircuit/10); } //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::Position(Int_t idCircuit) const { -// returns position of idCircuit in correcponding Module - return TMath::Abs(idCircuit)-TMath::Abs(Module(idCircuit))*10; +void AliMUONTriggerCircuit::Init(Int_t iCircuit, const AliMUONTriggerCrateStore& crates) +{ +/// initialize circuit characteristics + fILocalBoard=iCircuit+1;//AliMUONTriggerConstants::CircuitId(iCircuit); + + LoadXPos(crates); + LoadYPos(crates); + } -//---------------------------------------------------------------------- -void AliMUONTriggerCircuit::LoadX2() { -// initialize fX2m, fX2ud and fOrMud +//--------------------------------------------------------------------- +void AliMUONTriggerCircuit::LoadYPos(const AliMUONTriggerCrateStore& crates) +{ +/// fill fYpos11 and fYpos21 -> y position of X declusterized strips - 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 + const AliMUONLocalTriggerBoard* localBoard = crates.LocalBoard(fILocalBoard); -// 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; - } - -// 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; - } - } - - } else { // upper strips within same module - fOrMud[1]=0; + if (!localBoard) + { + AliError(Form("Did not get localboard %d",fILocalBoard)); + return; } -} - -//---------------------------------------------------------------------- -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++; - } - -// 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; istripPrint("CONF");); + + Int_t ichamber = 0; + Int_t icathode = 0; + + const AliMpVSegmentation* seg; + + Int_t zeroDown = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kZeroDown); + Int_t zeroUp = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kZeroUp); + + //--- first plane + ichamber = 10; + + char side; + Int_t iline, icol; + DecodeBoardName(localBoard->GetName(),side,iline,icol); + + Int_t detElemId = DetElemId(ichamber,side,iline); + seg = AliMpSegmentation::Instance() + ->GetMpSegmentation(detElemId, AliMp::GetCathodType(icathode)); + + Int_t iFirstStrip = FirstStrip(localBoard->GetName()); + Int_t iLastStrip = iFirstStrip + 16; + Int_t iStripCircuit = 0; + FillXstrips(seg,detElemId,icol, + iFirstStrip,iLastStrip,iStripCircuit,fYpos11); + + //--- second plane + ichamber = 12; + + detElemId = DetElemId(ichamber,side,iline); + seg = AliMpSegmentation::Instance() + ->GetMpSegmentation(detElemId, AliMp::GetCathodType(icathode)); + + // second plane middle part + Int_t iFirstStripMiddle = FirstStrip(localBoard->GetName()); + Int_t iLastStripMiddle = iFirstStrip + 16; + iStripCircuit=8; + FillXstrips(seg,detElemId,icol, + iFirstStripMiddle,iLastStripMiddle,iStripCircuit,fYpos21); + + // second plane upper part + if (zeroUp == 0) { // something up + Int_t iFirstStripUp; + Int_t iLastStripUp; + Int_t icolUp=icol; + // check if we need to move to another detElemId + AliMpPad pad = seg->PadByIndices(AliMpIntPair(icol-1,iLastStripMiddle+1),kFALSE); + if (pad.IsValid()) { // upper strips within same detElemId + iFirstStripUp = iLastStripMiddle; + iLastStripUp = iFirstStripUp + 8; + // icolUp = icol; + } else { // upper strips in another detElemId + detElemId = DetElemId(ichamber,side,iline+1); // get detElemId + seg = AliMpSegmentation::Instance() + ->GetMpSegmentation(detElemId, 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(seg,detElemId,icolUp, + iFirstStripUp,iLastStripUp,iStripCircuit,fYpos21); + + // fill strip between middle and upper part + fYpos21[47]=(fYpos21[46]+fYpos21[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 + detElemId = DetElemId(ichamber,side,iline); + seg = AliMpSegmentation::Instance() + ->GetMpSegmentation(detElemId, 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 = seg->PadByIndices(AliMpIntPair(icol-1,iFirstStripMiddle-1),kFALSE); + if (pad.IsValid()) { // lower strips within same detElemId + iFirstStripDo = iFirstStripMiddle - 8; + iLastStripDo = iFirstStripDo + 8; + // icolDo = icol; + } else { // lower strips in another detElemId + detElemId = DetElemId(ichamber,side,iline-1); // get detElemId + seg = AliMpSegmentation::Instance() + ->GetMpSegmentation(detElemId, AliMp::GetCathodType(icathode)); + + // get iFirstStrip in this module + const AliMpTriggerSegmentation* trig = (AliMpTriggerSegmentation*)(seg); + const AliMpTrigger* t = trig->Slat(); + 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++; - } - } -} - -//---------------------------------------------------------------------- -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) + iStripCircuit=0; + FillXstrips(seg,detElemId,icolDo, + iFirstStripDo,iLastStripDo,iStripCircuit,fYpos21); - 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(const AliMUONTriggerCrateStore& crates) +{ +/// 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) - 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); + const AliMUONLocalTriggerBoard* localBoard = crates.LocalBoard(fILocalBoard); - 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.; - } - } -} - -//---------------------------------------------------------------------- -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 + if (!localBoard) + { + AliError(Form("Did not get localboard %d",fILocalBoard)); + return; + } + StdoutToAliDebug(1,localBoard->Print("CONF");); - 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); + Int_t ichamber = 10; + Int_t icathode = 1; - 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 + Int_t x2u = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kX2u); + Int_t x2m = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kX2m); + Int_t x2d = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kX2d); + Int_t zeroAllYLSB = localBoard->GetSwitch(AliMUONLocalTriggerBoard::kZeroAllYLSB); + Int_t iStripCircuit = 0; + Int_t iFirstStrip = 0; + Int_t iLastStrip = 0; + Bool_t doubling = kFALSE; -// 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 AliMpVSegmentation* seg; -//---------------------------------------------------------------------- -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] - - // Int_t jdev = idev - 15; // jdev in [-15+15] - Int_t istripX2=istripX+idev+1; // find istripX2 using istripX and idev + char side; + Int_t iline, icol; + + DecodeBoardName(localBoard->GetName(),side,iline,icol); + + Int_t detElemId=DetElemId(ichamber,side,iline); // get detElemId + seg = AliMpSegmentation::Instance() + ->GetMpSegmentation(detElemId, AliMp::GetCathodType(icathode)); - Float_t yPosX1=fYpos11[istripX]; - Float_t yPosX2=fYpos21[istripX2]; - Float_t xPosY1=fXpos11[istripY]; + // check if one needs a strip doubling or not + if ( (x2u == 1 || x2m == 1 || x2d == 1) && x2m == 1) doubling = kTRUE; - Float_t zf=975.; - Float_t z1=AliMUONConstants::DefaultChamberZ(10); - Float_t z2=AliMUONConstants::DefaultChamberZ(12); - 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; + // check if one starts at strip = 0 or 8 (boards 26-29 and 143-146) + if (zeroAllYLSB == 1) iStripCircuit = 8; + + // get iFirstStrip in this module + const AliMpTriggerSegmentation* trig = (AliMpTriggerSegmentation*)(seg); + const AliMpTrigger* t = trig->Slat(); + const AliMpSlat* slat = t->GetLayer(0); + + const AliMpPCB* pcb = slat->GetPCB(icol-1); + iFirstStrip = pcb->Ixmin(); + + if (doubling) iLastStrip = iFirstStrip + 8; + else iLastStrip = iFirstStrip + 16; + + FillYstrips(seg,detElemId, + iFirstStrip,iLastStrip,iStripCircuit,doubling); } //---------------------------------------------------------------------- -//--- methods which return member data related info -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdCircuit() const { -// returns circuit Id - return fIdCircuit; -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdModule() const { -// returns module Id - return Module(fIdCircuit); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetNstripX() const { -// returns the number of X strips in the module where the circuit is sitting - return AliMUONTriggerConstants::NstripX(ModuleNumber(Module(fIdCircuit))); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetNstripY() const { -// returns the number of Y strips in the module where the circuit is sitting - return AliMUONTriggerConstants::NstripY(ModuleNumber(Module(fIdCircuit))); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetPosCircuit() const { -// returns the position of the circuit in its module - return Position(fIdCircuit); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdCircuitD() const { -// 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); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetICircuitD() const { -// 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); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetIdCircuitU() const { -// 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() const { -// 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); -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetX2m() const { -// returns fX2m - return fX2m; -} -//---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetX2ud() const { -// returns fX2ud - return fX2ud; -} -//---------------------------------------------------------------------- -void AliMUONTriggerCircuit::GetOrMud(Int_t orMud[2]) const { -// returns fOrMud - orMud[0]=fOrMud[0]; - orMud[1]=fOrMud[1]; +void AliMUONTriggerCircuit::FillYstrips( + const AliMpVSegmentation* seg, + const Int_t detElemId, + const Int_t iFirstStrip, const Int_t iLastStrip, Int_t liStripCircuit, + const Bool_t doubling) +{ +/// fill + Double_t xyGlobal[4]={0.,0.,0.,0.}; + for (Int_t istrip=iFirstStrip; istripPadByIndices(AliMpIntPair(istrip,0),kTRUE); + if ( !pad.IsValid() ) + { + StdoutToAliError(cout << "Pad not found in seg " << endl; + seg->Print(); + cout << " ix,iy=" << istrip << "," << 0 << endl; + ); + } + XYGlobal(detElemId,pad,xyGlobal); + + if (!doubling) { + fXpos11[liStripCircuit]=xyGlobal[0]; + } else if (doubling) { + fXpos11[2*liStripCircuit]=TMath::Sign(1.,xyGlobal[0]) * + (TMath::Abs(xyGlobal[0]) - xyGlobal[2]/2.); + fXpos11[2*liStripCircuit+1]=TMath::Sign(1.,xyGlobal[0]) * + (TMath::Abs(fXpos11[2*liStripCircuit]) + xyGlobal[2]); + } + liStripCircuit++; + } } + //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetXcode(Int_t chamber, Int_t istrip) const { -// returns X code of circuit/chamber/istrip (warning : chamber in [0,3]) - return fXcode[chamber][istrip]; +void AliMUONTriggerCircuit::FillXstrips( + const AliMpVSegmentation* seg, + const Int_t detElemId, const Int_t icol, + const Int_t iFirstStrip, const Int_t iLastStrip, + Int_t liStripCircuit, Float_t *tab) +{ +/// fill + Double_t xyGlobal[4]={0.,0.,0.,0.}; + for (Int_t istrip=iFirstStrip; istripPadByIndices(AliMpIntPair(icol-1,istrip),kTRUE); + if ( !pad.IsValid() ) + { + StdoutToAliError(cout << "Pad not found in seg " << endl; + seg->Print(); + cout << " ix,iy=" << icol-1 << "," << istrip << endl; + ); + } + + XYGlobal(detElemId,pad,xyGlobal); + + tab[2*liStripCircuit]=xyGlobal[1]; + if (istrip!=(iLastStrip-1)) tab[2*liStripCircuit+1]=xyGlobal[1]+xyGlobal[3]; + liStripCircuit++; + } } + //---------------------------------------------------------------------- -Int_t AliMUONTriggerCircuit::GetYcode(Int_t chamber, Int_t istrip) const { -// returns Y code of circuit/chamber/istrip (warning : chamber in [0,3]) - return fYcode[chamber][istrip]; -} +//--- methods which return member data related info //---------------------------------------------------------------------- Float_t AliMUONTriggerCircuit::GetY11Pos(Int_t istrip) const { -// returns Y position of X strip istrip in MC11 +/// returns Y position of X strip istrip in MC11 return fYpos11[istrip]; } //---------------------------------------------------------------------- Float_t AliMUONTriggerCircuit::GetY21Pos(Int_t istrip) const { -// returns Y position of X strip istrip in MC21 +/// returns Y position of X strip istrip in MC21 return fYpos21[istrip]; } //---------------------------------------------------------------------- Float_t AliMUONTriggerCircuit::GetX11Pos(Int_t istrip) const { -// returns X position of Y strip istrip in MC11 +/// returns X position of Y strip istrip in MC11 return fXpos11[istrip]; } //---------------------------------------------------------------------- //--- end of methods which return member data related info //---------------------------------------------------------------------- +/* removed tmp +void AliMUONTriggerCircuit::dump(const char* what, const Float_t* array, Int_t size) +{ + cout << what << " " << endl; + for ( Int_t i = 0; i < size; ++i ) + { + cout << array[i] << " , "; + } + cout << endl; +} +void AliMUONTriggerCircuit::dump(const char* what, const Int_t* array, Int_t size) +{ + cout << what << " " << endl; + for ( Int_t i = 0; i < size; ++i ) + { + cout << array[i] << " , "; + } + cout << endl; +} +//_____________________________________________________________________________ +void AliMUONTriggerCircuit::Print(Option_t* ) const +{ + cout << "IdCircuit " << fILocalBoard << " X2m,X2ud=" << fX2m << "," + << fX2ud; + for ( Int_t i = 0; i < 2; ++i ) + { + cout << " OrMud[" << i << "]=" << fOrMud[i]; + } + cout << endl; + dump("Xpos11",fXpos11,16); + dump("Ypos11",fYpos11,31); + dump("Ypos21",fYpos21,63); + for ( Int_t i = 0; i < 4; ++i ) + { + char s[80]; + sprintf(s,"Xcode[%d]",i); + dump(s,fXcode[i],32); + sprintf(s,"Ycode[%d]",i); + dump(s,fYcode[i],32); + } + // Int_t fILocalBoard; // circuit Id number + // Int_t fX2m; // internal info needed by TriggerDecision + // Int_t fX2ud; // internal info needed by TriggerDecision + // Int_t fOrMud[2]; // internal info needed by TriggerDecision + // Int_t fXcode[4][32]; // code of X strips + // Int_t fYcode[4][32]; // code of Y strips + // Float_t fXpos11[16]; // X position of Y strips in MC11 + // Float_t fYpos11[31]; // Y position of X strips in MC11 + // Float_t fYpos21[63]; // Y position of X strips in MC21 + +} +removed tmp*/ +//---------------------------------------------------------------------- +Int_t AliMUONTriggerCircuit::DetElemId(Int_t ichamber, char side, Int_t iline) +{ +/// returns detection element Id for chamber iChamber, side side and line iline + Int_t itmp=0; + if ( side == 'R' ) { // right side + switch (iline) // (from 1 to 9, from bottom to top) + { + case 1: + itmp = 14; + break; + case 2: + itmp = 15; + break; + case 3: + itmp = 16; + break; + case 4: + itmp = 17; + break; + case 5: + itmp = 0; + break; + case 6: + itmp = 1; + break; + case 7: + itmp = 2; + break; + case 8: + itmp = 3; + break; + case 9: + itmp = 4; + break; + } + } else if ( side == 'L' ) { // left side + switch (iline) // (from 1 to 9, from bottom to top) + { + case 1: + itmp = 13; + break; + case 2: + itmp = 12; + break; + case 3: + itmp = 11; + break; + case 4: + itmp = 10; + break; + case 5: + itmp = 9; + break; + case 6: + itmp = 8; + break; + case 7: + itmp = 7; + break; + case 8: + itmp = 6; + break; + case 9: + itmp = 5; + break; + } + } + return ((ichamber+1)*100)+itmp; +} +//---------------------------------------------------------------------- +Int_t +AliMUONTriggerCircuit::DetElemId(Int_t iChamber, const char* boardName) +{ +/// returns detection element Id for chamber iChamber and board boardName + char side; + Int_t iline; + Int_t icol; + DecodeBoardName(boardName, side, iline, icol); + + return DetElemId(iChamber,side,iline); +} + +//---------------------------------------------------------------------- +void +AliMUONTriggerCircuit::DecodeBoardName(const char* boardName, + char& side, + Int_t& iLine, + Int_t& iCol) +{ +/// get side, line and col from board boardName +/// note: icol = icol -1 for iline = 5 w.r.t other ilines + side = boardName[0]; + iLine = boardName[4] - '0'; + iCol = boardName[2] - '0'; + if ( iLine == 5 ) --iCol; +} +//---------------------------------------------------------------------- +Int_t +AliMUONTriggerCircuit::FirstStrip(const char* boardName) +{ +/// 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(boardName+6); + char side; + Int_t iline, icol; + DecodeBoardName(boardName,side,iline,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 = iFirstStrip + 16; // special case + return iFirstStrip; +} +//---------------------------------------------------------------------- +void AliMUONTriggerCircuit::XYGlobal( + Int_t detElemId, const AliMpPad& pad, + Double_t xyGlobal[4]) +{ +/// 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.Position().X(); + Double_t yl1 = pad.Position().Y(); + xyGlobal[2] = pad.Dimensions().X(); // half size! + xyGlobal[3] = pad.Dimensions().Y(); // half size! + Double_t zg1=0; + + // positions from local to global + fTransformer->Local2Global(detElemId, xl1, yl1, 0, + xyGlobal[0], xyGlobal[1], zg1); +}