**************************************************************************/
// $Id$
-// $MpId: AliMpPCB.cxx,v 1.5 2006/01/11 10:14:17 ivana Exp $
+// $MpId: AliMpPCB.cxx,v 1.8 2006/05/24 13:58:50 ivana Exp $
#include "AliMpPCB.h"
#include "AliMpMotif.h"
+#include "AliMpSlatMotifMap.h"
#include "AliMpMotifPosition.h"
#include "AliMpMotifSpecial.h"
#include "AliMpMotifType.h"
#include "AliLog.h"
#include "Riostream.h"
+#include "TList.h"
+#include "TObjString.h"
#include "TMath.h"
#include <sstream>
+///
+/// \class AliMpPCB
+///
+/// A PCB for station 3,4 or 5
+///
+/// A PCB is a group of pads having the same size
+/// Pads are grouped in motifs, where 1 motif = 1 MANU
+///
+/// The notion of PCB enveloppe is due to the fact that not all PCBs are
+/// "full" of pads, e.g. the rounded or short ones miss some pads,
+/// but the enveloppe is a virtual size that should be constant
+/// across the slats, and is 400x400 mm.
+/// It's a usefull notion to compute e.g. slat center in a uniform way,
+/// considering that a slat is N PCBs, of the same "virtual" size, that of
+/// the enveloppe.
+///
+/// \author L. Aphecetche
+
+/// \cond CLASSIMP
ClassImp(AliMpPCB)
+/// \endcond
//_____________________________________________________________________________
AliMpPCB::AliMpPCB()
- : TObject(), fId(""),
- fPadSizeX(0), fPadSizeY(0),
- fEnveloppeSizeX(0), fEnveloppeSizeY(0),
- fXoffset(0),
- fActiveXmin(0), fActiveXmax(0),
- fIxmin(99999), fIxmax(0), fIymin(99999), fIymax(0)
+: TObject(),
+ fId(""),
+ fPadSizeX(0),
+ fPadSizeY(0),
+ fEnveloppeSizeX(0),
+ fEnveloppeSizeY(0),
+ fXoffset(0),
+ fActiveXmin(0),
+ fActiveXmax(0),
+ fIxmin(99999),
+ fIxmax(0),
+ fIymin(99999),
+ fIymax(0),
+ fMotifPositions(),
+ fNofPads(0),
+ fMotifMap(0)
{
//
// Default ctor.
//
+#ifdef WITH_ROOT
+ fMotifPositions.SetOwner(kTRUE);
+#endif
+ AliDebug(1,Form("this=%p",this));
}
//_____________________________________________________________________________
-AliMpPCB::AliMpPCB(const char* id, Double_t padSizeX, Double_t padSizeY,
+AliMpPCB::AliMpPCB(AliMpSlatMotifMap* motifMap, const char* id, Double_t padSizeX, Double_t padSizeY,
Double_t enveloppeSizeX, Double_t enveloppeSizeY)
- : TObject(), fId(id),
- fPadSizeX(padSizeX), fPadSizeY(padSizeY),
- fEnveloppeSizeX(enveloppeSizeX), fEnveloppeSizeY(enveloppeSizeY),
- fXoffset(0),
- fActiveXmin(0), fActiveXmax(0),
- fIxmin(99999), fIxmax(0), fIymin(99999), fIymax(0)
+: TObject(),
+ fId(id),
+ fPadSizeX(padSizeX),
+ fPadSizeY(padSizeY),
+ fEnveloppeSizeX(enveloppeSizeX),
+ fEnveloppeSizeY(enveloppeSizeY),
+ fXoffset(0),
+ fActiveXmin(0),
+ fActiveXmax(0),
+ fIxmin(99999),
+ fIxmax(0),
+ fIymin(99999),
+ fIymax(0),
+ fMotifPositions(),
+ fNofPads(0),
+ fMotifMap(motifMap)
{
//
// Normal ctor. Must be fed with the PCB's name (id), the pad dimensions
// and the global dimension of the virtual enveloppe of the PCB
// (usually 400x400 mm)
+#ifdef WITH_ROOT
+ fMotifPositions.SetOwner(kTRUE);
+#endif
+ AliDebug(1,Form("this=%p id=%s",this,id));
}
//_____________________________________________________________________________
AliMpPCB::AliMpPCB(const AliMpPCB& o)
- : TObject(o),
- fPadSizeX(0), fPadSizeY(0),
- fEnveloppeSizeX(0), fEnveloppeSizeY(0),
- fXoffset(0),
- fActiveXmin(0), fActiveXmax(0),
- fIxmin(99999), fIxmax(0), fIymin(99999), fIymax(0)
-{
+: TObject(o),
+ fId(0),
+ fPadSizeX(0),
+ fPadSizeY(0),
+ fEnveloppeSizeX(0),
+ fEnveloppeSizeY(0),
+ fXoffset(0),
+ fActiveXmin(0),
+ fActiveXmax(0),
+ fIxmin(99999),
+ fIxmax(0),
+ fIymin(99999),
+ fIymax(0),
+ fMotifPositions(),
+ fNofPads(0),
+ fMotifMap(0x0)
+{
+ //
+ // Copy constructor
+
+#ifdef WITH_ROOT
+ fMotifPositions.SetOwner(kTRUE);
+#endif
+ AliDebug(1,Form("this=%p (copy ctor) : begin",this));
o.Copy(*this);
+ AliDebug(1,Form("this=%p (copy ctor) : end",this));
}
//_____________________________________________________________________________
AliMpPCB::AliMpPCB(const char* id, AliMpMotifSpecial* ms)
-: TObject(), fId(id), fPadSizeX(-1.0), fPadSizeY(-1.0),
- fXoffset(0)
+: TObject(),
+ fId(id),
+ fPadSizeX(-1.0),
+ fPadSizeY(-1.0),
+ fEnveloppeSizeX(ms->Dimensions().X()*2.0),
+ fEnveloppeSizeY(ms->Dimensions().Y()*2.0),
+ fXoffset(0.0),
+ fActiveXmin(0.0),
+ fActiveXmax(fEnveloppeSizeX),
+ fIxmin(0),
+ fIxmax(ms->GetMotifType()->GetNofPadsX()-1),
+ fIymin(0),
+ fIymax(ms->GetMotifType()->GetNofPadsY()-1),
+ fMotifPositions(),
+ fNofPads(ms->GetMotifType()->GetNofPads()),
+ fMotifMap(0x0)
{
//
// Very special ctor to be used by trigger stations only (and for a very
// This limitation might not be justified, except that it's all we need
// so far ;-)
//
- fXoffset = 0.0;
- fEnveloppeSizeX = ms->Dimensions().X()*2.0;
- fEnveloppeSizeY = ms->Dimensions().Y()*2.0;
- fActiveXmin = 0.0;
- fActiveXmax = fEnveloppeSizeX;
- fIxmin = fIymin = 0;
- fIxmax = ms->GetMotifType()->GetNofPadsX()-1;
- fIymax = ms->GetMotifType()->GetNofPadsY()-1;
+
+ AliDebug(1,Form("this=%p (ctor special motif)",this));
+
+#ifdef WITH_ROOT
+ fMotifPositions.SetOwner(kTRUE);
+#endif
TVector2 position(ms->Dimensions());
AliMpMotifPosition* mp = new AliMpMotifPosition(-1,ms,position);
mp->SetLowIndicesLimit(AliMpIntPair(fIxmin,fIymin));
mp->SetHighIndicesLimit(AliMpIntPair(fIxmax,fIymax));
#ifdef WITH_ROOT
- fMotifs.AddLast(mp);
+ fMotifPositions.AddLast(mp);
#else
- fMotifs.push_back(mp);
+ fMotifPositions.push_back(mp);
#endif
}
AliMpPCB&
AliMpPCB::operator=(const AliMpPCB& o)
{
+ // Assignment operator
+
+ AliDebug(1,Form("this=%p (assignment op) : begin",this));
o.Copy(*this);
+ AliDebug(1,Form("this=%p (assignment op) : end",this));
return *this;
}
//
// Dtor.
//
+ AliDebug(1,Form("this=%p",this));
+#ifndef WITH_ROOT
+ for ( size_t i = 0; i < fMotifPositions.size(); ++i )
+ {
+ delete fMotifPositions[i];
+ }
+#endif
+
}
//_____________________________________________________________________________
// (ix<0,iy<0) : top-right corner
// (ix>0,iy<0) : top-left corner.
- AliMpVMotif* motif =
- new AliMpMotif(mt->GetID(),mt,TVector2(PadSizeX()/2.0,PadSizeY()/2.0));
+ TString id(Form("%s-%e-%e",mt->GetID().Data(),PadSizeX(),PadSizeY()));
+
+ AliMpVMotif* motif = fMotifMap->FindMotif(id);
+
+ if (!motif)
+ {
+ motif = new AliMpMotif(id,mt,TVector2(PadSizeX()/2.0,PadSizeY()/2.0));
+ AliDebug(1,Form("Adding motif %s to motifMap",id.Data()));
+ fMotifMap->AddMotif(motif);
+ }
+ else
+ {
+ AliDebug(1,Form("Got motif %s from motifMap",id.Data()));
+ }
+
TVector2 position;
Int_t ixmin(-1);
Int_t iymin(-1);
mp->SetHighIndicesLimit(AliMpIntPair(ixmax,iymax));
#ifdef WITH_ROOT
- fMotifs.AddLast(mp);
+ fMotifPositions.AddLast(mp);
#else
- fMotifs.push_back(mp);
+ fMotifPositions.push_back(mp);
#endif
- fIxmin = std::min(fIxmin,ixmin);
- fIxmax = std::max(fIxmax,ixmax);
- fIymin = std::min(fIymin,iymin);
- fIymax = std::max(fIymax,iymax);
+ fIxmin = TMath::Min(fIxmin,ixmin);
+ fIxmax = TMath::Max(fIxmax,ixmax);
+ fIymin = TMath::Min(fIymin,iymin);
+ fIymax = TMath::Max(fIymax,iymax);
fActiveXmin = fIxmin*PadSizeX();
fActiveXmax = (fIxmax+1)*PadSizeX();
+ fNofPads += mt->GetNofPads();
+}
+
+//_____________________________________________________________________________
+AliMpArea
+AliMpPCB::Area() const
+{
+ return AliMpArea(TVector2( (Xmin()+Xmax())/2.0,DY()),
+ TVector2( DX(), DY() ) );
}
//_____________________________________________________________________________
//
// Return a full copy of this object.
//
- return new AliMpPCB(*this);
+ AliDebug(1,"begin");
+ TObject* object = new AliMpPCB(*this);
+ AliDebug(1,"end");
+ return object;
}
//_____________________________________________________________________________
// pcb of a slat will start at ix=0 (only relevant for rounded pcbs).
//
+ AliDebug(1,"begin");
+
// First get a full clone.
AliMpPCB* pcb = static_cast<AliMpPCB*>(Clone());
return 0;
}
- AliMpIntPair shift(-fIxmin,0);
+ AliMpIntPair shift(-fIxmin+ixOffset,0);
// Then change the internal MotifPositions wrt manu id
// and position (offset in x).
{
AliMpMotifPosition* mp = pcb->GetMotifPosition(i);
mp->SetID(manuids[i]);
- mp->SetPosition(mp->Position() + TVector2(xOffset,0));
- mp->SetLowIndicesLimit(mp->GetLowIndicesLimit()+
- shift+
- AliMpIntPair(ixOffset,0));
- mp->SetHighIndicesLimit(mp->GetHighIndicesLimit()+
- shift+
- AliMpIntPair(ixOffset,0));
+ TVector2 pos(mp->Position());
+ pos += TVector2(xOffset,0);
+ mp->SetPosition(pos);
+ AliMpIntPair offset(ixOffset,0);
+ AliMpIntPair low(mp->GetLowIndicesLimit());
+ low += shift;
+ mp->SetLowIndicesLimit(low);
+ AliMpIntPair high(mp->GetHighIndicesLimit());
+ high += shift;
+ mp->SetHighIndicesLimit(high);
}
- pcb->fIxmin += ixOffset + shift.GetFirst();
- pcb->fIxmax += ixOffset + shift.GetFirst();
+ pcb->fIxmin += shift.GetFirst();
+ pcb->fIxmax += shift.GetFirst();
pcb->fXoffset = xOffset;
pcb->fActiveXmin += xOffset;
pcb->fActiveXmax += xOffset;
+ AliDebug(1,"end");
+
return pcb;
}
void
AliMpPCB::Copy(TObject& o) const
{
+ // Copy *this into o
+
+ AliDebug(1,"begin");
+
TObject::Copy(o);
AliMpPCB& pcb = static_cast<AliMpPCB&>(o);
pcb.fId = fId;
pcb.fActiveXmax = fActiveXmax;
#ifdef WITH_ROOT
- pcb.fMotifs.Clear();
+ AliDebug(1,"Deleting pcb.fMotifPositions");
+ pcb.fMotifPositions.Delete();
+ AliDebug(1,"Deleting pcb.fMotifPositions : done");
#else
- pcb.fMotifs.clear();
+ for ( Size_t i = 0; i < pcb.fMotifPositions.size(); ++i )
+ {
+ delete pcb.fMotifPositions[i];
+ }
#endif
#ifdef WITH_ROOT
- for ( Size_t i = 0; i < fMotifs.GetEntriesFast(); ++i )
+ for ( Size_t i = 0; i < fMotifPositions.GetEntriesFast(); ++i )
#else
- for ( Size_t i = 0; i < fMotifs.size(); ++i )
+ for ( Size_t i = 0; i < fMotifPositions.size(); ++i )
#endif
{
- AliMpMotifPosition* pos = (AliMpMotifPosition*)fMotifs[i];
- AliMpMotifPosition* pcbpos =
- new AliMpMotifPosition(pos->GetID(),pos->GetMotif(),pos->Position());
+ AliMpMotifPosition* pos = (AliMpMotifPosition*)fMotifPositions[i];
+ AliMpMotifPosition* pcbpos = new AliMpMotifPosition(pos->GetID(),
+ pos->GetMotif(),
+ pos->Position());
pcbpos->SetLowIndicesLimit(pos->GetLowIndicesLimit());
pcbpos->SetHighIndicesLimit(pos->GetHighIndicesLimit());
#ifdef WITH_ROOT
- pcb.fMotifs.AddLast(pcbpos);
+ pcb.fMotifPositions.AddLast(pcbpos);
#else
- pcb.fMotifs.push_back(pcbpos);
+ pcb.fMotifPositions.push_back(pcbpos);
#endif
}
+
+ pcb.fNofPads = fNofPads;
+
+ pcb.fMotifMap = fMotifMap; // warning : we do share the motifmap.
+
+ AliDebug(1,"end");
}
//_____________________________________________________________________________
//
#ifdef WITH_ROOT
- for (Size_t i = 0; i < fMotifs.GetEntriesFast(); ++i )
+ for (Size_t i = 0; i < fMotifPositions.GetEntriesFast(); ++i )
#else
- for (Size_t i = 0; i < fMotifs.size(); ++i )
+ for (Size_t i = 0; i < fMotifPositions.size(); ++i )
#endif
{
- AliMpMotifPosition* mp = (AliMpMotifPosition*)fMotifs[i];
+ AliMpMotifPosition* mp = (AliMpMotifPosition*)fMotifPositions[i];
if ( mp->HasPad(AliMpIntPair(ix,iy)) )
{
return mp;
//
#ifdef WITH_ROOT
- for (Size_t i = 0; i < fMotifs.GetEntriesFast(); ++i )
+ for (Size_t i = 0; i < fMotifPositions.GetEntriesFast(); ++i )
#else
- for (Size_t i = 0; i < fMotifs.size(); ++i )
+ for (Size_t i = 0; i < fMotifPositions.size(); ++i )
#endif
{
- AliMpMotifPosition* mp = (AliMpMotifPosition*)fMotifs[i];
+ AliMpMotifPosition* mp = (AliMpMotifPosition*)fMotifPositions[i];
TVector2 localPos( TVector2(x,y) - mp->Position() );
//
#ifdef WITH_ROOT
- if ( i >= fMotifs.GetEntriesFast() ) return 0;
+ if ( i >= fMotifPositions.GetEntriesFast() ) return 0;
#else
- if ( i >= fMotifs.size() ) return 0;
+ if ( i >= fMotifPositions.size() ) return 0;
#endif
- return (AliMpMotifPosition*)fMotifs[i];
+ return (AliMpMotifPosition*)fMotifPositions[i];
}
//_____________________________________________________________________________
//
#ifdef WITH_ROOT
- return fMotifs.GetEntriesFast();
+ return fMotifPositions.GetEntriesFast();
#else
- return fMotifs.size();
+ return fMotifPositions.size();
#endif
}
if ( option && option[0] == 'M' )
{
#ifdef WITH_ROOT
- for ( Size_t i = 0; i < fMotifs.GetEntriesFast(); ++i )
+ for ( Size_t i = 0; i < fMotifPositions.GetEntriesFast(); ++i )
#else
- for ( Size_t i = 0; i < fMotifs.size(); ++i )
+ for ( Size_t i = 0; i < fMotifPositions.size(); ++i )
#endif
{
if (option)
{
- fMotifs[i]->Print(option+1);
+ fMotifPositions[i]->Print(option+1);
}
else
{
- fMotifs[i]->Print();
+ fMotifPositions[i]->Print();
}
}
}
}
+//_____________________________________________________________________________
+void
+AliMpPCB::Save() const
+{
+ //
+ // Save this PCB in the file
+
+ TString fileName(fId);
+ fileName += ".pcb";
+ TList lines;
+ lines.SetOwner(kTRUE);
+
+#ifdef WITH_ROOT
+ for ( Size_t i = 0; i < fMotifPositions.GetEntriesFast(); ++i )
+#else
+ for ( Size_t i = 0; i < fMotifPositions.size(); ++i )
+#endif
+ {
+ AliMpMotifPosition* pos = GetMotifPosition(i);
+ AliMpVMotif* motif = pos->GetMotif();
+ TVector2 lowerLeft(pos->Position()-pos->Dimensions());
+ TString id(motif->GetID());
+ // id is supposed to be of the form %s-%e-%e, and we're only
+ // interested in the %s part of it
+ Ssiz_t index = id.Index("-");
+ if ( index < 1 )
+ {
+ AliError(Form("id=%s does not meet expectations",id.Data()));
+ return;
+ }
+ TString motifName(id(0,index));
+ lines.Add(new TObjString(Form("MOTIF %s %d %d",
+ motifName.Data(),
+ TMath::Nint(lowerLeft.X()/fPadSizeX),
+ TMath::Nint(lowerLeft.Y()/fPadSizeY))));
+ }
+
+ ofstream out(fileName.Data());
+ out.precision(9);
+ out << "SIZES " << fPadSizeX << " " << fPadSizeY
+ << " " << fEnveloppeSizeX << " " << fEnveloppeSizeY
+ << endl;
+
+ TIter next(&lines);
+ TObjString* s;
+ while ( ( s = (TObjString*)next() ) )
+ {
+ out << s->String().Data() << endl;
+ }
+ out.close();
+}
+
//_____________________________________________________________________________
Double_t
AliMpPCB::X() const
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff