1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 //____________________________________________________________________
20 // Forward Multiplicity Detector based on Silicon wafers. This class
21 // contains the base procedures for the Forward Multiplicity detector
22 // Detector consists of 5 Si volumes covered pseudorapidity interval
25 // This class contains the detailed version of the FMD - that is, hits
26 // are produced during simulation.
28 // The actual code is done by various separate classes. Below is
29 // diagram showing the relationship between the various FMD classes
30 // that handles the geometry
33 // +----------+ +----------+
34 // | AliFMDv1 | | AliFMDv1 |
35 // +----------+ +----------+
37 // +----+--------------+
39 // | +------------+ 1 +---------------+
40 // | +- | AliFMDRing |<>--| AliFMDPolygon |
41 // V 2 | +------------+ +---------------+
44 // +--------+<>--+ V 1..2
45 // 3 | +-------------------+
46 // +-| AliFMDSubDetector |
47 // +-------------------+
50 // +-------------+-------------+
52 // +---------+ +---------+ +---------+
53 // | AliFMD1 | | AliFMD2 | | AliFMD3 |
54 // +---------+ +---------+ +---------+
57 // See also the class AliFMD for a more detailed explanation of the
59 #include <TVirtualMC.h> // ROOT_TVirtualMC
60 #include <AliRun.h> // ALIRUN_H
61 #include <AliMC.h> // ALIMC_H
62 #include <AliLog.h> // ALILOG_H
63 #include "AliFMDv1.h" // ALIFMDV1_H
65 //____________________________________________________________________
69 //____________________________________________________________________
71 AliFMDv1::StepManager()
74 // Called for every step in the Forward Multiplicity Detector
76 // The procedure is as follows:
78 // - IF NOT track is alive THEN RETURN ENDIF
79 // - IF NOT particle is charged THEN RETURN ENDIF
80 // - IF NOT volume name is "STRI" or "STRO" THEN RETURN ENDIF
81 // - Get strip number (volume copy # minus 1)
82 // - Get phi division number (mother volume copy #)
83 // - Get module number (grand-mother volume copy #)
84 // - section # = 2 * module # + phi division # - 1
85 // - Get ring Id from volume name
86 // - Get detector # from grand-grand-grand-mother volume name
87 // - Get pointer to sub-detector object.
88 // - Get track position
89 // - IF track is entering volume AND track is inside real shape THEN
90 // - Reset energy deposited
91 // - Get track momentum
92 // - Get particle ID #
94 // - IF track is inside volume AND inside real shape THEN
95 /// - Update energy deposited
97 // - IF track is inside real shape AND (track is leaving volume,
98 // or it died, or it is stopped THEN
103 AliDebug(10, Form("Is inside %s", gMC->CurrentVolName()));
105 // If the track is gone, return
106 if (!gMC->IsTrackAlive()) return;
108 // Only process charged particles
109 if(TMath::Abs(gMC->TrackCharge()) <= 0) return;
111 // Only do stuff is the track is in one of the strips.
112 // TString vol(gMC->CurrentVolName());
113 // if (!vol.Contains("STR")) return;
115 Int_t volumeId = gMC->CurrentVolID(copy);
116 // The ring ID is encoded in the volume name
118 if (volumeId == fInner->GetStripId()) ring = 'I';
119 else if (volumeId == fOuter->GetStripId()) ring = 'O';
122 // Get the strip number. Note, that GEANT numbers divisions from 1,
123 // so we subtract one
124 Int_t strip = copy - 1;
126 // Get the phi division of the module
127 Int_t phiDiv; // * The phi division number (1 or 2)
128 gMC->CurrentVolOffID(1, phiDiv); // in the module
130 // Active volume number - not used.
132 // gMC->CurrentVolOffID(2, active);
134 // Get the module number in the ring.
136 gMC->CurrentVolOffID(3, module);
138 // Ring copy number - the same as the detector number - not used
139 // Int_t ringCopy; // * Ring copy number
140 // gMC->CurrentVolOffID(4, ringCopy); // Same as detector number
142 // Get the detector number from the path name
143 Int_t detector = Int_t((gMC->CurrentVolOffName(5)[3]) - 48);
145 // The sector number, calculated from module and phi division #
146 Int_t sector = 2 * module + phiDiv - 1;
149 // Get a pointer to the sub detector structure
150 AliFMDSubDetector* det = 0;
152 case 1: det = fFMD1; break;
153 case 2: det = fFMD2; break;
154 case 3: det = fFMD3; break;
158 // Get the current track position
160 gMC->TrackPosition(v);
161 // Check that the track is actually within the active area
162 Bool_t isWithin = det->CheckHit(ring, module, v.X(), v.Y());
163 Bool_t entering = gMC->IsTrackEntering() && isWithin;
164 Bool_t inside = gMC->IsTrackInside() && isWithin;
165 Bool_t out = (gMC->IsTrackExiting()
166 || gMC->IsTrackDisappeared()
167 || gMC->IsTrackStop()
170 AliDebug(2, Form("Is inside FMD%d%c[%02d,%03d]: particle is %s",
171 detector, ring, sector, strip,
172 (entering ? "entering" :
176 // Reset the energy deposition for this track, and update some of
181 // Get production vertex and momentum of the track
183 gMC->TrackMomentum(fCurrentP);
184 fCurrentPdg = gMC->IdFromPDG(gMC->TrackPid());
187 // fAnalyser->Update(detector, ring, isWithin, v.X(), v.Y());
190 // If the track is inside, then update the energy deposition
191 if (inside && fCurrentDeltaE >= 0)
192 fCurrentDeltaE += 1000 * gMC->Edep();
194 // The track exits the volume, or it disappeared in the volume, or
195 // the track is stopped because it no longer fulfills the cuts
196 // defined, then we create a hit.
197 if (out && fCurrentDeltaE >= 0) {
198 fCurrentDeltaE += 1000 * gMC->Edep();
200 AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(),
201 detector, ring, sector, strip,
202 fCurrentV.X(), fCurrentV.Y(), fCurrentV.Z(),
203 fCurrentP.X(), fCurrentP.Y(), fCurrentP.Z(),
204 fCurrentDeltaE, fCurrentPdg, fCurrentV.T());
208 //___________________________________________________________________