/*
$Log$
+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 <iostream.h>
+#include <iomanip.h>
+#include <stdlib.h>
+#include <TShape.h>
+#include <TBRIK.h>
+
#include "AliITSgeomSSD.h"
ClassImp(AliITSgeomSSD)
+
+
AliITSgeomSSD::AliITSgeomSSD(){
- //
- // default constructor
- //
- fShapeSSD = new TBRIK("ActiveSSD","Active volume of SSD","SI",
- 3.0E-2/2.,1.73/2.,4.0/2.);
+// 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;i<fNp;i++) fLowEdgeP[i] = p[i];
+ for(i=0;i<fNn;i++) fLowEdgeN[i] = n[i];
+}
+//______________________________________________________________________
+AliITSgeomSSD::~AliITSgeomSSD(){
+ // Destructor.
+
+ if(fLowEdgeP) delete [] fLowEdgeP; fLowEdgeP = 0;
+ if(fLowEdgeN) delete [] fLowEdgeN; fLowEdgeN = 0;
+ if(fShapeSSD) delete fShapeSSD; fShapeSSD = 0;
+ fNp = 0;
+ fNn = 0;
+ fAngleP = 0.0;
+ fAngleN = 0.0;
+}
+AliITSgeomSSD::AliITSgeomSSD(const AliITSgeomSSD &source){
+////////////////////////////////////////////////////////////////////////
+// copy constructor
+////////////////////////////////////////////////////////////////////////
+ Int_t i;
+
+ if(this == &source) return;
+ 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;i<fNp;i++) this->fLowEdgeP[i] = source.fLowEdgeP[i];
+ for(i=0;i<fNn;i++) this->fLowEdgeN[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;i<fNp;i++) this->fLowEdgeP[i] = source.fLowEdgeP[i];
+ for(i=0;i<fNn;i++) this->fLowEdgeN[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<fLowEdgeP[0]) i=-1;
+ else for(i=0;i<fNp;i++){
+ if(fLowEdgeP[i]<d) break;
+ } // end for i
+ a = i;
+ if(b<fLowEdgeN[0]) i=-1;
+ else for(i=0;i<fNn;i++){
+ if(fLowEdgeN[i]<b) break;
+ } // end for i
+ c = i;
+ return;
+}
+//______________________________________________________________________
+void AliITSgeomSSD::Det2Local(Int_t a,Int_t c,Float_t &x,Float_t &z){
+// Float_t d,b;
+// Int_t i;
+
+ return;
+}
+//______________________________________________________________________
+void AliITSgeomSSD::Print(ostream *os) const {
+////////////////////////////////////////////////////////////////////////
+// Standard output format for this class.
+////////////////////////////////////////////////////////////////////////
+ Int_t i;
+#if defined __GNUC__
+#if __GNUC__ > 2
+ ios::fmtflags fmt;
+#else
+ Int_t fmt;
+#endif
+#else
+#if defined __ICC
+ 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;i<fNp;i++) *os << setprecision(16) << fLowEdgeP[i] << " ";
+ for(i=0;i<fNn;i++) *os << setprecision(16) << fLowEdgeN[i] << " ";
+ *os << endl;
+ os->flags(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<fNp;i++) *is >> fLowEdgeP[i];
+ for(i=0;i<fNn;i++) *is >> 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.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<kNstrips;i++){
+ leA[i+1] = leA[i] + kpitch;
+ leC[i+1] = leC[i] - kpitch;
+ } // end for i
+ leA[kNstrips] = kDxyz[0];
+ leC[kNstrips] = -kDxyz[0];
+// cout << "AliITSgeomSSD175 default creator called: start" << endl;
+ AliITSgeomSSD::ResetSSD(kDxyz,kangle,-kangle,
+ kNstrips+1,leA,kNstrips+1,leC);
+ delete leA;
+ delete leC;
+// cout << "AliITSgeomSSD175 default creator called: end" << endl;
+}
+//________________________________________________________________________
+ostream &operator<<(ostream &os,AliITSgeomSSD175 &p){
+////////////////////////////////////////////////////////////////////////
+// Standard output streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ p.Print(&os);
+ return os;
+}
+//----------------------------------------------------------------------
+istream &operator>>(istream &is,AliITSgeomSSD175 &r){
+////////////////////////////////////////////////////////////////////////
+// Standard input streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ r.Read(&is);
+ return is;
+}
+//======================================================================
+/*
+$Log$
+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<kNstrips;i++){
+ leA[i+1] = leA[i] + kpitch;
+ leC[i+1] = leC[i] - kpitch;
+ } // end for i
+ leA[kNstrips] = kDxyz[0];
+ leC[kNstrips] = -kDxyz[0];
+// cout << "AliITSgeomSSD275and75 default creator called: start" << endl;
+ AliITSgeomSSD::ResetSSD(par,kangleA,kangleC,
+ kNstrips+1,leA,kNstrips+1,leC);
+ delete leA;
+ delete leC;
+// cout << "AliITSgeomSSD275and75 default creator called: end" << endl;
+}
+//________________________________________________________________________
+ostream &operator<<(ostream &os,AliITSgeomSSD275and75 &p){
+////////////////////////////////////////////////////////////////////////
+// Standard output streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ p.Print(&os);
+ return os;
+}
+//----------------------------------------------------------------------
+istream &operator>>(istream &is,AliITSgeomSSD275and75 &r){
+////////////////////////////////////////////////////////////////////////
+// Standard input streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ r.Read(&is);
+ return is;
+}
+//======================================================================
+/*
+$Log$
+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<kNstrips;i++){
+ leA[i+1] = leA[i] + kpitch;
+ leC[i+1] = leC[i] - kpitch;
+ } // end for i
+ leA[kNstrips] = kDxyz[0];
+ leC[kNstrips] = -kDxyz[0];
+// cout << "AliITSgeomSSD275and75 default creator called: start" << endl;
+ AliITSgeomSSD::ResetSSD(par,kangleA,kangleC,
+ kNstrips+1,leA,kNstrips+1,leC);
+ delete leA;
+ delete leC;
+// cout << "AliITSgeomSSD275and75 default creator called: end" << endl;
+}
+//________________________________________________________________________
+ostream &operator<<(ostream &os,AliITSgeomSSD75and275 &p){
+////////////////////////////////////////////////////////////////////////
+// Standard output streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ p.Print(&os);
+ return os;
+}
+//----------------------------------------------------------------------
+istream &operator>>(istream &is,AliITSgeomSSD75and275 &r){
+////////////////////////////////////////////////////////////////////////
+// Standard input streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ r.Read(&is);
+ return is;
}
+//======================================================================
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff