Add Reaction Plane

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

diff --git a/ITS/AliITSPlaneEffSPD.cxx b/ITS/AliITSPlaneEffSPD.cxx
index eeb6637ac85c9bcb3a2f1201cc6b99d8262d3197..6de1135ccf966acdf49ebefd4834f837f923f50b 100644 (file)
--- a/ITS/AliITSPlaneEffSPD.cxx
+++ b/ITS/AliITSPlaneEffSPD.cxx
@@ -341,11 +341,17 @@ Int_t nf=-1;
 Int_t nt=-1;
 if(fo) {
  AliWarning("PlaneEff: you asked for FO efficiency");
- nf=fFound[GetKey(im,ic,fo,bcm4)];
- nt=fTried[GetKey(im,ic,fo,bcm4)];
+ UInt_t key=GetKey(im,ic,fo,bcm4);
+ if(key<kNModule*kNChip*(kNClockPhase+1)) {
+   nf=fFound[key];
+   nt=fTried[key];
+ }
 } else {
- nf=fFound[GetKey(im,ic)];
- nt=fTried[GetKey(im,ic)];
+ UInt_t key=GetKey(im,ic);
+ if (key<kNModule*kNChip) {
+  nf=fFound[key];
+  nt=fTried[key];
+ }
 }
 return AliITSPlaneEff::PlaneEff(nf,nt);
 }
@@ -366,11 +372,17 @@ Int_t nf=-1;
 Int_t nt=-1;
 if(fo) {
  AliWarning("ErrPlaneEff: you asked for FO efficiency");
- nf=fFound[GetKey(im,ic,fo,bcm4)];
- nt=fTried[GetKey(im,ic,fo,bcm4)];
+ UInt_t key=GetKey(im,ic,fo,bcm4);
+ if(key<kNModule*kNChip*(kNClockPhase+1)) {
+   nf=fFound[key];
+   nt=fTried[key];
+ }
 } else {
- nf=fFound[GetKey(im,ic)];
- nt=fTried[GetKey(im,ic)];
+ UInt_t key=GetKey(im,ic);
+ if (key<kNModule*kNChip) {
+   nf=fFound[key];
+   nt=fTried[key];
+ }
 }
 return AliITSPlaneEff::ErrPlaneEff(nf,nt);
 } 
@@ -383,14 +395,20 @@ if (im>=kNModule || ic>=kNChip)
 if(fo && bcm4>=kNClockPhase)
  {AliError("UpDatePlaneEff: you asked for Fast Or in a wrong phase"); return kFALSE;}
 if (!fo) {
- fTried[GetKey(im,ic)]++;
- if(Kfound) fFound[GetKey(im,ic)]++;
- return kTRUE;
+ UInt_t key=GetKey(im,ic);
+ if(key<kNModule*kNChip) {
+   fTried[key]++;
+   if(Kfound) fFound[key]++;
+   return kTRUE;
+ }
 }
 else {
- fTried[GetKey(im,ic,fo,bcm4)]++;
- if(Kfound) fFound[GetKey(im,ic,fo,bcm4)]++;
- return kTRUE;
+ UInt_t key=GetKey(im,ic,fo,bcm4);
+ if(key<kNModule*kNChip*(kNClockPhase+1)) {
+   fTried[key]++;
+   if(Kfound) fFound[key]++;
+   return kTRUE;
+ }
 }
 return kFALSE;
 }