X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSdigit.cxx;h=1689dc8367ab87947f6ee7bdcef32f09531a6bc9;hb=af730ce7691c942cd93bc1e8110fcfc9e72a3b50;hp=80c9dbfcccd39b7416eaa7e1b10f35f60bb5e026;hpb=ac74f489d065b3e59af328b0b5ffb403fa6de594;p=u%2Fmrichter%2FAliRoot.git

diff --git a/ITS/AliITSdigit.cxx b/ITS/AliITSdigit.cxx
index 80c9dbfcccd..1689dc8367a 100644
--- a/ITS/AliITSdigit.cxx
+++ b/ITS/AliITSdigit.cxx
@@ -1,5 +1,5 @@
 /**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * Copyright(c) 2004-2006, ALICE Experiment at CERN, All rights reserved. *
  *                                                                        *
  * Author: The ALICE Off-line Project.                                    *
  * Contributors are mentioned in the code where appropriate.              *
@@ -17,416 +17,62 @@
 
 ////////////////////////////////////////////////
 //  Digits classes for all ITS detectors      //
+//                                            //
+//                                            //
 ////////////////////////////////////////////////
-#include <TObjArray.h>
-#include <TArrayI.h>
-#include <TArrayF.h>
-#include <TMath.h>
-#include "AliITSdigit.h"
-
-//______________________________________________________________________
-ClassImp(AliITSdigit)
-AliITSdigit::AliITSdigit(const Int_t *digits) {
-  // Creates a real data digit object
-
-  fCoord1       = digits[0];
-  fCoord2       = digits[1];
-  fSignal       = digits[2];
-}
-//______________________________________________________________________
-void AliITSdigit::Print(ostream *os){
-    //Standard output format for this class
-
-    *os << fCoord1 <<","<< fCoord2 <<","<< fSignal;
-}
-//______________________________________________________________________
-void AliITSdigit::Read(istream *os){
-    //Standard input for this class
-
-    *os >> fCoord1 >> fCoord2 >> fSignal;
-}
-//______________________________________________________________________
-ostream &operator<<(ostream &os,AliITSdigit &source){
-    // Standard output streaming function.
-
-    source.Print(&os);
-    return os;
-}
-//______________________________________________________________________
-istream &operator>>(istream &os,AliITSdigit &source){
-    // Standard output streaming function.
-
-    source.Read(&os);
-    return os;
-}
-//______________________________________________________________________
-ClassImp(AliITSdigitSPD)
-AliITSdigitSPD::AliITSdigitSPD():AliITSdigit(){
-    // default constructor, zero coordinates and set array
-    // elements to clearly unphysical values. A value of 0 may
-    // be a valide track of hit number.
-    Int_t i;
-
-    for(i=0;i<fkSspd;i++) fTracks[i]  = -3;
-    for(i=0;i<fkSspd;i++) fHits[i]    = -1;
-}
-//______________________________________________________________________
-AliITSdigitSPD::AliITSdigitSPD(const Int_t *digits){
-    // Creates a SPD digit object
-    Int_t i;
-
-    for(i=0;i<fkSspd;i++) fTracks[i]  = -3;
-    for(i=0;i<fkSspd;i++) fHits[i]    = -1;
-    fCoord1       = digits[0];
-    fCoord2       = digits[1];
-    fSignal       = 1;
-    fSignalSPD    = digits[2];
-}
-//______________________________________________________________________
-AliITSdigitSPD::AliITSdigitSPD(const Int_t *digits,const Int_t *tracks,
-			       const Int_t *hits){
-    // Creates a simulated SPD digit object
-
-    for(Int_t i=0; i<fkSspd; i++) {
-	fTracks[i] = tracks[i];
-	fHits[i]   = hits[i];
-    } // end for i
-    fCoord1       = digits[0];
-    fCoord2       = digits[1];
-    fSignal       = 1;
-    fSignalSPD    = digits[2];
-}
-//______________________________________________________________________
-Int_t AliITSdigitSPD::GetListOfTracks(TArrayI &t){
-    // Fills the TArrayI t with the tracks found in fTracks removing
-    // duplicated tracks, but otherwise in the same order. It will return
-    // the number of tracks and fill the remaining elements to the array
-    // t with -1.
-    // Inputs:
-    //   TArrayI  &t Reference to a TArrayI to contain the list of
-    //               nonduplicated track numbers.
-    // Output:
-    //   TArrayI  &t The input array filled with the nonduplicated track
-    //               numbers.
-    // Return:
-    //   Int_t The number of none -1 entries in the TArrayI t.
-    Int_t nt = t.GetSize();
-    Int_t nth = this->GetNTracks();
-    Int_t n = 0,i,j;
-    Bool_t inlist = kFALSE;
-
-    t.Reset(-1); // -1 array.
-    for(i=0;i<nth;i++) {
-	if(this->GetTrack(i) == -1) continue;
-	inlist = kFALSE;
-	for(j=0;j<n;j++)if(this->GetTrack(i) == t.At(j)) inlist = kTRUE;
-	if(!inlist){ // add to end of list
-	    t.AddAt(this->GetTrack(i),n);
-	    if(n<nt) n++;
-	} // end if
-    } // end for i
-    return n;
-}
-//______________________________________________________________________
-void AliITSdigitSPD::Print(ostream *os){
-    //Standard output format for this class
-    Int_t i;
-
-    AliITSdigit::Print(os);
-    for(i=0;i<fkSspd;i++) *os <<","<< fTracks[i];
-    for(i=0;i<fkSspd;i++) *os <<","<< fHits[i];
-    *os << "," << fSignalSPD;
-}
-//______________________________________________________________________
-void AliITSdigitSPD::Read(istream *os){
-    //Standard input for this class
-    Int_t i;
-
-    AliITSdigit::Read(os);
-    for(i=0;i<fkSspd;i++) *os >> fTracks[i];
-    for(i=0;i<fkSspd;i++) *os >> fHits[i];
-    *os >> fSignalSPD;
-}
-//______________________________________________________________________
-ostream &operator<<(ostream &os,AliITSdigitSPD &source){
-    // Standard output streaming function.
-
-    source.Print(&os);
-    return os;
-}
-//______________________________________________________________________
-istream &operator>>(istream &os,AliITSdigitSPD &source){
-    // Standard output streaming function.
-
-    source.Read(&os);
-    return os;
-}
-//______________________________________________________________________
-ClassImp(AliITSdigitSDD)
-AliITSdigitSDD::AliITSdigitSDD():AliITSdigit(){
-    // default constructor, zero coordinates and set array
-    // elements to clearly unphysical values. A value of 0 may
-    // be a valide track of hit number.
-    Int_t i;
 
-    for(i=0;i<fkSsdd;i++) fTracks[i] = -3;
-    for(i=0;i<fkSsdd;i++) fHits[i]   = -1;
-    fPhysics = 0;
-    for(i=0;i<fkSsdd;i++) fTcharges[i] = 0;
-}
-//________________________________________________________________________
-AliITSdigitSDD::AliITSdigitSDD(Float_t phys,const Int_t *digits):
-    AliITSdigit(digits){
-    // Creates a simulated SDD digit object to be updated
+#include "AliITSdigit.h"
 
-    fPhysics = phys;
-}
-//_____________________________________________________________________________
-AliITSdigitSDD::AliITSdigitSDD(Float_t phys,const Int_t *digits,
-			       const Int_t *tracks,const Int_t *hits,
-			       const Float_t *charges):
-    AliITSdigit(digits){
-    // Creates a simulated SDD digit object
 
-    fPhysics = phys;
-    for(Int_t i=0; i<fkSsdd; i++) {
-	fTcharges[i] = charges[i];
-	fTracks[i]   = tracks[i];
-	fHits[i]     = hits[i];
-    } // end for i
-}
 //______________________________________________________________________
-Int_t AliITSdigitSDD::GetListOfTracks(TArrayI &t,TArrayF &c){
-    // Fills the TArrayI t with the tracks found in fTracks removing
-    // duplicated tracks, summing up their charge, and ordering the tracks
-    // by the charge contributed to this digit. It will return
-    // the number of tracks and fill the remaining elements to the array
-    // t with -1.
-    // Inputs:
-    //   TArrayI  &t Reference to a TArrayI to contain the list of
-    //               nonduplicated track numbers.
-    //   TArrayF  &c Reference to a TArrayF to contain the summed charge
-    //               contributed by each track.
-    // Output:
-    //   TArrayI  &t The input array filled with the nonduplicated track
-    //               numbers.
-    //   TArrayF  &c The input array filled with the summed charge 
-    //               contributed by the corresponding track in the array t.
-    // Return:
-    //   Int_t The number of none -1 entries in the TArrayI t.
-    Int_t nt = t.GetSize();
-    nt = TMath::Min(nt,c.GetSize());
-    Int_t nth = this->GetNTracks();
-    Int_t n = 0,i,j;
-    Bool_t inlist = kFALSE;
-
-    t.Reset(-1); // -1 array.
-    c.Reset(0.0); // zero array.
-    for(i=0;i<nth;i++) {
-	if(this->GetTrack(i) == -1) continue;
-	inlist = kFALSE;
-	for(j=0;j<n;j++)if(this->GetTrack(i) == t.At(j)){
-	    inlist = kTRUE;
-	    c.AddAt(this->GetCharge(i)+c.At(j),j);
-	} // end for j/end if
-	if(!inlist){ // add to end of list
-	    t.AddAt(this->GetTrack(i),n);
-	    c.AddAt(this->GetCharge(i),n);
-	    if(n<nt) n++;
-	} // end if
-    } // end for i
+ClassImp(AliITSdigit)
 
-    // Now lets sort the TArrays according to the charge. This algorithm
-    // is based on the method from Chapter 8 section 1 Straight Insertion
-    // sort. Wiliam H. Press, Saul A. Teukolsky, William T. Vetterling
-    // and Brian P. Flannery, "Numerical Recipeis in C, The Art of Scientific
-    // Computing", second Edition page 330 (1997).
-    Int_t   tr;
-    Float_t ch;
-    for(i=0;i<n;i++){
-	tr = t.At(i);
-	ch = c.At(i);
-	j = i-1;
-	while(j>-1 && c.At(j)>ch){
-	    t.AddAt(t.At(j+1),j);
-	    c.AddAt(c.At(j+1),j);
-	    j--;
-	} // end while
-	t.AddAt(tr,j+1);
-	c.AddAt(ch,j+1);
-    } // end for i
-    //
-    return n;
+AliITSdigit::AliITSdigit():
+fCoord1(0),
+fCoord2(0),
+fSignal(0){
+  //default constructor. zero all values.
+  for(Int_t i=0;i<fgkSize;i++) fTracks[i] = 0;
+  for(Int_t i=0;i<fgkSize;i++) fHits[i]   = 0;
+  
 }
-//______________________________________________________________________
-void AliITSdigitSDD::Print(ostream *os){
-    //Standard output format for this class
-    Int_t i;
 
-    AliITSdigit::Print(os);
-    *os <<","<< fPhysics;
-    for(i=0; i<fkSsdd; i++) *os <<","<< fTcharges[i];
-    for(i=0; i<fkSsdd; i++) *os <<","<< fTracks[i];
-    for(i=0; i<fkSsdd; i++) *os <<","<< fHits[i];
-}
-//______________________________________________________________________
-void AliITSdigitSDD::Read(istream *os){
-    //Standard input for this class
-    Int_t i;
 
-    AliITSdigit::Read(os);
-    *os >>fPhysics;
-    for(i=0; i<fkSsdd; i++) *os >> fTcharges[i];
-    for(i=0; i<fkSsdd; i++) *os >> fTracks[i];
-    for(i=0; i<fkSsdd; i++) *os >> fHits[i];
-}
-//______________________________________________________________________
-ostream &operator<<(ostream &os,AliITSdigitSDD &source){
-    // Standard output streaming function.
+AliITSdigit::AliITSdigit(const Int_t *digits):
+fCoord1(digits[0]),
+fCoord2(digits[1]),
+fSignal(digits[2]){
+  // Creates a real data digit object
+  for(Int_t i=0;i<fgkSize;i++) fTracks[i] = 0;
+  for(Int_t i=0;i<fgkSize;i++) fHits[i]   = 0;
 
-    source.Print(&os);
-    return os;
-}
-//______________________________________________________________________
-istream &operator>>(istream &os,AliITSdigitSDD &source){
-    // Standard output streaming function.
 
-    source.Read(&os);
-    return os;
 }
 //______________________________________________________________________
-ClassImp(AliITSTransientDigit)
-AliITSTransientDigit::AliITSTransientDigit(Float_t phys,const Int_t *digits): 
-    AliITSdigitSDD(phys,digits) {
-    // Creates a digit object in a list of digits to be updated
-
-    fTrackList   = new TObjArray;  
-}
-//__________________________________________________________________________
-AliITSTransientDigit::AliITSTransientDigit(const AliITSTransientDigit &source):
- AliITSdigitSDD(source){
-    // Copy Constructor 
-    if(&source == this) return;
-    this->fTrackList = source.fTrackList;
-    return;
-}
-//_________________________________________________________________________
-AliITSTransientDigit& AliITSTransientDigit::operator=(
-    const AliITSTransientDigit &source) {
-    // Assignment operator
-    if(&source == this) return *this;
-    this->fTrackList = source.fTrackList;
-    return *this;
-}
-//______________________________________________________________________
-void AliITSTransientDigit::Print(ostream *os){
+void AliITSdigit::Print(ostream *os) {
     //Standard output format for this class
 
-    AliITSdigitSDD::Print(os);
+    *os << fCoord1 <<","<< fCoord2 <<","<< fSignal;
 }
 //______________________________________________________________________
-void AliITSTransientDigit::Read(istream *os){
+void AliITSdigit::Read(istream *os) {
     //Standard input for this class
 
-    AliITSdigitSDD::Read(os);
+    *os >> fCoord1 >> fCoord2 >> fSignal;
 }
 //______________________________________________________________________
-ostream &operator<<(ostream &os,AliITSTransientDigit &source){
+ostream &operator<<(ostream &os,AliITSdigit &source){
     // Standard output streaming function.
 
     source.Print(&os);
     return os;
 }
 //______________________________________________________________________
-istream &operator>>(istream &os,AliITSTransientDigit &source){
+istream &operator>>(istream &os,AliITSdigit &source){
     // Standard output streaming function.
 
     source.Read(&os);
     return os;
 }
-//______________________________________________________________________
-ClassImp(AliITSdigitSSD)
-AliITSdigitSSD::AliITSdigitSSD():AliITSdigit(){
-    // default constructor
-    Int_t i;
-
-    for(i=0; i<fkSssd; i++) fTracks[i] = -3;
-    for(i=0; i<fkSssd; i++) fHits[i] = -1;
-}
-//__________________________________________________________________________
-AliITSdigitSSD::AliITSdigitSSD(const Int_t *digits):AliITSdigit(digits){
-    // Creates a real SSD digit object
-}
-//_____________________________________________________________________________
-AliITSdigitSSD::AliITSdigitSSD(const Int_t *digits,const Int_t *tracks,
-			       const Int_t *hits):AliITSdigit(digits){
-    // Creates a simulated SSD digit object
-
-    for(Int_t i=0; i<fkSssd; i++) {
-	fTracks[i] = tracks[i];
-	fHits[i]   = hits[i];
-    } // end for i
-}
-//______________________________________________________________________
-Int_t AliITSdigitSSD::GetListOfTracks(TArrayI &t){
-    // Fills the TArrayI t with the tracks found in fTracks removing
-    // duplicated tracks, but otherwise in the same order. It will return
-    // the number of tracks and fill the remaining elements to the array
-    // t with -1.
-    // Inputs:
-    //   TArrayI  &t Reference to a TArrayI to contain the list of
-    //               nonduplicated track numbers.
-    // Output:
-    //   TArrayI  &t The input array filled with the nonduplicated track
-    //               numbers.
-    // Return:
-    //   Int_t The number of none -1 entries in the TArrayI t.
-    Int_t nt = t.GetSize();
-    Int_t nth = this->GetNTracks();
-    Int_t n = 0,i,j;
-    Bool_t inlist = kFALSE;
-
-    t.Reset(-1); // -1 array.
-    for(i=0;i<nth;i++) {
-	if(this->GetTrack(i) == -1) continue;
-	inlist = kFALSE;
-	for(j=0;j<n;j++)if(this->GetTrack(i) == t.At(j)) inlist = kTRUE;
-	if(!inlist){ // add to end of list
-	    t.AddAt(this->GetTrack(i),n);
-	    if(n<nt) n++;
-	} // end if
-    } // end for i
-    return n;
-}
-//______________________________________________________________________
-void AliITSdigitSSD::Print(ostream *os){
-    //Standard output format for this class
-    Int_t i;
-
-    AliITSdigit::Print(os);
-    for(i=0; i<fkSssd; i++) *os <<","<< fTracks[i];
-    for(i=0; i<fkSssd; i++) *os <<","<< fHits[i];
-}
-//______________________________________________________________________
-void AliITSdigitSSD::Read(istream *os){
-    //Standard input for this class
-    Int_t i;
 
-    AliITSdigit::Read(os);
-    for(i=0; i<fkSssd; i++) *os >> fTracks[i];
-    for(i=0; i<fkSssd; i++) *os >> fHits[i];
-}
-//______________________________________________________________________
-ostream &operator<<(ostream &os,AliITSdigitSSD &source){
-    // Standard output streaming function.
 
-    source.Print(&os);
-    return os;
-}
-//______________________________________________________________________
-istream &operator>>(istream &os,AliITSdigitSSD &source){
-    // Standard output streaming function.
-
-    source.Read(&os);
-    return os;
-}