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

/**************************************************************************
 * Copyright(c) 2007-2009, 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.                  *
 **************************************************************************/

/* $Id$ */

///////////////////////////////////////////////////////////////////
//                                                               //
// Implementation of the class with array of AliITSDriftSpeedSDD //
// Origin: F.Prino, Torino, prino@to.infn.it                     //
//                                                               //
///////////////////////////////////////////////////////////////////

#include <TClonesArray.h>
#include "AliITSDriftSpeedArraySDD.h"
#include "AliITSDriftSpeedSDD.h"
#include "AliLog.h"

ClassImp(AliITSDriftSpeedArraySDD)
//______________________________________________________________________
AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD():
TObject(),
fNEvents(0),
fDriftSpeedSDD(0){
  // default constructor
  fDriftSpeedSDD=new TClonesArray("AliITSDriftSpeedSDD",100);
}
//______________________________________________________________________
AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD(Int_t numEv):
TObject(),
fNEvents(0),
fDriftSpeedSDD(0){
  // standard constructor
  fDriftSpeedSDD=new TClonesArray("AliITSDriftSpeedSDD",numEv);
}
//______________________________________________________________________
AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD(const AliITSDriftSpeedArraySDD& array):
TObject(),
fNEvents(array.fNEvents),
fDriftSpeedSDD(array.fDriftSpeedSDD){
  // copy constructor
}
//______________________________________________________________________
AliITSDriftSpeedArraySDD& AliITSDriftSpeedArraySDD::operator=(const AliITSDriftSpeedArraySDD& array){
  // assignment operator
  this->~AliITSDriftSpeedArraySDD();
  new(this) AliITSDriftSpeedArraySDD(array);
  return *this;
}

//______________________________________________________________________
AliITSDriftSpeedArraySDD::~AliITSDriftSpeedArraySDD(){
  // destructor
  if(fDriftSpeedSDD){
    fDriftSpeedSDD->Delete();
    delete fDriftSpeedSDD;
  }
}
//______________________________________________________________________
void AliITSDriftSpeedArraySDD::AddDriftSpeed(AliITSDriftSpeedSDD* drSpeed){
  // adds an AliITSDriftSpeedSDD object in the array
  TClonesArray &arr = *fDriftSpeedSDD;
  new(arr[fNEvents]) AliITSDriftSpeedSDD(*drSpeed);
  fNEvents++;
  
}
//______________________________________________________________________
void AliITSDriftSpeedArraySDD::PrintAll() const{
  // print drift speed parameters for all elements in the array
  printf("Array Size=%d\n",fDriftSpeedSDD->GetSize());
  printf("Array Elements =%d\n",fNEvents);
  for(Int_t i=0;i<fNEvents; i++){
    printf("     ====== Array el. #%d ======\n",i);
    AliITSDriftSpeedSDD *d=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(i);
    if(d) d->PrintDriftSpeedParameters();
  }
}
//______________________________________________________________________
Float_t AliITSDriftSpeedArraySDD::GetDriftSpeed(Int_t iEvent, Float_t iAnode) const{
  // returns drift speed for given event number and anode
  if(!fDriftSpeedSDD->IsSorted()) fDriftSpeedSDD->Sort();
  if(fNEvents==1){
    AliITSDriftSpeedSDD *d=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(0);
    return d->GetDriftSpeedAtAnode(iAnode);
  }else{
    Int_t nInjEv1=-1;
    Int_t nInjEv2=-1;
    AliITSDriftSpeedSDD *d1=NULL;
    AliITSDriftSpeedSDD *d2=NULL;
    for(Int_t i=0;i<fNEvents; i++){
      d1=d2;
      d2=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(i);
      nInjEv1=nInjEv2;
      if(d2!=0){
	nInjEv2=d2->GetEventNumber();
	if(nInjEv2>=iEvent){
	  if(i==0) d1=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(i+1);
	  nInjEv1=d1->GetEventNumber();
	  break;
	}
      }
    }
    if(nInjEv1>=0 && nInjEv2>=0 && nInjEv1!=nInjEv2){
      Float_t v1=d1->GetDriftSpeedAtAnode(iAnode);
      Float_t v2=d2->GetDriftSpeedAtAnode(iAnode);
      Float_t vdrift=(v2-v1)*(iEvent-nInjEv1)/(nInjEv2-nInjEv1)+v1;
      return vdrift;
    }
  }
  AliWarning("Vdrift interpolation error\n");
  return -999.;
}
Commit	Line	Data
18da6e54	1	/**************************************************************************
	2	* Copyright(c) 2007-2009, 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	/* $Id$ */
	17
	18	///////////////////////////////////////////////////////////////////
	19	// //
	20	// Implementation of the class with array of AliITSDriftSpeedSDD //
	21	// Origin: F.Prino, Torino, prino@to.infn.it //
	22	// //
	23	///////////////////////////////////////////////////////////////////
	24
	25	#include <TClonesArray.h>
	26	#include "AliITSDriftSpeedArraySDD.h"
	27	#include "AliITSDriftSpeedSDD.h"
	28	#include "AliLog.h"
	29
	30	ClassImp(AliITSDriftSpeedArraySDD)
	31	//______________________________________________________________________
	32	AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD():
	33	TObject(),
	34	fNEvents(0),
	35	fDriftSpeedSDD(0){
	36	// default constructor
	37	fDriftSpeedSDD=new TClonesArray("AliITSDriftSpeedSDD",100);
	38	}
	39	//______________________________________________________________________
	40	AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD(Int_t numEv):
	41	TObject(),
	42	fNEvents(0),
	43	fDriftSpeedSDD(0){
	44	// standard constructor
	45	fDriftSpeedSDD=new TClonesArray("AliITSDriftSpeedSDD",numEv);
	46	}
	47	//______________________________________________________________________
	48	AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD(const AliITSDriftSpeedArraySDD& array):
	49	TObject(),
	50	fNEvents(array.fNEvents),
	51	fDriftSpeedSDD(array.fDriftSpeedSDD){
	52	// copy constructor
	53	}
	54	//______________________________________________________________________
	55	AliITSDriftSpeedArraySDD& AliITSDriftSpeedArraySDD::operator=(const AliITSDriftSpeedArraySDD& array){
	56	// assignment operator
	57	this->~AliITSDriftSpeedArraySDD();
	58	new(this) AliITSDriftSpeedArraySDD(array);
	59	return *this;
	60	}
	61
	62	//______________________________________________________________________
	63	AliITSDriftSpeedArraySDD::~AliITSDriftSpeedArraySDD(){
	64	// destructor
65	if(fDriftSpeedSDD){
66	fDriftSpeedSDD->Delete();
67	delete fDriftSpeedSDD;
68	}
69	}
70	//______________________________________________________________________
71	void AliITSDriftSpeedArraySDD::AddDriftSpeed(AliITSDriftSpeedSDD* drSpeed){
72	// adds an AliITSDriftSpeedSDD object in the array
73	TClonesArray &arr = *fDriftSpeedSDD;
74	new(arr[fNEvents]) AliITSDriftSpeedSDD(*drSpeed);
75	fNEvents++;
76
77	}
78	//______________________________________________________________________
79	void AliITSDriftSpeedArraySDD::PrintAll() const{
80	// print drift speed parameters for all elements in the array
81	printf("Array Size=%d\n",fDriftSpeedSDD->GetSize());
82	printf("Array Elements =%d\n",fNEvents);
83	for(Int_t i=0;i<fNEvents; i++){
84	printf(" ====== Array el. #%d ======\n",i);
85	AliITSDriftSpeedSDD d=(AliITSDriftSpeedSDD)fDriftSpeedSDD->At(i);
86	if(d) d->PrintDriftSpeedParameters();
87	}
88	}
89	//______________________________________________________________________
90	Float_t AliITSDriftSpeedArraySDD::GetDriftSpeed(Int_t iEvent, Float_t iAnode) const{
91	// returns drift speed for given event number and anode
92	if(!fDriftSpeedSDD->IsSorted()) fDriftSpeedSDD->Sort();
93	if(fNEvents==1){
94	AliITSDriftSpeedSDD d=(AliITSDriftSpeedSDD)fDriftSpeedSDD->At(0);
95	return d->GetDriftSpeedAtAnode(iAnode);
96	}else{
97	Int_t nInjEv1=-1;
98	Int_t nInjEv2=-1;
99	AliITSDriftSpeedSDD *d1=NULL;
100	AliITSDriftSpeedSDD *d2=NULL;
101	for(Int_t i=0;i<fNEvents; i++){
102	d1=d2;
103	d2=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(i);
104	nInjEv1=nInjEv2;
105	if(d2!=0){
106	nInjEv2=d2->GetEventNumber();
107	if(nInjEv2>=iEvent){
108	if(i==0) d1=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(i+1);
109	nInjEv1=d1->GetEventNumber();
110	break;
111	}
112	}
113	}
114	if(nInjEv1>=0 && nInjEv2>=0 && nInjEv1!=nInjEv2){
115	Float_t v1=d1->GetDriftSpeedAtAnode(iAnode);
116	Float_t v2=d2->GetDriftSpeedAtAnode(iAnode);
117	Float_t vdrift=(v2-v1)*(iEvent-nInjEv1)/(nInjEv2-nInjEv1)+v1;
118	return vdrift;
119	}
120	}
121	AliWarning("Vdrift interpolation error\n");
122	return -999.;
123	}