]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MUON/AliMUONTriggerGeometryBuilder.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Replaced with new AliMUONSegFactory and AliMpSegFactory classes
[u/mrichter/AliRoot.git] / MUON / AliMUONTriggerGeometryBuilder.cxx
diff --git a/MUON/AliMUONTriggerGeometryBuilder.cxx b/MUON/AliMUONTriggerGeometryBuilder.cxx
index 9421e8f009662e01a3057018adc42d15741fab11..e073536db7b508b31b4116667966612997f2059c 100644 (file)
--- a/MUON/AliMUONTriggerGeometryBuilder.cxx
+++ b/MUON/AliMUONTriggerGeometryBuilder.cxx
@@ -28,7 +28,7 @@
 
 #include "AliMUONTriggerGeometryBuilder.h"
 #include "AliMUON.h"
-#include "AliMUONChamber.h"
+#include "AliMUONConstants.h"
 #include "AliMUONGeometryModule.h"
 #include "AliMUONGeometryEnvelopeStore.h"
 
@@ -36,11 +36,7 @@ ClassImp(AliMUONTriggerGeometryBuilder)
 
 //______________________________________________________________________________
 AliMUONTriggerGeometryBuilder::AliMUONTriggerGeometryBuilder(AliMUON* muon)
- : AliMUONVGeometryBuilder("trigger.dat",
-                           muon->Chamber(10).GetGeometry(), 
-                          muon->Chamber(11).GetGeometry(),
-                           muon->Chamber(12).GetGeometry(),
-                          muon->Chamber(13).GetGeometry()),
+ : AliMUONVGeometryBuilder(10, 11, 12, 13),
    fMUON(muon)
 {
 // Standard constructor
@@ -122,16 +118,18 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
     const Float_t kYMIN=34.;                              
     const Float_t kYMAX=51.;                              
 // inner/outer radius of flange between beam shield. and chambers (1/station)
-    const Float_t kRMIN[2]={50.,50.};
-    const Float_t kRMAX[2]={64.,68.};
+//    const Float_t kRMIN[2]={50.,50.};
+//    const Float_t kRMAX[2]={64.,68.};
 // z position of the middle of the gas gap in mother vol 
     const Float_t kZm=-3.6;
-    const Float_t kZp=+3.6;     
-
-    AliMUONChamber *iChamber, *iChamber1;
-    iChamber1 = &fMUON->Chamber(10);
-    Float_t zpos1=-iChamber1->Z(); 
+    const Float_t kZp=+3.6;
     
+    Float_t zpos1= AliMUONConstants::DefaultChamberZ(10); 
+    Double_t dstation =  ( ( - AliMUONConstants::DefaultChamberZ(11)) - 
+                           ( - AliMUONConstants::DefaultChamberZ(10)) ) /2.1;
+    Float_t par[3];
+    par[2] = dstation;
+
 // ratio of zpos1m/zpos1p and inverse for first plane
     Float_t zmp=(zpos1-3.6)/(zpos1+3.6);
     Float_t zpm=1./zmp;
@@ -142,12 +140,17 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
        for (Int_t iplane=0; iplane<2; iplane++) { // loop on detection planes
            
            Int_t iVolNum=1; // counter Volume Number
-           icount = Int_t(iplane*TMath::Power(2,0))+
-               Int_t(istation*TMath::Power(2,1));
+           icount = Int_t(iplane<<0)+Int_t(istation<<1);
            
-           iChamber = &fMUON->Chamber(10+icount);
-           Float_t zpos = - iChamber->Z();          
+           par[0] = AliMUONConstants::Rmin(5+istation); 
+           par[1] = AliMUONConstants::Rmax(5+istation);
+           Char_t volName[6];
+           sprintf(volName,"%s%d", "SC",11+icount);
+           gMC->Gsvolu(volName,"TUBE", idAir, par, 3);
+           GetGeometry(10+icount)->SetVolume(volName);
+           Float_t zpos =  AliMUONConstants::DefaultChamberZ(10+icount);            
            
+/* removed 03/18/05
 // Flange between beam shielding and RPC 
            tpar[0]= kRMIN[istation];
            tpar[1]= kRMAX[istation];
@@ -159,6 +162,7 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
            //gMC->Gspos(volFlange,1,"ALIC",0.,0.,zpos,0,"MANY");
            iChamber->GetGeometry()->GetEnvelopeStore()
              ->AddEnvelope(volFlange, 0, false, "MANY");
+*/
            
 // scaling factor
            Float_t zRatio = zpos / zpos1;
@@ -184,13 +188,12 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
 // chamber prototype
            tpar[0]= 0.;
            tpar[1]= 0.;
-           tpar[2]= 0.;            
-           char volAlu[5];     // Alu 
+           tpar[2]= 0.;            char volAlu[5];     // Alu 
            char volBak[5];     // Bakelite
            char volGaz[5];     // Gas streamer     
            sprintf(volAlu,"SC%dA",icount+1);
            sprintf(volBak,"SB%dA",icount+1);
-           sprintf(volGaz,"SG%dA",icount+1);
+           sprintf(volGaz,"S%dG",icount+11);
            gMC->Gsvolu(volAlu,"BOX",idAlu1,tpar,0);         // Al
            gMC->Gsvolu(volBak,"BOX",idtmed[1107],tpar,0);   // Bakelite
            gMC->Gsvolu(volGaz,"BOX",idtmed[1106],tpar,0);   // Gas streamer
@@ -213,10 +216,10 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
            dpar[0] = ((kXMAX-kXMED)/2.)*zRatio;
            dpar[1] = kYMIN * zRatio;
 
-           detElemId = (10+icount+1)*100+4;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[4], detElemId, true, TGeoTranslation(xEnv,yEnvM,kZm));
-           detElemId = (10+icount+1)*100+50+4;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[13], detElemId, true, TGeoTranslation(-xEnv,yEnvP,kZp));
+           detElemId = (10+icount+1)*100;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[4], detElemId, true, TGeoTranslation(xEnv,yEnvP,kZp));
+           detElemId = (10+icount+1)*100+9;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[13], detElemId, true, TGeoTranslation(-xEnv,yEnvM,kZm),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
 
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[4],iVolNum++,3, dpar);
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[13],iVolNum++,3, dpar);      
@@ -232,19 +235,19 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
            xpos = kXMIN/2. * zRatio;
            ypos = (kYMIN - kYMIN/4.) * zRatio;
 
-           detElemId = (10+icount+1)*100+3;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[3], detElemId, true, TGeoTranslation( xEnv,-yEnvP,kZp));
-           detElemId = (10+icount+1)*100+5;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[5], detElemId, true, TGeoTranslation( xEnv, yEnvP,kZp));
-           detElemId = (10+icount+1)*100+50+3;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[12], detElemId, true, TGeoTranslation(-xEnv,-yEnvM,kZm));
-           detElemId = (10+icount+1)*100+50+5;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[14], detElemId, true, TGeoTranslation(-xEnv, yEnvM,kZm));
+           detElemId = (10+icount+1)*100+17;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[3], detElemId, true, TGeoTranslation( xEnv,-yEnvM,kZm));            
+           detElemId = (10+icount+1)*100+1;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[5], detElemId, true, TGeoTranslation( xEnv, yEnvM,kZm));
+           detElemId = (10+icount+1)*100+10;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[12], detElemId, true, TGeoTranslation(-xEnv,-yEnvP,kZp),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
+           detElemId = (10+icount+1)*100+8;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[14], detElemId, true, TGeoTranslation(-xEnv, yEnvP,kZp),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
 
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
-           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
-           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);
+           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
+           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
 
 // chamber type C (note: same Z than type B)
            dpar[0] = (kXMAX/2)*zRatio;
@@ -254,9 +257,9 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
 
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[3],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[5],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
-           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(-xpos,-ypos,0.),3,dpar);
-           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(-xpos, ypos,0.),3,dpar);
-    
+           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[12],iVolNum++,TGeoTranslation(xpos,-ypos,0.),3,dpar);
+           GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[14],iVolNum++,TGeoTranslation(xpos, ypos,0.),3,dpar);
+
 // chamber type D, E and F (same size)
 // D       
            yEnvPsave = yEnvP;
@@ -266,14 +269,14 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
            dpar[0] = (kXMAX/2.)*zRatio;
            dpar[1] =  kYMIN*zRatio;
 
+           detElemId = (10+icount+1)*100+16;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[2], detElemId, true, TGeoTranslation(xEnv,-yEnvP,kZp));
            detElemId = (10+icount+1)*100+2;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[2], detElemId, true, TGeoTranslation(xEnv,-yEnvM,kZm));
-           detElemId = (10+icount+1)*100+6;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[6], detElemId, true, TGeoTranslation(xEnv, yEnvM,kZm));
-           detElemId = (10+icount+1)*100+50+2;
-            GetEnvelopes(10+icount)->AddEnvelope(volEnv[11], detElemId, true, TGeoTranslation(-xEnv,-yEnvP,kZp));
-           detElemId = (10+icount+1)*100+50+6;
-            GetEnvelopes(10+icount)->AddEnvelope(volEnv[15], detElemId, true, TGeoTranslation(-xEnv, yEnvP,kZp));
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[6], detElemId, true, TGeoTranslation(xEnv, yEnvP,kZp));
+           detElemId = (10+icount+1)*100+11;
+            GetEnvelopes(10+icount)->AddEnvelope(volEnv[11], detElemId, true, TGeoTranslation(-xEnv,-yEnvM,kZm),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
+           detElemId = (10+icount+1)*100+7;
+            GetEnvelopes(10+icount)->AddEnvelope(volEnv[15], detElemId, true, TGeoTranslation(-xEnv, yEnvM,kZm),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
 
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[2],iVolNum++,3, dpar);
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[6],iVolNum++,3, dpar);
@@ -286,14 +289,14 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
            yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
            yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
 
-           detElemId = (10+icount+1)*100+1;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[1], detElemId, true, TGeoTranslation(xEnv,-yEnvP,kZp));
-           detElemId = (10+icount+1)*100+7;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[7], detElemId, true, TGeoTranslation(xEnv, yEnvP,kZp));
-           detElemId = (10+icount+1)*100+50+1;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[10], detElemId, true, TGeoTranslation(-xEnv,-yEnvM,kZm));
-           detElemId = (10+icount+1)*100+50+7;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[16], detElemId, true, TGeoTranslation(-xEnv, yEnvM,kZm));
+           detElemId = (10+icount+1)*100+15;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[1], detElemId, true, TGeoTranslation(xEnv,-yEnvM,kZm));
+           detElemId = (10+icount+1)*100+3;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[7], detElemId, true, TGeoTranslation(xEnv, yEnvM,kZm));
+           detElemId = (10+icount+1)*100+12;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[10], detElemId, true, TGeoTranslation(-xEnv,-yEnvP,kZp),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
+           detElemId = (10+icount+1)*100+6;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[16], detElemId, true, TGeoTranslation(-xEnv, yEnvP,kZp),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
 
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[1],iVolNum++,3,dpar);
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[7],iVolNum++,3,dpar);
@@ -307,14 +310,14 @@ void AliMUONTriggerGeometryBuilder::CreateGeometry()
            yEnvP = (yEnvMsave + kYMIN * zRatio ) * zpm + kYMIN * zRatio;
            yEnvM = (yEnvPsave + kYMIN * zRatio ) * zmp + kYMIN * zRatio;
 
-           detElemId = (10+icount+1)*100;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[0], detElemId, true, TGeoTranslation(xEnv,-yEnvM,kZm));
-           detElemId = (10+icount+1)*100+8;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[8], detElemId, true, TGeoTranslation(xEnv, yEnvM,kZm));
-           detElemId = (10+icount+1)*100+50;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[9], detElemId, true, TGeoTranslation(-xEnv,-yEnvP,kZp));
-           detElemId = (10+icount+1)*100+50+8;
-           GetEnvelopes(10+icount)->AddEnvelope(volEnv[17], detElemId, true, TGeoTranslation(-xEnv, yEnvP,kZp));
+           detElemId = (10+icount+1)*100+14;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[0], detElemId, true, TGeoTranslation(xEnv,-yEnvP,kZp));
+           detElemId = (10+icount+1)*100+4;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[8], detElemId, true, TGeoTranslation(xEnv, yEnvP,kZp));
+           detElemId = (10+icount+1)*100+13;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[9], detElemId, true, TGeoTranslation(-xEnv,-yEnvM,kZm),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
+           detElemId = (10+icount+1)*100+5;
+           GetEnvelopes(10+icount)->AddEnvelope(volEnv[17], detElemId, true, TGeoTranslation(-xEnv, yEnvM,kZm),TGeoRotation("rot1",90.,180.,90.,90.,180.,0.));
            
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[0],iVolNum++,3,dpar);
            GetEnvelopes(10+icount)->AddEnvelopeConstituentParam(volAlu,volEnv[8],iVolNum++,3,dpar);
@@ -330,28 +333,19 @@ void AliMUONTriggerGeometryBuilder::SetTransformations()
 {
 // Defines the transformations for the trigger chambers.
 // ---
-    Double_t zpos1, zpos2;    
-    AliMUONChamber *iChamber1, *iChamber2;
 
-    iChamber1 = &fMUON->Chamber(10);
-    zpos1= - iChamber1->Z(); 
-    iChamber1->GetGeometry()
-       ->SetTranslation(TGeoTranslation(0., 0., zpos1));
-    
-    iChamber2 = &fMUON->Chamber(11);
-    zpos2 = - iChamber2->Z(); 
-    iChamber2->GetGeometry()
-       ->SetTranslation(TGeoTranslation(0., 0., zpos2));
-
-    iChamber1 = &fMUON->Chamber(12);
-    zpos1 = - iChamber1->Z(); 
-    iChamber1->GetGeometry()
-       ->SetTranslation(TGeoTranslation(0., 0., zpos1));
+    Double_t zpos1= AliMUONConstants::DefaultChamberZ(10); 
+    SetTranslation(10, TGeoTranslation(0., 0., zpos1));
     
-    iChamber2 = &fMUON->Chamber(13);
-    zpos2 = - iChamber2->Z(); 
-    iChamber2->GetGeometry()
-       ->SetTranslation(TGeoTranslation(0., 0., zpos2));
+    zpos1= AliMUONConstants::DefaultChamberZ(11); 
+    SetTranslation(11, TGeoTranslation(0., 0., zpos1));
+
+    zpos1= AliMUONConstants::DefaultChamberZ(12); 
+    SetTranslation(12, TGeoTranslation(0., 0., zpos1));
+
+    zpos1= AliMUONConstants::DefaultChamberZ(13); 
+    SetTranslation(13, TGeoTranslation(0., 0., zpos1));
+
 }
 
 //______________________________________________________________________________
@@ -360,9 +354,9 @@ void AliMUONTriggerGeometryBuilder::SetSensitiveVolumes()
 // Defines the sensitive volumes for trigger station chambers.
 // ---
 
-  GetGeometry(10)->SetSensitiveVolume("SG1A");
-  GetGeometry(11)->SetSensitiveVolume("SG2A");
-  GetGeometry(12)->SetSensitiveVolume("SG3A");
-  GetGeometry(13)->SetSensitiveVolume("SG4A");
+  GetGeometry(10)->SetSensitiveVolume("S11G");
+  GetGeometry(11)->SetSensitiveVolume("S12G");
+  GetGeometry(12)->SetSensitiveVolume("S13G");
+  GetGeometry(13)->SetSensitiveVolume("S14G");
 }