Bug fix (selection of RawReaderDate from Yuri was not correct)

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

diff --git a/TPC/AliTPCROC.cxx b/TPC/AliTPCROC.cxx
index 1b8c39db4a624d346ca49aa1a724b954b5b33f8f..e72058edd77e564a21a3dddd6cd5311e568d1d83 100644 (file)
--- a/TPC/AliTPCROC.cxx
+++ b/TPC/AliTPCROC.cxx
@@ -197,13 +197,92 @@ void AliTPCROC::SetGeometry()
   //   SetOuterWWPitch(kOuterWWPitch);
   //   SetROuterFirstWire(kROuterFirstWire);
   //   SetROuterLastWire(kROuterLastWire);  
+
+  UInt_t i=0;
+  Float_t firstrow = fInnerRadiusLow + 1.575;   
+  for( i= 0;i<fNRowLow;i++)
+    {
+      Float_t x = firstrow + fInnerPadPitchLength*(Float_t)i;  
+      fPadRowLow[i]=x;
+      fYInner[i+1]  = x*TMath::Tan(fInnerAngle/2.)-fInnerWireMount;
+      fNPadsLow[i] = GetNPads(0,i) ;     // ROC implement     
+    }
+  // cross talk rows
+  fYInner[0]=(fPadRowLow[0]-fInnerPadPitchLength)*TMath::Tan(fInnerAngle/2.)-fInnerWireMount;
+  fYInner[fNRowLow+1]=(fPadRowLow[fNRowLow-1]+fInnerPadPitchLength)*TMath::Tan(fInnerAngle/2.)-fInnerWireMount; 
+  firstrow = fOuterRadiusLow + 1.6;
+  for(i=0;i<fNRowUp;i++)
+    {
+      if(i<fNRowUp1){
+       Float_t x = firstrow + fOuter1PadPitchLength*(Float_t)i; 
+       fPadRowUp[i]=x;
+       fYOuter[i+1]= x*TMath::Tan(fOuterAngle/2.)-fOuterWireMount;
+       fNPadsUp[i] =  GetNPads(36,i) ;     // ROC implement      
+       if(i==fNRowUp1-1) {
+         fLastWireUp1=fPadRowUp[i] +0.625;
+         firstrow = fPadRowUp[i] + 0.5*(fOuter1PadPitchLength+fOuter2PadPitchLength);
+       }
+      }
+      else
+       {
+         Float_t x = firstrow + fOuter2PadPitchLength*(Float_t)(i-64);
+         fPadRowUp[i]=x;
+         fNPadsUp[i] =  GetNPads(36,i) ;     // ROC implement
+       }
+      fYOuter[i+1]  = fPadRowUp[i]*TMath::Tan(fOuterAngle/2.)-fOuterWireMount;
+    }
+  
+
+
 } 
 
 
 
 
 //_____________________________________________________________________________
-AliTPCROC::AliTPCROC():TObject(), fNSectorsAll(0)
+AliTPCROC::AliTPCROC()
+          :TObject(), 
+           fNSectorsAll(0),
+          fInnerRadiusLow(0.),
+          fInnerRadiusUp(0.),
+          fOuterRadiusUp(0.),
+          fOuterRadiusLow(0.),
+          fInnerFrameSpace(0.),
+          fOuterFrameSpace(0.),
+          fInnerWireMount(0.),
+          fOuterWireMount(0.),
+          fZLength(0.),
+          fInnerAngle(0.),
+          fOuterAngle(0.),
+          fNInnerWiresPerPad(0),
+          fInnerWWPitch(0.),
+          fInnerDummyWire(0),
+          fInnerOffWire(0.),
+          fRInnerFirstWire(0.),
+          fRInnerLastWire(0.),
+          fLastWireUp1(0.),
+          fNOuter1WiresPerPad(0),
+          fNOuter2WiresPerPad(0),
+          fOuterWWPitch(0.),
+          fOuterDummyWire(0),
+          fOuterOffWire(0),
+          fROuterFirstWire(0.),
+          fROuterLastWire(0),
+          fInnerPadPitchLength(0.),
+          fInnerPadPitchWidth(0.),
+          fInnerPadLength(0.),
+          fInnerPadWidth(0.),
+          fOuter1PadPitchLength(0.),
+          fOuter2PadPitchLength(0),
+          fOuterPadPitchWidth(0),
+          fOuter1PadLength(0.),
+          fOuter2PadLength(0),
+          fOuterPadWidth(0),
+          fNRowLow(0),
+          fNRowUp1(0),
+          fNRowUp2(0),
+          fNRowUp(0),
+          fNtRows(0)
 {
   //
   // Default constructor
@@ -218,7 +297,50 @@ AliTPCROC::AliTPCROC():TObject(), fNSectorsAll(0)
 
 
 //_____________________________________________________________________________
-AliTPCROC::AliTPCROC(const AliTPCROC &roc):TObject(roc)
+AliTPCROC::AliTPCROC(const AliTPCROC &roc)
+          :TObject(roc),
+           fNSectorsAll(0),
+          fInnerRadiusLow(0.),
+          fInnerRadiusUp(0.),
+          fOuterRadiusUp(0.),
+          fOuterRadiusLow(0.),
+          fInnerFrameSpace(0.),
+          fOuterFrameSpace(0.),
+          fInnerWireMount(0.),
+          fOuterWireMount(0.),
+          fZLength(0.),
+          fInnerAngle(0.),
+          fOuterAngle(0.),
+          fNInnerWiresPerPad(0),
+          fInnerWWPitch(0.),
+          fInnerDummyWire(0),
+          fInnerOffWire(0.),
+          fRInnerFirstWire(0.),
+          fRInnerLastWire(0.),
+          fLastWireUp1(0.),
+          fNOuter1WiresPerPad(0),
+          fNOuter2WiresPerPad(0),
+          fOuterWWPitch(0.),
+          fOuterDummyWire(0),
+          fOuterOffWire(0),
+          fROuterFirstWire(0.),
+          fROuterLastWire(0),
+          fInnerPadPitchLength(0.),
+          fInnerPadPitchWidth(0.),
+          fInnerPadLength(0.),
+          fInnerPadWidth(0.),
+          fOuter1PadPitchLength(0.),
+          fOuter2PadPitchLength(0),
+          fOuterPadPitchWidth(0),
+          fOuter1PadLength(0.),
+          fOuter2PadLength(0),
+          fOuterPadWidth(0),
+          fNRowLow(0),
+          fNRowUp1(0),
+          fNRowUp2(0),
+          fNRowUp(0),
+          fNtRows(0)
+
 {
   //
   // AliTPCROC copy constructor
@@ -249,7 +371,15 @@ AliTPCROC::AliTPCROC(const AliTPCROC &roc):TObject(roc)
     fRowPosIndex[1][irow] = roc.fRowPosIndex[1][irow];
   }
 }
-
+//____________________________________________________________________________
+AliTPCROC & AliTPCROC::operator =(const AliTPCROC & roc)
+{
+  //
+  // assignment operator - dummy
+  //
+  fZLength = roc.fZLength;
+  return (*this);
+}
 //_____________________________________________________________________________
 AliTPCROC::~AliTPCROC()
 {
@@ -264,3 +394,36 @@ AliTPCROC::~AliTPCROC()
   
 }
 
+
+
+
+void AliTPCROC::GetPositionLocal(UInt_t sector, UInt_t row, UInt_t pad, Float_t *pos){
+  //
+  // get position of center of pad - ideal frame used
+  //
+  pos[2]=fZLength;
+  if (sector<36){
+   pos[0] = fPadRowLow[row];
+   pos[1] = fInnerPadPitchWidth*(Int_t(pad)-Int_t(fNPads[0][row])/2);
+  }else{
+    pos[0] = fPadRowUp[row];
+    pos[1] = fOuterPadPitchWidth*(Int_t(pad)-Int_t(fNPads[1][row])/2);    
+  }
+  if ((sector%36)>=18){
+    pos[2] *= -1.;
+    pos[1] *= -1.;
+  }
+}
+
+
+void AliTPCROC::GetPositionGlobal(UInt_t sector, UInt_t row, UInt_t pad, Float_t *pos){
+  //
+  // get position of center of pad - ideal frame used 
+  //
+  GetPositionLocal(sector,row,pad,pos);
+  Double_t alpha = TMath::DegToRad()*(10.+20.*(sector%18));
+  Float_t gx = pos[0]*TMath::Cos(alpha)-pos[1]*TMath::Sin(alpha);
+  Float_t gy = pos[1]*TMath::Cos(alpha)+pos[0]*TMath::Sin(alpha);
+  pos[0] = gx;
+  pos[1] = gy;
+}