/************************************************************************** * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * * * Author: The ALICE Off-line Project. * * Contributors are mentioned in the code where appropriate. * * * * Permission to use, copy, modify and distribute this software and its * * documentation strictly for non-commercial purposes is hereby granted * * without fee, provided that the above copyright notice appears in all * * copies and that both the copyright notice and this permission notice * * appear in the supporting documentation. The authors make no claims * * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ /* $Log$ Revision 1.18 2002/12/03 09:03:06 hristov Changes needed on Itanium (F.Carminati) Revision 1.17 2002/10/22 14:45:41 alibrary Introducing Riostream.h Revision 1.16 2002/10/14 14:57:00 hristov Merging the VirtualMC branch to the main development branch (HEAD) Revision 1.14.6.1 2002/06/10 17:51:15 hristov Merged with v3-08-02 Revision 1.15 2002/05/19 18:17:03 hristov Changes needed by ICC/IFC compiler (Intel) Revision 1.14 2001/11/19 16:17:03 nilsen Applyed fixes to bugs found by Rene Brun. With many thanks. Some additonal bugs found by Rene require more work to fix. Will be fixed soon. Revision 1.13 2001/10/12 22:07:20 nilsen A patch for C++ io manipulation functions so that they will work both with GNU gcc 2.96 and GNU gcc 3.01 compilers. Needs to be tested with other platforms. Revision 1.12 2001/08/24 21:06:37 nilsen Added more documentation, fixed up some coding violations, and some forward declorations. Revision 1.11 2001/05/16 08:17:49 hristov Bug fixed in the StepManager to account for the difference in the geometry tree for the ITS pixels. This fixes both the funny distribution of pixel coordinates and the missing hits/digits/points in many sectors of the ITS pixel barrel. Also included is a patch to properly get and use the detector dimensions through out the ITS code. (B.Nilsen) Revision 1.10 2001/02/09 00:00:57 nilsen Fixed compatibility problem with HP unix {ios::fmtflags -> Int_t}. Fixed bugs in iostream based streamers used to read and write .det files. Fixed some detector sizes. Fixed bugs in some default-special constructors. Revision 1.9 2001/02/03 00:00:30 nilsen New version of AliITSgeom and related files. Now uses automatic streamers, set up for new formatted .det file which includes detector information. Additional smaller modifications are still to come. Revision 1.8 2000/10/02 16:32:43 barbera Forward declaration added Revision 1.2.4.8 2000/10/02 15:53:49 barbera Forward declaration added Revision 1.7 2000/07/10 16:07:18 fca Release version of ITS code Revision 1.2.4.2 2000/03/04 23:55:59 nilsen Fixed up the comments/documentation Revision 1.2.4.1 2000/01/12 19:03:32 nilsen This is the version of the files after the merging done in December 1999. See the ReadMe110100.txt file for details Revision 1.2 1999/09/29 09:24:20 fca Introduction of the Copyright and cvs Log */ //////////////////////////////////////////////////////////////////////// // This class is for the Silicon Strip Detector, SSD, specific geometry. // It is being replaced by AliITSsegmentationSSD class. This file also // constains classes derived from AliITSgeomSSD which do nothing but // initilize this one with predefined values. //////////////////////////////////////////////////////////////////////// #include #include #include #include #include "AliITSgeomSSD.h" ClassImp(AliITSgeomSSD) AliITSgeomSSD::AliITSgeomSSD(){ // Default constructor fShapeSSD = 0; fNp = 0; fNn = 0; fLowEdgeP = 0; fLowEdgeN = 0; fAngleP = 0.0; fAngleN = 0.0; } //---------------------------------------------------------------------- AliITSgeomSSD::AliITSgeomSSD(const Float_t *box,Float_t ap,Float_t an, Int_t np,Float_t *p,Int_t nn,Float_t *n){ //////////////////////////////////////////////////////////////////////// // Standard Constructor. *box={dx,dy,dz}, ap=anode angle, an=cathode angle, // nn= number of cathodes+1,*n= array of cathode low edges+highest edge, // np= number of anodes+1, *p= array of anode low edges+lighest edge. /////////////////////////////////////////////////////////////////////// fShapeSSD = 0; fNp = 0; fNn = 0; fLowEdgeP = 0; fLowEdgeN = 0; fAngleP = 0.0; fAngleN = 0.0; ResetSSD(box,ap,an,np,p,nn,n); } //---------------------------------------------------------------------- void AliITSgeomSSD::ResetSSD(const Float_t *box,Float_t ap,Float_t an, Int_t np,Float_t *p,Int_t nn,Float_t *n){ //////////////////////////////////////////////////////////////////////// // Standard Filler. *box={dx,dy,dz}, ap=anode angle, an=cathode angle, // nn= number of cathodes+1,*n= array of cathode low edges+highest edge, // np= number of anodes+1, *p= array of anode low edges+lighest edge. /////////////////////////////////////////////////////////////////////// Int_t i; fShapeSSD = new TBRIK("ActiveSSD","Active volume of SSD","SSD SI DET", box[0],box[1],box[2]); // if(fLowEdgeP!=0) delete fLowEdgeP; // if(fLowEdgeN!=0) delete fLowEdgeN; fNp = np; fNn = nn; fAngleP = ap; fAngleN = an; fLowEdgeP = new Float_t[fNp]; fLowEdgeN = new Float_t[fNn]; for(i=0;ifShapeSSD = new TBRIK(*(source.fShapeSSD)); this->fNp = source.fNp; this->fNn = source.fNn; delete fLowEdgeP; delete fLowEdgeN; this->fAngleP = source.fAngleP; this->fAngleN = source.fAngleN; fLowEdgeP = new Float_t[fNp]; fLowEdgeN = new Float_t[fNn]; for(i=0;ifLowEdgeP[i] = source.fLowEdgeP[i]; for(i=0;ifLowEdgeN[i] = source.fLowEdgeN[i]; return; } AliITSgeomSSD& AliITSgeomSSD::operator=(const AliITSgeomSSD &source) { //////////////////////////////////////////////////////////////////////// // assignment operator //////////////////////////////////////////////////////////////////////// Int_t i; if(this == &source) return *this; this->fShapeSSD = new TBRIK(*(source.fShapeSSD)); this->fNp = source.fNp; this->fNn = source.fNn; delete fLowEdgeP; delete fLowEdgeN; this->fAngleP = source.fAngleP; this->fAngleN = source.fAngleN; fLowEdgeP = new Float_t[fNp]; fLowEdgeN = new Float_t[fNn]; for(i=0;ifLowEdgeP[i] = source.fLowEdgeP[i]; for(i=0;ifLowEdgeN[i] = source.fLowEdgeN[i]; return *this; } //______________________________________________________________________ void AliITSgeomSSD::Local2Det(Float_t x,Float_t z,Int_t &a,Int_t &c){ // Given a GEANT detector local coordinate, cm, this function returns // the detector specific P and N side strip numbers. // Inputs are: // Float_t x Geant detector local x coordinate in cm // Float_t z Geant detector local z coordinate in cm // outputs are: // Int_t &a Detector anode strip number (P side) // Int_t &c Detector cathode strip number (N side) Float_t d,b; Int_t i; // project on to bonding edges. d = x*TMath::Cos(fAngleP)+z*TMath::Sin(fAngleP); b = x*TMath::Cos(fAngleN)+z*TMath::Sin(fAngleN); if(d 2 ios::fmtflags fmt; #else Int_t fmt; #endif #else #if defined __ICC || defined __ECC ios::fmtflags fmt; #else Int_t fmt; #endif #endif fmt = os->setf(ios::scientific); // set scientific floating point output *os << "TBRIK" << " "; *os << setprecision(16) << GetDx() << " "; *os << setprecision(16) << GetDy() << " "; *os << setprecision(16) << GetDz() << " "; *os << fNp << " " << fNn << " "; *os << setprecision(16) << fAngleP << " "; *os << setprecision(16) << fAngleN << " "; for(i=0;iflags(fmt); // reset back to old formating. return; } //______________________________________________________________________ void AliITSgeomSSD::Read(istream *is){ //////////////////////////////////////////////////////////////////////// // Standard input format for this class. //////////////////////////////////////////////////////////////////////// Float_t dx,dy,dz; Int_t i; char shp[20]; *is >> shp; *is >> dx >> dy >> dz; if(fShapeSSD!=0) delete fShapeSSD; fShapeSSD = new TBRIK("ActiveSSD","Active volume of SSD","SSD SI DET", dx,dy,dz); *is >> fNp >> fNn; *is >> fAngleP >> fAngleN; if(fLowEdgeP !=0) delete fLowEdgeP; if(fLowEdgeN !=0) delete fLowEdgeN; fLowEdgeP = new Float_t[fNp]; fLowEdgeN = new Float_t[fNn]; for(i=0;i> fLowEdgeP[i]; for(i=0;i> fLowEdgeN[i]; return; } //---------------------------------------------------------------------- ostream &operator<<(ostream &os,AliITSgeomSSD &p){ //////////////////////////////////////////////////////////////////////// // Standard output streaming function. //////////////////////////////////////////////////////////////////////// p.Print(&os); return os; } //---------------------------------------------------------------------- istream &operator>>(istream &is,AliITSgeomSSD &r){ //////////////////////////////////////////////////////////////////////// // Standard input streaming function. //////////////////////////////////////////////////////////////////////// r.Read(&is); return is; } //====================================================================== /* $Log$ Revision 1.18 2002/12/03 09:03:06 hristov Changes needed on Itanium (F.Carminati) Revision 1.17 2002/10/22 14:45:41 alibrary Introducing Riostream.h Revision 1.16 2002/10/14 14:57:00 hristov Merging the VirtualMC branch to the main development branch (HEAD) Revision 1.14.6.1 2002/06/10 17:51:15 hristov Merged with v3-08-02 Revision 1.15 2002/05/19 18:17:03 hristov Changes needed by ICC/IFC compiler (Intel) Revision 1.14 2001/11/19 16:17:03 nilsen Applyed fixes to bugs found by Rene Brun. With many thanks. Some additonal bugs found by Rene require more work to fix. Will be fixed soon. Revision 1.13 2001/10/12 22:07:20 nilsen A patch for C++ io manipulation functions so that they will work both with GNU gcc 2.96 and GNU gcc 3.01 compilers. Needs to be tested with other platforms. Revision 1.12 2001/08/24 21:06:37 nilsen Added more documentation, fixed up some coding violations, and some forward declorations. Revision 1.11 2001/05/16 08:17:49 hristov Bug fixed in the StepManager to account for the difference in the geometry tree for the ITS pixels. This fixes both the funny distribution of pixel coordinates and the missing hits/digits/points in many sectors of the ITS pixel barrel. Also included is a patch to properly get and use the detector dimensions through out the ITS code. (B.Nilsen) Revision 1.10 2001/02/09 00:00:57 nilsen Fixed compatibility problem with HP unix {ios::fmtflags -> Int_t}. Fixed bugs in iostream based streamers used to read and write .det files. Fixed some detector sizes. Fixed bugs in some default-special constructors. Revision 1.9 2001/02/03 00:00:30 nilsen New version of AliITSgeom and related files. Now uses automatic streamers, set up for new formatted .det file which includes detector information. Additional smaller modifications are still to come. */ //#include "AliITSgeomSSD175.h" ClassImp(AliITSgeomSSD175) AliITSgeomSSD175::AliITSgeomSSD175() : AliITSgeomSSD(){ //////////////////////////////////////////////////////////////////////// // default constructor //////////////////////////////////////////////////////////////////////// const Float_t kDxyz[] ={3.6500,0.0150,2.000};//cm. (Geant 3.12 units) // Size of sensitive detector area x,y(thickness),z const Float_t kangle = 0.0175; // angle in rad. of anode and cathodes const Float_t kpitch = 0.0095;// cm anode separation. const Int_t kNstrips = 768; // number of anode or cathode strips. Float_t *leA,*leC; // array of low edges anode and cathorde. Int_t i; leA = new Float_t[kNstrips+1]; leC = new Float_t[kNstrips+1]; leA[0] = -kDxyz[0]; leA[1] = -kpitch*(0.5*kNstrips-1); leC[0] = kDxyz[0]; leC[1] = kpitch*(0.5*kNstrips-1); for(i=1;i>(istream &is,AliITSgeomSSD175 &r){ //////////////////////////////////////////////////////////////////////// // Standard input streaming function. //////////////////////////////////////////////////////////////////////// r.Read(&is); return is; } //====================================================================== /* $Log$ Revision 1.18 2002/12/03 09:03:06 hristov Changes needed on Itanium (F.Carminati) Revision 1.17 2002/10/22 14:45:41 alibrary Introducing Riostream.h Revision 1.16 2002/10/14 14:57:00 hristov Merging the VirtualMC branch to the main development branch (HEAD) Revision 1.14.6.1 2002/06/10 17:51:15 hristov Merged with v3-08-02 Revision 1.15 2002/05/19 18:17:03 hristov Changes needed by ICC/IFC compiler (Intel) Revision 1.14 2001/11/19 16:17:03 nilsen Applyed fixes to bugs found by Rene Brun. With many thanks. Some additonal bugs found by Rene require more work to fix. Will be fixed soon. Revision 1.13 2001/10/12 22:07:20 nilsen A patch for C++ io manipulation functions so that they will work both with GNU gcc 2.96 and GNU gcc 3.01 compilers. Needs to be tested with other platforms. Revision 1.12 2001/08/24 21:06:37 nilsen Added more documentation, fixed up some coding violations, and some forward declorations. Revision 1.11 2001/05/16 08:17:49 hristov Bug fixed in the StepManager to account for the difference in the geometry tree for the ITS pixels. This fixes both the funny distribution of pixel coordinates and the missing hits/digits/points in many sectors of the ITS pixel barrel. Also included is a patch to properly get and use the detector dimensions through out the ITS code. (B.Nilsen) Revision 1.10 2001/02/09 00:00:57 nilsen Fixed compatibility problem with HP unix {ios::fmtflags -> Int_t}. Fixed bugs in iostream based streamers used to read and write .det files. Fixed some detector sizes. Fixed bugs in some default-special constructors. Revision 1.9 2001/02/03 00:00:30 nilsen New version of AliITSgeom and related files. Now uses automatic streamers, set up for new formatted .det file which includes detector information. Additional smaller modifications are still to come. */ //#include "AliITSgeomSSD275and75.h" ClassImp(AliITSgeomSSD275and75) AliITSgeomSSD275and75::AliITSgeomSSD275and75() : AliITSgeomSSD(){ //////////////////////////////////////////////////////////////////////// // default constructor //////////////////////////////////////////////////////////////////////// } //---------------------------------------------------------------------- AliITSgeomSSD275and75::AliITSgeomSSD275and75(Int_t npar,Float_t *par) : AliITSgeomSSD(){ // Default constructor for AliITSgeomSSD with strip angles of // 275 miliradians and 75 miliradians. This constructor initlizes // AliITSgeomSSD with the correct values. This is the miror image // of the AliITSgeomSSD75and275 class. const Float_t kDxyz[] ={3.6500,0.0150,2.000};//cm. (Geant 3.12 units) // Size of sensitive detector area x,y(thickness),z const Float_t kangleA = 0.0275; // angle in rad. of anode and cathodes const Float_t kangleC = 0.0075; // angle in rad. of anode and cathodes const Float_t kpitch = 0.0095;// cm anode separation. const Int_t kNstrips = 768; // number of anode or cathode strips. Float_t *leA,*leC; // array of low edges anode and cathorde. Int_t i; leA = new Float_t[kNstrips+1]; leC = new Float_t[kNstrips+1]; leA[0] = -kDxyz[0]; leA[1] = -kpitch*(0.5*kNstrips-1); leC[0] = kDxyz[0]; leC[1] = kpitch*(0.5*kNstrips-1); for(i=1;i>(istream &is,AliITSgeomSSD275and75 &r){ //////////////////////////////////////////////////////////////////////// // Standard input streaming function. //////////////////////////////////////////////////////////////////////// r.Read(&is); return is; } //====================================================================== /* $Log$ Revision 1.18 2002/12/03 09:03:06 hristov Changes needed on Itanium (F.Carminati) Revision 1.17 2002/10/22 14:45:41 alibrary Introducing Riostream.h Revision 1.16 2002/10/14 14:57:00 hristov Merging the VirtualMC branch to the main development branch (HEAD) Revision 1.14.6.1 2002/06/10 17:51:15 hristov Merged with v3-08-02 Revision 1.15 2002/05/19 18:17:03 hristov Changes needed by ICC/IFC compiler (Intel) Revision 1.14 2001/11/19 16:17:03 nilsen Applyed fixes to bugs found by Rene Brun. With many thanks. Some additonal bugs found by Rene require more work to fix. Will be fixed soon. Revision 1.13 2001/10/12 22:07:20 nilsen A patch for C++ io manipulation functions so that they will work both with GNU gcc 2.96 and GNU gcc 3.01 compilers. Needs to be tested with other platforms. Revision 1.12 2001/08/24 21:06:37 nilsen Added more documentation, fixed up some coding violations, and some forward declorations. Revision 1.11 2001/05/16 08:17:49 hristov Bug fixed in the StepManager to account for the difference in the geometry tree for the ITS pixels. This fixes both the funny distribution of pixel coordinates and the missing hits/digits/points in many sectors of the ITS pixel barrel. Also included is a patch to properly get and use the detector dimensions through out the ITS code. (B.Nilsen) Revision 1.10 2001/02/09 00:00:57 nilsen Fixed compatibility problem with HP unix {ios::fmtflags -> Int_t}. Fixed bugs in iostream based streamers used to read and write .det files. Fixed some detector sizes. Fixed bugs in some default-special constructors. Revision 1.9 2001/02/03 00:00:30 nilsen New version of AliITSgeom and related files. Now uses automatic streamers, set up for new formatted .det file which includes detector information. Additional smaller modifications are still to come. */ //#include "AliITSgeomSSD75and275.h" ClassImp(AliITSgeomSSD75and275) AliITSgeomSSD75and275::AliITSgeomSSD75and275() : AliITSgeomSSD(){ //////////////////////////////////////////////////////////////////////// // default constructor //////////////////////////////////////////////////////////////////////// } AliITSgeomSSD75and275::AliITSgeomSSD75and275(Int_t npar,Float_t *par) : AliITSgeomSSD(){ // Default constructor for AliITSgeomSSD with strip angles of // 75 miliradians and 275 miliradians. This constructor initlizes // AliITSgeomSSD with the correct values. This is the miror image // of the AliITSgeomSSD275and75 class. const Float_t kDxyz[] ={3.6500,0.0150,2.000};//cm. (Geant 3.12 units) // Size of sensitive detector area x,y(thickness),z const Float_t kangleA = 0.0075; // angle in rad. of anode and cathodes const Float_t kangleC = 0.0275; // angle in rad. of anode and cathodes const Float_t kpitch = 0.0095;// cm anode separation. const Int_t kNstrips = 768; // number of anode or cathode strips. Float_t *leA,*leC; // array of low edges anode and cathorde. Int_t i; leA = new Float_t[kNstrips+1]; leC = new Float_t[kNstrips+1]; leA[0] = -kDxyz[0]; leA[1] = -kpitch*(0.5*kNstrips-1); leC[0] = kDxyz[0]; leC[1] = kpitch*(0.5*kNstrips-1); for(i=1;i>(istream &is,AliITSgeomSSD75and275 &r){ //////////////////////////////////////////////////////////////////////// // Standard input streaming function. //////////////////////////////////////////////////////////////////////// r.Read(&is); return is; } //======================================================================