- gMC->Gsvolu("S05P","BOX",pcbMaterial,pcbpar,3);
- gMC->Gsvolu("S06P","BOX",pcbMaterial,pcbpar,3);
-
- // create the sensitive volumes,
- gMC->Gsvolu("S05G","BOX",sensMaterial,0,0);
- gMC->Gsvolu("S06G","BOX",sensMaterial,0,0);
-
-
- // create the vertical frame volume
-
- gMC->Gsvolu("S05V","BOX",vFrameMaterial,vFramepar,3);
- gMC->Gsvolu("S06V","BOX",vFrameMaterial,vFramepar,3);
-
- // create the horizontal frame volume
-
- gMC->Gsvolu("S05H","BOX",hFrameMaterial,hFramepar,3);
- gMC->Gsvolu("S06H","BOX",hFrameMaterial,hFramepar,3);
-
- // create the horizontal border volume
-
- gMC->Gsvolu("S05B","BOX",bFrameMaterial,bFramepar,3);
- gMC->Gsvolu("S06B","BOX",bFrameMaterial,bFramepar,3);
-
- index=0;
- for (i = 0; i<nSlats3; i++){
- sprintf(volNam5,"S05%d",i);
- sprintf(volNam6,"S06%d",i);
- Float_t xvFrame = (slatLength3[i] - vFrameLength)/2.;
- // position the vertical frames
- if (i!=1) {
- gMC->Gspos("S05V",2*i-1,volNam5, xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S05V",2*i ,volNam5,-xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S06V",2*i-1,volNam6, xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S06V",2*i ,volNam6,-xvFrame, 0., 0. , 0, "ONLY");
- }
- // position the panels and the insulating material
- for (j=0; j<nPCB3[i]; j++){
- index++;
- Float_t xx = sensLength * (-nPCB3[i]/2.+j+.5);
-
- Float_t zPanel = spar[2] - panelpar[2];
- gMC->Gspos("S05C",2*index-1,volNam5, xx, 0., zPanel , 0, "ONLY");
- gMC->Gspos("S05C",2*index ,volNam5, xx, 0.,-zPanel , 0, "ONLY");
- gMC->Gspos("S06C",2*index-1,volNam6, xx, 0., zPanel , 0, "ONLY");
- gMC->Gspos("S06C",2*index ,volNam6, xx, 0.,-zPanel , 0, "ONLY");
-
- gMC->Gspos("S05I",index,volNam5, xx, 0., 0 , 0, "ONLY");
- gMC->Gspos("S06I",index,volNam6, xx, 0., 0 , 0, "ONLY");
- }
- }
-
- // position the rohacell volume inside the panel volume
- gMC->Gspos("S05R",1,"S05C",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S06R",1,"S06C",0.,0.,0.,0,"ONLY");
-
- // position the PCB volume inside the insulating material volume
- gMC->Gspos("S05P",1,"S05I",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S06P",1,"S06I",0.,0.,0.,0,"ONLY");
- // position the horizontal frame volume inside the PCB volume
- gMC->Gspos("S05H",1,"S05P",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S06H",1,"S06P",0.,0.,0.,0,"ONLY");
- // position the sensitive volume inside the horizontal frame volume
- gMC->Gsposp("S05G",1,"S05H",0.,0.,0.,0,"ONLY",senspar,3);
- gMC->Gsposp("S06G",1,"S06H",0.,0.,0.,0,"ONLY",senspar,3);
- // position the border volumes inside the PCB volume
- Float_t yborder = ( pcbHeight - bFrameHeight ) / 2.;
- gMC->Gspos("S05B",1,"S05P",0., yborder,0.,0,"ONLY");
- gMC->Gspos("S05B",2,"S05P",0.,-yborder,0.,0,"ONLY");
- gMC->Gspos("S06B",1,"S06P",0., yborder,0.,0,"ONLY");
- gMC->Gspos("S06B",2,"S06P",0.,-yborder,0.,0,"ONLY");
-
- // create the NULOC volume and position it in the horizontal frame
-
- gMC->Gsvolu("S05N","BOX",nulocMaterial,nulocpar,3);
- gMC->Gsvolu("S06N","BOX",nulocMaterial,nulocpar,3);
- index = 0;
- for (xx = -xxmax; xx<=xxmax; xx+=3*nulocLength) {
- index++;
- gMC->Gspos("S05N",2*index-1,"S05B", xx, 0.,-bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S05N",2*index ,"S05B", xx, 0., bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S06N",2*index-1,"S06B", xx, 0.,-bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S06N",2*index ,"S06B", xx, 0., bFrameWidth/4., 0, "ONLY");
- }
-
- // position the volumes approximating the circular section of the pipe
- Float_t yoffs = sensHeight/2. - yOverlap;
- Float_t epsilon = 0.001;
- Int_t ndiv=6;
- Float_t divpar[3];
- Double_t dydiv= sensHeight/ndiv;
- Double_t ydiv = yoffs -dydiv - yOverlap/2.;
- Int_t imax=0;
- // for (Int_t islat=0; islat<nSlats3; islat++) imax += nPCB3[islat];
- imax = 1;
- Float_t rmin = 35.;
- Float_t z1 = spar[2], z2=2*spar[2]*1.01;
- for (Int_t idiv=0;idiv<ndiv; idiv++){
- ydiv+= dydiv;
- Float_t xdiv = 0.;
- if (ydiv<rmin) xdiv= rmin * TMath::Sin( TMath::ACos(ydiv/rmin) );
- divpar[0] = (pcbLength-xdiv)/2.;
- divpar[1] = dydiv/2. - epsilon;
- divpar[2] = sensWidth/2.;
- Float_t xvol=(pcbLength+xdiv)/2.+1.999;
- Float_t yvol=ydiv + dydiv/2.;
- gMC->Gsposp("S05G",imax+4*idiv+1,"C05M", xvol, yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S06G",imax+4*idiv+1,"C06M", xvol, yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S05G",imax+4*idiv+2,"C05M", xvol,-yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S06G",imax+4*idiv+2,"C06M", xvol,-yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S05G",imax+4*idiv+3,"C05M",-xvol, yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S06G",imax+4*idiv+3,"C06M",-xvol, yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S05G",imax+4*idiv+4,"C05M",-xvol,-yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S06G",imax+4*idiv+4,"C06M",-xvol,-yvol, z1-z2, 0, "ONLY",divpar,3);
- }
- }
-
-
- if (stations[3]) {
-
-//********************************************************************
-// Station 4 **
-//********************************************************************
- // indices 1 and 2 for first and second chambers in the station
- // iChamber (first chamber) kept for other quanties than Z,
- // assumed to be the same in both chambers
- iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[6];
- iChamber2 =(AliMUONChamber*) (*fChambers)[7];
- zpos1=iChamber1->Z();
- zpos2=iChamber2->Z();
- dstation = zpos2 - zpos1;
-// zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; // not used any more
-
-//
-// Mother volume
- tpar[0] = iChamber->RInner()-dframep;
- tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
- tpar[2] = 3.252;
-
- gMC->Gsvolu("C07M", "TUBE", idAir, tpar, 3);
- gMC->Gsvolu("C08M", "TUBE", idAir, tpar, 3);
- gMC->Gspos("C07M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
- gMC->Gspos("C08M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
-
-
- const Int_t nSlats4 = 6; // number of slats per quadrant
- const Int_t nPCB4[nSlats4] = {4,4,5,5,4,3}; // n PCB per slat
- const Float_t xpos4[nSlats4] = {37.5, 40., 0., 0., 0., 0.};
- Float_t slatLength4[nSlats4];
-
- // create and position the slat (mother) volumes
-
- char volNam7[5];
- char volNam8[5];
- Float_t xSlat4;
- Float_t ySlat4;
-
- for (i = 0; i<nSlats4; i++){
- slatLength4[i] = pcbLength * nPCB4[i] + 2. * dSlatLength;
- xSlat4 = slatLength4[i]/2. - vFrameLength/2. + xpos4[i];
- if (i==1) slatLength4[i] -= 2. *dSlatLength; // frame out in PCB with circular border
- ySlat4 = sensHeight * i - yOverlap *i;
-
- spar[0] = slatLength4[i]/2.;
- spar[1] = slatHeight/2.;
- spar[2] = slatWidth/2.*1.01;
- Float_t dzCh4=spar[2]*1.01;
- // zSlat to be checked (odd downstream or upstream?)
- Float_t zSlat = (i%2 ==0)? spar[2] : -spar[2];
- sprintf(volNam7,"S07%d",i);
- gMC->Gsvolu(volNam7,"BOX",slatMaterial,spar,3);
- gMC->Gspos(volNam7, i*4+1,"C07M", xSlat4, ySlat4, zSlat+2.*dzCh4, 0, "ONLY");
- gMC->Gspos(volNam7, i*4+2,"C07M",-xSlat4, ySlat4, zSlat-2.*dzCh4, 0, "ONLY");
- if (i>0) {
- gMC->Gspos(volNam7, i*4+3,"C07M", xSlat4,-ySlat4, zSlat+2.*dzCh4, 0, "ONLY");
- gMC->Gspos(volNam7, i*4+4,"C07M",-xSlat4,-ySlat4, zSlat-2.*dzCh4, 0, "ONLY");
- }
- sprintf(volNam8,"S08%d",i);
- gMC->Gsvolu(volNam8,"BOX",slatMaterial,spar,3);
- gMC->Gspos(volNam8, i*4+1,"C08M", xSlat4, ySlat4, zSlat+2.*dzCh4, 0, "ONLY");
- gMC->Gspos(volNam8, i*4+2,"C08M",-xSlat4, ySlat4, zSlat-2.*dzCh4, 0, "ONLY");
- if (i>0) {
- gMC->Gspos(volNam8, i*4+3,"C08M", xSlat4,-ySlat4, zSlat+2.*dzCh4, 0, "ONLY");
- gMC->Gspos(volNam8, i*4+4,"C08M",-xSlat4,-ySlat4, zSlat-2.*dzCh4, 0, "ONLY");
- }
- }
-
-
- // create the panel volume
-
- gMC->Gsvolu("S07C","BOX",panelMaterial,panelpar,3);
- gMC->Gsvolu("S08C","BOX",panelMaterial,panelpar,3);
-
- // create the rohacell volume
-
- gMC->Gsvolu("S07R","BOX",rohaMaterial,rohapar,3);
- gMC->Gsvolu("S08R","BOX",rohaMaterial,rohapar,3);
-
- // create the insulating material volume
-
- gMC->Gsvolu("S07I","BOX",insuMaterial,insupar,3);
- gMC->Gsvolu("S08I","BOX",insuMaterial,insupar,3);
-
- // create the PCB volume
-
- gMC->Gsvolu("S07P","BOX",pcbMaterial,pcbpar,3);
- gMC->Gsvolu("S08P","BOX",pcbMaterial,pcbpar,3);
-
- // create the sensitive volumes,
-
- gMC->Gsvolu("S07G","BOX",sensMaterial,0,0);
- gMC->Gsvolu("S08G","BOX",sensMaterial,0,0);
-
- // create the vertical frame volume
-
- gMC->Gsvolu("S07V","BOX",vFrameMaterial,vFramepar,3);
- gMC->Gsvolu("S08V","BOX",vFrameMaterial,vFramepar,3);
-
- // create the horizontal frame volume
-
- gMC->Gsvolu("S07H","BOX",hFrameMaterial,hFramepar,3);
- gMC->Gsvolu("S08H","BOX",hFrameMaterial,hFramepar,3);
-
- // create the horizontal border volume
-
- gMC->Gsvolu("S07B","BOX",bFrameMaterial,bFramepar,3);
- gMC->Gsvolu("S08B","BOX",bFrameMaterial,bFramepar,3);
-
- index=0;
- for (i = 0; i<nSlats4; i++){
- sprintf(volNam7,"S07%d",i);
- sprintf(volNam8,"S08%d",i);
- Float_t xvFrame = (slatLength4[i] - vFrameLength)/2.;
- // position the vertical frames
- if (i!=1) {
- gMC->Gspos("S07V",2*i-1,volNam7, xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S07V",2*i ,volNam7,-xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S08V",2*i-1,volNam8, xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S08V",2*i ,volNam8,-xvFrame, 0., 0. , 0, "ONLY");
- }
- // position the panels and the insulating material
- for (j=0; j<nPCB4[i]; j++){
- index++;
- Float_t xx = sensLength * (-nPCB4[i]/2.+j+.5);
-
- Float_t zPanel = spar[2] - panelpar[2];
- gMC->Gspos("S07C",2*index-1,volNam7, xx, 0., zPanel , 0, "ONLY");
- gMC->Gspos("S07C",2*index ,volNam7, xx, 0.,-zPanel , 0, "ONLY");
- gMC->Gspos("S08C",2*index-1,volNam8, xx, 0., zPanel , 0, "ONLY");
- gMC->Gspos("S08C",2*index ,volNam8, xx, 0.,-zPanel , 0, "ONLY");
-
- gMC->Gspos("S07I",index,volNam7, xx, 0., 0 , 0, "ONLY");
- gMC->Gspos("S08I",index,volNam8, xx, 0., 0 , 0, "ONLY");
- }
- }
-
- // position the rohacell volume inside the panel volume
- gMC->Gspos("S07R",1,"S07C",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S08R",1,"S08C",0.,0.,0.,0,"ONLY");
-
- // position the PCB volume inside the insulating material volume
- gMC->Gspos("S07P",1,"S07I",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S08P",1,"S08I",0.,0.,0.,0,"ONLY");
- // position the horizontal frame volume inside the PCB volume
- gMC->Gspos("S07H",1,"S07P",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S08H",1,"S08P",0.,0.,0.,0,"ONLY");
- // position the sensitive volume inside the horizontal frame volume
- gMC->Gsposp("S07G",1,"S07H",0.,0.,0.,0,"ONLY",senspar,3);
- gMC->Gsposp("S08G",1,"S08H",0.,0.,0.,0,"ONLY",senspar,3);
- // position the border volumes inside the PCB volume
- Float_t yborder = ( pcbHeight - bFrameHeight ) / 2.;
- gMC->Gspos("S07B",1,"S07P",0., yborder,0.,0,"ONLY");
- gMC->Gspos("S07B",2,"S07P",0.,-yborder,0.,0,"ONLY");
- gMC->Gspos("S08B",1,"S08P",0., yborder,0.,0,"ONLY");
- gMC->Gspos("S08B",2,"S08P",0.,-yborder,0.,0,"ONLY");
-
- // create the NULOC volume and position it in the horizontal frame
-
- gMC->Gsvolu("S07N","BOX",nulocMaterial,nulocpar,3);
- gMC->Gsvolu("S08N","BOX",nulocMaterial,nulocpar,3);
- index = 0;
- for (xx = -xxmax; xx<=xxmax; xx+=3*nulocLength) {
- index++;
- gMC->Gspos("S07N",2*index-1,"S07B", xx, 0.,-bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S07N",2*index ,"S07B", xx, 0., bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S08N",2*index-1,"S08B", xx, 0.,-bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S08N",2*index ,"S08B", xx, 0., bFrameWidth/4., 0, "ONLY");
- }
-
- // position the volumes approximating the circular section of the pipe
- Float_t yoffs = sensHeight/2. - yOverlap/2.;
- Float_t epsilon = 0.001;
- Int_t ndiv=6;
- Float_t divpar[3];
- Double_t dydiv= sensHeight/ndiv;
- Double_t ydiv = yoffs -dydiv - yOverlap/2.;
- Int_t imax=0;
- // for (Int_t islat=0; islat<nSlats3; islat++) imax += nPCB3[islat];
- imax = 1;
- Float_t rmin = 40.;
- Float_t z1 = -spar[2], z2=2*spar[2]*1.01;
- for (Int_t idiv=0;idiv<ndiv; idiv++){
- ydiv+= dydiv;
- Float_t xdiv = 0.;
- if (ydiv<rmin) xdiv= rmin * TMath::Sin( TMath::ACos(ydiv/rmin) );
- divpar[0] = (pcbLength-xdiv)/2.;
- divpar[1] = dydiv/2. - epsilon;
- divpar[2] = sensWidth/2.;
- Float_t xvol=(pcbLength+xdiv)/2.+1.999;
- Float_t yvol=ydiv + dydiv/2.;
- gMC->Gsposp("S07G",imax+4*idiv+1,"C07M", xvol, yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S08G",imax+4*idiv+1,"C08M", xvol, yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S07G",imax+4*idiv+2,"C07M", xvol,-yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S08G",imax+4*idiv+2,"C08M", xvol,-yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S07G",imax+4*idiv+3,"C07M",-xvol, yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S08G",imax+4*idiv+3,"C08M",-xvol, yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S07G",imax+4*idiv+4,"C07M",-xvol,-yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S08G",imax+4*idiv+4,"C08M",-xvol,-yvol, z1-z2, 0, "ONLY",divpar,3);
- }
-
-
-
-
-
- }
-
- if (stations[4]) {
-
-
-//********************************************************************
-// Station 5 **
-//********************************************************************
- // indices 1 and 2 for first and second chambers in the station
- // iChamber (first chamber) kept for other quanties than Z,
- // assumed to be the same in both chambers
- iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[8];
- iChamber2 =(AliMUONChamber*) (*fChambers)[9];
- zpos1=iChamber1->Z();
- zpos2=iChamber2->Z();
- dstation = zpos2 - zpos1;
-// zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; // not used any more
-
-//
-// Mother volume
- tpar[0] = iChamber->RInner()-dframep;
- tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi);
- tpar[2] = dstation/5.;
-
- gMC->Gsvolu("C09M", "TUBE", idAir, tpar, 3);
- gMC->Gsvolu("C10M", "TUBE", idAir, tpar, 3);
- gMC->Gspos("C09M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY");
- gMC->Gspos("C10M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY");
-
-
- const Int_t nSlats5 = 7; // number of slats per quadrant
- const Int_t nPCB5[nSlats5] = {5,5,6,6,5,4,3}; // n PCB per slat
- const Float_t xpos5[nSlats5] = {37.5, 40., 0., 0., 0., 0., 0.};
- Float_t slatLength5[nSlats5];
- char volNam9[5];
- char volNam10[5];
- Float_t xSlat5;
- Float_t ySlat5;
-
- for (i = 0; i<nSlats5; i++){
- slatLength5[i] = pcbLength * nPCB5[i] + 2. * dSlatLength;
- xSlat5 = slatLength5[i]/2. - vFrameLength/2. +xpos5[i];
- if (i==1) slatLength5[i] -= 2. *dSlatLength; // frame out in PCB with circular border
- ySlat5 = sensHeight * i - yOverlap * i;
- spar[0] = slatLength5[i]/2.;
- spar[1] = slatHeight/2.;
- spar[2] = slatWidth/2. * 1.01;
- Float_t dzCh5=spar[2]*1.01;
- // zSlat to be checked (odd downstream or upstream?)
- Float_t zSlat = (i%2 ==0)? -spar[2] : spar[2];
- sprintf(volNam9,"S09%d",i);
- gMC->Gsvolu(volNam9,"BOX",slatMaterial,spar,3);
- gMC->Gspos(volNam9, i*4+1,"C09M", xSlat5, ySlat5, zSlat+2.*dzCh5, 0, "ONLY");
- gMC->Gspos(volNam9, i*4+2,"C09M",-xSlat5, ySlat5, zSlat-2.*dzCh5, 0, "ONLY");
- if (i>0) {
- gMC->Gspos(volNam9, i*4+3,"C09M", xSlat5,-ySlat5, zSlat+2.*dzCh5, 0, "ONLY");
- gMC->Gspos(volNam9, i*4+4,"C09M",-xSlat5,-ySlat5, zSlat-2.*dzCh5, 0, "ONLY");
- }
- sprintf(volNam10,"S10%d",i);
- gMC->Gsvolu(volNam10,"BOX",slatMaterial,spar,3);
- gMC->Gspos(volNam10, i*4+1,"C10M", xSlat5, ySlat5, zSlat+2.*dzCh5, 0, "ONLY");
- gMC->Gspos(volNam10, i*4+2,"C10M",-xSlat5, ySlat5, zSlat-2.*dzCh5, 0, "ONLY");
- if (i>0) {
- gMC->Gspos(volNam10, i*4+3,"C10M", xSlat5,-ySlat5, zSlat+2.*dzCh5, 0, "ONLY");
- gMC->Gspos(volNam10, i*4+4,"C10M",-xSlat5,-ySlat5, zSlat-2.*dzCh5, 0, "ONLY");
- }
- }
-
- // create the panel volume
-
- gMC->Gsvolu("S09C","BOX",panelMaterial,panelpar,3);
- gMC->Gsvolu("S10C","BOX",panelMaterial,panelpar,3);
-
- // create the rohacell volume
-
- gMC->Gsvolu("S09R","BOX",rohaMaterial,rohapar,3);
- gMC->Gsvolu("S10R","BOX",rohaMaterial,rohapar,3);
-
- // create the insulating material volume
-
- gMC->Gsvolu("S09I","BOX",insuMaterial,insupar,3);
- gMC->Gsvolu("S10I","BOX",insuMaterial,insupar,3);
-
- // create the PCB volume
-
- gMC->Gsvolu("S09P","BOX",pcbMaterial,pcbpar,3);
- gMC->Gsvolu("S10P","BOX",pcbMaterial,pcbpar,3);
-
- // create the sensitive volumes,
-
- gMC->Gsvolu("S09G","BOX",sensMaterial,0,0);
- gMC->Gsvolu("S10G","BOX",sensMaterial,0,0);
-
- // create the vertical frame volume
-
- gMC->Gsvolu("S09V","BOX",vFrameMaterial,vFramepar,3);
- gMC->Gsvolu("S10V","BOX",vFrameMaterial,vFramepar,3);
-
- // create the horizontal frame volume
-
- gMC->Gsvolu("S09H","BOX",hFrameMaterial,hFramepar,3);
- gMC->Gsvolu("S10H","BOX",hFrameMaterial,hFramepar,3);
-
- // create the horizontal border volume
-
- gMC->Gsvolu("S09B","BOX",bFrameMaterial,bFramepar,3);
- gMC->Gsvolu("S10B","BOX",bFrameMaterial,bFramepar,3);
-
- index=0;
- for (i = 0; i<nSlats5; i++){
- sprintf(volNam9,"S09%d",i);
- sprintf(volNam10,"S10%d",i);
- Float_t xvFrame = (slatLength5[i] - vFrameLength)/2.;
- // position the vertical frames
- if (i!=1) {
- gMC->Gspos("S09V",2*i-1,volNam9, xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S09V",2*i ,volNam9,-xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S10V",2*i-1,volNam10, xvFrame, 0., 0. , 0, "ONLY");
- gMC->Gspos("S10V",2*i ,volNam10,-xvFrame, 0., 0. , 0, "ONLY");
- }
-
- // position the panels and the insulating material
- for (j=0; j<nPCB5[i]; j++){
- index++;
- Float_t xx = sensLength * (-nPCB5[i]/2.+j+.5);
-
- Float_t zPanel = spar[2] - panelpar[2];
- gMC->Gspos("S09C",2*index-1,volNam9, xx, 0., zPanel , 0, "ONLY");
- gMC->Gspos("S09C",2*index ,volNam9, xx, 0.,-zPanel , 0, "ONLY");
- gMC->Gspos("S10C",2*index-1,volNam10, xx, 0., zPanel , 0, "ONLY");
- gMC->Gspos("S10C",2*index ,volNam10, xx, 0.,-zPanel , 0, "ONLY");
-
- gMC->Gspos("S09I",index,volNam9, xx, 0., 0 , 0, "ONLY");
- gMC->Gspos("S10I",index,volNam10, xx, 0., 0 , 0, "ONLY");
- }
- }
-
- // position the rohacell volume inside the panel volume
- gMC->Gspos("S09R",1,"S09C",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S10R",1,"S10C",0.,0.,0.,0,"ONLY");
-
- // position the PCB volume inside the insulating material volume
- gMC->Gspos("S09P",1,"S09I",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S10P",1,"S10I",0.,0.,0.,0,"ONLY");
- // position the horizontal frame volume inside the PCB volume
- gMC->Gspos("S09H",1,"S09P",0.,0.,0.,0,"ONLY");
- gMC->Gspos("S10H",1,"S10P",0.,0.,0.,0,"ONLY");
- // position the sensitive volume inside the horizontal frame volume
- gMC->Gsposp("S09G",1,"S09H",0.,0.,0.,0,"ONLY",senspar,3);
- gMC->Gsposp("S10G",1,"S10H",0.,0.,0.,0,"ONLY",senspar,3);
- // position the border volumes inside the PCB volume
- Float_t yborder = ( pcbHeight - bFrameHeight ) / 2.;
- gMC->Gspos("S09B",1,"S09P",0., yborder,0.,0,"ONLY");
- gMC->Gspos("S09B",2,"S09P",0.,-yborder,0.,0,"ONLY");
- gMC->Gspos("S10B",1,"S10P",0., yborder,0.,0,"ONLY");
- gMC->Gspos("S10B",2,"S10P",0.,-yborder,0.,0,"ONLY");
-
- // create the NULOC volume and position it in the horizontal frame
-
- gMC->Gsvolu("S09N","BOX",nulocMaterial,nulocpar,3);
- gMC->Gsvolu("S10N","BOX",nulocMaterial,nulocpar,3);
- index = 0;
- for (xx = -xxmax; xx<=xxmax; xx+=3*nulocLength) {
- index++;
- gMC->Gspos("S09N",2*index-1,"S09B", xx, 0.,-bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S09N",2*index ,"S09B", xx, 0., bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S10N",2*index-1,"S10B", xx, 0.,-bFrameWidth/4., 0, "ONLY");
- gMC->Gspos("S10N",2*index ,"S10B", xx, 0., bFrameWidth/4., 0, "ONLY");
- }
- // position the volumes approximating the circular section of the pipe
- Float_t yoffs = sensHeight/2. - yOverlap/2.;
- Float_t epsilon = 0.001;
- Int_t ndiv=6;
- Float_t divpar[3];
- Double_t dydiv= sensHeight/ndiv;
- Double_t ydiv = yoffs -dydiv - yOverlap/2.;
- Int_t imax=0;
- // for (Int_t islat=0; islat<nSlats3; islat++) imax += nPCB3[islat];
- imax = 1;
- Float_t rmin = 40.;
- Float_t z1 = spar[2], z2=2*spar[2]*1.01;
- for (Int_t idiv=0;idiv<ndiv; idiv++){
- ydiv+= dydiv;
- Float_t xdiv = 0.;
- if (ydiv<rmin) xdiv= rmin * TMath::Sin( TMath::ACos(ydiv/rmin) );
- divpar[0] = (pcbLength-xdiv)/2.;
- divpar[1] = dydiv/2. - epsilon;
- divpar[2] = sensWidth/2.;
- Float_t xvol=(pcbLength+xdiv)/2. + 1.999;
- Float_t yvol=ydiv + dydiv/2.;
- gMC->Gsposp("S09G",imax+4*idiv+1,"C09M", xvol, yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S10G",imax+4*idiv+1,"C10M", xvol, yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S09G",imax+4*idiv+2,"C09M", xvol,-yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S10G",imax+4*idiv+2,"C10M", xvol,-yvol, z1+z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S09G",imax+4*idiv+3,"C09M",-xvol, yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S10G",imax+4*idiv+3,"C10M",-xvol, yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S09G",imax+4*idiv+4,"C09M",-xvol,-yvol, z1-z2, 0, "ONLY",divpar,3);
- gMC->Gsposp("S10G",imax+4*idiv+4,"C10M",-xvol,-yvol, z1-z2, 0, "ONLY",divpar,3);
- }
-
- }
-
-
-///////////////////////////////////////
-// GEOMETRY FOR THE TRIGGER CHAMBERS //
-///////////////////////////////////////
-
-// 03/00 P. Dupieux : introduce a slighly more realistic
-// geom. of the trigger readout planes with
-// 2 Zpos per trigger plane (alternate
-// between left and right of the trigger)
-
-// Parameters of the Trigger Chambers
-
-
- const Float_t kXMC1MIN=34.;
- const Float_t kXMC1MED=51.;
- const Float_t kXMC1MAX=272.;
- const Float_t kYMC1MIN=34.;
- const Float_t kYMC1MAX=51.;
- const Float_t kRMIN1=50.;
- const Float_t kRMAX1=62.;
- const Float_t kRMIN2=50.;
- const Float_t kRMAX2=66.;
-
-// zposition of the middle of the gas gap in mother vol
- const Float_t kZMCm=-3.6;
- const Float_t kZMCp=+3.6;
-
-
-// TRIGGER STATION 1 - TRIGGER STATION 1 - TRIGGER STATION 1
-
- // iChamber 1 and 2 for first and second chambers in the station
- // iChamber (first chamber) kept for other quanties than Z,
- // assumed to be the same in both chambers
- iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[10];
- iChamber2 =(AliMUONChamber*) (*fChambers)[11];
-
- // 03/00
- // zpos1 and zpos2 are now the middle of the first and second
- // plane of station 1 :
- // zpos1=(16075+15995)/2=16035 mm, thick/2=40 mm
- // zpos2=(16225+16145)/2=16185 mm, thick/2=40 mm
- //
- // zpos1m=15999 mm , zpos1p=16071 mm (middles of gas gaps)
- // zpos2m=16149 mm , zpos2p=16221 mm (middles of gas gaps)
- // rem : the total thickness accounts for 1 mm of al on both
- // side of the RPCs (see zpos1 and zpos2), as previously
-
- zpos1=iChamber1->Z();
- zpos2=iChamber2->Z();
-
-
-// Mother volume definition
- tpar[0] = iChamber->RInner();
- tpar[1] = iChamber->ROuter();
- tpar[2] = 4.0;
- gMC->Gsvolu("CM11", "TUBE", idAir, tpar, 3);
- gMC->Gsvolu("CM12", "TUBE", idAir, tpar, 3);
-
-// Definition of the flange between the beam shielding and the RPC
- tpar[0]= kRMIN1;
- tpar[1]= kRMAX1;
- tpar[2]= 4.0;
-
- gMC->Gsvolu("CF1A", "TUBE", idAlu1, tpar, 3); //Al
- gMC->Gspos("CF1A", 1, "CM11", 0., 0., 0., 0, "MANY");
- gMC->Gspos("CF1A", 2, "CM12", 0., 0., 0., 0, "MANY");
-
-
-// FIRST PLANE OF STATION 1
-
-// ratios of zpos1m/zpos1p and inverse for first plane
- Float_t zmp=(zpos1-3.6)/(zpos1+3.6);
- Float_t zpm=1./zmp;
-
-
-// Definition of prototype for chambers in the first plane
-
- tpar[0]= 0.;
- tpar[1]= 0.;
- tpar[2]= 0.;
-
- gMC->Gsvolu("CC1A", "BOX ", idAlu1, tpar, 0); //Al
- gMC->Gsvolu("CB1A", "BOX ", idtmed[1107], tpar, 0); //Bakelite
- gMC->Gsvolu("CG1A", "BOX ", idtmed[1106], tpar, 0); //Gas streamer
-
-// chamber type A
- tpar[0] = -1.;
- tpar[1] = -1.;
-
- const Float_t kXMC1A=kXMC1MED+(kXMC1MAX-kXMC1MED)/2.;
- const Float_t kYMC1Am=0.;
- const Float_t kYMC1Ap=0.;
-
- tpar[2] = 0.1;
- gMC->Gsposp("CG1A", 1, "CB1A", 0., 0., 0., 0, "ONLY",tpar,3);
- tpar[2] = 0.3;
- gMC->Gsposp("CB1A", 1, "CC1A", 0., 0., 0., 0, "ONLY",tpar,3);
-
- tpar[2] = 0.4;
- tpar[0] = (kXMC1MAX-kXMC1MED)/2.;
- tpar[1] = kYMC1MIN;
-
- gMC->Gsposp("CC1A", 1, "CM11",kXMC1A,kYMC1Am,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 2, "CM11",-kXMC1A,kYMC1Ap,kZMCp, 0, "ONLY", tpar, 3);
-
-// chamber type B
- Float_t tpar1save=tpar[1];
- Float_t y1msave=kYMC1Am;
- Float_t y1psave=kYMC1Ap;
-
- tpar[0] = (kXMC1MAX-kXMC1MIN)/2.;
- tpar[1] = (kYMC1MAX-kYMC1MIN)/2.;
-
- const Float_t kXMC1B=kXMC1MIN+tpar[0];
- const Float_t kYMC1Bp=(y1msave+tpar1save)*zpm+tpar[1];
- const Float_t kYMC1Bm=(y1psave+tpar1save)*zmp+tpar[1];
-
- gMC->Gsposp("CC1A", 3, "CM11",kXMC1B,kYMC1Bp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 4, "CM11",-kXMC1B,kYMC1Bm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 5, "CM11",kXMC1B,-kYMC1Bp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 6, "CM11",-kXMC1B,-kYMC1Bm,kZMCm, 0, "ONLY", tpar, 3);
-
-// chamber type C (end of type B !!)
- tpar1save=tpar[1];
- y1msave=kYMC1Bm;
- y1psave=kYMC1Bp;
-
- tpar[0] = kXMC1MAX/2;
- tpar[1] = kYMC1MAX/2;
-
- const Float_t kXMC1C=tpar[0];
-// warning : same Z than type B
- const Float_t kYMC1Cp=(y1psave+tpar1save)*1.+tpar[1];
- const Float_t kYMC1Cm=(y1msave+tpar1save)*1.+tpar[1];
-
- gMC->Gsposp("CC1A", 7, "CM11",kXMC1C,kYMC1Cp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 8, "CM11",-kXMC1C,kYMC1Cm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 9, "CM11",kXMC1C,-kYMC1Cp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 10, "CM11",-kXMC1C,-kYMC1Cm,kZMCm, 0, "ONLY", tpar, 3);
-
-// chamber type D, E and F (same size)
- tpar1save=tpar[1];
- y1msave=kYMC1Cm;
- y1psave=kYMC1Cp;
-
- tpar[0] = kXMC1MAX/2.;
- tpar[1] = kYMC1MIN;
-
- const Float_t kXMC1D=tpar[0];
- const Float_t kYMC1Dp=(y1msave+tpar1save)*zpm+tpar[1];
- const Float_t kYMC1Dm=(y1psave+tpar1save)*zmp+tpar[1];
-
- gMC->Gsposp("CC1A", 11, "CM11",kXMC1D,kYMC1Dm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 12, "CM11",-kXMC1D,kYMC1Dp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 13, "CM11",kXMC1D,-kYMC1Dm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 14, "CM11",-kXMC1D,-kYMC1Dp,kZMCp, 0, "ONLY", tpar, 3);
-
-
- tpar1save=tpar[1];
- y1msave=kYMC1Dm;
- y1psave=kYMC1Dp;
- const Float_t kYMC1Ep=(y1msave+tpar1save)*zpm+tpar[1];
- const Float_t kYMC1Em=(y1psave+tpar1save)*zmp+tpar[1];
-
- gMC->Gsposp("CC1A", 15, "CM11",kXMC1D,kYMC1Ep,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 16, "CM11",-kXMC1D,kYMC1Em,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 17, "CM11",kXMC1D,-kYMC1Ep,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 18, "CM11",-kXMC1D,-kYMC1Em,kZMCm, 0, "ONLY", tpar, 3);
-
- tpar1save=tpar[1];
- y1msave=kYMC1Em;
- y1psave=kYMC1Ep;
- const Float_t kYMC1Fp=(y1msave+tpar1save)*zpm+tpar[1];
- const Float_t kYMC1Fm=(y1psave+tpar1save)*zmp+tpar[1];
-
- gMC->Gsposp("CC1A", 19, "CM11",kXMC1D,kYMC1Fm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 20, "CM11",-kXMC1D,kYMC1Fp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 21, "CM11",kXMC1D,-kYMC1Fm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC1A", 22, "CM11",-kXMC1D,-kYMC1Fp,kZMCp, 0, "ONLY", tpar, 3);
-
-// Positioning first plane in ALICE
- gMC->Gspos("CM11", 1, "ALIC", 0., 0., zpos1, 0, "ONLY");
-
-// End of geometry definition for the first plane of station 1
-
-
-
-// SECOND PLANE OF STATION 1 : proj ratio = zpos2/zpos1
-
- const Float_t kZ12=zpos2/zpos1;
-
-// Definition of prototype for chambers in the second plane of station 1
-
- tpar[0]= 0.;
- tpar[1]= 0.;
- tpar[2]= 0.;
-
- gMC->Gsvolu("CC2A", "BOX ", idAlu1, tpar, 0); //Al
- gMC->Gsvolu("CB2A", "BOX ", idtmed[1107], tpar, 0); //Bakelite
- gMC->Gsvolu("CG2A", "BOX ", idtmed[1106], tpar, 0); //Gas streamer
-
-// chamber type A
- tpar[0] = -1.;
- tpar[1] = -1.;
-
- const Float_t kXMC2A=kXMC1A*kZ12;
- const Float_t kYMC2Am=0.;
- const Float_t kYMC2Ap=0.;
-
- tpar[2] = 0.1;
- gMC->Gsposp("CG2A", 1, "CB2A", 0., 0., 0., 0, "ONLY",tpar,3);
- tpar[2] = 0.3;
- gMC->Gsposp("CB2A", 1, "CC2A", 0., 0., 0., 0, "ONLY",tpar,3);
-
- tpar[2] = 0.4;
- tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ12;
- tpar[1] = kYMC1MIN*kZ12;
-
- gMC->Gsposp("CC2A", 1, "CM12",kXMC2A,kYMC2Am,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 2, "CM12",-kXMC2A,kYMC2Ap,kZMCp, 0, "ONLY", tpar, 3);
-
-
-// chamber type B
-
- tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ12;
- tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ12;
-
- const Float_t kXMC2B=kXMC1B*kZ12;
- const Float_t kYMC2Bp=kYMC1Bp*kZ12;
- const Float_t kYMC2Bm=kYMC1Bm*kZ12;
- gMC->Gsposp("CC2A", 3, "CM12",kXMC2B,kYMC2Bp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 4, "CM12",-kXMC2B,kYMC2Bm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 5, "CM12",kXMC2B,-kYMC2Bp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 6, "CM12",-kXMC2B,-kYMC2Bm,kZMCm, 0, "ONLY", tpar, 3);
-
-
-// chamber type C (end of type B !!)
-
- tpar[0] = (kXMC1MAX/2)*kZ12;
- tpar[1] = (kYMC1MAX/2)*kZ12;
-
- const Float_t kXMC2C=kXMC1C*kZ12;
- const Float_t kYMC2Cp=kYMC1Cp*kZ12;
- const Float_t kYMC2Cm=kYMC1Cm*kZ12;
- gMC->Gsposp("CC2A", 7, "CM12",kXMC2C,kYMC2Cp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 8, "CM12",-kXMC2C,kYMC2Cm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 9, "CM12",kXMC2C,-kYMC2Cp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 10, "CM12",-kXMC2C,-kYMC2Cm,kZMCm, 0, "ONLY", tpar, 3);
-
-// chamber type D, E and F (same size)
-
- tpar[0] = (kXMC1MAX/2.)*kZ12;
- tpar[1] = kYMC1MIN*kZ12;
-
- const Float_t kXMC2D=kXMC1D*kZ12;
- const Float_t kYMC2Dp=kYMC1Dp*kZ12;
- const Float_t kYMC2Dm=kYMC1Dm*kZ12;
- gMC->Gsposp("CC2A", 11, "CM12",kXMC2D,kYMC2Dm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 12, "CM12",-kXMC2D,kYMC2Dp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 13, "CM12",kXMC2D,-kYMC2Dm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 14, "CM12",-kXMC2D,-kYMC2Dp,kZMCp, 0, "ONLY", tpar, 3);
-
- const Float_t kYMC2Ep=kYMC1Ep*kZ12;
- const Float_t kYMC2Em=kYMC1Em*kZ12;
- gMC->Gsposp("CC2A", 15, "CM12",kXMC2D,kYMC2Ep,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 16, "CM12",-kXMC2D,kYMC2Em,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 17, "CM12",kXMC2D,-kYMC2Ep,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 18, "CM12",-kXMC2D,-kYMC2Em,kZMCm, 0, "ONLY", tpar, 3);
-
-
- const Float_t kYMC2Fp=kYMC1Fp*kZ12;
- const Float_t kYMC2Fm=kYMC1Fm*kZ12;
- gMC->Gsposp("CC2A", 19, "CM12",kXMC2D,kYMC2Fm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 20, "CM12",-kXMC2D,kYMC2Fp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 21, "CM12",kXMC2D,-kYMC2Fm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC2A", 22, "CM12",-kXMC2D,-kYMC2Fp,kZMCp, 0, "ONLY", tpar, 3);
-
-// Positioning second plane of station 1 in ALICE
-
- gMC->Gspos("CM12", 1, "ALIC", 0., 0., zpos2, 0, "ONLY");
-
-// End of geometry definition for the second plane of station 1
-
-
-
-// TRIGGER STATION 2 - TRIGGER STATION 2 - TRIGGER STATION 2
-
- // 03/00
- // zpos3 and zpos4 are now the middle of the first and second
- // plane of station 2 :
- // zpos3=(17075+16995)/2=17035 mm, thick/2=40 mm
- // zpos4=(17225+17145)/2=17185 mm, thick/2=40 mm
- //
- // zpos3m=16999 mm , zpos3p=17071 mm (middles of gas gaps)
- // zpos4m=17149 mm , zpos4p=17221 mm (middles of gas gaps)
- // rem : the total thickness accounts for 1 mm of al on both
- // side of the RPCs (see zpos3 and zpos4), as previously
- iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[12];
- iChamber2 =(AliMUONChamber*) (*fChambers)[13];
- Float_t zpos3=iChamber1->Z();
- Float_t zpos4=iChamber2->Z();
-
-
-// Mother volume definition
- tpar[0] = iChamber->RInner();
- tpar[1] = iChamber->ROuter();
- tpar[2] = 4.0;
-
- gMC->Gsvolu("CM21", "TUBE", idAir, tpar, 3);
- gMC->Gsvolu("CM22", "TUBE", idAir, tpar, 3);
-
-// Definition of the flange between the beam shielding and the RPC
-// ???? interface shielding
-
- tpar[0]= kRMIN2;
- tpar[1]= kRMAX2;
- tpar[2]= 4.0;
-
- gMC->Gsvolu("CF2A", "TUBE", idAlu1, tpar, 3); //Al
- gMC->Gspos("CF2A", 1, "CM21", 0., 0., 0., 0, "MANY");
- gMC->Gspos("CF2A", 2, "CM22", 0., 0., 0., 0, "MANY");
-
-
-
-// FIRST PLANE OF STATION 2 : proj ratio = zpos3/zpos1
-
- const Float_t kZ13=zpos3/zpos1;
-
-// Definition of prototype for chambers in the first plane of station 2
- tpar[0]= 0.;
- tpar[1]= 0.;
- tpar[2]= 0.;
-
- gMC->Gsvolu("CC3A", "BOX ", idAlu1, tpar, 0); //Al
- gMC->Gsvolu("CB3A", "BOX ", idtmed[1107], tpar, 0); //Bakelite
- gMC->Gsvolu("CG3A", "BOX ", idtmed[1106], tpar, 0); //Gas streamer
-
-
-// chamber type A
- tpar[0] = -1.;
- tpar[1] = -1.;
-
- const Float_t kXMC3A=kXMC1A*kZ13;
- const Float_t kYMC3Am=0.;
- const Float_t kYMC3Ap=0.;
-
- tpar[2] = 0.1;
- gMC->Gsposp("CG3A", 1, "CB3A", 0., 0., 0., 0, "ONLY",tpar,3);
- tpar[2] = 0.3;
- gMC->Gsposp("CB3A", 1, "CC3A", 0., 0., 0., 0, "ONLY",tpar,3);
-
- tpar[2] = 0.4;
- tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ13;
- tpar[1] = kYMC1MIN*kZ13;
- gMC->Gsposp("CC3A", 1, "CM21",kXMC3A,kYMC3Am,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 2, "CM21",-kXMC3A,kYMC3Ap,kZMCp, 0, "ONLY", tpar, 3);
-
-
-// chamber type B
- tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ13;
- tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ13;
-
- const Float_t kXMC3B=kXMC1B*kZ13;
- const Float_t kYMC3Bp=kYMC1Bp*kZ13;
- const Float_t kYMC3Bm=kYMC1Bm*kZ13;
- gMC->Gsposp("CC3A", 3, "CM21",kXMC3B,kYMC3Bp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 4, "CM21",-kXMC3B,kYMC3Bm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 5, "CM21",kXMC3B,-kYMC3Bp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 6, "CM21",-kXMC3B,-kYMC3Bm,kZMCm, 0, "ONLY", tpar, 3);
-
-
-// chamber type C (end of type B !!)
- tpar[0] = (kXMC1MAX/2)*kZ13;
- tpar[1] = (kYMC1MAX/2)*kZ13;
-
- const Float_t kXMC3C=kXMC1C*kZ13;
- const Float_t kYMC3Cp=kYMC1Cp*kZ13;
- const Float_t kYMC3Cm=kYMC1Cm*kZ13;
- gMC->Gsposp("CC3A", 7, "CM21",kXMC3C,kYMC3Cp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 8, "CM21",-kXMC3C,kYMC3Cm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 9, "CM21",kXMC3C,-kYMC3Cp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 10, "CM21",-kXMC3C,-kYMC3Cm,kZMCm, 0, "ONLY", tpar, 3);
-
-
-// chamber type D, E and F (same size)
-
- tpar[0] = (kXMC1MAX/2.)*kZ13;
- tpar[1] = kYMC1MIN*kZ13;
-
- const Float_t kXMC3D=kXMC1D*kZ13;
- const Float_t kYMC3Dp=kYMC1Dp*kZ13;
- const Float_t kYMC3Dm=kYMC1Dm*kZ13;
- gMC->Gsposp("CC3A", 11, "CM21",kXMC3D,kYMC3Dm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 12, "CM21",-kXMC3D,kYMC3Dp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 13, "CM21",kXMC3D,-kYMC3Dm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 14, "CM21",-kXMC3D,-kYMC3Dp,kZMCp, 0, "ONLY", tpar, 3);
-
- const Float_t kYMC3Ep=kYMC1Ep*kZ13;
- const Float_t kYMC3Em=kYMC1Em*kZ13;
- gMC->Gsposp("CC3A", 15, "CM21",kXMC3D,kYMC3Ep,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 16, "CM21",-kXMC3D,kYMC3Em,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 17, "CM21",kXMC3D,-kYMC3Ep,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 18, "CM21",-kXMC3D,-kYMC3Em,kZMCm, 0, "ONLY", tpar, 3);
-
- const Float_t kYMC3Fp=kYMC1Fp*kZ13;
- const Float_t kYMC3Fm=kYMC1Fm*kZ13;
- gMC->Gsposp("CC3A", 19, "CM21",kXMC3D,kYMC3Fm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 20, "CM21",-kXMC3D,kYMC3Fp,kZMCp, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 21, "CM21",kXMC3D,-kYMC3Fm,kZMCm, 0, "ONLY", tpar, 3);
- gMC->Gsposp("CC3A", 22, "CM21",-kXMC3D,-kYMC3Fp,kZMCp, 0, "ONLY", tpar, 3);