* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.8 2000/10/02 16:32:43 barbera
-Forward declaration added
+/* $Id$ */
-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 <Riostream.h>
#include <TBRIK.h>
+#include <TMath.h>
+#include <TShape.h>
#include "AliITSgeomSSD.h"
ClassImp(AliITSgeomSSD)
+
+
+AliITSgeomSSD::AliITSgeomSSD():
+TObject(),
+fName(),
+fTitle(),
+fMat(),
+fDx(0.0),
+fDy(0.0),
+fDz(0.0),
+fNp(0),
+fNn(0),
+fLowEdgeP(0),
+fLowEdgeN(0),
+fAngleP(0.0),
+fAngleN(0.0){
+// Default constructor
+ 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){
+ Int_t np,Float_t *p,Int_t nn,Float_t *n):
+TObject(),
+fName(),
+fTitle(),
+fMat(),
+fDx(0.0),
+fDy(0.0),
+fDz(0.0),
+fNp(0),
+fNn(0),
+fLowEdgeP(0),
+fLowEdgeN(0),
+fAngleP(0.0),
+fAngleN(0.0){
////////////////////////////////////////////////////////////////////////
// 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.
+///////////////////////////////////////////////////////////////////////
+ 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]);
+ fName = "ActiveSSD";
+ fTitle = "Active volume of SSD";
+ fMat = "SSD Si Det";
+ fDx = box[0];
+ fDy = box[1];
+ fDz = box[2];
if(fLowEdgeP!=0) delete fLowEdgeP;
if(fLowEdgeN!=0) delete fLowEdgeN;
fNp = np;
AliITSgeomSSD::~AliITSgeomSSD(){
// Destructor.
- delete fLowEdgeP; fLowEdgeP = 0;
- delete fLowEdgeN; fLowEdgeN = 0;
- delete fShapeSSD; fShapeSSD = 0;
+ if(fLowEdgeP) delete [] fLowEdgeP; fLowEdgeP = 0;
+ if(fLowEdgeN) delete [] fLowEdgeN; fLowEdgeN = 0;
fNp = 0;
fNn = 0;
fAngleP = 0.0;
fAngleN = 0.0;
}
-AliITSgeomSSD::AliITSgeomSSD(const AliITSgeomSSD &source){
+//______________________________________________________________________
+AliITSgeomSSD::AliITSgeomSSD(const AliITSgeomSSD &source) : TObject(source),
+fName(source.fName),fTitle(source.fTitle),fMat(source.fMat),fDx(source.fDx),fDy(source.fDy),fDz(source.fDz),fNp(source.fNp),fNn(source.fNn),fLowEdgeP(0),fLowEdgeN(0),fAngleP(source.fAngleP),fAngleN(source.fAngleN){
////////////////////////////////////////////////////////////////////////
// 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];
////////////////////////////////////////////////////////////////////////
// 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;
+ this->~AliITSgeomSSD();
+ new(this) AliITSgeomSSD(source);
+ 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;
void AliITSgeomSSD::Det2Local(Int_t a,Int_t c,Float_t &x,Float_t &z){
// Float_t d,b;
// Int_t i;
+ // use AliITSsegmentationSSD.
+ x=a;
+ z=c;
+ Error("Det2Locat","Use AliITSsegmentationSSD");
return;
}
//______________________________________________________________________
-void AliITSgeomSSD::Print(ostream *os){
+void AliITSgeomSSD::Print(ostream *os) const {
////////////////////////////////////////////////////////////////////////
// Standard output format for this class.
////////////////////////////////////////////////////////////////////////
- ios::fmtflags fmt;
Int_t i;
+#if defined __GNUC__
+#if __GNUC__ > 2
+ ios::fmtflags fmt;
+#else
+ Int_t fmt;
+#endif
+#else
+#if defined __ICC || defined __ECC || defined __xlC__
+ ios::fmtflags fmt;
+#else
+ Int_t fmt;
+#endif
+#endif
fmt = os->setf(ios::scientific); // set scientific floating point output
*os << "TBRIK" << " ";
////////////////////////////////////////////////////////////////////////
// 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 >> fDx >> fDy >> fDz;
*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;0<fNp;i++) *is >> fLowEdgeP[i];
+ for(i=0;i<fNp;i++) *is >> fLowEdgeP[i];
for(i=0;i<fNn;i++) *is >> fLowEdgeN[i];
return;
}
return is;
}
//======================================================================
-/*
-$Log$
-*/
-
-//#include "AliITSgeomSSD175.h"
ClassImp(AliITSgeomSSD175)
} // end for i
leA[kNstrips] = kDxyz[0];
leC[kNstrips] = -kDxyz[0];
-// cout << "AliITSgeomSSD175 default creator called: start" << endl;
- AliITSgeomSSD::AliITSgeomSSD(kDxyz,kangle,-kangle,
+ 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){
r.Read(&is);
return is;
}
-//======================================================================
-/*
-$Log$
-*/
+AliITSgeomSSD175& AliITSgeomSSD175::operator=(const AliITSgeomSSD175 &source) {
+////////////////////////////////////////////////////////////////////////
+// assignment operator
+////////////////////////////////////////////////////////////////////////
+
-//#include "AliITSgeomSSD275and75.h"
+ if(this == &source) return *this;
+ Error("AliITSgeomSSD175","Not allowed to make a = with "
+ "AliITSgeomSSD175 Using default creater instead");
+
+ return *this;
+}
+//======================================================================
ClassImp(AliITSgeomSSD275and75)
////////////////////////////////////////////////////////////////////////
// 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
Float_t *leA,*leC; // array of low edges anode and cathorde.
Int_t i;
+ if(npar<3){
+ Error("AliITSgeomSSD275and75",
+ "npar=%d<3. array par must be [3] or larger.",npar);
+ return;
+ } // end if
leA = new Float_t[kNstrips+1];
leC = new Float_t[kNstrips+1];
leA[0] = -kDxyz[0];
} // end for i
leA[kNstrips] = kDxyz[0];
leC[kNstrips] = -kDxyz[0];
-// cout << "AliITSgeomSSD275and75 default creator called: start" << endl;
- AliITSgeomSSD::AliITSgeomSSD(kDxyz,kangleA,kangleC,
+ AliITSgeomSSD::ResetSSD(par,kangleA,kangleC,
kNstrips+1,leA,kNstrips+1,leC);
- delete leA;
- delete leC;
-// cout << "AliITSgeomSSD275and75 default creator called: end" << endl;
+ delete [] leA;
+ delete [] leC;
}
+
//________________________________________________________________________
ostream &operator<<(ostream &os,AliITSgeomSSD275and75 &p){
////////////////////////////////////////////////////////////////////////
r.Read(&is);
return is;
}
+AliITSgeomSSD275and75& AliITSgeomSSD275and75::operator=(const AliITSgeomSSD275and75 &source) {
+////////////////////////////////////////////////////////////////////////
+// assignment operator
+////////////////////////////////////////////////////////////////////////
+
+
+ if(this == &source) return *this;
+ Error("AliITSgeomSSD275and75","Not allowed to make a = with "
+ "AliITSgeomSSD275and75 Using default creater instead");
+
+ return *this;
+}
+
//======================================================================
-/*
-$Log$
-*/
-//#include "AliITSgeomSSD75and275.h"
ClassImp(AliITSgeomSSD75and275)
////////////////////////////////////////////////////////////////////////
// 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
Float_t *leA,*leC; // array of low edges anode and cathorde.
Int_t i;
+ if(npar<3){
+ Error("AliITSgeomSSD75and275",
+ "npar=%d<3. array par must be [3] or larger.",npar);
+ return;
+ } // end if
leA = new Float_t[kNstrips+1];
leC = new Float_t[kNstrips+1];
leA[0] = -kDxyz[0];
} // end for i
leA[kNstrips] = kDxyz[0];
leC[kNstrips] = -kDxyz[0];
-// cout << "AliITSgeomSSD275and75 default creator called: start" << endl;
- AliITSgeomSSD::AliITSgeomSSD(kDxyz,kangleA,kangleC,
+ 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){
r.Read(&is);
return is;
}
+AliITSgeomSSD75and275& AliITSgeomSSD75and275::operator=(const AliITSgeomSSD75and275 &source) {
+////////////////////////////////////////////////////////////////////////
+// assignment operator
+////////////////////////////////////////////////////////////////////////
+
+
+ if(this == &source) return *this;
+ Error("AliITSgeomSSD75and275","Not allowed to make a = with "
+ "AliITSgeomSSD75and275 Using default creater instead");
+
+ return *this;
+}
+
//======================================================================