Adding anti-Lambda and anti-Xi

[u/mrichter/AliRoot.git] / TPC / AliSimDigits.cxx

diff --git a/TPC/AliSimDigits.cxx b/TPC/AliSimDigits.cxx
index c57d6662c438098b391925d97da5167a7a17cb25..0ebe639593771336a2f62adb41801a08079d8c22 100644 (file)
--- a/TPC/AliSimDigits.cxx
+++ b/TPC/AliSimDigits.cxx
@@ -30,8 +30,7 @@
 #include "TError.h"
 #include "AliSegmentID.h"
 #include "AliH2F.h"
-#include "AliArrayI.h"
-#include "AliArrayS.h"
+#include "TArrayI.h"
 #include "AliDigits.h"
 #include "AliSimDigits.h"
 #include "AliTPCdigit.h"
@@ -45,12 +44,30 @@
 ClassImp(AliSimDigits)
 
 AliSimDigits::AliSimDigits()
+             :AliDigits(),
+             fTracks(0),
+             fTrIndex(0),
+             fNlevel(0),
+             fTrBufType(0)  
 {
   //  AliDigits::Invalite();
-  fTracks = 0;
-  fTrIndex = 0;  
+ 
   InvalidateTrack();
 }
+//
+AliSimDigits::AliSimDigits(const AliSimDigits &param)
+             :AliDigits(),
+             fTracks(0),
+             fTrIndex(0),
+             fNlevel(0),
+             fTrBufType(0) 
+{
+  //
+  // dummy
+  //
+  fTrIndex = param.fTrIndex;
+}
+//
 AliSimDigits::~AliSimDigits()
 {
 
@@ -62,15 +79,26 @@ AliSimDigits::~AliSimDigits()
   } 
 
 }
+AliSimDigits & AliSimDigits::operator =(const AliSimDigits & param)
+{
+  //
+  // assignment operator - dummy
+  //
+  if(this!=&param){
+    fTrIndex=param.fTrIndex;
+  }
+  return (*this);
+}
 
+//__________________________________________________________________
 void AliSimDigits::InvalidateTrack() 
 { 
   //
   //set default (invalid parameters)
   if ( fTracks != 0) delete fTracks;
-  fTracks = new AliArrayI;
+  fTracks = new TArrayI;
   if ( fTrIndex  != 0) delete fTrIndex;
-  fTrIndex = new AliArrayI;
+  fTrIndex = new TArrayI;
 
   for (Int_t i = 0; i<3; i++){
     fTracks->Set(0);
@@ -187,7 +215,7 @@ void  AliSimDigits::ExpandTrackBuffer1()
   Int_t all   = fNrows*fNcols;  //total number of digits
   Int_t elems = all*fNlevel;  //length of the buffer
 
-  AliArrayI * buf = new AliArrayI;
+  TArrayI * buf = new TArrayI;
   buf->Set(elems);
   fTrIndex->Set(0);
   //
@@ -228,9 +256,9 @@ void  AliSimDigits::CompresTrackBuffer1()
   //
   fTrBufType = 1;  
 
-  AliArrayI *  buf = new AliArrayI;   //create  new buffer 
+  TArrayI *  buf = new TArrayI;   //create  new buffer 
   buf->Set(fNrows*fNcols*fNlevel); //lets have the nearly the "worst case"
-  AliArrayI *  index = new AliArrayI;
+  TArrayI *  index = new TArrayI;
   index->Set(fNcols*fNlevel);
   //  Int_t * pindex = 
 
@@ -255,7 +283,7 @@ void  AliSimDigits::CompresTrackBuffer1()
        if (id <= 0) {
          if ( inum> 0 ) { //if we have some tracks in buffer
            icurrent++;
-           if ((icurrent+1)>=buf->fN) buf->Expand(icurrent*2+1); //MI change - allocate +1
+           if ((icurrent+1)>=buf->fN) buf->Set(icurrent*2+1); //MI change - allocate +1
            (*buf)[icurrent] = inum;
            icurrent++;
            (*buf)[icurrent] = lastID;  
@@ -269,7 +297,7 @@ void  AliSimDigits::CompresTrackBuffer1()
            if ( izero > 0 ) { 
              //if we have currently izero count of non tracks digits
              icurrent++;         
-             if (icurrent>=buf->fN) buf->Expand(icurrent*2+1);
+             if (icurrent>=buf->fN) buf->Set(icurrent*2+1);
              (*buf)[icurrent]= -izero;  //write how many under zero
              inum++;
              izero = 0;             
@@ -278,7 +306,7 @@ void  AliSimDigits::CompresTrackBuffer1()
            else{ 
              //if we change track id from another track id         
              icurrent++;         
-             if ((icurrent+1)>=buf->fN) buf->Expand(icurrent*2+1);
+             if ((icurrent+1)>=buf->fN) buf->Set(icurrent*2+1);
              (*buf)[icurrent] = inum;
              icurrent++;
              (*buf)[icurrent] = lastID;        
@@ -294,19 +322,19 @@ void  AliSimDigits::CompresTrackBuffer1()
       if ( izero > 0 ) { 
        //if we have currently izero count of non tracks digits
        icurrent++;       
-       if (icurrent>=buf->fN) buf->Expand(icurrent*2);
+       if (icurrent>=buf->fN) buf->Set(icurrent*2);
        (*buf)[icurrent]= -izero;  //write how many under zero  
       }
       if ( inum> 0 ) { //if we have some tracks in buffer
        icurrent++;
-       if ((icurrent+1)>=buf->fN) buf->Expand(icurrent*2);
+       if ((icurrent+1)>=buf->fN) buf->Set(icurrent*2);
        (*buf)[icurrent] = inum;
        icurrent++;
        (*buf)[icurrent] = id;  
       }      
     }//end of loop over columns
   }//end of loop over differnet track level  
-  buf->Expand(icurrent+1);
+  buf->Set(icurrent+1);
   delete fTracks;
   fTracks = buf;
   delete fTrIndex;
@@ -346,7 +374,8 @@ AliH2F *  AliSimDigits::DrawTracks( const char *option,Int_t level,
   //  
   //make digits histo 
   char ch[30];
-  sprintf(ch,"Track Segment_%d level %d ",GetID(),level );
+  //sprintf(ch,"Track Segment_%d level %d ",GetID(),level );
+  snprintf(ch,30,"Track Segment_%d level %d ",GetID(),level );
   if ( (fNrows<1)|| (fNcols<1)) {
     return 0;
   }
@@ -395,4 +424,37 @@ TClonesArray *  AliSimDigits::GenerTPCClonesArray(TClonesArray * arr)
     }    
   return digits;
 }
+void AliSimDigits::GlitchFilter(){
+  //
+  //  glitch filter, optionally
+  //
+  
+  for (Int_t i=0;i<fNcols;i++){ //pads
+    for(Int_t j=1;j<fNrows-1;j++){ //time bins
+      // first and last time bins are checked separately
+      if(GetDigitFast(j,i)){// nonzero digit
+        if (!GetDigitFast(j-1,i) && !GetDigitFast(j+1,i)) {
+          SetDigitFast(0,j,i);
+          SetTrackIDFast(-2,j,i,0);
+          SetTrackIDFast(-2,j,i,1);
+          SetTrackIDFast(-2,j,i,2);
+       }
+      }
+    }//time
+   
+    if(GetDigitFast(0,i) && !GetDigitFast(1,i)) {
+        SetDigitFast(0,0,i);
+        SetTrackIDFast(-2,0,i,0);
+        SetTrackIDFast(-2,0,i,1);
+        SetTrackIDFast(-2,0,i,2);
+    }
+    if(GetDigitFast(fNrows-1,i) && !GetDigitFast(fNrows-2,i)){ 
+       SetDigitFast(0,fNrows-1,i);
+       SetTrackIDFast(-2,fNrows-1,i,0);
+       SetTrackIDFast(-2,fNrows-1,i,1);
+       SetTrackIDFast(-2,fNrows-1,i,2);    
+    }
+  }//pads
+ 
+}