1 /**************************************************************************
2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////
20 // Implementation of the class with array of AliITSDriftSpeedSDD //
21 // Origin: F.Prino, Torino, prino@to.infn.it //
23 ///////////////////////////////////////////////////////////////////
25 #include <TClonesArray.h>
26 #include "AliITSDriftSpeedArraySDD.h"
27 #include "AliITSDriftSpeedSDD.h"
30 ClassImp(AliITSDriftSpeedArraySDD)
31 //______________________________________________________________________
32 AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD():
36 // default constructor
37 fDriftSpeedSDD=new TClonesArray("AliITSDriftSpeedSDD",100);
39 //______________________________________________________________________
40 AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD(Int_t numEv):
44 // standard constructor
45 fDriftSpeedSDD=new TClonesArray("AliITSDriftSpeedSDD",numEv);
47 //______________________________________________________________________
48 AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD(const AliITSDriftSpeedArraySDD& array):
50 fNEvents(array.fNEvents),
51 fDriftSpeedSDD(array.fDriftSpeedSDD){
54 //______________________________________________________________________
55 AliITSDriftSpeedArraySDD& AliITSDriftSpeedArraySDD::operator=(const AliITSDriftSpeedArraySDD& array){
56 // assignment operator
57 this->~AliITSDriftSpeedArraySDD();
58 new(this) AliITSDriftSpeedArraySDD(array);
62 //______________________________________________________________________
63 AliITSDriftSpeedArraySDD::~AliITSDriftSpeedArraySDD(){
66 fDriftSpeedSDD->Delete();
67 delete fDriftSpeedSDD;
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);
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();
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();
94 AliITSDriftSpeedSDD *d=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(0);
95 return d->GetDriftSpeedAtAnode(iAnode);
99 AliITSDriftSpeedSDD *d1=NULL;
100 AliITSDriftSpeedSDD *d2=NULL;
101 for(Int_t i=0;i<fNEvents; i++){
103 d2=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(i);
106 nInjEv2=d2->GetEventNumber();
108 if(i==0) d1=(AliITSDriftSpeedSDD*)fDriftSpeedSDD->At(i+1);
109 nInjEv1=d1->GetEventNumber();
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;
121 AliWarning("Vdrift interpolation error\n");