Go from pointer to ifstream to ifstream.

[u/mrichter/AliRoot.git] / MUON / mapping / AliMpRowSegment.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

// $Id$
// $MpId: AliMpRowSegment.cxx,v 1.10 2006/05/24 13:58:46 ivana Exp $
// Category: sector

//-----------------------------------------------------------------------------
// Class AliMpRowSegment
// ---------------------
// Class describing a row segment composed of the 
// the identic motifs.
// Included in AliRoot: 2003/05/02
// Authors: David Guez, Ivana Hrivnacova; IPN Orsay
//-----------------------------------------------------------------------------

#include "AliMpRowSegment.h"
#include "AliMpRow.h"
#include "AliMpVMotif.h"
#include "AliMpMotifType.h"
#include "AliMpMotifTypePadIterator.h"
#include "AliMpMotifMap.h"
#include "AliMpMotifPosition.h"
#include "AliMpConstants.h"
#include "AliMpEncodePair.h"

#include "AliLog.h"

#include <TMath.h>
#include <Riostream.h>

using std::endl;
/// \cond CLASSIMP
ClassImp(AliMpRowSegment)
/// \endcond

//_____________________________________________________________________________
AliMpRowSegment::AliMpRowSegment(AliMpRow* row, AliMpVMotif* motif, 
                                 Int_t padOffsetX, Int_t padOffsetY,
                                 Int_t nofMotifs,
                                 Int_t motifPositionId, Int_t motifPositionDId)
  : AliMpVRowSegment(),
    fNofMotifs(nofMotifs),
    fLPadOffset(AliMp::Pair(padOffsetX,padOffsetY)),
    fOffsetX(0.),
    fOffsetY(0.),
    fRow(row),
    fMotif(motif),
    fMotifPositionId(motifPositionId),
    fMotifPositionDId(motifPositionDId)
{
/// Standard constructor
 
  // Keep pad offset in the low indices limits
  SetLowIndicesLimit(fLPadOffset);
}

//_____________________________________________________________________________
AliMpRowSegment::AliMpRowSegment() 
  : AliMpVRowSegment(),
    fNofMotifs(0),
    fLPadOffset(0),
    fOffsetX(0.),
    fOffsetY(0.),
    fRow(0),
    fMotif(0),
    fMotifPositionId(0),
    fMotifPositionDId(0)
{
/// Default constructor
}

//_____________________________________________________________________________
AliMpRowSegment::~AliMpRowSegment() 
{
/// Destructor  
}

//
// private methods  
//

//_____________________________________________________________________________
Double_t AliMpRowSegment::FirstMotifCenterX() const
{
/// Return the x coordinate of the first motif center
/// in the global coordinate system.

  return fOffsetX;
}  

//_____________________________________________________________________________
Double_t AliMpRowSegment::LastMotifCenterX() const
{
/// Return the x coordinate of the last motif center
/// in the global coordinate system.

  return fOffsetX + 2.*(fNofMotifs-1)*fMotif->DimensionX();
}

//_____________________________________________________________________________
Double_t AliMpRowSegment::MotifCenterX(Int_t motifPositionId) const
{
/// Return the x coordinate of the motif specified with
/// the given position identifier.

  // Check if x is in the row segment range
  if (! HasMotifPosition(motifPositionId)) {
    AliErrorStream() << "Outside row segment region" << endl;
    return 0;
  }
  
  // Find the position number in the segment  
  Int_t num = (motifPositionId - fMotifPositionId) *  fMotifPositionDId;

  return fOffsetX + num*(fMotif->DimensionX() * 2.0);
}

//_____________________________________________________________________________
Double_t AliMpRowSegment::MotifCenterY(Int_t motifPositionId) const
{
/// Return the y coordinate of the motif specified with
/// the given position identifier.

  // Check if x is in the row segment range
  if (! HasMotifPosition(motifPositionId)) {
    AliErrorStream() << "Outside row segment region" << endl;
    return 0;
  }
  
  return GetRow()->GetPositionY() + fOffsetY;
}

//_____________________________________________________________________________
Bool_t AliMpRowSegment::IsInside(Double_t x, Double_t y, Bool_t warn) const
{
/// Check if the position is inside some motif of this row segment.

  Double_t minY = GetRow()->GetPositionY() + fOffsetY - fMotif->DimensionY();
  Double_t maxY = GetRow()->GetPositionY() + fOffsetY + fMotif->DimensionY();

  if ( x < LeftBorderX() || x > RightBorderX() ||
       y < minY || y > maxY ) {

    if (warn)
      AliWarningStream() << "Outside row segment region" << endl;
    return false;
  }
  else
    return true;
}    

//
// public methods  
//

//_____________________________________________________________________________
Double_t  AliMpRowSegment::LeftBorderX() const
{
/// Return the x coordinate of the left row segment border
/// in the global coordinate system.

  return FirstMotifCenterX() - fMotif->DimensionX();
}

//_____________________________________________________________________________
Double_t  AliMpRowSegment::RightBorderX() const
{
/// Return the x coordinate of the right row segment border
/// in the global coordinate system.

  return LastMotifCenterX() + fMotif->DimensionX();
}

//_____________________________________________________________________________
Double_t  AliMpRowSegment::HalfSizeY() const
{
/// Return the size in y of this row segment.

  return fMotif->DimensionY() + fOffsetY;
}

//_____________________________________________________________________________
AliMpVMotif*  AliMpRowSegment::FindMotif(Double_t x, Double_t y) const
{
/// Return the motif of this row; 

  if ( IsInside(x, y, false) )
    return fMotif;
  else  
    return 0;
}  

//_____________________________________________________________________________
Int_t AliMpRowSegment::FindMotifPositionId(Double_t x, Double_t y) const
{
/// Return the motif position identified for the given
/// geometric position.

  if ( ! IsInside(x, y, false) ) return 0;

  // Find the position number in the segment  
  Int_t num 
    = Int_t((x - LeftBorderX()) / (fMotif->DimensionX() * 2.0));

  // Calculate the position Id
  return fMotifPositionId + num*fMotifPositionDId;  
}

//_____________________________________________________________________________
Bool_t AliMpRowSegment::HasMotifPosition(Int_t motifPositionId) const
{
/// Return true if the motif specified with the given position identifier
/// is in this segment.

  Int_t minId = TMath::Min(fMotifPositionId, 
                    fMotifPositionId + (fNofMotifs-1)*fMotifPositionDId);
  Int_t maxId = TMath::Max(fMotifPositionId, 
                    fMotifPositionId + (fNofMotifs-1)*fMotifPositionDId);

  if (motifPositionId >= minId && motifPositionId <= maxId) {
    return true;
  }
  else 
    return false;
}

//_____________________________________________________________________________
void AliMpRowSegment::MotifCenter(Int_t motifPositionId,
                                  Double_t& x, Double_t& y) const
{
/// Return the coordinates of the motif specified with
/// the given position identifier.

  x = MotifCenterX(motifPositionId);
  y = MotifCenterY(motifPositionId);
}

//_____________________________________________________________________________
Double_t AliMpRowSegment::GetPositionX() const
{
/// Return the x position of the row segment centre.

  return (LeftBorderX() + RightBorderX())/2.;		    
}

//_____________________________________________________________________________
Double_t AliMpRowSegment::GetPositionY() const
{
/// Return the y position of the row segment centre.

  return GetRow()->GetPositionY();  
}

//_____________________________________________________________________________
Double_t AliMpRowSegment::GetDimensionX() const
{
/// Return the halflengths of the row segment in x, y.
// ---

  return (RightBorderX() - LeftBorderX())/2.;		    
}

//_____________________________________________________________________________
Double_t AliMpRowSegment::GetDimensionY() const
{
/// Return the halflengths of the row segment in x, y.
// ---

  return GetRow()->GetDimensionY();  
}

//_____________________________________________________________________________
void   AliMpRowSegment::SetOffset(Double_t x, Double_t y)
{
/// Calculate offset from given offset and 
/// stored offset in pads.

  AliMpMotifTypePadIterator iter(fMotif->GetMotifType());
  iter.First();

  Int_t ix = iter.CurrentItem().GetIx();
  Int_t iy = iter.CurrentItem().GetIy();
  
  Double_t dx, dy;
  fMotif->GetPadDimensionsByIndices(ix, iy, dx, dy);  

  fOffsetX 
     = x + 2.*AliMp::PairFirst(fLPadOffset) * dx + fMotif->DimensionX(); 

  fOffsetY
    = y + AliMp::PairSecond(fLPadOffset) * dy; 
}

//_____________________________________________________________________________
void AliMpRowSegment::SetGlobalIndices(AliMpRow* /*rowBefore*/)
{
/// Set global indices limits.

  // The low/high indices limits has to be taken as the highest/lowest from all 
  // motif positions
  Int_t ixl = 9999;
  Int_t iyl = 9999;
  Int_t ixh = AliMpConstants::StartPadIndex();
  Int_t iyh = AliMpConstants::StartPadIndex();

  for (Int_t i=0; i<GetNofMotifs(); i++) {
     
     AliMpMotifPosition* mPos 
       = GetRow()->GetMotifMap()->FindMotifPosition(GetMotifPositionId(i));
       
     // Check if the motif positions has the limits set
     if ( !mPos->HasValidIndices() )
       Fatal("SetGlobalIndices", 
             "Indices of motif positions have to be set first.");
            
     if ( mPos->GetLowLimitIx() < ixl ) 
       ixl = mPos->GetLowLimitIx();
       
     if ( mPos->GetLowLimitIy() < iyl ) 
       iyl = mPos->GetLowLimitIy();

     if ( mPos->GetHighLimitIx() > ixh ) 
       ixh = mPos->GetHighLimitIx();
       
     if ( mPos->GetHighLimitIy() > iyh ) 
       iyh = mPos->GetHighLimitIy();
  }     

  SetLowIndicesLimit(ixl, iyl);
  SetHighIndicesLimit(ixh, iyh);
}  

//_____________________________________________________________________________
Int_t AliMpRowSegment::SetIndicesToMotifPosition(Int_t i, MpPair_t indices)
{
/// Set global indices to i-th motif position and returns next index
/// in x.

  // Get motif position
  AliMpMotifPosition* motifPosition
    = GetRow()->GetMotifMap()->FindMotifPosition(GetMotifPositionId(i));

  // Low limit
  MpPair_t low = indices + AliMp::Pair(0, GetLowLimitIy());
  motifPosition->SetLowIndicesLimit(low);
	  
  // High limit
  AliMpMotifType* motifType = motifPosition->GetMotif()->GetMotifType();    
  MpPair_t high
    = motifPosition->GetLowIndicesLimit()
      + AliMp::Pair(motifType->GetNofPadsX()-1, motifType->GetNofPadsY()-1);
  motifPosition->SetHighIndicesLimit(high);

  // Return next index in x
  return AliMp::PairFirst(high)+1;
  // return motifType->GetNofPadsX();
}


//_____________________________________________________________________________
AliMpRow*  AliMpRowSegment::GetRow() const
{
/// Return the row.which this row segment belongs to.

  return fRow;
}  

//_____________________________________________________________________________
Int_t  AliMpRowSegment::GetNofMotifs() const
{
/// Return number of motifs in this this row segment.

  return fNofMotifs;
}  

//_____________________________________________________________________________
Int_t  AliMpRowSegment::GetMotifPositionId(Int_t i) const
{
/// Return number of motifs in this this row segment.

  return fMotifPositionId + i*fMotifPositionDId;
}  

//_____________________________________________________________________________
AliMpVMotif*  AliMpRowSegment::GetMotif(Int_t /*i*/) const
{
/// Return the motif of this row segment.

  return fMotif;
}