TRD nSigma OADB related new codes and modifications and OADB root file -- Xianguo Lu
[u/mrichter/AliRoot.git] / ITS / AliITSPlaneEff.cxx
CommitLineData
4a66240a 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// Base Plane Efficiency class for ITS
17// Specific subdetector implementation is done in
18// AliITSPlaneEffSPD
19// AliITSPlaneEffSDD
20// AliITSPlaneEffSSD
21//
22// Author: G.E. Bruno
23// giuseppe.bruno@ba.infn.it
24//
25///////////////////////////////////////////////////////////////////////////
26
27/* $Id$ */
28
29#include <TMath.h>
30#include "AliITSPlaneEff.h"
31#include "AliLog.h"
32#include "AliCDBManager.h"
33#include "AliCDBStorage.h"
34
35ClassImp(AliITSPlaneEff)
36//______________________________________________________________________
7167ae53 37AliITSPlaneEff::AliITSPlaneEff(): AliPlaneEff(),
4a66240a 38fRunNumber(0),
39fCDBUri(""),
5fbd4fd6 40fInitCDBCalled(kFALSE),
41fHis(kFALSE)
4a66240a 42{
43 // Default constructor
44 // Inputs:
45 // none.
46 // Outputs:
47 // none.
48 // Return:
49 // a default constructed AliITSPlaneEff class
b0b80d33 50 //InitCDB();
4a66240a 51}
52//______________________________________________________________________
7167ae53 53AliITSPlaneEff::AliITSPlaneEff(const AliITSPlaneEff &s) : AliPlaneEff(s),
4a66240a 54fRunNumber(s.fRunNumber),
55fCDBUri(s.fCDBUri),
5fbd4fd6 56fInitCDBCalled(s.fInitCDBCalled),
57fHis(s.fHis)
4a66240a 58{
59 // Copy Constructor
60 // Inputs:
61 // const AliITSPlaneEff &s simulation class to copy from
62 // Outputs:
63 // none.
64 // Return:
65 // a standard constructed AliITSPlaneEff class with values the same
66 // as that of s.
67
68}
69//_________________________________________________________________________
70//void AliITSPlaneEff::operator+=(const AliITSPlaneEff &add){
71 // Add to me operator
72 // Inputs:
73 // const AliITSPlaneEff &add simulation class to be added
74 // Outputs:
75 // none.
76 // Return:
77 // none
78
79// return;
80//}
81//_________________________________________________________________________
82AliITSPlaneEff& AliITSPlaneEff::operator=(const AliITSPlaneEff &source){
83 // Assignment operator
84 // Inputs:
85 // const AliITSPlaneEff &source simulation class to copy from
86 // Outputs:
87 // none.
88 // Return:
89 // a standard constructed AliITSPlaneEff class with values the same
90 // as that of s.
91 if(this != &source){
92 source.Copy(*this);
93 }
94 return *this;
4a66240a 95}
7167ae53 96//_________________________________________________________________________
97/*
98AliPlaneEff& AliITSPlaneEff::operator=(const
99 AliPlaneEff &s){
100 // Assignment operator
101 // Inputs:
102 // AliPlaneEff &s The original class for which
103 // this class is a copy of
104 // Outputs:
105 // none.
106 // Return:
107
108 if(&s == this) return *this;
109 AliWarning("AliITSPlaneEff Not allowed to make a = Using default creator instead");
110 return *this;
111}
112*/
113//_________________________________________________________________________
4a66240a 114void AliITSPlaneEff::Copy(TObject &obj) const {
115 // copy this to obj
116 ((AliITSPlaneEff& ) obj).fRunNumber = fRunNumber;
117 ((AliITSPlaneEff& ) obj).fCDBUri = fCDBUri;
118 ((AliITSPlaneEff& ) obj).fInitCDBCalled = fInitCDBCalled;
5fbd4fd6 119 ((AliITSPlaneEff& ) obj).fHis = fHis;
4a66240a 120}
121//_________________________________________________________________________
122Double_t AliITSPlaneEff::PlaneEff(Int_t nf,Int_t nt) const {
123 // Compute the efficiency for a basic block,
124 // Inputs:
125 // number of associated cluslters (nf)
126 // number of used tracks (nt)
127 // Outputs:
128 // none.
129 // Return:
130 // the efficiency
131if(nf<0 || nt<=0 || nt<nf) {
132 AliInfo(Form("AliITSPlaneEff::PlaneEff: nfound= %i, ntried= %i",nf,nt));
133 return -1.;}
134 Double_t eff=nf;
135 return eff/=nt;
136}
137//_________________________________________________________________________
138Double_t AliITSPlaneEff::ErrPlaneEff(Int_t nf,Int_t nt) const{
139 // Compute the statistical error on efficiency for a basic block,
140 // using binomial statistics
141 // Inputs:
142 // number of associated cluslters (nf)
143 // number of used tracks (nt)
144 // Outputs:
145 // none.
146 // Return:
147 // the error on the efficiency
148if(nf<0 || nt<=0 || nt<nf) {
149 AliInfo(Form("AliITSPlaneEff::ErrPlaneEff: nfound= %i, ntried= %i",nf,nt));
150 return -1.;}
1cc5cedc 151 Double_t err=TMath::Sqrt((Double_t)nf*(1.-(Double_t)nf/(Double_t)nt));
152 return err/=(Double_t)nt;
4a66240a 153}
154//______________________________________________________________________
155Int_t AliITSPlaneEff::GetNTracksForGivenEff(Double_t eff, Double_t RelErr) const {
156 // Estimate of the number of tracks needed for measuring efficiency
157 // with a given relative error using binomial statistics
158 // Inputs:
159 // exspected efficiency eff (e.g. from previous measurements)
160 // wished relative error RelErr
161 // Outputs:
162 // none.
163 // Return: number of tracks given as the nearest integer
164if(eff<=0 || eff>1 || RelErr<=0 ) {
165 AliInfo(Form("AliITSPlaneEff::GetNTracksForGivenEff: eff= %f, RelErr= %f",eff,RelErr));
166 return -1;}
167return TMath::Nint((1-eff)/(eff*RelErr*RelErr));
168}
169//________________________________________________________________________
170void AliITSPlaneEff::InitCDB()
171{
172// activate a default CDB storage
173// First check if we have any CDB storage set, because it is used
174// to retrieve the calibration and alignment constants
175
b0b80d33 176 //if (fInitCDBCalled) return;
4a66240a 177 fInitCDBCalled = kTRUE;
178
179 AliCDBManager* man = AliCDBManager::Instance();
180 if (man->IsDefaultStorageSet())
181 {
182 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
183 AliWarning("Default CDB storage has been already set !");
184 AliWarning(Form("Ignoring the default storage declared in AliITSPlaneEff: %s",fCDBUri.Data()));
185 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
186 fCDBUri = man->GetDefaultStorage()->GetURI();
187 }
188 else {
189 if (fCDBUri.Length() > 0)
190 {
191 AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
192 AliDebug(2, Form("Default CDB storage is set to: %s", fCDBUri.Data()));
193 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
194 } else {
162637e4 195 fCDBUri="local://$ALICE_ROOT/OCDB";
4a66240a 196 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
197 AliWarning("Default CDB storage not yet set !!!!");
198 AliWarning(Form("Setting it now to: %s", fCDBUri.Data()));
199 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
200
201 }
202 man->SetDefaultStorage(fCDBUri);
b0b80d33 203 man->SetRun(fRunNumber);
4a66240a 204 }
b0b80d33 205 return;
4a66240a 206}
207//_____________________________________________________________________________
208void AliITSPlaneEff::SetDefaultStorage(const char* uri) {
209// Store the desired default CDB storage location
210// Activate it later within the Run() method
211
212 fCDBUri = uri;
213
214}
215