[u/mrichter/AliRoot.git] / ITS / AliITSgeomSPD425.cxx

/**************************************************************************
 * 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.3  2000/06/10 10:42:49  nilsen
Fixed bug in copy and operator =.


*/

#include <TShape.h>

#include "AliITSgeomSPD425.h"

ClassImp(AliITSgeomSPD425)

AliITSgeomSPD425::AliITSgeomSPD425(){
////////////////////////////////////////////////////////////////////////
//    default constructor, for ITS post TDR geometry.
////////////////////////////////////////////////////////////////////////

const Float_t kdx=0.6400,kdy=0.0075,kdz=4.23625;// cm; Standard pixel detector
                                                // size is 2dx wide, 2dz long,
                                                // and 2dy thick. Geant 3.12
                                                // units.
const Float_t kbinx0 = 0.0050; // cm; Standard pixel size in x direction.
const Int_t   knbinx = 256;    // number of pixels along x direction.
const Float_t kbinz0 = 0.0425; // cm; Standard pixel size in z direction.
const Float_t kbinz1 = 0.0625; // cm; Special pixel size in z direction.
const Int_t   knbinz = 197;    // number of pixels along z direction.
    Int_t i;

    fdx = kdx;  // default value.
    fdy = kdy;  // default value.
    fdz = kdz;  // default value.
    fNbinx = knbinx; // default number of bins in x.
    fNbinz = knbinz; // default number of bins in z.

    fBinSizeX = new Float_t[fNbinx]; // array of bin sizes along x.
    for(i=0;i<fNbinx;i++) fBinSizeX[i] = kbinx0; // default x bin size.
    fBinSizeZ = new Float_t[fNbinz]; // array of bin sizes along z.
    for(i=0;i<fNbinz;i++) fBinSizeZ[i] = kbinz0; // default z bin size.
    fBinSizeZ[ 31] = kbinz1;
    fBinSizeZ[ 32] = kbinz1;

    fBinSizeZ[ 64] = kbinz1;
    fBinSizeZ[ 65] = kbinz1;

    fBinSizeZ[ 97] = kbinz1;
    fBinSizeZ[ 98] = kbinz1;

    fBinSizeZ[130] = kbinz1;
    fBinSizeZ[131] = kbinz1;

    fBinSizeZ[163] = kbinz1;
    fBinSizeZ[164] = kbinz1;

    // correct detector size for bin size.
    fdx = 0.0;
    for(i=0;i<fNbinx;i++) fdx +=fBinSizeX[i];
    fdx *= 0.5;
    fdz = 0.0;
    for(i=0;i<fNbinz;i++) fdz +=fBinSizeZ[i];
    fdz *= 0.5;

    fShapeSPD = new TBRIK("ActiveSPD","Active volume of SPD","SPD SI DET",
			  fdx,fdy,fdz);
}
//______________________________________________________________________
AliITSgeomSPD425::AliITSgeomSPD425(Float_t dy,Int_t nx,Float_t *bx,
				                Int_t nz,Float_t *bz){
////////////////////////////////////////////////////////////////////////
//    default constructor, for a User modified TDR based geometry.
////////////////////////////////////////////////////////////////////////
    Int_t i;
    fdx = 0.0;
    fdy =  dy;
    fdz = 0.0;
    fNbinx = nx; // new number of bins in x.
    fNbinz = nz; // new number of bins in z.

    fBinSizeX = new Float_t[fNbinx]; // array of bin sizes along x.
    for(i=0;i<fNbinx;i++) fBinSizeX[i] = bx[i]; // new x bin size.
    fBinSizeZ = new Float_t[fNbinz]; // array of bin sizes along z.
    for(i=0;i<fNbinz;i++) fBinSizeZ[i] = bz[i]; // new z bin size.

    // correct detector size for bin size.
    for(i=0;i<fNbinx;i++) fdx +=fBinSizeX[i];
    fdx *= 0.5;
    for(i=0;i<fNbinz;i++) fdz +=fBinSizeZ[i];
    fdz *= 0.5;

    fShapeSPD = new TBRIK("ActiveSPD","Active volume of SPD","SPD SI DET",
			  fdx,fdy,fdz);
}
//______________________________________________________________________
AliITSgeomSPD425::AliITSgeomSPD425(AliITSgeomSPD425 &source){
  // copy constructor
    Int_t i;
    if(&source == this) return;
    this->fShapeSPD = new TBRIK(*(source.fShapeSPD));
    this->fdx = source.fdx;
    this->fdy = source.fdy;
    this->fdz = source.fdz;
    if(this->fBinSizeX) delete[] this->fBinSizeX; 
    if(this->fBinSizeX) delete[] this->fBinSizeZ;
    this->fNbinx = source.fNbinx;
    this->fBinSizeX = new Float_t[this->fNbinx];
    this->fNbinz = source.fNbinz;
    this->fBinSizeZ = new Float_t[this->fNbinz];
    for(i=0;i<fNbinx;i++) this->fBinSizeX[i] = source.fBinSizeX[i];
    for(i=0;i<fNbinz;i++) this->fBinSizeZ[i] = source.fBinSizeZ[i];
}
//______________________________________________________________________
AliITSgeomSPD425& AliITSgeomSPD425::operator=(AliITSgeomSPD425 &source){
  // = operator
    Int_t i;
    if(&source == this) return *this;
    this->fShapeSPD = new TBRIK(*(source.fShapeSPD));
    this->fdx = source.fdx;
    this->fdy = source.fdy;
    this->fdz = source.fdz;
    if(this->fBinSizeX) delete[] this->fBinSizeX; 
    if(this->fBinSizeX) delete[] this->fBinSizeZ;
    this->fNbinx = source.fNbinx;
    this->fBinSizeX = new Float_t[this->fNbinx];
    this->fNbinz = source.fNbinz;
    this->fBinSizeZ = new Float_t[this->fNbinz];
    for(i=0;i<fNbinx;i++) this->fBinSizeX[i] = source.fBinSizeX[i];
    for(i=0;i<fNbinz;i++) this->fBinSizeZ[i] = source.fBinSizeZ[i];
    return *this;
}
//______________________________________________________________________
AliITSgeomSPD425::~AliITSgeomSPD425(){
  // destructor
    delete[] fBinSizeX;
    delete[] fBinSizeZ;
    delete   fShapeSPD;
}
//______________________________________________________________________
void AliITSgeomSPD425::ReSetBins(Float_t dy,Int_t nx,Float_t *bx,
                                        Int_t nz,Float_t *bz){
////////////////////////////////////////////////////////////////////////
//    default constructor, for a User modified TDR based geometry.
////////////////////////////////////////////////////////////////////////
    Int_t i;
    fdx = 0.0;
    fdy =  dy;
    fdz = 0.0;
    fNbinx = nx; // new number of bins in x.
    fNbinz = nz; // new number of bins in z.

    if(fBinSizeX!=0) delete[] fBinSizeX;
    fBinSizeX = new Float_t[fNbinx]; // array of bin sizes along x.
    for(i=0;i<fNbinx;i++) fBinSizeX[i] = bx[i]; // new x bin size.
    if(fBinSizeZ!=0) delete[] fBinSizeZ;
    fBinSizeZ = new Float_t[fNbinz]; // array of bin sizes along z.
    for(i=0;i<fNbinz;i++) fBinSizeZ[i] = bz[i]; // new z bin size.

    // correct detector size for bin size.
    for(i=0;i<fNbinx;i++) fdx +=fBinSizeX[i];
    fdx *= 0.5;
    for(i=0;i<fNbinz;i++) fdz +=fBinSizeZ[i];
    fdz *= 0.5;

    fShapeSPD = new TBRIK("ActiveSPD","Active volume of SPD","SPD SI DET",
			  fdx,fdy,fdz);
}
//----------------------------------------------------------------------
void AliITSgeomSPD425::Streamer(TBuffer &R__b){
    // Streamer function for the class AliITSgeomSPD425.
    UInt_t R__s, R__c;

    if(R__b.IsReading()){
	Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
	if(R__v==1){
	    TObject::Streamer(R__b);
	    fShapeSPD->Streamer(R__b);
	    R__b >> fdx;
	    R__b >> fdy;
	    R__b >> fdz;
	}else if (R__v==2){
	    AliITSgeomSPD::Streamer(R__b);
	    fShapeSPD->Streamer(R__b);
	    R__b >> fdx;
	    R__b >> fdy;
	    R__b >> fdz;
	    R__b >> fNbinx;
	    if(fBinSizeX!=0) delete[] fBinSizeX;
	    fBinSizeX = new Float_t[fNbinx];
	    for(Int_t i=0;i<fNbinx;i++) R__b >> fBinSizeX[i];
	    R__b >> fNbinz;
	    if(fBinSizeZ!=0) delete[] fBinSizeZ;
	    fBinSizeZ = new Float_t[fNbinz];
	    for(Int_t i=0;i<fNbinz;i++) R__b >> fBinSizeZ[i];
	    R__b.CheckByteCount(R__s, R__c, AliITSgeomSPD425::IsA());
	} // end if R__v==1
    } else { // IsWriting.
	R__c = R__b.WriteVersion(AliITSgeomSPD425::IsA(), kTRUE);
	AliITSgeomSPD::Streamer(R__b);
	fShapeSPD->Streamer(R__b);
	R__b << fdx;
	R__b << fdy;
	R__b << fdz;
	R__b << fNbinx;
	for(Int_t i=0;i<fNbinx;i++) R__b << fBinSizeX[i];
	R__b << fNbinz;
	for(Int_t i=0;i<fNbinz;i++) R__b << fBinSizeZ[i];
	R__b.SetByteCount(R__c, kTRUE);
    } // end if R__b.IsReading()
}
//----------------------------------------------------------------------
Commit	Line	Data
d53869b0	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/*
	17	$Log$
b3acdf28	18	Revision 1.3 2000/06/10 10:42:49 nilsen
	19	Fixed bug in copy and operator =.
	20
d53869b0	21
	22	*/
	23
	24	#include <TShape.h>
	25
	26	#include "AliITSgeomSPD425.h"
	27
	28	ClassImp(AliITSgeomSPD425)
	29
	30	AliITSgeomSPD425::AliITSgeomSPD425(){
	31	////////////////////////////////////////////////////////////////////////
	32	// default constructor, for ITS post TDR geometry.
	33	////////////////////////////////////////////////////////////////////////
	34
	35	const Float_t kdx=0.6400,kdy=0.0075,kdz=4.23625;// cm; Standard pixel detector
	36	// size is 2dx wide, 2dz long,
	37	// and 2dy thick. Geant 3.12
	38	// units.
	39	const Float_t kbinx0 = 0.0050; // cm; Standard pixel size in x direction.
	40	const Int_t knbinx = 256; // number of pixels along x direction.
	41	const Float_t kbinz0 = 0.0425; // cm; Standard pixel size in z direction.
	42	const Float_t kbinz1 = 0.0625; // cm; Special pixel size in z direction.
	43	const Int_t knbinz = 197; // number of pixels along z direction.
b3acdf28	44	Int_t i;
d53869b0	45
	46	fdx = kdx; // default value.
	47	fdy = kdy; // default value.
	48	fdz = kdz; // default value.
	49	fNbinx = knbinx; // default number of bins in x.
	50	fNbinz = knbinz; // default number of bins in z.
	51
	52	fBinSizeX = new Float_t[fNbinx]; // array of bin sizes along x.
b3acdf28	53	for(i=0;i<fNbinx;i++) fBinSizeX[i] = kbinx0; // default x bin size.
d53869b0	54	fBinSizeZ = new Float_t[fNbinz]; // array of bin sizes along z.
b3acdf28	55	for(i=0;i<fNbinz;i++) fBinSizeZ[i] = kbinz0; // default z bin size.
d53869b0	56	fBinSizeZ[ 31] = kbinz1;
	57	fBinSizeZ[ 32] = kbinz1;
	58
	59	fBinSizeZ[ 64] = kbinz1;
	60	fBinSizeZ[ 65] = kbinz1;
	61
	62	fBinSizeZ[ 97] = kbinz1;
	63	fBinSizeZ[ 98] = kbinz1;
	64
	65	fBinSizeZ[130] = kbinz1;
	66	fBinSizeZ[131] = kbinz1;
	67
	68	fBinSizeZ[163] = kbinz1;
	69	fBinSizeZ[164] = kbinz1;
	70
	71	// correct detector size for bin size.
	72	fdx = 0.0;
b3acdf28	73	for(i=0;i<fNbinx;i++) fdx +=fBinSizeX[i];
d53869b0	74	fdx *= 0.5;
d53869b0	75	fdz = 0.0;
b3acdf28	76	for(i=0;i<fNbinz;i++) fdz +=fBinSizeZ[i];
d53869b0	77	fdz *= 0.5;
	78
	79	fShapeSPD = new TBRIK("ActiveSPD","Active volume of SPD","SPD SI DET",
	80	fdx,fdy,fdz);
	81	}
	82	//______________________________________________________________________
	83	AliITSgeomSPD425::AliITSgeomSPD425(Float_t dy,Int_t nx,Float_t *bx,
	84	Int_t nz,Float_t *bz){
	85	////////////////////////////////////////////////////////////////////////
	86	// default constructor, for a User modified TDR based geometry.
	87	////////////////////////////////////////////////////////////////////////
b3acdf28	88	Int_t i;
d53869b0	89	fdx = 0.0;
	90	fdy = dy;
	91	fdz = 0.0;
	92	fNbinx = nx; // new number of bins in x.
	93	fNbinz = nz; // new number of bins in z.
	94
	95	fBinSizeX = new Float_t[fNbinx]; // array of bin sizes along x.
b3acdf28	96	for(i=0;i<fNbinx;i++) fBinSizeX[i] = bx[i]; // new x bin size.
d53869b0	97	fBinSizeZ = new Float_t[fNbinz]; // array of bin sizes along z.
b3acdf28	98	for(i=0;i<fNbinz;i++) fBinSizeZ[i] = bz[i]; // new z bin size.
d53869b0	99
d53869b0	100	// correct detector size for bin size.
b3acdf28	101	for(i=0;i<fNbinx;i++) fdx +=fBinSizeX[i];
d53869b0	102	fdx *= 0.5;
b3acdf28	103	for(i=0;i<fNbinz;i++) fdz +=fBinSizeZ[i];
d53869b0	104	fdz *= 0.5;
	105
	106	fShapeSPD = new TBRIK("ActiveSPD","Active volume of SPD","SPD SI DET",
	107	fdx,fdy,fdz);
	108	}
	109	//______________________________________________________________________
	110	AliITSgeomSPD425::AliITSgeomSPD425(AliITSgeomSPD425 &source){
	111	// copy constructor
b3acdf28	112	Int_t i;
d53869b0	113	if(&source == this) return;
93a31784	114	this->fShapeSPD = new TBRIK(*(source.fShapeSPD));
d53869b0	115	this->fdx = source.fdx;
	116	this->fdy = source.fdy;
	117	this->fdz = source.fdz;
	118	if(this->fBinSizeX) delete[] this->fBinSizeX;
	119	if(this->fBinSizeX) delete[] this->fBinSizeZ;
	120	this->fNbinx = source.fNbinx;
	121	this->fBinSizeX = new Float_t[this->fNbinx];
	122	this->fNbinz = source.fNbinz;
	123	this->fBinSizeZ = new Float_t[this->fNbinz];
b3acdf28	124	for(i=0;i<fNbinx;i++) this->fBinSizeX[i] = source.fBinSizeX[i];
b3acdf28	125	for(i=0;i<fNbinz;i++) this->fBinSizeZ[i] = source.fBinSizeZ[i];
d53869b0	126	}
	127	//______________________________________________________________________
	128	AliITSgeomSPD425& AliITSgeomSPD425::operator=(AliITSgeomSPD425 &source){
	129	// = operator
b3acdf28	130	Int_t i;
d53869b0	131	if(&source == this) return *this;
93a31784	132	this->fShapeSPD = new TBRIK(*(source.fShapeSPD));
d53869b0	133	this->fdx = source.fdx;
	134	this->fdy = source.fdy;
	135	this->fdz = source.fdz;
	136	if(this->fBinSizeX) delete[] this->fBinSizeX;
	137	if(this->fBinSizeX) delete[] this->fBinSizeZ;
	138	this->fNbinx = source.fNbinx;
	139	this->fBinSizeX = new Float_t[this->fNbinx];
	140	this->fNbinz = source.fNbinz;
	141	this->fBinSizeZ = new Float_t[this->fNbinz];
b3acdf28	142	for(i=0;i<fNbinx;i++) this->fBinSizeX[i] = source.fBinSizeX[i];
b3acdf28	143	for(i=0;i<fNbinz;i++) this->fBinSizeZ[i] = source.fBinSizeZ[i];
d53869b0	144	return *this;
	145	}
	146	//______________________________________________________________________
	147	AliITSgeomSPD425::~AliITSgeomSPD425(){
	148	// destructor
	149	delete[] fBinSizeX;
	150	delete[] fBinSizeZ;
	151	delete fShapeSPD;
	152	}
	153	//______________________________________________________________________
	154	void AliITSgeomSPD425::ReSetBins(Float_t dy,Int_t nx,Float_t *bx,
	155	Int_t nz,Float_t *bz){
	156	////////////////////////////////////////////////////////////////////////
	157	// default constructor, for a User modified TDR based geometry.
	158	////////////////////////////////////////////////////////////////////////
b3acdf28	159	Int_t i;
d53869b0	160	fdx = 0.0;
	161	fdy = dy;
	162	fdz = 0.0;
	163	fNbinx = nx; // new number of bins in x.
	164	fNbinz = nz; // new number of bins in z.
	165
	166	if(fBinSizeX!=0) delete[] fBinSizeX;
	167	fBinSizeX = new Float_t[fNbinx]; // array of bin sizes along x.
b3acdf28	168	for(i=0;i<fNbinx;i++) fBinSizeX[i] = bx[i]; // new x bin size.
d53869b0	169	if(fBinSizeZ!=0) delete[] fBinSizeZ;
d53869b0	170	fBinSizeZ = new Float_t[fNbinz]; // array of bin sizes along z.
b3acdf28	171	for(i=0;i<fNbinz;i++) fBinSizeZ[i] = bz[i]; // new z bin size.
d53869b0	172
d53869b0	173	// correct detector size for bin size.
b3acdf28	174	for(i=0;i<fNbinx;i++) fdx +=fBinSizeX[i];
d53869b0	175	fdx *= 0.5;
b3acdf28	176	for(i=0;i<fNbinz;i++) fdz +=fBinSizeZ[i];
d53869b0	177	fdz *= 0.5;
	178
	179	fShapeSPD = new TBRIK("ActiveSPD","Active volume of SPD","SPD SI DET",
	180	fdx,fdy,fdz);
	181	}
	182	//----------------------------------------------------------------------
	183	void AliITSgeomSPD425::Streamer(TBuffer &R__b){
	184	// Streamer function for the class AliITSgeomSPD425.
	185	UInt_t R__s, R__c;
	186
	187	if(R__b.IsReading()){
	188	Version_t R__v = R__b.ReadVersion(&R__s, &R__c);
	189	if(R__v==1){
	190	TObject::Streamer(R__b);
	191	fShapeSPD->Streamer(R__b);
	192	R__b >> fdx;
	193	R__b >> fdy;
	194	R__b >> fdz;
	195	}else if (R__v==2){
	196	AliITSgeomSPD::Streamer(R__b);
	197	fShapeSPD->Streamer(R__b);
	198	R__b >> fdx;
	199	R__b >> fdy;
	200	R__b >> fdz;
	201	R__b >> fNbinx;
	202	if(fBinSizeX!=0) delete[] fBinSizeX;
	203	fBinSizeX = new Float_t[fNbinx];
	204	for(Int_t i=0;i<fNbinx;i++) R__b >> fBinSizeX[i];
	205	R__b >> fNbinz;
	206	if(fBinSizeZ!=0) delete[] fBinSizeZ;
	207	fBinSizeZ = new Float_t[fNbinz];
	208	for(Int_t i=0;i<fNbinz;i++) R__b >> fBinSizeZ[i];
	209	R__b.CheckByteCount(R__s, R__c, AliITSgeomSPD425::IsA());
	210	} // end if R__v==1
	211	} else { // IsWriting.
	212	R__c = R__b.WriteVersion(AliITSgeomSPD425::IsA(), kTRUE);
	213	AliITSgeomSPD::Streamer(R__b);
	214	fShapeSPD->Streamer(R__b);
	215	R__b << fdx;
	216	R__b << fdy;
	217	R__b << fdz;
	218	R__b << fNbinx;
	219	for(Int_t i=0;i<fNbinx;i++) R__b << fBinSizeX[i];
	220	R__b << fNbinz;
	221	for(Int_t i=0;i<fNbinz;i++) R__b << fBinSizeZ[i];
	222	R__b.SetByteCount(R__c, kTRUE);
	223	} // end if R__b.IsReading()
	224	}
	225	//----------------------------------------------------------------------