add AliMUONTriggerSegmentation

diff --git a/MUON/AliMUONTest.cxx b/MUON/AliMUONTest.cxx
index d3cdc4e94defcdf26acc7c39d81f5456299b5b21..c79c7a056671d39038f4cd3e61c23ae4b5ca1092 100644 (file)
--- a/MUON/AliMUONTest.cxx
+++ b/MUON/AliMUONTest.cxx
@@ -39,6 +39,8 @@
 #include "AliMUONGeometrySegmentation.h"
 #include "AliMUONSt12QuadrantSegmentation.h"
 #include "AliMUONSt345SlatSegmentation.h"
+#include "AliMUONTriggerSegmentation.h"
+#include "AliMUONTriggerConstants.h"
 
 ClassImp(AliMUONTest)
 
@@ -548,6 +550,280 @@ AliMUONTest::CreateSlatSegmentation(Int_t chamberId, Int_t cathod)
   return chamberSeg;
   
 } 
+//_____________________________________________________________________________
+AliMUONGeometrySegmentation*     
+AliMUONTest::CreateTriggerSegmentation(Int_t chamberId, Int_t cathod)
+{
+// Create Trigger geometry segmentation for given chamber and cathod
+
+    printf("in CreateTriggerSegmentation chamber=%d cathode=%d\n",
+          chamberId,cathod);
+  AliMUON* muon = (AliMUON*)gAlice->GetModule("MUON");
+  if (!muon) {
+    AliFatal("MUON detector not defined.");
+    return 0;
+  }  
+
+  AliMUONGeometrySegmentation *chamberSeg = new AliMUONGeometrySegmentation(muon->Chamber(chamberId).GetGeometry());
+
+   //Trigger Segmentations
+  AliMUONTriggerSegmentation *trigSegX[9]; 
+  AliMUONTriggerSegmentation *trigSegY[9]; 
+  for(Int_t i=0; i<9; i++) {
+    trigSegX[i] = new AliMUONTriggerSegmentation(0);
+    trigSegY[i] = new AliMUONTriggerSegmentation(1);
+  }
+
+  AliMUONChamber *iChamber, *iChamber1;
+  iChamber1 = &muon->Chamber(10);
+  iChamber  = &muon->Chamber(chamberId);
+  Float_t zpos1= - iChamber1->Z();  
+  Float_t zpos = - iChamber->Z();           
+  Float_t zRatio = zpos / zpos1;
+
+// init
+  Int_t nStrip[7]={0,0,0,0,0,0,0};       
+  Float_t stripYsize[7]={0.,0.,0.,0.,0.,0.,0.};
+  Float_t stripXsize[7]={0.,0.,0.,0.,0.,0.,0.};
+
+// chamber 8 cathode 0
+  for (Int_t i=0; i<7; i++) nStrip[i]=16;
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio; 
+  trigSegX[8]->Init(0,1,nStrip,stripYsize,stripXsize,0.);  
+ 
+// chamber 8 cathode 1
+  for (Int_t i=0; i<6; i++) nStrip[i]=8;
+  nStrip[6]=16;
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;
+  for (Int_t i=0; i<7; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[8]->Init(0,1,nStrip,stripYsize,stripXsize,0.);  
+ 
+// chamber 7 cathode 0
+  for (Int_t i=0; i<6; i++) nStrip[i]=32;
+  nStrip[7]=16;  
+  for (Int_t i=0; i<6; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripYsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  trigSegX[7]->Init(0,2,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 7 cathode 1
+  for (Int_t i=0; i<6; i++) nStrip[i]=8;
+  nStrip[6]=16;
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;
+  for (Int_t i=0; i<7; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[7]->Init(0,2,nStrip,stripYsize,stripXsize,0.);  
+ 
+// chamber 6 cathode 0
+  for (Int_t i=0; i<6; i++) nStrip[i]=32;
+  nStrip[7]=16;  
+  for (Int_t i=0; i<6; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripYsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  trigSegX[6]->Init(0,3,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 6 cathode 1
+  for (Int_t i=0; i<5; i++) nStrip[i]=16;
+  nStrip[5]=8;
+  nStrip[6]=16;  
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;  
+  for (Int_t i=0; i<5; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripXsize[5]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[6]->Init(0,3,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 5 cathode 0
+  nStrip[0]=48;
+  for (Int_t i=1; i<3; i++) nStrip[i]=64;  
+  for (Int_t i=3; i<6; i++) nStrip[i]=32;
+  nStrip[6]=16;
+  for (Int_t i=0; i<3; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(0)*zRatio;
+  for (Int_t i=3; i<6; i++) 
+  stripYsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripYsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  trigSegX[5]->Init(0,4,nStrip,stripYsize,stripXsize,AliMUONTriggerConstants::StripLength(0)*zRatio);  
+
+// chamber 5 cathode 1
+  for (Int_t i=0; i<5; i++) nStrip[i]=16;
+  nStrip[5]=8;  
+  nStrip[6]=16;
+  stripYsize[0]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  for (Int_t i=1; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;
+  for (Int_t i=0; i<5; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripXsize[5]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[5]->Init(0,4,nStrip,stripYsize,stripXsize,AliMUONTriggerConstants::StripLength(0)*zRatio);  
+
+// chamber 4 cathode 0
+  nStrip[0]=0;
+  for (Int_t i=1; i<3; i++) nStrip[i]=64;  
+  for (Int_t i=3; i<6; i++) nStrip[i]=32;  
+  nStrip[6]=16;
+  stripYsize[0]=0.;
+  for (Int_t i=1; i<3; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(0)*zRatio;
+  for (Int_t i=3; i<6; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripYsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  stripXsize[0]=0;  
+  stripXsize[1]=AliMUONTriggerConstants::StripLength(0)*zRatio;
+  for (Int_t i=2; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  trigSegX[4]->Init(0,5,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 4 cathode 1
+  nStrip[0]=0;  
+  nStrip[1]=8;  
+  for (Int_t i=2; i<5; i++) nStrip[i]=16;
+  nStrip[5]=8;
+  nStrip[6]=16;  
+  stripYsize[0]=0.;  
+  for (Int_t i=1; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;
+  stripXsize[0]=0.;
+  for (Int_t i=1; i<5; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripXsize[5]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[4]->Init(0,5,nStrip,stripYsize,stripXsize,0.);
+
+// chamber 3 cathode 0
+  nStrip[0]=48;
+  for (Int_t i=1; i<3; i++) nStrip[i]=64;  
+  for (Int_t i=3; i<6; i++) nStrip[i]=32;
+  nStrip[6]=16;
+  for (Int_t i=0; i<3; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(0)*zRatio;
+  for (Int_t i=3; i<6; i++) 
+  stripYsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripYsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  trigSegX[3]->Init(0,6,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 3 cathode 1
+  for (Int_t i=0; i<5; i++) nStrip[i]=16;
+  nStrip[5]=8;  
+  nStrip[6]=16;
+  stripYsize[0]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  for (Int_t i=1; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;
+  for (Int_t i=0; i<5; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripXsize[5]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[3]->Init(0,6,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 2 cathode 0
+  for (Int_t i=0; i<6; i++) nStrip[i]=32;
+  nStrip[7]=16;  
+  for (Int_t i=0; i<6; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripYsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  trigSegX[2]->Init(0,7,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 2 cathode 1
+  for (Int_t i=0; i<5; i++) nStrip[i]=16;
+  nStrip[5]=8;
+  nStrip[6]=16;  
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;  
+  for (Int_t i=0; i<5; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripXsize[5]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[2]->Init(0,7,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 1 cathode 0
+  for (Int_t i=0; i<6; i++) nStrip[i]=32;
+  nStrip[7]=16;  
+  for (Int_t i=0; i<6; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(1)*zRatio;
+  stripYsize[6]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio;
+  trigSegX[1]->Init(0,8,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 1 cathode 1
+  for (Int_t i=0; i<6; i++) nStrip[i]=8;
+  nStrip[6]=16;
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;
+  for (Int_t i=0; i<7; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[1]->Init(0,8,nStrip,stripYsize,stripXsize,0.);  
+
+// chamber 0 cathode 0
+  for (Int_t i=0; i<7; i++) nStrip[i]=16;
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  for (Int_t i=0; i<6; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripLength(1)*zRatio;
+  stripXsize[6]=AliMUONTriggerConstants::StripLength(2)*zRatio; 
+  trigSegX[0]->Init(0,9,nStrip,stripYsize,stripXsize,0.);  
+ 
+// chamber 0 cathode 1
+  for (Int_t i=0; i<6; i++) nStrip[i]=8;
+  nStrip[6]=16;
+  for (Int_t i=0; i<7; i++) 
+      stripYsize[i]=AliMUONTriggerConstants::StripLength(3)*zRatio;
+  for (Int_t i=0; i<7; i++) 
+      stripXsize[i]=AliMUONTriggerConstants::StripWidth(2)*zRatio;
+  trigSegY[0]->Init(0,9,nStrip,stripYsize,stripXsize,0.);
+
+  Int_t icount=chamberId-10;  // chamber counter (0 1 2 3)
+  Int_t id0=(10+icount+1)*100;
+
+  printf("in CreateTriggerSegmentation here 0 id0=%i \n",id0);  
+
+  for (Int_t i=0; i<9; i++) {      
+      if (cathod==0) {   
+         chamberSeg->Add(id0+i,     trigSegX[i]);
+         chamberSeg->Add(id0+50+i,  trigSegX[i]);
+      } else if (cathod==1) {       
+         chamberSeg->Add(id0+i,     trigSegY[i]);
+         chamberSeg->Add(id0+50+i,  trigSegY[i]);
+      }
+  }
+
+  printf("in CreateTriggerSegmentation here 1\n");  
+
+  if (!id0) {
+      AliWarning(Form("Segmentation for chamber %d , cathod %d is not yet defined",chamberId, cathod));
+      return 0;      
+  }
+
+  DrawSegmentation(chamberSeg);  
+  return chamberSeg;
+}
+
 
 //
 // public methods
@@ -589,9 +865,7 @@ AliMUONTest::CreateSegmentation(Int_t chamberId, Int_t cath)
     case 11: 
     case 12: 
     case 13:
-        //return CreateTriggerSegmentation(chamberId, cath);
-        AliWarning("Not yet implemented");
-       return 0;
+        return CreateTriggerSegmentation(chamberId, cath);
        break;
 
     default:
@@ -600,7 +874,6 @@ AliMUONTest::CreateSegmentation(Int_t chamberId, Int_t cath)
        break;
   }    
 }              
-
 //______________________________________________________________________________
 void  AliMUONTest::DetElemTransforms()
 {
@@ -709,7 +982,8 @@ void AliMUONTest::ForWhole(AliMUONTests testCase)
   timer.Start();  
 
   // Loop over chambers
-  for (Int_t iChamber=0; iChamber<AliMUONConstants::NCh(); iChamber++) {
+//  for (Int_t iChamber=0; iChamber<AliMUONConstants::NCh(); iChamber++) {
+  for (Int_t iChamber=10; iChamber<14; iChamber++) {
 
     // Loop over cathods
     //for (Int_t cath=0; cath<2; cath++) {
@@ -933,3 +1207,5 @@ void AliMUONTest::DrawSegmentation(AliMUONGeometrySegmentation *seg)
   }
 }
 
+
+

diff --git a/MUON/AliMUONTest.h b/MUON/AliMUONTest.h
index db742b412ed82a0ab2e2bc4cd12be66cc23b80db..3566c1d311b25b240b86a5ac684b0ad0d59a03d1 100644 (file)
--- a/MUON/AliMUONTest.h
+++ b/MUON/AliMUONTest.h
@@ -75,6 +75,8 @@ class AliMUONTest : public  TObject
                                        Int_t chamberId, Int_t cathod);
     AliMUONGeometrySegmentation* CreateSlatSegmentation(
                                        Int_t chamberId, Int_t cathod);
+        AliMUONGeometrySegmentation* CreateTriggerSegmentation(
+                                       Int_t chamberId, Int_t cathod);
                                       
     // data members
     TCanvas* fCanvas; // The canvas for drawing                                       

diff --git a/MUON/MUONbaseLinkDef.h b/MUON/MUONbaseLinkDef.h
index c5bd2b780d55ceda88782dff097b7ad9690231d8..17d4edc62f44bab878c897b3b5781d1dedb953c5 100644 (file)
--- a/MUON/MUONbaseLinkDef.h
+++ b/MUON/MUONbaseLinkDef.h
@@ -22,6 +22,7 @@
 #pragma link C++ class AliMUONSegmentManuIndex+;
 #pragma link C++ class AliMUONSegmentPosition+;
 #pragma link C++ class AliMUONSt345SlatSegmentation+;
+#pragma link C++ class AliMUONTriggerSegmentation+;
 // geometry & response
 #pragma link C++ class AliMUONResponse+; 
 #pragma link C++ class AliMUONResponseV0+;

diff --git a/MUON/libMUONbase.pkg b/MUON/libMUONbase.pkg
index b3e793114f812601c4a288a58d28e9e948ae93cb..c7277ed56fae2e16d0f7ba8db9d4a15862aa5844 100644 (file)
--- a/MUON/libMUONbase.pkg
+++ b/MUON/libMUONbase.pkg
@@ -46,7 +46,8 @@ SRCS:=  AliMUON.cxx AliMUONv1.cxx \
         AliMUONDisplay.cxx \
         AliMUONRecoDisplay.cxx \
         AliMUONRecoCheck.cxx \
-        AliMUONSt345SlatSegmentation.cxx
+        AliMUONSt345SlatSegmentation.cxx \
+        AliMUONTriggerSegmentation.cxx
 
 HDRS:= $(SRCS:.cxx=.h)