#include "AliRun.h"
#include "AliMUON.h"
#include "AliMUONTriggerConstants.h"
-#include "AliSegmentation.h"
#include "AliMUONGeometrySegmentation.h"
-#include "AliMUONChamber.h"
+#include "AliMUONSegmentation.h"
#include "AliMUONConstants.h"
#include "AliLog.h"
fIdCircuit(0),
fX2m(0),
fX2ud(0)
- // fSegmentationType(1)
{
// Constructor
void AliMUONTriggerCircuit::Init(Int_t iCircuit) {
// initialize circuit characteristics
fIdCircuit=AliMUONTriggerConstants::CircuitId(iCircuit);
- AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON");
LoadX2();
LoadXCode();
LoadYCode();
- if (pMUON->WhichSegmentation() == 1) {
- LoadXPos();
- LoadYPos();
- } else {
- LoadXPos2();
- LoadYPos2();
- }
+ LoadXPos2();
+ LoadYPos2();
+
}
//----------------------------------------------------------------------
}
}
-//----------------------------------------------------------------------
-void AliMUONTriggerCircuit::LoadYPos(){
-// fill fYpos11 and fYpos21 -> 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);
-
- 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);
-
- 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
-
- 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 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
-
-// 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; istrip<nStrY; istrip++) {
- segmentation->GetPadC(idModule,istrip,x,y,z);
- fXpos11[istrip]=x;
- }
-// second case : mixing 8 and 16 strips within same circuit
- } else {
- for (istrip=0; istrip<nStrY; istrip++) {
- if (nStrY!=8) { printf(" bug in LoadXpos \n");}
- segmentation->GetPadC(idModule,istrip,x,y,z);
- fXpos11[2*istrip]=x-width/4.;
- fXpos11[2*istrip+1]=fXpos11[2*istrip]+width/2.;
- }
- }
-}
-
//----------------------------------------------------------------------
Float_t AliMUONTriggerCircuit::PtCal(Int_t istripX, Int_t idev, Int_t istripY){
// returns calculated pt for circuit/istripX/idev/istripY according
Int_t istrip, idDE;
AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON");
- AliMUONChamber* iChamber;
AliMUONGeometrySegmentation* segmentation;
// first plane (11)
chamber=11;
cathode=1;
- iChamber = &(pMUON->Chamber(chamber-1));
- segmentation=iChamber->SegmentationModel2(cathode);
+ segmentation
+ = pMUON->GetSegmentation()->GetModuleSegmentation(chamber-1, cathode-1);
if (!segmentation) {
AliWarning("Segmentation not defined.");
// second plane (21)
chamber=13;
cathode=1;
- iChamber = &(pMUON->Chamber(chamber-1));
- segmentation=iChamber->SegmentationModel2(cathode);
+ segmentation
+ = pMUON->GetSegmentation()->GetModuleSegmentation(chamber-1, cathode-1);
for (istrip=0; istrip<32; istrip++) {
code=fXcode[2][istrip]; // decode current strip
Int_t chamber=11;
Int_t cathode=2;
AliMUON *pMUON = (AliMUON*)gAlice->GetModule("MUON");
- AliMUONChamber* iChamber;
AliMUONGeometrySegmentation* segmentation;
- iChamber = &(pMUON->Chamber(chamber-1));
- segmentation=iChamber->SegmentationModel2(cathode);
-
+ segmentation
+ = pMUON->GetSegmentation()->GetModuleSegmentation(chamber-1, cathode-1);
+
if (!segmentation) {
AliWarning("Segmentation not defined.");
return;
Int_t idModule=Module(fIdCircuit); // corresponding module Id.
// number of Y strips
idDE = DetElemId(chamber, idModule);
+
Int_t nStrY=AliMUONTriggerConstants::NstripY(ModuleNumber(idModule));
Int_t idSector=segmentation->Sector(idDE, idModule,0); // corresp. sector
Float_t width=segmentation->Dpx(idDE, idSector); // corresponding strip width
if ((nStrY==16)||(nStrY==8&&fX2m==0&&fX2ud==0)) {
for (istrip=0; istrip<nStrY; istrip++) {
segmentation->GetPadC(idDE, idModule,istrip,x,y,z);
+ AliDebug(1,Form("idDE %d idModule %d istrip %d x,y,z=%e,%e,%e\n",
+ idDE,idModule,istrip,x,y,z));
fXpos11[istrip]=x;
}
// second case : mixing 8 and 16 strips within same circuit
//--- end of methods which return member data related info
//----------------------------------------------------------------------
+void 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 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 " << fIdCircuit << " 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 fIdCircuit; // 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
+
+}
Int_t AliMUONTriggerCircuit::DetElemId(Int_t ichamber, Int_t idModule)
{
+// adpated to official numbering (09/20/05)
// returns the detection element Id for given chamber and module
-// ichamber (from 11 to 14), idModule (from 11 to 97)
- Int_t itmp = (idModule > 0) ? 0 : 50; // right=0, left=50
- return (ichamber*100)+itmp+(9-Int_t(TMath::Abs(idModule)/10));
+// ichamber (from 11 to 14), idModule (from -97 to 97)
+//
+ Int_t itmp=0;
+ Int_t linenumber=Int_t(idModule/10);
+ switch (linenumber) // (from 1 to 9, from top to bottom)
+ {
+ case 1:
+ itmp = 4;
+ break;
+ case 2:
+ itmp = 3;
+ break;
+ case 3:
+ itmp = 2;
+ break;
+ case 4:
+ itmp = 1;
+ break;
+ case 5:
+ itmp = 0;
+ break;
+ case 6:
+ itmp = 17;
+ break;
+ case 7:
+ itmp = 16;
+ break;
+ case 8:
+ itmp = 15;
+ break;
+ case 9:
+ itmp = 14;
+ break;
+// left
+ case -1:
+ itmp = 5;
+ break;
+ case -2:
+ itmp = 6;
+ break;
+ case -3:
+ itmp = 7;
+ break;
+ case -4:
+ itmp = 8;
+ break;
+ case -5:
+ itmp = 9;
+ break;
+ case -6:
+ itmp = 10;
+ break;
+ case -7:
+ itmp = 11;
+ break;
+ case -8:
+ itmp = 12;
+ break;
+ case -9:
+ itmp = 13;
+ break;
+ }
+ return (ichamber*100)+itmp;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff