fChargeDispersion(25.e-4),
fSingleStepForChargeDispersion(0),
fNStepForChargeDispersion(4),
- fDensitySupportOverSi(0.15),
+ fDensitySupportOverSi(0.05),
fFileNameForUnderyingEvent(0),
fFileNameForPileUpEvents(0),
fNPileUpEvents(0),
fChargeDispersion(25.e-4),
fSingleStepForChargeDispersion(0),
fNStepForChargeDispersion(4),
- fDensitySupportOverSi(0.15),
+ fDensitySupportOverSi(0.05),
fFileNameForUnderyingEvent(0),
fFileNameForPileUpEvents(0),
fNPileUpEvents(0),
fChargeDispersion(25.e-4),
fSingleStepForChargeDispersion(0),
fNStepForChargeDispersion(4),
- fDensitySupportOverSi(0.15),
+ fDensitySupportOverSi(0.05),
fFileNameForUnderyingEvent(0),
fFileNameForPileUpEvents(0),
fNPileUpEvents(0),
AliDebug(2, Form("Entering StepManager: TVirtualMC::GetMC()->CurrentVolName() = %s", TVirtualMC::GetMC()->CurrentVolName()));
- if (!fSegmentation) AliFatal("No segmentation available"); // DO WE HAVE A SEGMENTATION???
+ if (!fSegmentation) AliFatal("No segmentation available");
if (!(this->IsActive())) return;
if (!(TVirtualMC::GetMC()->TrackCharge())) return;
TString planeNumber = TVirtualMC::GetMC()->CurrentVolName();
TString detElemNumber = TVirtualMC::GetMC()->CurrentVolName();
- if (planeNumber.Contains("support")) return;
- if (planeNumber.Contains("readout")) return;
+ if (!(planeNumber.Contains("active"))) return;
planeNumber.Remove(0,9);
- detElemNumber.Remove(0,19);
+ detElemNumber.Remove(0,18);
planeNumber.Remove(2);
detElemNumber.Remove(3);
Int_t detElemID = fSegmentation->GetDetElemGlobalID(planeNumber.Atoi(), detElemNumber.Atoi());
TGeoMedium *silicon = gGeoManager->GetMedium("MFT_Si");
TGeoMedium *readout = gGeoManager->GetMedium("MFT_Readout");
TGeoMedium *support = gGeoManager->GetMedium("MFT_Support");
- TGeoMedium *carbon = gGeoManager->GetMedium("MFT_Carbon");
- TGeoMedium *alu = gGeoManager->GetMedium("MFT_Alu");
- TGeoMedium *water = gGeoManager->GetMedium("MFT_Water");
- TGeoMedium *si02 = gGeoManager->GetMedium("MFT_SiO2");
- TGeoMedium *inox = gGeoManager->GetMedium("MFT_Inox");
+ TGeoMedium *carbon = gGeoManager->GetMedium("MFT_Carbon");
+ TGeoMedium *alu = gGeoManager->GetMedium("MFT_Alu");
+ TGeoMedium *water = gGeoManager->GetMedium("MFT_Water");
+ TGeoMedium *si02 = gGeoManager->GetMedium("MFT_SiO2");
+ TGeoMedium *inox = gGeoManager->GetMedium("MFT_Inox");
// ---- Cage & Services Description --------------------------------------------
// R. Tieulent - 17/01/2014 - Basic description for ITS/TPC matching studies
void SetMFTGeometry() {
- const Int_t nPlanes = 6;
+ const Int_t nPlanes = 5;
- const Float_t zCenter[nPlanes] = { -45.5, -49.3, -53.1, -56.9, -68.8, -76.8 }; // expressed in cm
+ const Float_t zCenter[nPlanes] = { -46.0, -49.3, -53.1, -68.7, -76.8 }; // expressed in cm
- const Float_t rMin[nPlanes] = { 2.50, 2.50, 2.50, 3.00, 3.50, 3.90 }; // expressed in cm
- const Float_t rMax[nPlanes] = { 10.80, 11.50, 12.30, 13.00, 15.30, 16.90 }; // expressed in cm
+ const Float_t rMin[nPlanes] = { 2.30, 2.30, 2.30, 3.30, 3.60 }; // expressed in cm
+ const Float_t rMax[nPlanes] = { 10.00, 10.50, 11.10, 13.90, 15.30 }; // expressed in cm
- const Float_t pixelSizeX[nPlanes] = { 20.e-4, 20.e-4, 20.e-4, 20.e-4, 20.e-4, 20.e-4 }; // expressed in cm
- const Float_t pixelSizeY[nPlanes] = { 20.e-4, 20.e-4, 20.e-4, 20.e-4, 20.e-4, 20.e-4 }; // expressed in cm
+ const Float_t pixelSizeX[nPlanes] = { 28.e-4, 28.e-4, 28.e-4, 28.e-4, 28.e-4 }; // expressed in cm
+ const Float_t pixelSizeY[nPlanes] = { 28.e-4, 28.e-4, 28.e-4, 28.e-4, 28.e-4 }; // expressed in cm
- const Float_t thicknessActive[nPlanes] = { 50.e-4, 50.e-4, 50.e-4, 50.e-4, 50.e-4, 50.e-4 }; // expressed in cm
- const Float_t thicknessSupport[nPlanes] = {2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4, 2000.e-4 }; // expressed in cm
- const Float_t thicknessReadout[nPlanes] = { 50.e-4, 50.e-4, 50.e-4, 50.e-4, 50.e-4, 50.e-4 }; // expressed in cm
+ const Float_t thicknessActive[nPlanes] = { 50.e-4, 50.e-4, 50.e-4, 50.e-4, 50.e-4 }; // expressed in cm
+ const Float_t thicknessSupport[nPlanes] = { 1. , 1. , 1. , 1. , 1. }; // expressed in cm
+ const Float_t thicknessReadout[nPlanes] = { 50.e-4, 50.e-4, 50.e-4, 50.e-4, 50.e-4 }; // expressed in cm
- const Float_t equivalentSilicon[nPlanes] = { 400.e-4, 400.e-4, 400.e-4, 400.e-4, 400.e-4, 400.e-4 }; // expressed in cm
- const Float_t equivalentSiliconBeforeFront[nPlanes] = { 0.e-4, 0.e-4, 0.e-4, 0.e-4, 0.e-4, 0.e-4 }; // expressed in cm
- const Float_t equivalentSiliconBeforeBack[nPlanes] = { 350.e-4, 350.e-4, 350.e-4, 350.e-4, 350.e-4, 350.e-4 }; // expressed in cm
+ const Float_t equivalentSilicon[nPlanes] = { 600.e-4, 600.e-4, 600.e-4, 600.e-4, 600.e-4 }; // expressed in cm
+ const Float_t equivalentSiliconBeforeFront[nPlanes] = { 0.e-4, 0.e-4, 0.e-4, 0.e-4, 0.e-4 }; // expressed in cm
+ const Float_t equivalentSiliconBeforeBack[nPlanes] = { 550.e-4, 550.e-4, 550.e-4, 550.e-4, 550.e-4 }; // expressed in cm
- const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {0., 0., 0., 0., 0., 0.};
+ const Float_t hasPixelRectangularPatternAlongY[nPlanes] = {0., 0., 0., 0., 0.};
TNtuple *geomMFT = new TNtuple("AliMFTGeometry", "ALICE MFT Geometry", "zCenter:rMin:rMax:pixelSizeX:pixelSizeY:thicknessActive:thicknessSupport:thicknessReadout:equivalentSilicon:equivalentSiliconBeforeFront:equivalentSiliconBeforeBack:hasPixelRectangularPatternAlongY");