New FMD code from Alla + code cleaning
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 22 Dec 2000 11:17:39 +0000 (11:17 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 22 Dec 2000 11:17:39 +0000 (11:17 +0000)
FMD/AliFMD.cxx
FMD/AliFMD.h
FMD/AliFMDhit.cxx [new file with mode: 0644]
FMD/AliFMDhit.h [new file with mode: 0644]
FMD/AliFMDv0.cxx
FMD/AliFMDv0.h
FMD/AliFMDv1.cxx
FMD/AliFMDv1.h
FMD/Makefile
FMD/ViewFMD.C

index e3cc9ff..321c68f 100644 (file)
  * about the suitability of this software for any purpose. It is          *
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
-
-/*
-$Log$
-Revision 1.8  2000/10/02 21:28:07  fca
-Removal of useless dependecies via forward declarations
-
-Revision 1.7  2000/01/19 17:16:47  fca
-Introducing a list of lists of hits -- more hits allowed for detector now
-
-Revision 1.6  1999/09/29 09:24:14  fca
-Introduction of the Copyright and cvs Log
-
-*/
-
-///////////////////////////////////////////////////////////////////////////////
+ //////////////////////////////////////////////////////////////////////////////
 //                                                                           //
-//  Forward Multiplicity Detector                                            //
+//  Forward Multiplicity Detector based on Silicon plates                    //
 //  This class contains the base procedures for the Forward Multiplicity     //
 //  detector                                                                 //
+//  Detector consists of 6 Si volumes covered pseudorapidity interval         //
+//  from 1.6 to 6.0.                                                         //
 //                                                                           //
 //Begin_Html
 /*
-<img src="picts/AliFMDClass.gif">
+<img src="gif/AliFMDClass.gif">
 </pre>
 <br clear=left>
 <font size=+2 color=red>
 <p>The responsible person for this module is
-<a href="mailto:Valeri.Kondratiev@cern.ch">Valeri Kondratiev</a>.
+<a href="mailto:Alla.Maevskaia@cern.ch">Alla Maevskaia</a>.
 </font>
 <pre>
 */
@@ -51,12 +39,12 @@ Introduction of the Copyright and cvs Log
 #include <TTUBE.h>
 #include <TNode.h>
 #include <TGeometry.h>
-#include <TLorentzVector.h>
-
-#include "AliMC.h"
+#include <TTree.h>
 #include "AliRun.h"
+#include "AliMC.h"
 #include "AliFMD.h"
+#include "AliFMDhit.h"
+
 ClassImp(AliFMD)
  
 //_____________________________________________________________________________
@@ -78,80 +66,83 @@ AliFMD::AliFMD(const char *name, const char *title)
  
   //
   // Initialise Hit array
-  fHits   = new TClonesArray("AliFMDhit",  405);
-  gAlice->AddHitList(fHits);
+  fHits     = new TClonesArray("AliFMDhit", 1000);
   
   fIshunt     =  0;
+  fIdSens1    =  0;
+
+  SetMarkerColor(kRed);
 }
  
+
+AliFMD::~AliFMD()
+{
+  delete fHits;
+}
 //_____________________________________________________________________________
 void AliFMD::AddHit(Int_t track, Int_t *vol, Float_t *hits)
 {
   //
-  // Add a FMD hit
+  // Add a hit to the list
   //
   TClonesArray &lhits = *fHits;
   new(lhits[fNhits++]) AliFMDhit(fIshunt,track,vol,hits);
 }
 //_____________________________________________________________________________
 void AliFMD::BuildGeometry()
 {
   //
   // Build simple ROOT TNode geometry for event display
   //
-  TNode *Node, *Top;
+  TNode *node, *top;
   const int kColorFMD  = 7;
   //
-  Top=gAlice->GetGeometry()->GetNode("alice");
+  top=gAlice->GetGeometry()->GetNode("alice");
 
   // FMD define the different volumes
+  new TRotMatrix("rot901","rot901",  90, 0, 90, 90, 180, 0);
+
+  new TTUBE("S_FMD0","FMD  volume 0","void",4.73,17.7,1.5);
+  top->cd();
+  node = new TNode("FMD0","FMD0","S_FMD0",0,0,64,"");
+  node->SetLineColor(kColorFMD);
+  fNodes->Add(node);
+
+  new TTUBE("S_FMD1","FMD  volume 1","void",23.4,36.,1.5);
+  top->cd();
+  node = new TNode("FMD1","FMD1","S_FMD1",0,0,85,"");
+  node->SetLineColor(kColorFMD);
+  fNodes->Add(node);
   
-  new TTUBE("S_FMD1","FMD sensitive volume 1","void",4.5,10.5,3);
-  Top->cd();
-  Node = new TNode("FMD1","FMD1","S_FMD1",0,0,86.5,"");
-  Node->SetLineColor(kColorFMD);
-  fNodes->Add(Node);
-  
-  new TTUBE("S_FMD2","FMD sensitive volume 2","void",4.5,10.5,3);
-  Top->cd();
-  Node = new TNode("FMD2","FMD2","S_FMD2",0,0,-86.5,"");
-  Node->SetLineColor(kColorFMD);
-  fNodes->Add(Node);
-  
-  new TTUBE("S_FMD3","FMD sensitive volume 3","void",8,14,3);
-  Top->cd();
-  Node = new TNode("FMD3","FMD3","S_FMD3",0,0,71.2,"");
-  Node->SetLineColor(kColorFMD);
-  fNodes->Add(Node);
-  
-  new TTUBE("S_FMD4","FMD sensitive volume 4","void",8,14,3);
-  Top->cd();
-  Node = new TNode("FMD4","FMD4","S_FMD4",0,0,-71.2,"");
-  fNodes->Add(Node);
-  Node->SetLineColor(kColorFMD);
-  
-  new TTUBE("S_FMD5","FMD sensitive volume 5","void",8,17.5,3);
-  Top->cd();
-  Node = new TNode("FMD5","FMD5","S_FMD5",0,0,44,"");
-  Node->SetLineColor(kColorFMD);
-  fNodes->Add(Node);
+  new TTUBE("S_FMD2","FMD  volume 2","void",4.73,17.7,1.5);
+  top->cd();
+  node = new TNode("FMD2","FMD2","S_FMD2",0,0,-64,"");
+  node->SetLineColor(kColorFMD);
+  fNodes->Add(node);
+
+  new TTUBE("S_FMD3","FMD  volume 3","void",23.4,36.,1.5);
+  top->cd();
+  node = new TNode("FMD3","FMD3","S_FMD3",0,0,-85,"");
+  node->SetLineColor(kColorFMD);
+  fNodes->Add(node);
+
+  new TTUBE("S_FMD4","FMD  volume 4","void",5,15,0.015);
+  top->cd();
+  node = new TNode("FMD4","FMD4","S_FMD4",0,0,-270,"");
+  node->SetLineColor(kColorFMD);
+  fNodes->Add(node);
   
-  new TTUBE("S_FMD6","FMD sensitive volume 6","void",8,17.5,3);
-  Top->cd();
-  Node = new TNode("FMD6","FMD6","S_FMD6",0,0,-44,"");
-  Node->SetLineColor(kColorFMD);
-  fNodes->Add(Node);
   
-  new TTUBE("S_FMD7","FMD sensitive volume 7","void",4.2,13,3);
-  Top->cd();
-  Node = new TNode("FMD7","FMD7","S_FMD7",0,0,-231,"");
-  Node->SetLineColor(kColorFMD);
-  fNodes->Add(Node);
+  new TTUBE("S_FMD5","FMD  volume 5","void",5,14,0.015);
+  top->cd();
+  node = new TNode("FMD5","FMD5","S_FMD5",0,0,-630,"");
+  node->SetLineColor(kColorFMD);
+  fNodes->Add(node);
+
 }
  
 //_____________________________________________________________________________
-Int_t AliFMD::DistancetoPrimitive(Int_t , Int_t )
+Int_t AliFMD::DistanceToPrimitive(Int_t px, Int_t py)
 {
   //
   // Calculate the distance from the mouse to the FMD on the screen
@@ -161,30 +152,14 @@ Int_t AliFMD::DistancetoPrimitive(Int_t , Int_t )
 }
  
 //_____________________________________________________________________________
-void AliFMD::StepManager()
-{
-  //
-  // Called for each step in the FMD
-  //
-  
-  Float_t       hits[3];
-  Int_t         i,copy,vol[1];
-  TClonesArray &lhits = *fHits;
-  TLorentzVector p;
-  
-  gMC->CurrentVolID(copy);
-  vol[0] = copy;
-  gMC->TrackPosition(p);
-  for(i=0;i<3;++i) hits[i]=p[i];
-  new(lhits[fNhits++]) AliFMDhit(fIshunt,gAlice->CurrentTrack(),vol,hits);
-}
 
-//___________________________________________
+//-------------------------------------------------------------------------
 void AliFMD::Init()
 {
   //
   // Initialis the FMD after it has been built
   Int_t i;
+  AliMC* pMC = AliMC::GetMC();
   //
   printf("\n");
   for(i=0;i<35;i++) printf("*");
@@ -195,21 +170,41 @@ void AliFMD::Init()
   // Here the FMD initialisation code (if any!)
   for(i=0;i<80;i++) printf("*");
   printf("\n");
+ //
+ //
+  fIdSens1=pMC->VolId("GFSI"); //Si sensetive volume
+
 }
+//---------------------------------------------------------------------
+void AliFMD::MakeBranch(Option_t* option)
+{
+  // Create Tree branches for the FMD.
+  Int_t buffersize = 4000;
+  char branchname[10];
+  sprintf(branchname,"%s",GetName());
 
+  AliDetector::MakeBranch(option);
+
+  if (fDigits   && gAlice->TreeD()) {
+    gAlice->TreeD()->Branch(branchname,&fDigits, buffersize);
+    printf("Making Branch %s for digits\n",branchname);
+  }
+}
  
-ClassImp(AliFMDhit)
-//_____________________________________________________________________________
-AliFMDhit::AliFMDhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits):
-  AliHit(shunt, track)
+//---------------------------------------------------------------------
+
+void AliFMD::Eta2Radius(Float_t eta, Float_t zDisk, Float_t *radius)
 {
-  //
-  // Add a FMD hit
-  //
-  fVolume = vol[0];
-  fX=hits[0];
-  fY=hits[1];
-  fZ=hits[2];
+   Float_t expEta=TMath::Exp(-eta);
+   Float_t theta=TMath::ATan(expEta);
+   theta=2.*theta;
+   Float_t rad=zDisk*(TMath::Tan(theta));
+   *radius=rad;
+   
+   printf(" eta %f radius %f\n", eta, rad);
 }
  
+
+    
index 2b01c43..2895593 100644 (file)
@@ -1,49 +1,34 @@
-#ifndef FMD_H
-#define FMD_H
+#ifndef ALIFMD_H
+#define ALIFMD_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
-/* $Id$ */
-
 ////////////////////////////////////////////////
-//  Manager and hits classes for set:FMD     //
+//  Manager and hits classes for set:Si-FMD     //
 ////////////////////////////////////////////////
  
 #include "AliDetector.h"
-#include "AliHit.h"
  
 class AliFMD : public AliDetector {
  
 public:
   AliFMD();
   AliFMD(const char *name, const char *title);
-  virtual       ~AliFMD() {}
+  virtual       ~AliFMD(); 
   virtual void   AddHit(Int_t, Int_t*, Float_t*);
   virtual void   BuildGeometry();
-  virtual void   CreateGeometry() {}
-  virtual void   CreateMaterials() {}
-  Int_t          DistancetoPrimitive(Int_t, Int_t);
+  virtual void   CreateGeometry()=0;
+  virtual void   CreateMaterials()=0; 
+  virtual Int_t  DistanceToPrimitive(Int_t px, Int_t py);
   virtual Int_t  IsVersion() const =0;
   virtual void   Init();
-  virtual void   DrawModule()=0;
-  virtual void   StepManager();
+  virtual void   MakeBranch(Option_t *opt=" ");
+  virtual void   DrawDetector()=0;
+  virtual void   StepManager()=0;
+  void  Eta2Radius(Float_t, Float_t, Float_t*);
   
+ protected:
+  Int_t fIdSens1;     //Si sensetive volume
   ClassDef(AliFMD,1)  //Class for the FMD detector
 };
-
-//_____________________________________________________________________________
-class AliFMDhit : public AliHit {
-public:
-  Int_t      fVolume;  //Volume copy identifier
-  
-public:
-  AliFMDhit() {}
-  AliFMDhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits);
-  virtual ~AliFMDhit() {}
-  
-  ClassDef(AliFMDhit,1)  //Hits for detector FMD
-};
-
 #endif
diff --git a/FMD/AliFMDhit.cxx b/FMD/AliFMDhit.cxx
new file mode 100644 (file)
index 0000000..0381da4
--- /dev/null
@@ -0,0 +1,47 @@
+/**************************************************************************
+ * 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.                  *
+ **************************************************************************/
+#include "AliFMDhit.h"
+
+ClassImp(AliFMDhit)
+
+AliFMDhit::AliFMDhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits):
+  AliHit(shunt, track)
+{
+//Normal FMD  hit ctor
+  fVolume = vol[0];
+  fX=hits[0];
+  fY=hits[1];
+  fZ=hits[2];
+  fPx=hits[3];
+  fPy=hits[4];
+  fPz=hits[5];
+  fEdep=hits[6];
+  fParticle=(Int_t)hits[7];
+  fTime=hits[8];
+}
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/FMD/AliFMDhit.h b/FMD/AliFMDhit.h
new file mode 100644 (file)
index 0000000..938eef1
--- /dev/null
@@ -0,0 +1,91 @@
+#ifndef AliFMDhit_H
+#define AliFMDhit_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+////////////////////////////////////////////////
+//  Manager and hits classes for set:FMD     //
+////////////////////////////////////////////////
+#include "AliHit.h"
+class AliFMDhit : public AliHit {
+///////////////////////////////////////////////////////////////////////
+// AliFMDhit is the hit class for the FMD. Hits are the information
+// that comes from a Monte Carlo at each step as a particle mass through
+// sensitive detector elements as particles are transported through a
+// detector.
+//
+// Data members:
+//
+// Int_t fTrack
+//     See AliHit for a full description. The track number of the track
+// that made this hit.
+//
+// Float_t fX
+//     See AliHit for a full description. The global x position of the
+// hit (in the standard units of the Monte Carlo).
+//
+// Float_t fY
+//     See AliHit for a full description. The global y position of the
+// hit (in the standard units of the Monte Carlo).
+//
+// Float_t fZ
+//     See AliHit for a full description. The global z position of the
+// hit (in the standard units of the Monte Carlo).
+//
+// Int_t fStatus
+//     The track status flag. This flag indicates the track status
+// at the time of creating this hit. It is made up of the following 8
+// status bits from highest order to lowest order bits
+// 0           :  IsTrackAlive():    IsTrackStop():IsTrackDisappeared():
+// IsTrackOut():IsTrackExiting():IsTrackEntering():IsTrackInside()     .
+// See AliMC for a description of these functions. If the function is
+// true then the bit is set to one, otherwise it is zero.
+
+// Int_t fVolume
+//     The number of the FMD detector that contains this hit. 
+
+// Float_t fEdep
+//     The energy lost by the particle during the step ending in this
+// hit. The units are those determined by the Monte Carlo.
+//
+// Float_t fPx
+//     The x momentum, in global coordinates, of the particle that
+// "created" the hit at the time and position of the hit. The units
+// are those determined by the Monte Carlo.
+//
+// Float_t fPy
+//     The y momentum, in global coordinates, of the particle that
+// "created" the hit at the time and position of the hit. The units
+// are those determined by the Monte Carlo.
+//
+// Float_t fPz
+//     The z momentum, in global coordinates, of the particle that
+// "created" the hit at the time and position of the hit. The units
+// are those determined by the Monte Carlo.
+//
+///
+// Float_t fTime
+//     The time of flight associated with the particle ending in this
+// hit. The time is typically measured from the point of creation of the
+// original particle (if this particle is a daughter).  The units
+// are those determined by the Monte Carlo.
+
+private:
+  Int_t      fVolume;       //Volume copy identifier
+  Int_t      fParticle;     //Particle identificator
+  Float_t    fEdep;         //Energy deposition
+  Float_t    fPx;            // Particle's momentum X
+  Float_t    fPy;            // Particle's momentum Y
+  Float_t    fPz;            // Particle's momentum Z
+  Float_t    fTime;         // Particle's time of flight
+  
+public:
+  AliFMDhit() {}
+  AliFMDhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits);
+  virtual ~AliFMDhit() {}
+  
+  ClassDef(AliFMDhit,1)  //Hits for detector FMD
+};
+
+#endif
index f4580c2..f5fab4c 100644 (file)
  * about the suitability of this software for any purpose. It is          *
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
-
-/*
-$Log$
-Revision 1.7  1999/09/29 09:24:14  fca
-Introduction of the Copyright and cvs Log
-
-*/
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//  Forward Multiplicity Detector version 0                                  //
-//                                                                           //
-//Begin_Html
+ /////////////////////////////////////////////////////////////////////
+//                                                                 //
+// Forward Multiplicity detector based on Silicon version 0        //
+//
+//Begin Html       
 /*
-<img src="picts/AliFMDv0Class.gif">
-</pre>
-<br clear=left>
-<font size=+2 color=red>
-<p>The responsible person for this module is
-<a href="mailto:Valeri.Kondratiev@cern.ch">Valeri Kondratiev</a>.
-</font>
-<pre>
+<img src="gif/AliFMDv0Class.gif">
 */
-//End_Html
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
+//End Html
+//                                                                  //
+//                                                                  //
+//////////////////////////////////////////////////////////////////////
 
-#include "AliRun.h"
+#include <TMath.h>
+#include <TGeometry.h>
+#include <TTUBE.h>
+#include <TNode.h>
+#include <TLorentzVector.h>
 #include "AliFMDv0.h"
+#include "AliRun.h"
 #include "AliMC.h"
+#include <iostream.h>
+#include <fstream.h>
+
 #include "AliMagF.h"
+#include "AliFMDhit.h"
+#include <stdlib.h>
 
-#include "stdlib.h"
+//class TGeant3;
 ClassImp(AliFMDv0)
-//_____________________________________________________________________________
-AliFMDv0::AliFMDv0()
-{
-  //
-  // Default constructor for FMD version 0
-  //
-}
-//_____________________________________________________________________________
-AliFMDv0::AliFMDv0(const char *name, const char *title)
-  : AliFMD(name,title)
+
+//--------------------------------------------------------------------
+AliFMDv0::AliFMDv0(const char *name, const char *title):
+ AliFMD(name,title)
 {
   //
-  // Standard constructor for FMD version 0
+  // Standart constructor for Forward Multiplicity Detector version 0
   //
-  AliModule *start = gAlice->GetModule("START");
-  if(start) {
-    Error("ctor","This version of FMD is incompatible with START\n");
-    exit(1);
-  }
+  fIdSens1=0;
+//  setBufferSize(128000);
 }
-//___________________________________________
+//-------------------------------------------------------------------------
 void AliFMDv0::CreateGeometry()
 {
+ //
+  // Create the geometry of Forward Multiplicity Detector version 0
   //
-  // Creation of the geometry of the FMD version 0
-  //
-  //Begin_Html
-  /*
-    <img src="picts/AliFMDv0Tree.gif">
-  */
-  //End_Html
-  //Begin_Html
+  //Detector consists of 6 volumes: 
+  // 1st covered pseudorapidity interval from 3.3 to 2.0
+  // and placed on 65cm in Z-direction;
+  // 2nd - from 2.0 to 1.6 and Z=85 cm;
+  // 3d  - the same pseudorapidity interval as the 1st 
+  // but on the other side from the interaction point z=-65cm;
+  // 4th - simmetricaly with the 2nd : 
+  // pseudorapidity from 2.0 to 1.6, Z=-85cm   
+  // 5th - from 3.6 to 4.7, Z=-270cm
+  // 6th - from 4.5 to 5.5 , Z=-630cm.
+  // Each part has 400mkm Si (sensetive area, detector itself),
+  // 0.75cm of plastic simulated electronics material,
+  // Al support ring 2cm thickness and 1cm width placed on 
+  // the outer radius of each Si disk;
+  //    
+  // begin Html
   /*
-    <img src="picts/AliFMDv0.gif">
+   <img src="gif/AliFMDv0.gif">
   */
-  //End_Html
+  //
+
+  Int_t *idtmed = fIdtmed->GetArray();
+   
+  Int_t ifmd;
+  Int_t idrotm[999];
+  Float_t zfmd,par[3];
+  char name[5];
+
+  Float_t rin[6], rout[6],zpos;
+  Float_t etain[6]= {3.3, 2.0, 3.3, 2.0, 4.7, 5.5};
+  Float_t etaout[6]={2.0, 1.6, 2.0, 1.6, 3.6, 4.5};
+  Float_t z[6]={64., 85., -64., -85., -270., -630};
+  Float_t zDet=0.04;
+  Float_t zElectronic=0.75;
+  Float_t zSupport=1.;
+  Float_t zFMD=3.;
+//-------------------------------------------------------------------
+ //  FMD 
+ //------------------------------------------------------------------
+
+  AliMatrix(idrotm[901], 90, 0, 90, 90, 180, 0);
+
+  gMC->Gsvolu("GFSI","TUBE", idtmed[1], par, 0);
+  gMC->Gsvolu("GEL ","TUBE", idtmed[4], par, 0);
+  gMC->Gsvolu("GSUP","TUBE", idtmed[2], par, 0);
+
+  for (ifmd =0; ifmd < 6; ifmd++){
+
+    sprintf(name,"FMD%d",ifmd);
+    printf(name);
+    
+    zfmd=TMath::Abs(z[ifmd]);
+    printf("zfmd %f z[ifmd] %f",zfmd,z[ifmd]);
+    AliFMD::Eta2Radius(etain[ifmd],zfmd,&rin[ifmd]);
+    AliFMD::Eta2Radius(etaout[ifmd],zfmd,&rout[ifmd]);
+    
+    par[0]=rin[ifmd]; // pipe size
+    par[1]=rout[ifmd];
+    par[2]=zFMD/2;
+    gMC->Gsvolu(name,"TUBE", idtmed[3], par, 3);
+    
+    printf ("rin %f rout %f ZFMD %f\n",par[0],par[1],z[ifmd]);
+    if (z[ifmd] < 0){  
+      gMC->Gspos(name,1,"ALIC",0,0,z[ifmd],0, "ONLY");}
+    else { 
+      gMC->Gspos(name,1,"ALIC",0,0,z[ifmd],idrotm[901], "ONLY");}
+  //Silicon detector
+    par[2]=zDet/2;
+    zpos=zFMD/2 -par[2];
+    gMC->Gsposp("GFSI",ifmd+1,name,0,0,zpos,0, "ONLY",par,3);
+
+    //Plastic slice for electronics
+    par[2]=zElectronic/2;
+    zpos=zpos-zDet/2-par[2];
+    gMC->Gsposp("GEL ",ifmd+1,name,0,0,zpos,0, "ONLY",par,3);
+
+   //Simple Al support
+   par[1]=rout[ifmd];
+   par[0]=rout[ifmd]-2;
+   par[2]=zSupport/2;
+   zpos=zpos-zElectronic/2-par[2];
+   gMC->Gsposp("GSUP",ifmd+1,name,0,0,zpos,0, "ONLY",par,3);
+   
+
+  }  
+  //Silicon
 
   
-  Float_t rout, z;
-  Float_t par[3], rin;
-  
-  Int_t *idtmed = fIdtmed->GetArray()-899;    
-  
-  // ******************************************************** 
-  //       DEFINE RIGHT DISK#3  OF FMD 
-  // ******************************************************** 
-  
-  //       Define parameters 
-  
-  rin  = 4.5;
-  rout = 10.5;
-  z    = 85.;
-  //       Ring #1 
-  par[0] = rin;
-  par[1] = rout;
-  par[2] = 1.5;
-  gMC->Gsvolu("R1R3", "TUBE", idtmed[900], par, 3);
-  gMC->Gspos("R1R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
-  //       Ring #2 
-  par[0] = rout;
-  par[1] = rout + .65;
-  par[2] = 1.5;
-  gMC->Gsvolu("R2R3", "TUBE", idtmed[901], par, 3);
-  gMC->Gspos("R2R3", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
-  //       Ring #3 
-  par[0] = rout + .65;
-  par[1] = rout + 5.65;
-  par[2] = .025;
-  gMC->Gsvolu("R3R3", "TUBE", idtmed[902], par, 3);
-  gMC->Gspos("R3R3", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
-  // *********************************************************** 
-  //       DEFINE LEFT DISK#3 OF FMD 
-  // *********************************************************** 
-  
-  //       Define parameters 
-  
-  rin  = 4.5;
-  rout = 10.5;
-  z    = -85.;
-  //       Ring #1 
-  par[0] = rin;
-  par[1] = rout;
-  par[2] = 1.5;
-  gMC->Gsvolu("R1L3", "TUBE", idtmed[900], par, 3);
-  gMC->Gspos("R1L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #2 
-  par[0] = rout;
-  par[1] = rout + .65;
-  par[2] = 1.5;
-  gMC->Gsvolu("R2L3", "TUBE", idtmed[901], par, 3);
-  gMC->Gspos("R2L3", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #3 
-  par[0] = rout + .65;
-  par[1] = rout + 5.65;
-  par[2] = .025;
-  gMC->Gsvolu("R3L3", "TUBE", idtmed[902], par, 3);
-  gMC->Gspos("R3L3", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
-  // ******************************************************** 
-  //       DEFINE RIGHT DISK#2  OF FMD 
-  // ******************************************************** 
-  
-  //       Define parameters 
-  
-  rin  = 8.;
-  rout = 14.;
-  z    = 69.7;
-  //       Ring #1 
-  par[0] = rin;
-  par[1] = rout;
-  par[2] = 1.5;
-  gMC->Gsvolu("R1R2", "TUBE", idtmed[900], par, 3);
-  gMC->Gspos("R1R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
-  //       Ring #2 
-  par[0] = rout;
-  par[1] = rout + .65;
-  par[2] = 1.5;
-  gMC->Gsvolu("R2R2", "TUBE", idtmed[901], par, 3);
-  gMC->Gspos("R2R2", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
-  //       Ring #3 
-  par[0] = rout + .65;
-  par[1] = rout + 5.65;
-  par[2] = .025;
-  gMC->Gsvolu("R3R2", "TUBE", idtmed[902], par, 3);
-  gMC->Gspos("R3R2", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
-  // *********************************************************** 
-  //       DEFINE LEFT DISK#2 OF FMD 
-  // *********************************************************** 
-  
-  //       Define parameters 
-  
-  rin  = 8.;
-  rout = 14.;
-  z    = -69.7;
-  //       Ring #1 
-  par[0] = rin;
-  par[1] = rout;
-  par[2] = 1.5;
-  gMC->Gsvolu("R1L2", "TUBE", idtmed[900], par, 3);
-  gMC->Gspos("R1L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #2 
-  par[0] = rout;
-  par[1] = rout + .65;
-  par[2] = 1.5;
-  gMC->Gsvolu("R2L2", "TUBE", idtmed[901], par, 3);
-  gMC->Gspos("R2L2", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #3 
-  par[0] = rout + .65;
-  par[1] = rout + 5.65;
-  par[2] = .025;
-  gMC->Gsvolu("R3L2", "TUBE", idtmed[902], par, 3);
-  gMC->Gspos("R3L2", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
-  // ******************************************************** 
-  //       DEFINE RIGHT DISK#1  OF FMD 
-  // ******************************************************** 
-  
-  //       Define parameters 
-  
-  rin  = 8.;
-  rout = 17.5;
-  z    = 42.5;
-  //       Ring #1 
-  par[0] = rin;
-  par[1] = rout;
-  par[2] = 1.5;
-  gMC->Gsvolu("R1R1", "TUBE", idtmed[900], par, 3);
-  gMC->Gspos("R1R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
-  //       Ring #2 
-  par[0] = rout;
-  par[1] = rout + .65;
-  par[2] = 1.5;
-  gMC->Gsvolu("R2R1", "TUBE", idtmed[901], par, 3);
-  gMC->Gspos("R2R1", 1, "ALIC", 0., 0., z + 1.5, 0, "ONLY");
-  //       Ring #3 
-  par[0] = rout + .65;
-  par[1] = rout + 5.65;
-  par[2] = .025;
-  gMC->Gsvolu("R3R1", "TUBE", idtmed[902], par, 3);
-  gMC->Gspos("R3R1", 1, "ALIC", 0., 0., z + 1.525, 0, "ONLY");
-  // *********************************************************** 
-  //       DEFINE LEFT DISK#1 OF FMD 
-  // *********************************************************** 
-  
-  //       Define parameters 
-  
-  rin  = 8.;
-  rout = 17.5;
-  z    = -42.5;
-  //       Ring #1 
-  par[0] = rin;
-  par[1] = rout;
-  par[2] = 1.5;
-  gMC->Gsvolu("R1L1", "TUBE", idtmed[900], par, 3);
-  gMC->Gspos("R1L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #2 
-  par[0] = rout;
-  par[1] = rout + .65;
-  par[2] = 1.5;
-  gMC->Gsvolu("R2L1", "TUBE", idtmed[901], par, 3);
-  gMC->Gspos("R2L1", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #3 
-  par[0] = rout + .65;
-  par[1] = rout + 5.65;
-  par[2] = .025;
-  gMC->Gsvolu("R3L1", "TUBE", idtmed[902], par, 3);
-  gMC->Gspos("R3L1", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
-  // *********************************************************** 
-  //       DEFINE LEFT DISK#4 OF FMD 
-  // *********************************************************** 
+    
 
-  //       Define parameters 
+}    
+//------------------------------------------------------------------------
+void AliFMDv0::CreateMaterials()
+{
+ Int_t isxfld   = gAlice->Field()->Integ();
+ Float_t sxmgmx = gAlice->Field()->Max();
 
-  rin  = 4.2;
-  rout = 13.;
-  z    = -229.5;
-  //       Ring #1 
-  par[0] = rin;
-  par[1] = rout;
-  par[2] = 1.5;
-  gMC->Gsvolu("R1L4", "TUBE", idtmed[900], par, 3);
-  gMC->Gspos("R1L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #2 
-  par[0] = rout;
-  par[1] = rout + .65;
-  par[2] = 1.5;
-  gMC->Gsvolu("R2L4", "TUBE", idtmed[901], par, 3);
-  gMC->Gspos("R2L4", 1, "ALIC", 0., 0., z - 1.5, 0, "ONLY");
-  //       Ring #3 
-  par[0] = rout + .65;
-  par[1] = rout + 5.65;
-  par[2] = .025;
-  gMC->Gsvolu("R3L4", "TUBE", idtmed[902], par, 3);
-  gMC->Gspos("R3L4", 1, "ALIC", 0., 0., z - 1.525, 0, "ONLY");
-}
+ // Plastic CH
+ Float_t aPlastic[2]={1.01,12.01};
+ Float_t zPlastic[2]={1,6};
+ Float_t wPlastic[2]={1,1};
+ Float_t denPlastic=1.03;
+   //
+ //*** Definition Of avaible FMD materials ***
+ AliMaterial(0, "Si chip$", 28.0855,14.,2.33,9.36,999);
+ AliMaterial(1, "Al supprt$", 26.980,13.,2.70,8.9,999);
+ AliMaterial(2, "FMD Air$", 14.61, 7.3, .001205, 30423.,999); 
+ AliMixture( 5, "Plastic$",aPlastic,zPlastic,denPlastic,-2,wPlastic);
 
-//_____________________________________________________________________________
-void AliFMDv0::DrawModule()
+//**
+ AliMedium(1, "Si chip$", 0, 1, isxfld, sxmgmx, 1., .001, 1., .001, .001);
+ AliMedium(2, "Al support$", 1, 0, isxfld, sxmgmx, 1., .001, 1., .001, .001);
+ AliMedium(3, "FMD air$", 2, 0, isxfld, sxmgmx, 1., .001, 1., .001, .001);
+ AliMedium(4, "Plastic$", 5, 0,isxfld, sxmgmx,  10., .01, 1., .003, .003);
+
+
+}
+//---------------------------------------------------------------------
+void AliFMDv0::DrawDetector()
 {
-  //
-  // Draw a shaded view of the FMD version 0
-  //
+//
+// Draw a shaded view of the Forward multiplicity detector version 0
+//
 
-  // Set everything unseen
-  gMC->Gsatt("*", "seen", -1);
-  // 
-  // Set ALIC mother transparent
-  gMC->Gsatt("ALIC","SEEN",0);
-  //
-  // Set the volumes visible
-  gMC->Gsatt("R1R3","SEEN",1);
-  gMC->Gsatt("R2R3","SEEN",1);
-  gMC->Gsatt("R3R3","SEEN",1);
-  gMC->Gsatt("R1L3","SEEN",1);
-  gMC->Gsatt("R2L3","SEEN",1);
-  gMC->Gsatt("R3L3","SEEN",1);
-  gMC->Gsatt("R1R2","SEEN",1);
-  gMC->Gsatt("R2R2","SEEN",1);
-  gMC->Gsatt("R3R2","SEEN",1);
-  gMC->Gsatt("R1L2","SEEN",1);
-  gMC->Gsatt("R2L2","SEEN",1);
-  gMC->Gsatt("R3L2","SEEN",1);
-  gMC->Gsatt("R1R1","SEEN",1);
-  gMC->Gsatt("R2R1","SEEN",1);
-  gMC->Gsatt("R3R1","SEEN",1);
-  gMC->Gsatt("R1L1","SEEN",1);
-  gMC->Gsatt("R2L1","SEEN",1);
-  gMC->Gsatt("R3L1","SEEN",1);
-  gMC->Gsatt("R1L4","SEEN",1);
-  gMC->Gsatt("R2L4","SEEN",1);
-  gMC->Gsatt("R3L4","SEEN",1);
-  //
-  gMC->Gdopt("hide", "on");
-  gMC->Gdopt("shad", "on");
-  gMC->Gsatt("*", "fill", 7);
-  gMC->SetClipBox(".");
-  gMC->SetClipBox("*", 0, 1000, -1000, 1000, -1000, 1000);
-  gMC->DefaultRange();
-  gMC->Gdraw("alic", 40, 30, 0, 6, 9, .08, .08);
-  gMC->Gdhead(1111, "Forward Multiplicity Detector version 0");
-  gMC->Gdman(13, 9, "MAN");
+AliMC* pMC = AliMC::GetMC();
+
+//Set ALIC mother transparent
+pMC->Gsatt("ALIC","SEEN",0);
+//
+//Set volumes visible
+gMC->Gsatt("FMD0","SEEN",1);
+gMC->Gsatt("FMD1","SEEN",1);
+gMC->Gsatt("FMD2","SEEN",1);
+gMC->Gsatt("FMD3","SEEN",1);
+gMC->Gsatt("FMD4","SEEN",1);
+gMC->Gsatt("FMD5","SEEN",1);
+
+//
+gMC->Gdopt("hide","on");
+gMC->Gdopt("shad","on");
+gMC->SetClipBox(".");
+gMC->SetClipBox("*",0,1000,-1000,1000,-1000,1000);
+gMC->DefaultRange();
+gMC->Gdraw("alic",40,30,0,12,9.5,.2,0.2);
+gMC->Gdhead(1111,"Forward multiplicity detector");
+gMC->Gdopt("hide","off");
+}
+//-------------------------------------------------------------------
+void AliFMDv0::Init()
+{
+// Initialises version 0 of the Forward Multiplicity Detector
+//
+AliMC* gMC=AliMC::GetMC();
+AliFMD::Init();
+fIdSens1=gMC->VolId("GFSI");
+printf("*** FMD version 0 initialized ***\n");
 }
 
-//_____________________________________________________________________________
-void AliFMDv0::CreateMaterials()
+//-------------------------------------------------------------------
+
+void AliFMDv0::StepManager()
 {
   //
-  // Create Materials for version 0 of FMD
+  // Called for every step in the Forward Multiplicity Detector
   //
+  Int_t id,copy;
+  static Float_t hits[9];
+  static Int_t vol[1];
+  static Float_t de;
+  TLorentzVector pos;
+  TLorentzVector mom;
 
-  //     Material for ring #1 
-  Float_t ar1[8] = { 55.8,58.7,52.,47.9,16.,28.,207.2,27. };
-  Float_t zr1[8] = { 26.,28.,24.,22.,8.,14.,82.,13. };
-  Float_t wr1[8] = { .27,.081,.054,.045,.18,.25,.06,.06 };
-  //     Material for ring #2 
-  Float_t ar2[3] = { 55.8,27.,16. };
-  Float_t zr2[3] = { 26.,13.,8. };
-  Float_t wr2[3] = { .35,.34,.31 };
-  //     Material for ring #3 
-  Float_t ar3[3] = { 28.,27.,16. };
-  Float_t zr3[3] = { 14.,13.,8. };
-  Float_t wr3[3] = { .37,.33,.3 };
-  
-  Float_t epsil, stmin, deemax, tmaxfd, stemax;
-  
-  Int_t   ISXFLD = gAlice->Field()->Integ();
-  Float_t SXMGMX = gAlice->Field()->Max();
-  
-  //     Ring #1 
-  
-  AliMixture(1, "FMD_R1$", ar1, zr1, 2.69, 8, wr1);
-  
-  //     Ring #2 
-  
-  AliMixture(2, "FMD_R2$", ar2, zr2, 2.63, 3, wr2);
-  
-  //     Ring #3 
+
+  TClonesArray &lhits = *fHits;
+  AliMC* gMC=AliMC::GetMC();
+  if(!gMC->IsTrackAlive()) return; // particle has disappeared
+
+  Float_t charge = gMC->TrackCharge();
+  if(TMath::Abs(charge)<=0.) return; //take only charged particles
+
+  //  ((TGeant3*)gMC)->Gpcxyz();
+
+  id=gMC->CurrentVolID(copy);
   
-  AliMixture(3, "FMD_R3$", ar3, zr3, 3.15, 3, wr3);
-  // ******************************************************* 
-  //     Defines tracking media parameters. 
-  // ******************************************************* 
-  epsil  = .001; // Tracking precision, DLS 
-  stemax = -1.;  // Maximum displacement for multiple scattering 
-  tmaxfd = -20.; // Maximum angle due to field deflection 
-  deemax = -.3;  // Maximum fractional energy loss, DLS 
-  stmin  = -.8;
-  // ******************************************************** 
-  AliMedium(1, "FMD_R1_L3        ", 1, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(2, "FMD_R2_L3        ", 2, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(3, "FMD_R3_L3        ", 3, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-}
+// Check the sensetive volume
+   if(id==fIdSens1)
+     {
+       if(gMC->IsTrackEntering())
+        {
+          vol[0]=copy;
+          gMC->TrackPosition(pos);
+          hits[0]=pos[0];
+          hits[1]=pos[1];
+          hits[2]=pos[2];
+          
+          gMC->TrackMomentum(mom);
+          hits[3]=mom[0];
+          hits[4]=mom[1];
+          hits[5]=mom[2];
 
+          Int_t iPart= gMC->TrackPid();
+          Int_t partId=gMC->IdFromPDG(iPart);
+          hits[7]=partId;
+          hits[8]=1e9*gMC->TrackTime();
+          de=0.;
+        }
+       if(gMC->IsTrackInside()){
+          de=de+1000.*gMC->Edep();
+       }
+       
+       if(gMC->IsTrackExiting()
+         ||gMC->IsTrackDisappeared()||
+         gMC->IsTrackStop())
+        {
+            hits[6]=de+1000.*gMC->Edep();
+    new(lhits[fNhits++]) AliFMDhit(fIshunt,gAlice->CurrentTrack(),vol,hits);
+        } // IsTrackExiting()
+     }
+ }
 
index abdad2b..24f9b53 100644 (file)
@@ -1,11 +1,9 @@
-#ifndef FMDV0_H
-#define FMDV0_H
+#ifndef ALIFMDV0_H
+#define ALIFMDV0_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
-/* $Id$ */
-
-////////////////////////////////////////////////
+///////////////////////////////////////////////
 //  Manager and hits classes for set:FMD     //
 ////////////////////////////////////////////////
  
 class AliFMDv0 : public AliFMD {
   
 public:
-  AliFMDv0();
+  AliFMDv0() {};
   AliFMDv0(const char *name, const char *title);
   virtual       ~AliFMDv0() {}
   virtual void   CreateGeometry();
   virtual void   CreateMaterials();
-  virtual void   DrawModule();
+  virtual void   DrawDetector();
+  virtual void   Init();
   virtual Int_t  IsVersion() const {return 0;}
+  virtual void   StepManager();
+  
+protected:
+   Int_t fIdSens1; // Sensetive volume  in FMD
  
   ClassDef(AliFMDv0,1)  //Class for FMD version 0
 };
 
 #endif
+
+
index 1990b26..598fd69 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.10  2000/10/02 21:28:07  fca
+Removal of useless dependecies via forward declarations
+
 Revision 1.9  2000/05/10 21:56:07  fca
 Avoid clashes with ITS and add supports
 
@@ -820,8 +823,8 @@ void AliFMDv1::CreateMaterials()
   
   Float_t epsil, stmin, deemax, tmaxfd, stemax;
   
-  Int_t   ISXFLD = gAlice->Field()->Integ();
-  Float_t SXMGMX = gAlice->Field()->Max();
+  Int_t   isxfld = gAlice->Field()->Integ();
+  Float_t sxmgmx = gAlice->Field()->Max();
   
   //     Ring #1 
   
@@ -852,11 +855,11 @@ void AliFMDv1::CreateMaterials()
   deemax = -.3;  // Maximum fractional energy loss, DLS 
   stmin  = -.8;
   // ******************************************************** 
-  AliMedium(1, "FMD_R1_L3        ", 1, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(2, "FMD_R2_L3        ", 2, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(3, "FMD_R3_L3        ", 3, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(4, "FMD_BR_L3        ", 4, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(5, "FMD_CB_L3        ", 5, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(1, "FMD_R1_L3        ", 1, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(2, "FMD_R2_L3        ", 2, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(3, "FMD_R3_L3        ", 3, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(4, "FMD_BR_L3        ", 4, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(5, "FMD_CB_L3        ", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 }
 
 
index ebf3c82..c4b152d 100644 (file)
@@ -1,5 +1,5 @@
-#ifndef FMDV1_H
-#define FMDV1_H
+#ifndef ALIFMDV1_H
+#define ALIFMDV1_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
index 8be75d8..b99e269 100644 (file)
@@ -9,7 +9,7 @@ PACKAGE = FMD
 
 # C++ sources
 
-SRCS          = AliFMD.cxx AliFMDv0.cxx AliFMDv1.cxx 
+SRCS          = AliFMD.cxx AliFMDv0.cxx AliFMDv1.cxx AliFMDhit.cxx
 
 # C++ Headers
 
index 984359c..306545c 100644 (file)
@@ -1,94 +1,7 @@
 void ViewFMD()
 {
-   gMC->Gsatt("IWR3","seen",1);
-   gMC->Gsatt("FWR3","seen",1);
-   gMC->Gsatt("RWR3","seen",1);
-   gMC->Gsatt("MPR3","seen",1);
-   gMC->Gsatt("SPR3","seen",1);
-   gMC->Gsatt("SMR3","seen",1);
-   gMC->Gsatt("CPR3","seen",1);
-   gMC->Gsatt("C1R3","seen",1);
-   gMC->Gsatt("C2R3","seen",1);
-   gMC->Gsatt("GKR3","seen",1);
-   gMC->Gsatt("SCR3","seen",1);
-   gMC->Gsatt("SER3","seen",1);
-   gMC->Gsatt("CER3","seen",1);
-   gMC->Gsatt("IWL3","seen",1);
-   gMC->Gsatt("FWL3","seen",1);
-   gMC->Gsatt("RWL3","seen",1);
-   gMC->Gsatt("MPL3","seen",1);
-   gMC->Gsatt("SPL3","seen",1);
-   gMC->Gsatt("SML3","seen",1);
-   gMC->Gsatt("CPL3","seen",1);
-   gMC->Gsatt("C1L3","seen",1);
-   gMC->Gsatt("C2L3","seen",1);
-   gMC->Gsatt("GKL3","seen",1);
-   gMC->Gsatt("SCL3","seen",1);
-   gMC->Gsatt("SEL3","seen",1);
-   gMC->Gsatt("CEL3","seen",1);
-   gMC->Gsatt("IWR2","seen",1);
-   gMC->Gsatt("FWR2","seen",1);
-   gMC->Gsatt("RWR2","seen",1);
-   gMC->Gsatt("MPR2","seen",1);
-   gMC->Gsatt("SPR2","seen",1);
-   gMC->Gsatt("SMR2","seen",1);
-   gMC->Gsatt("CPR2","seen",1);
-   gMC->Gsatt("C1R2","seen",1);
-   gMC->Gsatt("C2R2","seen",1);
-   gMC->Gsatt("GKR2","seen",1);
-   gMC->Gsatt("SCR2","seen",1);
-   gMC->Gsatt("SER2","seen",1);
-   gMC->Gsatt("CER2","seen",1);
-   gMC->Gsatt("IWL2","seen",1);
-   gMC->Gsatt("FWL2","seen",1);
-   gMC->Gsatt("RWL2","seen",1);
-   gMC->Gsatt("MPL2","seen",1);
-   gMC->Gsatt("SPL2","seen",1);
-   gMC->Gsatt("SML2","seen",1);
-   gMC->Gsatt("CPL2","seen",1);
-   gMC->Gsatt("C1L2","seen",1);
-   gMC->Gsatt("C2L2","seen",1);
-   gMC->Gsatt("GKL2","seen",1);
-   gMC->Gsatt("SCL2","seen",1);
-   gMC->Gsatt("SEL2","seen",1);
-   gMC->Gsatt("CEL2","seen",1);
-   gMC->Gsatt("IWR1","seen",1);
-   gMC->Gsatt("FWR1","seen",1);
-   gMC->Gsatt("RWR1","seen",1);
-   gMC->Gsatt("MPR1","seen",1);
-   gMC->Gsatt("SPR1","seen",1);
-   gMC->Gsatt("SMR1","seen",1);
-   gMC->Gsatt("CPR1","seen",1);
-   gMC->Gsatt("C1R1","seen",1);
-   gMC->Gsatt("C2R1","seen",1);
-   gMC->Gsatt("GKR1","seen",1);
-   gMC->Gsatt("SCR1","seen",1);
-   gMC->Gsatt("SER1","seen",1);
-   gMC->Gsatt("CER1","seen",1);
-   gMC->Gsatt("IWL1","seen",1);
-   gMC->Gsatt("FWL1","seen",1);
-   gMC->Gsatt("RWL1","seen",1);
-   gMC->Gsatt("MPL1","seen",1);
-   gMC->Gsatt("SPL1","seen",1);
-   gMC->Gsatt("SML1","seen",1);
-   gMC->Gsatt("CPL1","seen",1);
-   gMC->Gsatt("C1L1","seen",1);
-   gMC->Gsatt("C2L1","seen",1);
-   gMC->Gsatt("GKL1","seen",1);
-   gMC->Gsatt("SCL1","seen",1);
-   gMC->Gsatt("SEL1","seen",1);
-   gMC->Gsatt("CEL1","seen",1);
-   gMC->Gsatt("IWL4","seen",1);
-   gMC->Gsatt("FWL4","seen",1);
-   gMC->Gsatt("RWL4","seen",1);
-   gMC->Gsatt("MPL4","seen",1);
-   gMC->Gsatt("SPL4","seen",1);
-   gMC->Gsatt("SML4","seen",1);
-   gMC->Gsatt("CPL4","seen",1);
-   gMC->Gsatt("C1L4","seen",1);
-   gMC->Gsatt("C2L4","seen",1);
-   gMC->Gsatt("GKL4","seen",1);
-   gMC->Gsatt("SCL4","seen",1);
-   gMC->Gsatt("SEL4","seen",1);
-   gMC->Gsatt("CEL4","seen",1);
+   gMC->Gsatt("FMD0","seen",1);
+   gMC->Gsatt("FMD1","seen",1);
+   gMC->Gsatt("FMD2","seen",1);
+   gMC->Gsatt("FMD3","seen",1);
 }