X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSv11.cxx;h=7b1cdbcafa70a433a774ff4abbaba8479af8491b;hb=e2e989f91c676946b600de9bbf5a88c4c8d481c0;hp=edb0eb64aab6575ea759bcd5078563c6d780d55c;hpb=dfefbaec6ca3cd3a201055c975a32c05746c28ef;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSv11.cxx b/ITS/AliITSv11.cxx index edb0eb64aab..7b1cdbcafa7 100644 --- a/ITS/AliITSv11.cxx +++ b/ITS/AliITSv11.cxx @@ -13,13 +13,7 @@ * provided "as is" without express or implied warranty. * **************************************************************************/ -/* -$Log$ -Revision 1.1 2003/01/20 23:32:49 nilsen -New ITS geometry. Only a Skeleton for now. - -$Id$ -*/ +/* $Id$ */ ////////////////////////////////////////////////////////////////////////////// // // @@ -37,1193 +31,363 @@ $Id$ // // ////////////////////////////////////////////////////////////////////////////// // See AliITSv11::StepManager(). -#include +// General C/C++ includes #include #include +// General Root includes +#include #include -#include -#include -#include -#include -#include -#include // only required for Tracking function? -#include +#include #include +#include #include #include -#include -#include -#include - - +// Root Geometry includes +#include +#include +#include +#include +#include // contaings TGeoTubeSeg +#include +#include +#include +#include +#include +#include +#include "AliITSv11GeometrySupport.h" +// General AliRoot includes #include "AliRun.h" #include "AliMagF.h" #include "AliConst.h" -#include "AliITSGeant3Geometry.h" -#include "AliITShit.h" -#include "AliITS.h" -#include "AliITSv11.h" -#include "AliITSgeom.h" +// ITS specific includes #include "AliITSgeomSPD.h" #include "AliITSgeomSDD.h" -#include "AliITSgeomSSD.h" #include "AliITSDetType.h" #include "AliITSresponseSPD.h" #include "AliITSresponseSDD.h" -#include "AliITSresponseSSD.h" #include "AliITSsegmentationSPD.h" #include "AliITSsegmentationSDD.h" -#include "AliITSsegmentationSSD.h" #include "AliITSsimulationSPD.h" #include "AliITSsimulationSDD.h" -#include "AliITSsimulationSSD.h" -#include "AliITSClusterFinderSPD.h" -#include "AliITSClusterFinderSDD.h" #include "AliITSClusterFinderSSD.h" - +#include "AliITSv11.h" +#include "AliITSv11GeometrySPD.h" +#include "AliITSv11GeometrySDD.h" +#include "AliITSv11GeometrySupport.h" ClassImp(AliITSv11) -//______________________________________________________________________ -AliITSv11::AliITSv11() : AliITS() { - //////////////////////////////////////////////////////////////////////// - // Standard default constructor for the ITS version 11. - //////////////////////////////////////////////////////////////////////// -} -//______________________________________________________________________ -AliITSv11::AliITSv11(const char *title) : AliITS("ITS", title){ - //////////////////////////////////////////////////////////////////////// - // Standard constructor for the ITS version 11. - //////////////////////////////////////////////////////////////////////// -} -//______________________________________________________________________ -AliITSv11::~AliITSv11() { - //////////////////////////////////////////////////////////////////////// - // Standard destructor for the ITS version 11. - //////////////////////////////////////////////////////////////////////// -} -//______________________________________________________________________ -void AliITSv11::Box(const char gnam[3],const TString &dis, - Double_t dx,Double_t dy,Double_t dz,Int_t med){ - // Interface to TMC->Gsvolu() for ITS bos geometries. Box with faces - // perpendicular to the axes. It has 3 paramters. See SetScale() for - // units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dx half-length of box in x-axis - // Double_t dy half-length of box in y-axis - // Double_t dz half-length of box in z-axis - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[3]; - - param[0] = fScale*dx; - param[1] = fScale*dy; - param[2] = fScale*dz; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"BOX ",fidmed[med],param,3); -} -//______________________________________________________________________ -void AliITSv11::Trapezoid1(const char gnam[3],const TString &dis, - Double_t dxn,Double_t dxp,Double_t dy,Double_t dz, - Int_t med){ - // Interface to TMC->Gsvolu() for ITS TRD1 geometries. Trapezoid with the - // x dimension varing along z. It has 4 parameters. See SetScale() for - // units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dxn half-length along x at the z surface positioned - // at -DZ - // Double_t dxp half-length along x at the z surface positioned - // at +DZ - // Double_t dy half-length along the y-axis - // Double_t dz half-length along the z-axis - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[4]; - - param[0] = fScale*dxn; - param[1] = fScale*dxp; - param[2] = fScale*dy; - param[3] = fScale*dz; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"TRD1",fidmed[med],param,4); -} -//______________________________________________________________________ -void AliITSv11::Trapezoid2(const char gnam[3],const TString &dis,Double_t dxn, - Double_t dxp,Double_t dyn,Double_t dyp,Double_t dz, - Int_t med){ - // Interface to TMC->Gsvolu() for ITS TRD2 geometries. Trapezoid with the - // x and y dimension varing along z. It has 5 parameters. See SetScale() - // for units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dxn half-length along x at the z surface positioned - // at -DZ - // Double_t dxp half-length along x at the z surface positioned - // at +DZ - // Double_t dyn half-length along x at the z surface positioned - // at -DZ - // Double_t dyp half-length along x at the z surface positioned - // at +DZ - // Double_t dz half-length along the z-axis - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[5]; - - param[0] = fScale*dxn; - param[1] = fScale*dxp; - param[2] = fScale*dyn; - param[3] = fScale*dyp; - param[4] = fScale*dz; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"TRD2",fidmed[med],param,5); -} -//______________________________________________________________________ -void AliITSv11::Trapezoid(const char gnam[3],const TString &dis,Double_t dz, - Double_t thet,Double_t phi,Double_t h1,Double_t bl1, - Double_t tl1,Double_t alp1,Double_t h2,Double_t bl2, - Double_t tl2,Double_t alp2,Int_t med){ - // Interface to TMC->Gsvolu() for ITS TRAP geometries. General Trapezoid, - // The faces perpendicular to z are trapezia and their centers are not - // necessarily on a line parallel to the z axis. This shape has 11 - // parameters, but only cosidering that the faces should be planar, only 9 - // are really independent. A check is performed on the user parameters and - // a message is printed in case of non-planar faces. Ignoring this warning - // may cause unpredictable effects at tracking time. See SetScale() - // for units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dz Half-length along the z-asix - // Double_t thet Polar angle of the line joing the center of the - // face at -dz to the center of the one at dz - // [degree]. - // Double_t phi aximuthal angle of the line joing the center of - // the face at -dz to the center of the one at +dz - // [degree]. - // Double_t h1 half-length along y of the face at -dz. - // Double_t bl1 half-length along x of the side at -h1 in y of - // the face at -dz in z. - // Double_t tl1 half-length along x of teh side at +h1 in y of - // the face at -dz in z. - // Double_t alp1 angle with respect to the y axis from the center - // of the side at -h1 in y to the cetner of the - // side at +h1 in y of the face at -dz in z - // [degree]. - // Double_t h2 half-length along y of the face at +dz - // Double_t bl2 half-length along x of the side at -h2 in y of - // the face at +dz in z. - // Double_t tl2 half-length along x of the side at _h2 in y of - // the face at +dz in z. - // Double_t alp2 angle with respect to the y axis from the center - // of the side at -h2 in y to the center of the - // side at +h2 in y of the face at +dz in z - // [degree]. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[11]; - - param[0] = fScale*dz; - param[1] = thet; - param[2] = phi; - param[3] = fScale*h1; - param[4] = fScale*bl1; - param[5] = fScale*tl1; - param[6] = alp1; - param[7] = fScale*h2; - param[8] = fScale*bl2; - param[9] = fScale*tl2; - param[10] = alp2; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"TRAP",fidmed[med],param,11); -} -//______________________________________________________________________ -void AliITSv11::Tube(const char gnam[3],const TString &dis,Double_t rmin, - Double_t rmax,Double_t dz,Int_t med){ - // Interface to TMC->Gsvolu() for ITS TUBE geometries. Simple Tube. It has - // 3 parameters. See SetScale() - // for units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t rmin Inside Radius. - // Double_t rmax Outside Radius. - // Double_t dz half-length along the z-axis - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[3]; - - param[0] = fScale*rmin; - param[1] = fScale*rmax; - param[2] = fScale*dz; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"TUBE",fidmed[med],param,3); -} -//______________________________________________________________________ -void AliITSv11::TubeSegment(const char gnam[3],const TString &dis, - Double_t rmin,Double_t rmax,Double_t dz, - Double_t phi1,Double_t phi2,Int_t med){ - // Interface to TMC->Gsvolu() for ITS TUBE geometries. Phi segment of a - // tube. It has 5 parameters. Phi1 should be smaller than phi2. If this is - // not the case, the system adds 360 degrees to phi2. See SetScale() - // for units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t rmin Inside Radius. - // Double_t rmax Outside Radius. - // Double_t dz half-length along the z-axis - // Double_t phi1 Starting angle of the segment [degree]. - // Double_t phi2 Ending angle of the segment [degree]. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[5]; - - param[0] = fScale*rmin; - param[1] = fScale*rmax; - param[2] = fScale*dz; - param[3] = phi1; - param[4] = phi2; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"TUBS",fidmed[med],param,5); -} -//______________________________________________________________________ -void AliITSv11::Cone(const char gnam[3],const TString &dis,Double_t dz, - Double_t rmin1,Double_t rmax1,Double_t rmin2, - Double_t rmax2,Int_t med){ - // Interface to TMC->Gsvolu() for ITS Cone geometries. Conical tube. It - // has 5 parameters. See SetScale() - // for units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dz half-length along the z-axis - // Double_t rmin1 Inside Radius at -dz. - // Double_t rmax1 Outside Radius at -dz. - // Double_t rmin2 inside radius at +dz. - // Double_t rmax2 outside radius at +dz. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[5]; - - param[0] = fScale*dz; - param[1] = fScale*rmin1; - param[2] = fScale*rmax1; - param[3] = fScale*rmin2; - param[4] = fScale*rmax2; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"CONS",fidmed[med],param,5); -} -//______________________________________________________________________ -void AliITSv11::ConeSegment(const char gnam[3],const TString &dis,Double_t dz, - Double_t rmin1,Double_t rmax1,Double_t rmin2, - Double_t rmax2,Double_t phi1,Double_t phi2, - Int_t med){ - // Interface to TMC->Gsvolu() for ITS ConS geometries. One segment of a - // conical tube. It has 7 parameters. Phi1 should be smaller than phi2. If - // this is not the case, the system adds 360 degrees to phi2. See - // SetScale() for units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dz half-length along the z-axis - // Double_t rmin1 Inside Radius at -dz. - // Double_t rmax1 Outside Radius at -dz. - // Double_t rmin2 inside radius at +dz. - // Double_t rmax2 outside radius at +dz. - // Double_t phi1 Starting angle of the segment [degree]. - // Double_t phi2 Ending angle of the segment [degree]. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[7]; - - param[0] = fScale*dz; - param[1] = fScale*rmin1; - param[2] = fScale*rmax1; - param[3] = fScale*rmin2; - param[4] = fScale*rmax2; - param[5] = phi1; - param[6] = phi2; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"CONS",fidmed[med],param,7); -} -//______________________________________________________________________ -void AliITSv11::Sphere(const char gnam[3],const TString &dis,Double_t rmin, - Double_t rmax,Double_t the1,Double_t the2,Double_t phi1, - Double_t phi2,Int_t med){ - // Interface to TMC->Gsvolu() for ITS SPHE geometries. Segment of a - // sphereical shell. It has 6 parameters. See SetScale() - // for units. Default units are geant 3 [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t rmin Inside Radius. - // Double_t rmax Outside Radius. - // Double_t the1 staring polar angle of the shell [degree]. - // Double_t the2 ending polar angle of the shell [degree]. - // Double_t phui staring asimuthal angle of the shell [degree]. - // Double_t phi2 ending asimuthal angle of the shell [degree]. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[6]; - - param[0] = fScale*rmin; - param[1] = fScale*rmax; - param[2] = the1; - param[3] = the2; - param[4] = phi1; - param[5] = phi2; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"SPHE",fidmed[med],param,6); -} -//______________________________________________________________________ -void AliITSv11::Parallelepiped(const char gnam[3],const TString &dis, - Double_t dx,Double_t dy,Double_t dz, - Double_t alph,Double_t thet,Double_t phi, - Int_t med){ - // Interface to TMC->Gsvolu() for ITS PARA geometries. Parallelepiped. It - // has 6 parameters. See SetScale() for units. Default units are geant 3 - // [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dx half-length allong x-axis - // Double_t dy half-length allong y-axis - // Double_t dz half-length allong z-axis - // Double_t alpha angle formed by the y axis and by the plane - // joining the center of teh faces parallel to the - // z-x plane at -dY and +dy [degree]. - // Double_t thet polar angle of the line joining the centers of - // the faces at -dz and +dz in z [degree]. - // Double_t phi azimuthal angle of teh line joing the centers of - // the faaces at -dz and +dz in z [degree]. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[6]; - - param[0] = fScale*dx; - param[1] = fScale*dy; - param[2] = fScale*dz; - param[3] = alpha; - param[4] = thet; - param[5] = phi; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"PARA",fidmed[med],param,6); -} -//______________________________________________________________________ -void AliITSv11::Polygon(const char gnam[3],const TString &dis,Double_t phi1, - Double_t dphi,Int_t npdv,Int_t nz,Double_t *z, - Double_t *rmin,Double_t *rmax,Double_t ,Int_t med){ - // Interface to TMC->Gsvolu() for ITS PGON geometry. Polygon It has 10 - // parameters or more. See SetScale() for units. Default units are geant 3 - // [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t phi1 the azimuthal angle at which the volume begins - // (angles are counted clouterclockwise) [degrees]. - // Double_t dphi opening angle of the volume, which extends from - // phi1 to phi1+dphi [degree]. - // Int_t npdv the number of sides of teh cross section between - // the given phi limits. - // Int_t nz number of planes perpendicular to the z axis - // where the dimension of the section is given - - // this number should be at least 2 and NP triples - // of number must follow. - // Double_t *z array [nz] of z coordiates of the sections.. - // Double_t *rmin array [nz] of radius of teh circle tangent to - // the sides of the inner polygon in teh - // cross-section. - // Double_t *rmax array [nz] of radius of the circle tangent to - // the sides of the outer polygon in the - // cross-section. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t *param; - Int_t n,i; - - n = 4+3*nz; - param = new Float_t[n] - param[0] = phi1; - param[1] = dphi; - param[2] = (Float_t)npdv; - param[3] = (Float_t)nz; - for(i=0;iGsvolu(name,"PGON",fidmed[med],param,n); - - delete[] param; -} -//______________________________________________________________________ -void AliITSv11::PolyCone(const char gnam[3],const TString &dis,Double_t phi1, - Double_t dphi,Int_t nz,Double_t *z,Double_t *rmin, - Double_t *rmax,Int_t med){ - // Interface to TMC->Gsvolu() for ITS PCON geometry. Poly-cone It has 9 - // parameters or more. See SetScale() for units. Default units are geant 3 - // [cm]. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t phi1 the azimuthal angle at which the volume begins - // (angles are counted clouterclockwise) [degrees]. - // Double_t dphi opening angle of the volume, which extends from - // phi1 to phi1+dphi [degree]. - // Int_t nz number of planes perpendicular to the z axis - // where the dimension of the section is given - - // this number should be at least 2 and NP triples - // of number must follow. - // Double_t *z Array [nz] of z coordinate of the section. - // Double_t *rmin Array [nz] of radius of teh inner circle in the - // cross-section. - // Double_t *rmax Array [nz] of radius of the outer circle in the - // cross-section. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t *param; - Int_t n,i; - - n = 3+3*nz; - param = new Float_t[n]; - param[0] = phi1; - param[1] = dphi; - param[2] = (Float_t) nz; - for(i=0;iGsvolu(name,"PCON",fidmed[med],param,n); - - delete[] param; -} -//______________________________________________________________________ -void AliITSv11::TubeElliptical(const char gnam[3],const TString &dis, - Double_t p1,Double_t p2,Double_t dz,Int_t med){ - // Interface to TMC->Gsvolu() for ITS ELTU geometries. Elliptical - // cross-section Tube. It has 3 parameters. See SetScale() - // for units. Default units are geant 3 [cm]. The equation of the surface - // is x^2 * p1^-2 + y^2 * p2^-2 = 1. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t p1 semi-axis of the elipse along x. - // Double_t p2 semi-axis of the elipse along y. - // Double_t dz half-length along the z-axis - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[3]; - - param[0] = fScale*p1; - param[1] = fScale*p2; - param[2] = fScale*dz; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"ELTU",fidmed[med],param,3); -} -//______________________________________________________________________ -void AliITSv11::HyperbolicTube(const char gnam[3],const TString &dis, - Double_t rmin,Double_t rmax,Double_t dz, - Double_t thet,Int_t med){ - // Interface to TMC->Gsvolu() for ITS HYPE geometries. Hyperbolic tube. - // Fore example the inner and outer surfaces are hyperboloids, as would be - // foumed by a system of cylinderical wires which were then rotated - // tangentially about their centers. It has 4 parameters. See SetScale() - // for units. Default units are geant 3 [cm]. The hyperbolic surfaces are - // given by r^2 = (ztan(thet)^2 + r(z=0)^2. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t rmin Inner radius at z=0 where tube is narrowest. - // Double_t rmax Outer radius at z=0 where tube is narrowest. - // Double_t dz half-length along the z-axis - // Double_t thet stero angel of rotation of the two faces - // [degrees]. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[4]; - - param[0] = fScale*rmin; - param[1] = fScale*rmax; - param[2] = fScale*dz; - param[3] = thet; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"HYPE",fidmed[med],param,4); -} -//______________________________________________________________________ -void AliITSv11::TwistedTrapezoid(const char gnam[3],const TString &dis, - Double_t dz,Double_t thet,Double_t phi, - Double_t twist,Double_t h1,Double_t bl1, - Double_t tl1,Double_t apl1,Double_t h2, - Double_t bl2,Double_t tl2,Double_t apl2, - Int_t med){ - // Interface to TMC->Gsvolu() for ITS GTRA geometries. General twisted - // trapazoid. The faces perpendicular to z are trapazia and their centers - // are not necessarily on a line parallel to the z axis as the TRAP. - // Additionally, the faces may be twisted so that none of their edges are - // parallel. It is a TRAP shape, exept that it is twisted in the x-y plane - // as a function of z. The parallel sides perpendicular to the x axis are - // rotated with respect to the x axis by an angle TWIST, which is one of - // the parameters. The shape is defined by the eight corners and is assumed - // to be constructed of straight lines joingin points on the boundry of the - // trapezoidal face at Z=-dz to the coresponding points on the face at - // z=+dz. Divisions are not allowed. It has 12 parameters. See SetScale() - // for units. Default units are geant 3 [cm]. Note: This shape suffers from - // the same limitations than the TRAP. The tracking routines assume that - // the faces are planar, but htis constraint is not easily expressed in - // terms of the 12 parameters. Additionally, no check on th efaces is - // performed in this case. Users should avoid to use this shape as much as - // possible, and if they have to do so, they should make sure that the - // faces are really planes. If this is not the case, the result of the - // trasport is unpredictable. To accelerat ethe computations necessary for - // trasport, 18 additioanl parameters are calculated for this shape are - // 1 DXODZ dx/dz of the line joing the centers of the faces at z=+_dz. - // 2 DYODZ dy/dz of the line joing the centers of the faces at z=+_dz. - // 3 XO1 x at z=0 for line joing the + on parallel side, perpendicular - // corners at z=+_dz. - // 4 YO1 y at z=0 for line joing the + on parallel side, + on - // perpendicular corners at z=+-dz. - // 5 DXDZ1 dx/dz for line joing the + on parallel side, + on - // perpendicular corners at z=+-dz. - // 6 DYDZ1 dy/dz for line joing the + on parallel side, + on - // perpendicular corners at z=+-dz. - // 7 X02 x at z=0 for line joing the - on parallel side, + on - // perpendicular corners at z=+-dz. - // 8 YO2 y at z=0 for line joing the - on parallel side, + on - // perpendicular corners at z=+-dz. - // 9 DXDZ2 dx/dz for line joing the - on parallel side, + on - // perpendicular corners at z=+-dz. - // 10 DYDZ2dy/dz for line joing the - on parallel side, + on - // perpendicular corners at z=+-dz. - // 11 XO3 x at z=0 for line joing the - on parallel side, - on - // perpendicular corners at z=+-dz. - // 12 YO3 y at z=0 for line joing the - on parallel side, - on - // perpendicular corners at z=+-dz. - // 13 DXDZ3 dx/dzfor line joing the - on parallel side, - on - // perpendicular corners at z=+-dz. - // 14 DYDZ3 dydz for line joing the - on parallel side, - on - // perpendicular corners at z=+-dz. - // 15 XO4 x at z=0 for line joing the + on parallel side, - on - // perpendicular corners at z=+-dz. - // 16 YO4 y at z=0 for line joing the + on parallel side, - on - // perpendicular corners at z=+-dz. - // 17 DXDZ4 dx/dz for line joing the + on parallel side, - on - // perpendicular corners at z=+-dz. - // 18 DYDZ4 dydz for line joing the + on parallel side, - on - // perpendicular corners at z=+-dz. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t dz half-length along the z axis. - // Double_t thet polar angle of the line joing the center of the - // face at -dz to the center of the one at +dz - // [degrees]. - // Double_t phi Azymuthal angle of teh line joing the centre of - // the face at -dz to the center of the one at +dz - // [degrees]. - // Double_t twist Twist angle of the faces parallel to the x-y - // plane at z=+-dz around an axis parallel to z - // passing through their centre [degrees]. - // Double_t h1 Half-length along y of the face at -dz. - // Double_t bl1 half-length along x of the side -h1 in y of the - // face at -dz in z. - // Double_t tl1 half-length along x of the side at +h1 in y of - // the face at -dz in z. - // Double_t apl1 Angle with respect to the y ais from the center - // of the side at -h1 in y to the centere of the - // side at +h1 in y of the face at -dz in z - // [degrees]. - // Double_t h2 half-length along the face at +dz. - // Double_t bl2 half-length along x of the side at -h2 in y of - // the face at -dz in z. - // Double_t tl2 half-length along x of the side at +h2 in y of - // the face at +dz in z. - // Double_t apl2 angle with respect to the y axis from the center - // of the side at -h2 in y to the center of the side - // at +h2 in y of the face at +dz in z [degrees]. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[12]; - - param[0] = fScale*dz; - param[1] = thet; - param[2] = phi; - param[3] = twist; - param[4] = fScale*h1; - param[5] = fScale*bl1; - param[6] = fScale*tl1; - param[7] = alp1; - param[8] = fScale*h2; - param[9] = fScale*bl2; - param[10] = fScale*tl2; - param[11] = alp2; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"GTRA",fidmed[med],param,12); -} -//______________________________________________________________________ -void AliITSv11::CutTube(const char gnam[3],const TString &dis,Double_t rmin, - Double_t rmax,Double_t dz,Double_t phi1,Double_t phi2, - Double_t lx,Double_t ly,Double_t lz,Double_t hx, - Double_t hy,Double_t hz,Int_t med){ - // Interface to TMC->Gsvolu() for ITS CTUB geometries. Cut tube. A tube cut - // at the extremities with planes not necessarily perpendicular tot he z - // axis. It has 11 parameters. See SetScale() for units. Default units are - // geant 3 [cm]. phi1 should be smaller than phi2. If this is not the case, - // the system adds 360 degrees to phi2. - // Inputs: - // const char gnam[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // TString &dis String containging part discription. - // Double_t rmin Inner radius at z=0 where tube is narrowest. - // Double_t rmax Outer radius at z=0 where tube is narrowest. - // Double_t dz half-length along the z-axis - // Double_t dz half-length along the z-axis - // Double_t phi1 Starting angle of the segment [degree]. - // Double_t phi2 Ending angle of the segment [degree]. - // Double_t lx x component of a unit vector perpendicular to - // the face at -dz. - // Double_t ly y component of a unit vector perpendicular to - // the face at -dz. - // Double_t lz z component of a unit vector perpendicular to - // the face at -dz. - // Double_t hx x component of a unit vector perpendicular to - // the face at +dz. - // Double_t hy y component of a unit vector perpendicular to - // the face at +dz. - // Double_t hz z component of a unit vector perpendicular to - // the face at +dz. - // Int_t med media index number. - // Output: - // none. - // Return. - // none. - char name[4]; - Float_t param[11]; +/* + Some temparary #define's used untill ROOT has addoppted the proper + Getter in it's classes. + These Below are for TGeoPcon functions. +*/ - param[0] = fScale*rmin; - param[1] = fScale*rmax; - param[2] = fScale*dz; - param[3] = phi1; - param[4] = phi2; - param[5] = lx; - param[6] = ly; - param[7] = lz; - param[8] = hx; - param[9] = hy; - param[10] = hz; - name[0] = 'I'; - for(Int_t i=0;i<3;i++) name[i+1] = gnam[i]; - gMC->Gsvolu(name,"CTUB",fidmed[med],param,11); -} //______________________________________________________________________ -void AliITSv11::Pos(const char vol[3],Int_t cn,const char moth[3],Double_t x, - Double_t y,Double_t z,Int_t irot){ - // Place a copy of a volume previously defined by a call to GSVOLU inside - // its mother volulme moth. +AliITSv11::AliITSv11() : + AliITS(), + fGeomDetOut(kFALSE), + fGeomDetIn(kFALSE), + fMajorVersion(11), + fMinorVersion(0), + fDet1(0.0), + fDet2(0.0), + fChip1(0.0), + fChip2(0.0), + fRails(0), + fFluid(1), + fSPDgeom(0), + fSupgeom(0){ + // Standard default constructor for the ITS version 11. // Inputs: - // const char vol[3] 3 character geant volume name. The letter "I" - // is appended to the front to indecate that this - // is an ITS volume. - // const char moth[3] 3 character geant volume name of the mother volume - // in which vol will be placed. The letter "I" is - // appended to the front to indecate that this is an - // ITS volume. - // Double_t x The x positon of the volume in the mother's - // reference system - // Double_t y The y positon of the volume in the mother's - // reference system - // Double_t z The z positon of the volume in the mother's - // reference system - // Int_t irot the index for the rotation matrix to be used. - // irot=-1 => unit rotation. + // none. // Outputs: - // none. - // Return: - // none. - char name[4],mother[4]; - Float_t param[3]; - Int_t r=0,i; - - param[0] = x; - param[1] = y; - param[2] = z; - name[0] = 'I'; - for(i=0;i<3;i++) name[i+1] = vol[i]; - mother[0] = 'I'; - for(i=0;i<3;i++) mother[i+1] = moth[i]; - if(irot>=0) r=fidrot[irot]; - fMC->Gspos(name,mother,param[0],param[1],param[2],r,"ONLY"); -} -//______________________________________________________________________ -void AliITSv11::Matrix(Int_t irot,Double_t thet1,Double_t phi1, - Double_t thet2,Double_t phi2, - Double_t thet3,Double_t phi3){ - // Defines a Geant rotation matrix. checks to see if it is the unit - // matrix. If so, then no additonal matrix is defined. Stores rotation - // matrix irot in the data structure JROTM. If the matrix is not - // orthonormal, it will be corrected by setting y' perpendicular to x' - // and z' = x' X y'. A warning message is printed in this case. + // none. + // Return + // A default constructed AliITSv11 class. + + //fITSV = 0; + //fcS = 0; +// fcD = 0; + + SetEUCLID(kFALSE); +} +//______________________________________________________________________ +AliITSv11::AliITSv11(Int_t debugITS,Int_t debugSPD,Int_t debugSDD, + Int_t debugSSD,Int_t debugSUP) : + AliITS(), + fGeomDetOut(kFALSE), + fGeomDetIn(kFALSE), + fMajorVersion(11), + fMinorVersion(0), + fDet1(0.0), + fDet2(0.0), + fChip1(0.0), + fChip2(0.0), + fRails(0), + fFluid(1), + fSPDgeom(0), + fSupgeom(0){ + // Standard default constructor for the ITS version 11. // Inputs: - // Int_t irot Intex specifing which rotation matrix. - // Double_t thet1 Polar angle for axisw x [degrees]. - // Double_t phi1 azimuthal angle for axis x [degrees]. - // Double_t thet12Polar angle for axisw y [degrees]. - // Double_t phi2 azimuthal angle for axis y [degrees]. - // Double_t thet3 Polar angle for axisw z [degrees]. - // Double_t phi3 azimuthal angle for axis z [degrees]. + // Int_t debugITS Debug flag for ITS code (required). + // Int_t debugSPD Debug flag for SPD geometry (default = 0). + // Int_t debugSDD Debug flag for SDD geometry (default = 0). + // Int_t debugSSD Debug flag for SSD geometry (default = 0). + // Int_t debugSUP Debug flag for SUPort geometry (default = 0). // Outputs: - // none. - // Return: - // none. - Float_t t1,p1,t2,p2,t3,p3; - - if(thet1==90.0&&phi1==0.0&&thet2==90.0&&phi2==90.0&&thet3==0.0&&phi3==0.0){ - fidrot[irot] = 0; // Unit matrix - }else{ - t1 = thet1; - p1 = phi1; - t2 = thet2; - p2 = phi2; - t3 = thet3; - p3 = phi3 - AliMatrix(fidrot[irot],t1,p1,t2,p2,t3,p3); - } // end if -} -//______________________________________________________________________ -void AliITSv11::Matrix(Int_t irot,Int_t axis,Double_t thet){ - // Defines a Geant rotation matrix. checks to see if it is the unit - // matrix. If so, then no additonal matrix is defined. Stores rotation - // matrix irot in the data structure JROTM. If the matrix is not - // orthonormal, it will be corrected by setting y' perpendicular to x' - // and z' = x' X y'. A warning message is printed in this case. + // none. + // Return + // A Constructor for ITS geometry version 11 useful for Geometry display. + + SetEUCLID(kFALSE); + debugSSD = debugSDD; // remove waring for unused variables. + SetDebug(debugITS); + fSPDgeom = new AliITSv11GeometrySPD(debugSPD); + fSDDgeom = new AliITSv11GeometrySDD(debugSDD); + fSupgeom = new AliITSv11GeometrySupport(debugSUP); +} +//______________________________________________________________________ +AliITSv11::AliITSv11(const char *title) : + AliITS("ITS", title), + fGeomDetOut(kFALSE), + fGeomDetIn(kFALSE), + fMajorVersion(11), + fMinorVersion(0), + fDet1(0.0), + fDet2(0.0), + fChip1(0.0), + fChip2(0.0), + fRails(0), + fFluid(1), + fSPDgeom(0), + fSupgeom(0){ + // Standard constructor for the ITS version 11. // Inputs: - // Int_t irot Intex specifing which rotation matrix. - // Int_t axis Axis about which rotation is to be done. - // Double_t thet Angle to rotate by [degrees]. + // const char *title The title of for this geometry. // Outputs: - // none. - // Return: - // none. - - if(thet==0.0){ - fidrot[irot] = 0; // Unit matrix - }else{ - switch (irot) { - case 0: //Rotate about x-axis, x-axis does not change. - AliMatrix(fidrot[irot],90.0,0.0,90.0+thet,90.0,thet,90.0); - break; - case 1: //Rotate about y-axis, y-axis does not change. - AliMatrix(fidrot[irot],-90.0-thet,0.0,90.0,90.0,thet,90.0); - break; - case 2: //Rotate about z-axis, z-axis does not change. - AliMatrix(fidrot[irot],90.0,thet,90.0,-thet-90.0,0.0,0.0); - break; - default: - Error("Matrix","axis must be either 0, 1, or 2. for matrix=%d", - irot); - break; - } // end switch - } // end if + // none. + // Return + // A Standard constructed AliITSv11 class. + SetEUCLID(kFALSE); + fSPDgeom = new AliITSv11GeometrySPD(GetDebug()); + fSupgeom = new AliITSv11GeometrySupport(GetDebug()); } //______________________________________________________________________ -void AliITSv11::Matrix(Int_t irot,Double_t rot[3][3]){ - // Defines a Geant rotation matrix. checks to see if it is the unit - // matrix. If so, then no additonal matrix is defined. Stores rotation - // matrix irot in the data structure JROTM. If the matrix is not - // orthonormal, it will be corrected by setting y' perpendicular to x' - // and z' = x' X y'. A warning message is printed in this case. +AliITSv11::~AliITSv11() { + // Standard destructor for the ITS version 11. // Inputs: - // Int_t irot Intex specifing which rotation matrix. - // Double_t rot[3][3] The 3 by 3 rotation matrix. + // none. // Outputs: - // none. - // Return: - // none. - - if(rot[0][0]==1.0&&rot[1][1]==1.0&&rot[2][2]==1.0&& - rot[0][1]==0.0&&rot[0][2]==0.0&&rot[1][0]==0.0&& - rot[1][2]==0.0&&rot[2][0]==0.0&&rot[2][1]==0.0){ - fidrot[irot] = 0; // Unit matrix - }else{ - Double_t si,c=180./TMath::Pi(); - Double_t ang[6]; - - ang[1] = TMath::ATan2(rot[0][1],rot[0][0]); - if(TMath::Cos(ang[1])!=0.0) si = rot[0][0]/TMath::Cos(ang[1]); - else si = rot[0][1]/TMath::Sin(ang[1]); - ang[0] = TMath::ATan2(si,rot[0][2]); - - ang[3] = TMath::ATan2(rot[1][1],rot[1][0]); - if(TMath::Cos(ang[3])!=0.0) si = rot[1][0]/TMath::Cos(ang[3]); - else si = rot[1][1]/TMath::Sin(ang[3]); - ang[2] = TMath::ATan2(si,rot[1][2]); - - ang[5] = TMath::ATan2(rot[2][1],rot[2][0]); - if(TMath::Cos(ang[5])!=0.0) si = rot[2][0]/TMath::Cos(ang[5]); - else si = rot[2][1]/TMath::Sin(ang[5]); - ang[4] = TMath::ATan2(si,rot[2][2]); + // none. + // Return + // none. - for(Int_t i=0;i<6;i++) {ang[i] *= c; if(ang[i]<0.0) ang[i] += 360.;} - AliMatrix(fidrot[irot],ang[0],ang[1],ang[2],ang[3],ang[4],ang[5]); - } // end if +// if(fITSV!=0) delete fITSV; +// if(fcS!=0) delete fcS; +// if(fcD!=0) delete fcD; } //______________________________________________________________________ -Float_t AliITSv11::GetA(Int_t z){ - // Returns the isotopicaly averaged atomic number. +AliITSv11::AliITSv11(const AliITSv11 &source) : AliITS(source){ + // Copy Constructor for ITS version 11. // Inputs: - // Int_t z Elemental number + // AliITSv11 &source class to be copied from. // Outputs: - // none. - // Return: - // The atomic mass number. - const Float_t A[]={ 1.00794 , 4.0026902, 6.941 , 9.012182 , 10.811 , - 12.01007 , 14.00674 , 15.9994 , 18.9984032, 20.1797 , - 22.98970 , 24.3050 , 26.981538, 28.0855 , 30.973761, - 32.066 , 35.4527 , 39.948 , 39.0983 , 40.078 , - 44.95591 , 47.867 , 50.9415 , 51.9961 , 54.938049, - 55.845 , 58.933200 , 58.6934 , 63.546 , 65.39 , - 69.723 , 72.61 , 74.92160 , 78.96 , 79.904 , - 83.80 , 85.4678 , 87.62 , 88.9085 , 91.224 , - 92.90638 , 95.94 , 97.907215, 101.07 ,102.90550 , - 106.42 ,107.8682 ,112.411 ,114.818 ,118.710 , - 121.760 ,127.60 ,126.90447 ,131.29 ,132.90545 , - 137.327 ,138.9055 ,140.116 ,140.90765 ,144.24 , - 144.912746,150.36 ,151.964 ,157.25 ,158.92534 , - 162.50 ,164.93032 ,167.26 ,168.93421 ,173.04 , - 174.967 ,178.49 ,180.9479 ,183.84 ,186.207 , - 190.23 ,192.217 ,195.078 ,196.96655 ,200.59 , - 204.3833 ,207.2 ,208.98038,208.982415 ,209.987131, - 222.017570 ,223.019731,226.025402,227.027747 ,232.0381 , - 231.03588 238.0289}; + // none. + // Return + // none. - if(z<1||z>92){ - Error("GetA","z must be 0BuildDisplayGeometry(); } //______________________________________________________________________ -Float_t AliITSv11::GetStandardEfraction(Int_t istd){ - // Returns one of a set of standard E fraction values. +void AliITSv11::CreateGeometry(){ + // This routine defines and Creates the geometry for version 11 of + // the ITS. The geometry is used by the particle trasport routines, + // and therefore, is very detailed. // Inputs: - // Int_t istd Index to indecate which standard. + // none. // Outputs: - // none. - // Return: - // The appropreate standard E fraction value [#]. - Float_t t[]={0.1, // default - 0.1, // Silicon detectors... - 0.1, // Air in central detectors region - 1.0 // Material in non-centeral region - }; - return t[istd]; + // none. + // Return + // none. + TVector3 t(0.0,0.0,0.0); + const Double_t kcm = 1.0; + + TGeoManager *mgr = gGeoManager; + TGeoVolume *vALIC = mgr->GetTopVolume(); + + TGeoPcon *sITS = new TGeoPcon("ITS Top Volume, Daughter of ALIC", + 0.0,360.0,2); + // DefineSection(section number, Z, Rmin, Rmax). + sITS->DefineSection(0,-300.0*kcm,0.01*kcm,50.0*kcm); + sITS->DefineSection(1,+300.0*kcm,0.01*kcm,50.0*kcm); + TGeoVolume *vITS = new TGeoVolume("ITSV",sITS,0); + mgr->AddVolume(vITS); + vITS->SetVisibility(kFALSE); + vALIC->AddNode(vITS,1,0); + // + fSPDgeom->CarbonFiberSector(vITS); + fSDDgeom->SetGeomParameters(); // needed + fSDDgeom->Layer3(vITS); + fSDDgeom->Layer4(vITS); + fSupgeom->SPDCone(vITS); + fSupgeom->SPDThermalSheald(vITS); + fSupgeom->SDDCone(vITS); + fSupgeom->SSDCone(vITS); + fSupgeom->ServicesCableSupport(vITS); } //______________________________________________________________________ -void AliITSv11::Element(Int_t imat,const char* name,Int_t z,Double_t dens, - Int_t istd){ - // Defines a Geant single element material and sets its Geant medium - // proporties. The average atomic A is assumed to be given by their - // natural abundances. Things like the radiation length are calculated - // for you. +void AliITSv11::CreateMaterials(){ + // Create ITS materials + // This function defines the default materials used in the Geant + // Monte Carlo simulations for the geometries AliITSv11. + // In general it is automatically replaced by // Inputs: - // Int_t imat Material number. - // const char* name Material name. No need to add a $ at the end. - // Int_t z The elemental number. - // Double_t dens The density of the material [g/cm^3]. - // Int_t istd Defines which standard set of transport parameters - // which should be used. - // Output: - // none. - // Return: - // none. - Float_t rad,Z,A=GetA(z),tmax,stemax,deemax,epsilon; - char *name2; - Int_t len; - - len = strlng(name)+1; - name2 = new char[len]; - strncpy(name2,name,len-1); - name2[len-1] = '\0'; - name2[len-2] = '$'; - Z = (Float_t)z; - rad = GetRadLength(z)/dens; - AliMaterial(imat,name2,A,Z,dens,rad,0.0,0,0); - tmax = GetStandardTheataMax(istd); // degree - stemax = GetStandardMaxStepSize(istd); // cm - deemax = GetStandardEfraction(istd); // # - epsilon = GetStandardEpsilon(istd); - AliMedium(imat,name2,imat,0,gAlice->Field()->Integ(), - gAlice->Field()->Max(),tmax,stemax,deemax,epsilon,0.0); - delete[] name2; -} -//______________________________________________________________________ -void AliITSv11::SSDCone(Double_t zShift){ - // Defines the volumes and materials for the ITS SSD Support cone. - // Based on drawings ALR-0767 and ALR-0767/3. Units are in mm. - // Inputs: - // Double_t zShift The z shift to be applied to the final volume. + // none. // Outputs: // none. - // Return: + // Return // none. - Double_t *za,*rmina,*rmaxa,phi0=0.0,dphi=360.0; - Int_t i,n,nz,nrad=0; - Double_t Cthick=1.5; //mm, Carbon finber thickness - Double_t r=15.0; // mm, Radius of curvature. - Double_t tc=51.0; // angle of SSD cone [degrees]. - Double_t t; // some general angle [degrees]. - Int_t SSDcf=; // SSD support cone Carbon Fiber materal number. - - SetScalemm(); - // Lets start with the upper left outer carbon fiber surface. - nz = 6; - za = new Double_t[nz]; - rmina = new Double_t[nz]; - rmaxa = new Double_t[nz]; - za[0] = 0.0; - rmaxa[0] = 985./2.; - rmina[0] = rmaxa[0] - Cthick; - za[1] = 13.5 - 5.0; // The original size of 13.5 (ALR-0767) is milled down - rmaxa[1] = rmaxa[0]; // by 5mm to give a well defined reference surface - rmina[1] = rmina[0]; // (ALR-0767/3). - // The curved section is given by the following fomula - // rmax = r*Sind(90.-t)+rmaxa[0]-r for 0<=t<=tc. - // rmin = (r-Cthick)*Sind(90.-t)+rmina[0]-(r-Cthick) for 0<=t<=tc. - // z = r*Cosd(90.-t)+za[1] for 0<=t<=tc. - za[2] = (r-Cthick)*Cosd(90.-tc)+za[1]; - rmina[2] = (r-Cthick)*Sind(90.0-tc)+rmaxa[0]-(r-Cthick); - t = 90.0 - ACosd((za[2]-za[1])/r); // angle for rmax - rmaxa[2] = r*Sind(90.-t)+rmaxa[0]-r; - rmaxa[3] = r*Sind(90.0-tc)+rmaxa[0]-r; - za[3] =r*Cosd(90.-tc)+za[1]; - // angled section. surface is given by the following equations - // R = -Tand(tc)*Z + rmaxa[3] or rmina[2] - rmina[3] = -Tand(tc)*za[3] + rmina[2]; - // Point of whole. Whole surface has fixed radius = 890.0/2 mm - rmina[4] = 890.0/2.+Cthick ; // Inner whole surface radius (ALR-0767) - za[4] = (rmina[4] - rmina[2])/(-Tand(tc)); - rmaxa[4] = -Tand(tc)*za[4]+rmaxa[3]; + //TGeoMaterial *C = new TGeoMaterial("ITSCarbon",12.0,6.0,2.265); + TGeoMaterial *matAl = new TGeoMaterial("ITSAluminum",26.981539,13.0,2.07); + TGeoMixture *matCfiber = new TGeoMixture("ITSCarbonFiber",6,1.930); + TGeoMixture *matRohacell = new TGeoMixture("ITSRohacell",6,1.930); + TGeoMixture *matStaselite = new TGeoMixture("ITSStaselite4411w",6,1.930); + TGeoMixture *matAir = new TGeoMixture("ITSAir",6,1.205*1.E-3); + TGeoMixture *matStainless = new TGeoMixture("ITSStainless",6,1.930); // - rmaxa[5] = rmina[4]; - rmina[5] = rmina[4]; - za[5] = (rmaxa[5] - rmaxa[3])/(-Tand(tc)); - PolyCone("SCAA","SSD Suport cone Carbon Fiber Surface outer left", - phi0,dphi,nz,*z,*rmin,*rmax,SSDcf); - Za[0] = 1.; Wa[0] = ; // Hydrogen Content - Za[1] = 6.; Wa[0] = ; // Carbon Content - MixtureByWeight(SSDcf,"Carbon Fiber for SSD support cone",Z,W,dens,3); - // Now for the Filler in the mounting regions Upper left first - zb[0] = za[0]; - rmaxb[0] = rmina[0]; - rminb[0] = 945./2.-Cthick; + Double_t medSPDcone[20]; + medSPDcone[0] = 1.0; // imat + medSPDcone[1] = 0.0; // isvol + medSPDcone[2] = gAlice->Field()->Integ(); // ifield + medSPDcone[3] = gAlice->Field()->Max(); // fieldm + medSPDcone[4] = 1.0; // tmaxfd [degrees] + medSPDcone[5] = 1.0; // stemax [cm] + medSPDcone[6] = 0.5; // deemax [fraction] + medSPDcone[7] = 1.0E-3; // epsil [cm] + medSPDcone[8] = 0.0; // stmin [cm] + new TGeoMedium("ITSspdCarbonFiber",1,matCfiber,medSPDcone); + medSPDcone[0] += 1.0; + new TGeoMedium("ITSspdStaselite4411w",2,matStaselite,medSPDcone); + medSPDcone[0] += 1.0; + new TGeoMedium("ITSspdRohacell50A",3,matRohacell,medSPDcone); + medSPDcone[0] += 1.0; + new TGeoMedium("ITSspdStainlesSteal",4,matStainless,medSPDcone); + medSPDcone[0] += 1.0; + new TGeoMedium("ITSspdAir",5,matAir,medSPDcone); + medSPDcone[0] += 1.0; + new TGeoMedium("ITSspdAl",6,matAl,medSPDcone); // - rmaxb[1] = rmina[1]; - - delete[] za;delete[] rmina;delete[] rmaxa; - // Set back to cm default scale before exiting. - SetScalecm(); - return; -} -//______________________________________________________________________ -void AliITSv11::CreateGeometry(){ - //////////////////////////////////////////////////////////////////////// - // This routine defines and Creates the geometry for version 11 of the ITS. - //////////////////////////////////////////////////////////////////////// -} -//______________________________________________________________________ -void AliITSv11::CreateMaterials(){ -//////////////////////////////////////////////////////////////////////// - // - // Create ITS materials - // This function defines the default materials used in the Geant - // Monte Carlo simulations for the geometries AliITSv1, AliITSv3, - // AliITSv11. - // In general it is automatically replaced by - // the CreatMaterials routine defined in AliITSv?. Should the function - // CreateMaterials not exist for the geometry version you are using this - // one is used. See the definition found in AliITSv5 or the other routine - // for a complete definition. - // + Double_t medSSDcone[20]; + medSSDcone[0] = 1.0; // imat + medSSDcone[1] = 0.0; // isvol + medSSDcone[2] = gAlice->Field()->Integ(); // ifield + medSSDcone[3] = gAlice->Field()->Max(); // fieldm + medSSDcone[4] = 1.0; // tmaxfd [degrees] + medSSDcone[5] = 1.0; // stemax [cm] + medSSDcone[6] = 0.5; // deemax [fraction] + medSSDcone[7] = 1.0E-3; // epsil [cm] + medSSDcone[8] = 0.0; // stmin [cm] + new TGeoMedium("ITSssdCarbonFiber",1,matCfiber,medSSDcone); + medSSDcone[0] += 1.0; + new TGeoMedium("ITSssdStaselite4411w",2,matStaselite,medSSDcone); + medSSDcone[0] += 1.0; + new TGeoMedium("ITSssdRohacell50A",3,matRohacell,medSSDcone); + medSSDcone[0] += 1.0; + new TGeoMedium("ITSssdStainlesSteal",4,matStainless,medSSDcone); + medSSDcone[0] += 1.0; + new TGeoMedium("ITSssdAir",5,matAir,medSSDcone); + medSSDcone[0] += 1.0; + new TGeoMedium("ITSssdAl",6,matAl,medSSDcone); } //______________________________________________________________________ void AliITSv11::InitAliITSgeom(){ - // Based on the geometry tree defined in Geant 3.21, this - // routine initilizes the Class AliITSgeom from the Geant 3.21 ITS geometry - // sturture. + // Based on the geometry tree defined in Geant 3.21, this + // routine initilizes the Class AliITSgeom from the Geant 3.21 ITS + // geometry sturture. + // Inputs: + // none. + // Outputs: + // none. + // Return + // none. } - //______________________________________________________________________ void AliITSv11::Init(){ - //////////////////////////////////////////////////////////////////////// - // Initialise the ITS after it has been created. - //////////////////////////////////////////////////////////////////////// + // Initialise the ITS after it has been created. + // Inputs: + // none. + // Outputs: + // none. + // Return + // none. } //______________________________________________________________________ void AliITSv11::SetDefaults(){ - // sets the default segmentation, response, digit and raw cluster classes + // Sets the default segmentation, response, digit and raw cluster + // classes to be used. These defaults can be overwritten in the + // macros that do these later steps. Defaults are give hear for the + // general user. + // Inputs: + // none. + // Outputs: + // none. + // Return + // none. } //______________________________________________________________________ void AliITSv11::DrawModule(){ - //////////////////////////////////////////////////////////////////////// - // Draw a shaded view of the FMD version 11. - //////////////////////////////////////////////////////////////////////// + // Draw a standard set of shaded view of the ITS version 11. + // Inputs: + // none. + // Outputs: + // none. + // Return + // none. } //______________________________________________________________________ void AliITSv11::StepManager(){ - //////////////////////////////////////////////////////////////////////// - // Called for every step in the ITS, then calles the AliITShit class + // Called for every step in the ITS, then calles the AliITShit class // creator with the information to be recoreded about that hit. - // The value of the macro ALIITSPRINTGEOM if set to 1 will allow the + // The value of the macro ALIITSPRINTGEOM if set to 1 will allow the // printing of information to a file which can be used to create a .det // file read in by the routine CreateGeometry(). If set to 0 or any other // value except 1, the default behavior, then no such file is created nor - // it the extra variables and the like used in the printing allocated. - //////////////////////////////////////////////////////////////////////// + // is the extra variables and the like used in the printing allocated. }