Changes related to the extraction of the V0 finder into a separate class (A. Dainese...
[u/mrichter/AliRoot.git] / ITS / AliITSDriftSpeedArraySDD.cxx
CommitLineData
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
18da6e54 25#include "AliITSDriftSpeedArraySDD.h"
26#include "AliITSDriftSpeedSDD.h"
27#include "AliLog.h"
28
29ClassImp(AliITSDriftSpeedArraySDD)
30//______________________________________________________________________
31AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD():
32TObject(),
33fNEvents(0),
7b208ef4 34fDriftSpeedSDD(10){
18da6e54 35 // default constructor
18da6e54 36}
37//______________________________________________________________________
38AliITSDriftSpeedArraySDD::AliITSDriftSpeedArraySDD(Int_t numEv):
39TObject(),
40fNEvents(0),
7b208ef4 41fDriftSpeedSDD(numEv){
18da6e54 42 // standard constructor
18da6e54 43}
44//______________________________________________________________________
45void AliITSDriftSpeedArraySDD::AddDriftSpeed(AliITSDriftSpeedSDD* drSpeed){
46 // adds an AliITSDriftSpeedSDD object in the array
4c1bcb5e 47 fDriftSpeedSDD.AddLast(drSpeed);
18da6e54 48 fNEvents++;
18da6e54 49}
50//______________________________________________________________________
51void AliITSDriftSpeedArraySDD::PrintAll() const{
52 // print drift speed parameters for all elements in the array
4c1bcb5e 53 printf("Array Size=%d\n",fDriftSpeedSDD.GetSize());
18da6e54 54 printf("Array Elements =%d\n",fNEvents);
55 for(Int_t i=0;i<fNEvents; i++){
56 printf(" ====== Array el. #%d ======\n",i);
4c1bcb5e 57 AliITSDriftSpeedSDD *d=(AliITSDriftSpeedSDD*)fDriftSpeedSDD.At(i);
18da6e54 58 if(d) d->PrintDriftSpeedParameters();
59 }
60}
61//______________________________________________________________________
4c1bcb5e 62Double_t AliITSDriftSpeedArraySDD::GetDriftSpeed(Int_t iEvent, Double_t iAnode){
18da6e54 63 // returns drift speed for given event number and anode
4c1bcb5e 64 if(!fDriftSpeedSDD.IsSorted()) fDriftSpeedSDD.Sort();
18da6e54 65 if(fNEvents==1){
4c1bcb5e 66 AliITSDriftSpeedSDD *d=(AliITSDriftSpeedSDD*)fDriftSpeedSDD.At(0);
18da6e54 67 return d->GetDriftSpeedAtAnode(iAnode);
68 }else{
69 Int_t nInjEv1=-1;
70 Int_t nInjEv2=-1;
71 AliITSDriftSpeedSDD *d1=NULL;
72 AliITSDriftSpeedSDD *d2=NULL;
73 for(Int_t i=0;i<fNEvents; i++){
74 d1=d2;
4c1bcb5e 75 d2=(AliITSDriftSpeedSDD*)fDriftSpeedSDD.At(i);
18da6e54 76 nInjEv1=nInjEv2;
77 if(d2!=0){
78 nInjEv2=d2->GetEventNumber();
79 if(nInjEv2>=iEvent){
4c1bcb5e 80 if(i==0) d1=(AliITSDriftSpeedSDD*)fDriftSpeedSDD.At(i+1);
18da6e54 81 nInjEv1=d1->GetEventNumber();
82 break;
83 }
84 }
85 }
86 if(nInjEv1>=0 && nInjEv2>=0 && nInjEv1!=nInjEv2){
4c1bcb5e 87 Double_t v1=d1->GetDriftSpeedAtAnode(iAnode);
88 Double_t v2=d2->GetDriftSpeedAtAnode(iAnode);
89 Double_t vdrift=(v2-v1)*(iEvent-nInjEv1)/(nInjEv2-nInjEv1)+v1;
18da6e54 90 return vdrift;
91 }
92 }
93 AliWarning("Vdrift interpolation error\n");
94 return -999.;
95}