]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - FMD/AliFMD.cxx
Fixed assignment operator
[u/mrichter/AliRoot.git] / FMD / AliFMD.cxx
index fd95c4dc377c832cde5bf630bed368734cd7816a..9be365dde77b8975cefeb6987ccdb802e30fbe2d 100644 (file)
@@ -84,7 +84,7 @@
 // These files are not in the same directory, so there's no reason to
 // ask the preprocessor to search in the current directory for these
 // files by including them with `#include "..."' 
-#include <math.h>               // __CMATH__
+#include <cmath>                // __CMATH__
 #include <TClonesArray.h>      // ROOT_TClonesArray
 #include <TGeometry.h>         // ROOT_TGeomtry
 #include <TNode.h>             // ROOT_TNode
 #include <TTUBE.h>             // ROOT_TTUBE
 #include <TTree.h>             // ROOT_TTree
 #include <TBrowser.h>          // ROOT_TBrowser
-#include <TMath.h>             // ROOT_TMath
-#include <TVirtualMC.h>                // ROOT_TVirtualMC
-#include <TVector2.h>
-#include <TVector3.h>
-#include <TMarker3DBox.h>
+// #include <TVirtualMC.h>     // ROOT_TVirtualMC
+#include <TVector2.h>           // ROOT_TVector2 
+#include <TGeoManager.h>        // ROOT_TGeoManager
 
 #include <AliRunDigitizer.h>   // ALIRUNDIGITIZER_H
 #include <AliLoader.h>         // ALILOADER_H
 #include <AliRun.h>            // ALIRUN_H
 #include <AliMC.h>             // ALIMC_H
-#include "AliMagF.h"           // ALIMAGF_H
+#include <AliMagF.h>           // ALIMAGF_H
 #include <AliLog.h>            // ALILOG_H
 #include "AliFMD.h"            // ALIFMD_H
-#include "AliFMDDigit.h"       // ALIFMDDIGIG_H
+#include "AliFMDDigit.h"       // ALIFMDDIGIT_H
+#include "AliFMDSDigit.h"      // ALIFMDSDIGIT_H
 #include "AliFMDHit.h"         // ALIFMDHIT_H
 #include "AliFMDGeometry.h"    // ALIFMDGEOMETRY_H
 #include "AliFMDDetector.h"    // ALIFMDDETECTOR_H
 #include "AliFMDRing.h"                // ALIFMDRING_H
 #include "AliFMDDigitizer.h"   // ALIFMDDIGITIZER_H
-#include "AliPoints.h"         // ALIPOINTS_H
-#include "AliFMDGeometryBuilder.h"
+#include "AliFMDSDigitizer.h"  // ALIFMDSDIGITIZER_H
+// #include "AliFMDGeometryBuilder.h"
 #include "AliFMDRawWriter.h"   // ALIFMDRAWWRITER_H
-
-
-class AliFMDPoints : public AliPoints 
-{
-public:
-  /** Constructor 
-      @param hit   Hit to draw
-      @param color Color of hit */
-  AliFMDPoints(AliFMDHit* hit, UInt_t color) 
-    : AliPoints(1), fMarker(0)
-  {
-    if (!hit) return;
-    Float_t  size  = TMath::Min(TMath::Max(hit->Edep() * .1, .1), 1.);
-    TVector3 p(hit->Px(), hit->Py(), hit->Pz());
-    fMarker = new TMarker3DBox(hit->X(), hit->Y(), hit->Z(), size, size, size,
-                              p.Theta(), p.Phi());
-    fMarker->SetLineColor(color);
-    fMarker->SetRefObject(this);
-    fP[0] = hit->X();
-    fP[1] = hit->Y();
-    fP[2] = hit->Z();
-  }
-  /** Destructor  */
-  ~AliFMDPoints() 
-  {
-    // if (fMarker) delete  fMarker;
-  }
-  void SetXYZ(Double_t x, Double_t y, Double_t z)
-  {
-    if (fMarker) fMarker->SetPosition(x, y, z);
-  }
-  Int_t DistancetoPrimitive(Int_t px, Int_t py) 
-  {
-    return fMarker->DistancetoPrimitive(px, py);
-  }
-  void Draw(Option_t* option) 
-  {
-    if (fMarker) fMarker->Draw(option);
-  }
-  void Paint(Option_t* option)
-  {
-    if (fMarker) fMarker->Paint(option);
-  }
-  void SetMarkerColor(Color_t colour) 
-  {
-    if (fMarker) fMarker->SetLineColor(colour);
-  }
-private:
-  TMarker3DBox* fMarker;
-};
-
+#include "AliFMDPoints.h"       // ALIFMDPOINTS_H
 
 //____________________________________________________________________
 ClassImp(AliFMD)
@@ -179,7 +128,9 @@ AliFMD::AliFMD()
     fSDigits(0), 
     fNsdigits(0),
     fDetailed(kTRUE),
-    fBad(0)
+    fUseOld(kFALSE),
+    fUseAssembly(kTRUE),
+    fBad(0) 
 {
   //
   // Default constructor for class AliFMD
@@ -188,38 +139,23 @@ AliFMD::AliFMD()
   fHits        = 0;
   fDigits      = 0;
   fIshunt      = 0;
-  fUseOld      = kFALSE;
-  fUseAssembly = kTRUE;
   fBad         = new TClonesArray("AliFMDHit");
 }
 
-//____________________________________________________________________
-AliFMD::AliFMD(const AliFMD& other)
-  : AliDetector(other),
-    fSDigits(other.fSDigits), 
-    fNsdigits(other.fNsdigits),
-    fDetailed(other.fDetailed),
-    fBad(other.fBad)
-{
-  // Copy constructor 
-  fUseOld      = other.fUseOld;
-  fUseAssembly = other.fUseAssembly;
-}
-
 //____________________________________________________________________
 AliFMD::AliFMD(const char *name, const char *title)
   : AliDetector (name, title),
     fSDigits(0),
     fNsdigits(0),
     fDetailed(kTRUE),
+    fUseOld(kFALSE),
+    fUseAssembly(kFALSE),
     fBad(0)
 {
   //
   // Standard constructor for Forward Multiplicity Detector
   //
   AliDebug(10, "\tStandard CTOR");
-  fUseOld      = kFALSE;
-  fUseAssembly = kFALSE;
   fBad         = new TClonesArray("AliFMDHit");
   
   // Initialise Hit array
@@ -232,8 +168,8 @@ AliFMD::AliFMD(const char *name, const char *title)
   
   // CHC: What is this?
   fIshunt = 0;
-  SetMarkerColor(kRed);
-  SetLineColor(kYellow);
+  //PH  SetMarkerColor(kRed);
+  //PH  SetLineColor(kYellow);
 }
 
 //____________________________________________________________________
@@ -262,18 +198,6 @@ AliFMD::~AliFMD ()
   }
 }
 
-//____________________________________________________________________
-AliFMD&
-AliFMD::operator=(const AliFMD& other)
-{
-  // Assignment operator
-  AliDetector::operator=(other);
-  fSDigits             = other.fSDigits; 
-  fNsdigits            = other.fNsdigits;
-  fDetailed            = other.fDetailed;
-  fBad                  = other.fBad;
-  return *this;
-}
 
 //====================================================================
 //
@@ -550,18 +474,18 @@ AliFMD::BuildGeometry()
        return;
       }
       Double_t    siThick  = r->GetSiThickness();
-      const Int_t nv       = r->GetNVerticies();
+      const Int_t knv      = r->GetNVerticies();
       Double_t    theta    = r->GetTheta();
       Int_t       nmod     = r->GetNModules();
       
-      TXTRU* shape = new TXTRU(r->GetName(), r->GetTitle(), "void", nv, 2);
-      for (Int_t j = 0; j < nv; j++) {
-       TVector2* vv = r->GetVertex(nv - 1 - j);
+      TXTRU* shape = new TXTRU(r->GetName(), r->GetTitle(), "void", knv, 2);
+      for (Int_t j = 0; j < knv; j++) {
+       TVector2* vv = r->GetVertex(knv - 1 - j);
        shape->DefineVertex(j, vv->X(), vv->Y());
       }
       shape->DefineSection(0, -siThick / 2, 1, 0, 0);
       shape->DefineSection(1, +siThick / 2, 1, 0, 0);
-      shape->SetLineColor(GetLineColor());
+      shape->SetLineColor(kYellow); //PH kYellow is the default line color in FMD
       
       TObjArray* rots = new TObjArray(nmod);
       for (Int_t j = 0; j < nmod; j++) {
@@ -636,7 +560,7 @@ AliFMD::BuildGeometry()
        TNode* mnod = new TNode(name.Data(), title.Data(), rshape, 
                                0, 0, rz - siThick / 2 
                                + TMath::Sign(offz,z), rot);
-       mnod->SetLineColor(GetLineColor());
+       mnod->SetLineColor(kYellow); //PH kYellow is the default line color in FMD
        fNodes->Add(mnod);
       } // for (Int_t k = 0 ; ...)
     } // for (Int_t j = 0 ; ...)
@@ -686,7 +610,7 @@ AliFMD::LoadPoints(Int_t /* track */)
     // than 1 keV - a MIP is 100 eV. 
     if (edep > absQ * absQ && poverm > 1) bad = kTRUE;
 
-    AliFMDPoints* p1 = new AliFMDPoints(hit, GetMarkerColor());
+    AliFMDPoints* p1 = new AliFMDPoints(hit, kRed); //PH kRed is the default marker color in FMD
     // AliPoints* p1 = new AliPoints();
     // p1->SetMarkerColor(GetMarkerColor());
     // p1->SetMarkerSize(GetMarkerSize());
@@ -731,7 +655,7 @@ AliFMD::DrawDetector()
 
 //____________________________________________________________________
 Int_t 
-AliFMD::DistanceToPrimitive(Int_t, Int_t)
+AliFMD::DistancetoPrimitive(Int_t, Int_t)
 {
   // Calculate the distance from the mouse to the FMD on the screen
   // Dummy routine.
@@ -1150,6 +1074,52 @@ AliFMD::Browse(TBrowser* b)
   b->Add(AliFMDGeometry::Instance());
 }
 
+//____________________________________________________________________ 
+void
+AliFMD::AddAlignableVolumes() const
+{
+  //
+  // Create entries for alignable volumes associating the symbolic volume
+  // name with the corresponding volume path. Needs to be syncronized with
+  // eventual changes in the geometry.
+  // 
+  // This code was made by Raffaele Grosso <rgrosso@mail.cern.ch>.  I
+  // (cholm) will probably want to change it.   For one, I think it
+  // should be the job of the geometry manager to deal with this. 
+  AliFMDGeometry::Instance()->SetAlignableVolumes();
+#if 0  
+  for(size_t f = 1; f <= 3; f++){ // Detector 1,2,3
+    for(size_t tb =  0; tb <2 ; tb++){ // Top/Bottom 
+      char     stb = tb == 0 ? 'T' : 'B';
+      unsigned min = tb == 0 ? 0   : 5;
+
+      TString halfVol(Form("/ALIC_1/F%dM%c_%d", f, stb, f));
+      TString halfSym(halfVol);
+      if(!gGeoManager->SetAlignableEntry(halfSym.Data(),halfVol.Data()))
+       AliFatal(Form("Alignable entry %s not created. "
+                     "Volume path %s not valid", 
+                     halfSym.Data(),halfVol.Data()));
+      for(size_t io = 0; io < 2; io++){ // inner, outer 
+       if (f==1 && io==1) continue; // Only one ring in FMD1 
+       if(tb == 1 && io==1) min=10;
+       char     sio = (io == 0 ? 'I' : 'O');
+       unsigned nio = (io == 0 ? 3   : 9);
+       unsigned max = (io == 0 ? 5   : 10) + min;
+       
+       for(size_t i = min; i < max; i++) { // Modules
+         TString modVol(Form("%s/F%c%cV_7%d/F%cSE_%d", halfVol.Data(), 
+                             sio, stb, nio, sio, i));
+         TString modSym(modVol);
+         if(!gGeoManager->SetAlignableEntry(modSym.Data(),modVol.Data()))
+           AliFatal(Form("Alignable entry %s not created. "
+                         "Volume path %s not valid", 
+                         modSym.Data(), modVol.Data()));
+       }
+      }
+    }
+  }
+#endif
+}
 //___________________________________________________________________
 //
 // EOF