- zcoor -= zSenStrip;
- j++;
- Int_t upDown = -1; // upDown=-1 -> Upper strip
- // upDown=+1 -> Lower strip
- do{
- ang = atan(zcoor/radius);
- ang *= kRaddeg;
- AliMatrix (idrotm[nrot], 90., 0.,90.-ang,90.,-ang, 90.);
- AliMatrix (idrotm[nrot+1],90.,180.,90.+ang,90., ang, 90.);
- ang /= kRaddeg;
- ycoor = -14.5+ kspace; //2 cm over front plate
- ycoor += (1-(upDown+1)/2)*gap;
- gMC->Gspos("FSTR",j ,"FLTA",0.,ycoor, zcoor,idrotm[nrot], "ONLY");
- gMC->Gspos("FSTR",j+1,"FLTA",0.,ycoor,-zcoor,idrotm[nrot+1],"ONLY");
- if(fDebug) {
- printf("%s: %f, St. %2i, Pl.3 ",ClassName(),ang*kRaddeg,i);
- printf("y = %f, z = %f, zpos = %f \n",ycoor,zcoor,zpos);
- }
- j += 2;
- upDown*= -1; // Alternate strips
- zcoor = zcoor-(zSenStrip/2)/TMath::Cos(ang)-
- upDown*gap*TMath::Tan(ang)-
- (zSenStrip/2)/TMath::Cos(ang);
- } while (zcoor-(stripWidth/2)*TMath::Cos(ang)>-t+zFLTC+zFLTB+db*2);
-
- zcoor = zcoor+(zSenStrip/2)/TMath::Cos(ang)+
- upDown*gap*TMath::Tan(ang)+
- (zSenStrip/2)/TMath::Cos(ang);
-
- gap = fGapB;
- zcoor = zcoor-(zSenStrip/2)/TMath::Cos(ang)-
- upDown*gap*TMath::Tan(ang)-
- (zSenStrip/2)/TMath::Cos(ang);
-
- ang = atan(zcoor/radius);
- ang *= kRaddeg;
- AliMatrix (idrotm[nrot], 90., 0.,90.-ang,90.,-ang, 90.);
- AliMatrix (idrotm[nrot+1],90.,180.,90.+ang,90., ang, 90.);
- ang /= kRaddeg;
-
- ycoor = -14.5+ kspace; //2 cm over front plate
- ycoor += (1-(upDown+1)/2)*gap;
- gMC->Gspos("FSTR",j ,"FLTA",0.,ycoor, zcoor,idrotm[nrot], "ONLY");
- gMC->Gspos("FSTR",j+1,"FLTA",0.,ycoor,-zcoor,idrotm[nrot+1],"ONLY");
- if(fDebug) {
- printf("%s: %f, St. %2i, Pl.3 ",ClassName(),ang*kRaddeg,i);
- printf("y = %f, z = %f, zpos = %f \n",ycoor,zcoor,zpos);
- }
- ycoor = -hTof/2.+ kspace;//2 cm over front plate
-
- // Plate B
-
- nrot = 0;
- i=1;
- upDown = 1;
- Float_t deadRegion = 1.0;//cm
-
- zpos = zcoor - (zSenStrip/2)/TMath::Cos(ang)-
- upDown*gap*TMath::Tan(ang)-
- (zSenStrip/2)/TMath::Cos(ang)-
- deadRegion/TMath::Cos(ang);
-
- ang = atan(zpos/radius);
- ang *= kRaddeg;
- AliMatrix (idrotm[nrot], 90., 0., 90.-ang,90.,ang, 270.);
- ang /= kRaddeg;
- ycoor = -hTof*0.5+ kspace ; //2 cm over front plate
- ycoor += (1-(upDown+1)/2)*gap;
- zcoor = zpos+(zFLTA*0.5+zFLTB*0.5+db); // Moves to the system of the modulus FLTB
- gMC->Gspos("FSTR",i, "FLTB", 0., ycoor, zcoor,idrotm[nrot], "ONLY");
- if(fDebug) {
- printf("%s: %f, St. %2i, Pl.4 ",ClassName(),ang*kRaddeg,i);
- printf("y = %f, z = %f, zpos = %f \n",ycoor,zcoor,zpos);
- }
- i++;
- upDown*=-1;
-
- do {
- zpos = zpos - (zSenStrip/2)/TMath::Cos(ang)-
- upDown*gap*TMath::Tan(ang)-
- (zSenStrip/2)/TMath::Cos(ang);
- ang = atan(zpos/radius);
- ang *= kRaddeg;
- AliMatrix (idrotm[nrot], 90., 0., 90.-ang,90.,ang, 270.);
- ang /= kRaddeg;
- ycoor = -hTof*0.5+ kspace ; //2 cm over front plate
- ycoor += (1-(upDown+1)/2)*gap;
- zcoor = zpos+(zFLTA*0.5+zFLTB*0.5+db); // Moves to the system of the modulus FLTB
- gMC->Gspos("FSTR",i, "FLTB", 0., ycoor, zcoor,idrotm[nrot], "ONLY");
- if(fDebug) {
- printf("%s: %f, St. %2i, Pl.4 ",ClassName(),ang*kRaddeg,i);
- printf("y = %f, z = %f, zpos = %f \n",ycoor,zcoor,zpos);
- }
- upDown*=-1;
- i++;
- } while (TMath::Abs(ang*kRaddeg)<22.5);
- //till we reach a tilting angle of 22.5 degrees
-
- ycoor = -hTof*0.5+ kspace ; //2 cm over front plate
- zpos = zpos - zSenStrip/TMath::Cos(ang);
-
- do {
- ang = atan(zpos/radius);
- ang *= kRaddeg;
- AliMatrix (idrotm[nrot], 90., 0., 90.-ang,90.,ang, 270.);
- ang /= kRaddeg;
- zcoor = zpos+(zFLTB/2+zFLTA/2+db);
- gMC->Gspos("FSTR",i, "FLTB", 0., ycoor, zcoor,idrotm[nrot], "ONLY");
- zpos = zpos - zSenStrip/TMath::Cos(ang);
- if(fDebug) {
- printf("%s: %f, St. %2i, Pl.4 ",ClassName(),ang*kRaddeg,i);
- printf("y = %f, z = %f, zpos = %f \n",ycoor,zcoor,zpos);
- }
- i++;
-
- } while (zpos-stripWidth*0.5/TMath::Cos(ang)>-t+zFLTC+db);
-
- // Plate C
-
- zpos = zpos + zSenStrip/TMath::Cos(ang);
-
- zpos = zpos - (zSenStrip/2)/TMath::Cos(ang)+
- gap*TMath::Tan(ang)-
- (zSenStrip/2)/TMath::Cos(ang);
-
- nrot = 0;
- i=0;
- ycoor= -hTof*0.5+kspace+gap;
-
- do {
- i++;
- ang = atan(zpos/radius);
- ang *= kRaddeg;
- AliMatrix (idrotm[nrot], 90., 0., 90.-ang,90.,ang, 270.);
- ang /= kRaddeg;
- zcoor = zpos+(zFLTC*0.5+zFLTB+zFLTA*0.5+db*2);
- gMC->Gspos("FSTR",i, "FLTC", 0., ycoor, zcoor,idrotm[nrot], "ONLY");
- if(fDebug) {
- printf("%s: %f, St. %2i, Pl.5 ",ClassName(),ang*kRaddeg,i);
- printf("y = %f, z = %f, zpos = %f \n",ycoor,zcoor,zpos);
- }
- zpos = zpos - zSenStrip/TMath::Cos(ang);
- } while (zpos-stripWidth*TMath::Cos(ang)*0.5>-t);
-
-
-////////// Layers after strips /////////////////
-// honeycomb (Polyethilene) Layer after (1.2cm)