// Copyright 2019-2020 CERN and copyright holders of ALICE O2.

// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.

// All rights not expressly granted are reserved.

//

// This software is distributed under the terms of the GNU General Public

// License v3 (GPL Version 3), copied verbatim in the file "COPYING".

//

// In applying this license CERN does not waive the privileges and immunities

// granted to it by virtue of its status as an Intergovernmental Organization

// or submit itself to any jurisdiction.

#include "TGeoManager.h" // for TGeoManager

#include "TMath.h"

#include "TString.h"

#include <fairlogger/Logger.h>

#include "FairVolume.h"

#include "FairRootManager.h"

#include "TOFSimulation/Detector.h"

#include <TVirtualMC.h> // for TVirtualMC, gMC

#include "DetectorsBase/GeometryManager.h"

#include "DetectorsBase/Stack.h"

using namespace o2::tof;

ClassImp(Detector);

Detector::Detector(Bool_t active)

: o2::base::DetImpl<Detector>("TOF", active), mEventNr(0), mTOFHoles(kTRUE), mHits(o2::utils::createSimVector<HitType>())

{

for (Int_t i = 0; i < Geo::NSECTORS; i++) {

mTOFSectors[i] = 1;

}

}

Detector::Detector(const Detector& rhs)

: o2::base::DetImpl<Detector>(rhs),

mEventNr(0),

mTOFHoles(rhs.mTOFHoles),

mHits(o2::utils::createSimVector<HitType>())

{

for (Int_t i = 0; i < Geo::NSECTORS; i++) {

mTOFSectors[i] = rhs.mTOFSectors[i];

}

}

Detector::~Detector()

{

o2::utils::freeSimVector(mHits);

}

void Detector::InitializeO2Detector()

{

TGeoVolume* v = gGeoManager->GetVolume("FPAD");

if (v == nullptr) {

printf("Sensitive volume FSEN not found!!!!!!!!");

} else {

AddSensitiveVolume(v);

}

}

Bool_t Detector::ProcessHits(FairVolume* v)

{

// This method is called from the MC stepping for the sensitive volume only

if (static_cast<int>(fMC->TrackCharge()) == 0) {

// set a very large step size for neutral particles

return kFALSE; // take only charged particles

}

float pos2x, pos2y, pos2z;

fMC->TrackPosition(pos2x, pos2y, pos2z);

Float_t radius = std::sqrt(pos2x * pos2x + pos2y * pos2y);

LOG(debug) << "Process hit in TOF volume ar R=" << radius << " - Z=" << pos2z;

Float_t enDep = fMC->Edep();

if (enDep < 1E-8) {

return kFALSE; // wo se need a threshold?

}

// ADD HIT

float posx, posy, posz;

fMC->TrackPosition(posx, posy, posz);

float time = fMC->TrackTime() * 1.0e09;

auto stack = static_cast<o2::data::Stack*>(fMC->GetStack());

int trackID = stack->GetCurrentTrackNumber();

int sensID = v->getMCid();

Int_t det[5];

Float_t pos[3] = {posx, posy, posz};

Float_t delta[3];

Geo::getPadDxDyDz(pos, det, delta);

auto channel = Geo::getIndex(det);

HitType newhit(posx, posy, posz, time, enDep, trackID, sensID);

if (channel != mLastChannelID || !isMergable(newhit, mHits->back())) {

mHits->push_back(newhit);

stack->addHit(GetDetId());

} else {

mHits->back().SetEnergyLoss(mHits->back().GetEnergyLoss() + newhit.GetEnergyLoss());

// LOG(info)<<"Merging hit "<<"\n";

// <<mHits->back().GetId()<<"with new hit "<<newhit.GetId()<<"\n";

}

mLastChannelID = channel;

return kTRUE;

}

void Detector::Register()

{

FairRootManager::Instance()->RegisterAny(addNameTo("Hit").data(), mHits, kTRUE);

}

void Detector::Reset()

{

// TODO: move this out of here

if (!o2::utils::ShmManager::Instance().isOperational()) {

mHits->clear();

}

mLastChannelID = -1;

}

void Detector::CreateMaterials()

{

Int_t isxfld = 2;

Float_t sxmgmx = 10.;

o2::base::Detector::initFieldTrackingParams(isxfld, sxmgmx);

//--- Quartz (SiO2) ---

Float_t aq[2] = {28.0855, 15.9994};

Float_t zq[2] = {14., 8.};

Float_t wq[2] = {1., 2.};

Float_t dq = 2.7; // (+5.9%)

Int_t nq = -2;

// --- Nomex (C14H22O2N2) ---

Float_t anox[4] = {12.011, 1.00794, 15.9994, 14.00674};

Float_t znox[4] = {6., 1., 8., 7.};

Float_t wnox[4] = {14., 22., 2., 2.};

// Float_t dnox = 0.048; //old value

Float_t dnox = 0.22; // (x 4.6)

Int_t nnox = -4;

// --- G10 {Si, O, C, H, O} ---

Float_t we[7], na[7];

Float_t ag10[5] = {28.0855, 15.9994, 12.011, 1.00794, 15.9994};

Float_t zg10[5] = {14., 8., 6., 1., 8.};

Float_t wmatg10[5];

Int_t nlmatg10 = 5;

na[0] = 1., na[1] = 2., na[2] = 0., na[3] = 0., na[4] = 0.;

MaterialMixer(we, ag10, na, 5);

wmatg10[0] = we[0] * 0.6;

wmatg10[1] = we[1] * 0.6;

na[0] = 0., na[1] = 0., na[2] = 14., na[3] = 20., na[4] = 3.;

MaterialMixer(we, ag10, na, 5);

wmatg10[2] = we[2] * 0.4;

wmatg10[3] = we[3] * 0.4;

wmatg10[4] = we[4] * 0.4;

// Float_t densg10 = 1.7; //old value

Float_t densg10 = 2.0; // (+17.8%)

// --- Water ---

Float_t awa[2] = {1.00794, 15.9994};

Float_t zwa[2] = {1., 8.};

Float_t wwa[2] = {2., 1.};

Float_t dwa = 1.0;

Int_t nwa = -2;

// --- Air ---

Float_t aAir[4] = {12.011, 14.00674, 15.9994, 39.948};

Float_t zAir[4] = {6., 7., 8., 18.};

Float_t wAir[4] = {0.000124, 0.755267, 0.231781, 0.012827};

Float_t dAir = 1.20479E-3;

// --- Fibre Glass ---

Float_t afg[4] = {28.0855, 15.9994, 12.011, 1.00794};

Float_t zfg[4] = {14., 8., 6., 1.};

Float_t wfg[4] = {0.12906, 0.29405, 0.51502, 0.06187};

// Float_t dfg = 1.111;

Float_t dfg = 2.05; // (x1.845)

Int_t nfg = 4;

// --- Freon C2F4H2 + SF6 ---

Float_t afre[4] = {12.011, 1.00794, 18.9984032, 32.0065};

Float_t zfre[4] = {6., 1., 9., 16.};

Float_t wfre[4] = {0.21250, 0.01787, 0.74827, 0.021355};

Float_t densfre = 0.00375;

Int_t nfre = 4;

// --- Cables and tubes {Al, Cu} ---

Float_t acbt[2] = {26.981539, 63.546};

Float_t zcbt[2] = {13., 29.};

Float_t wcbt[2] = {0.407, 0.593};

Float_t decbt = 0.68;

// --- Cable {CH2, Al, Cu} ---

Float_t asc[4] = {12.011, 1.00794, 26.981539, 63.546};

Float_t zsc[4] = {6., 1., 13., 29.};

Float_t wsc[4];

for (Int_t ii = 0; ii < 4; ii++) {

wsc[ii] = 0.;

}

Float_t wDummy[4], nDummy[4];

for (Int_t ii = 0; ii < 4; ii++) {

wDummy[ii] = 0.;

}

for (Int_t ii = 0; ii < 4; ii++) {

nDummy[ii] = 0.;

}

nDummy[0] = 1.;

nDummy[1] = 2.;

MaterialMixer(wDummy, asc, nDummy, 2);

wsc[0] = 0.4375 * wDummy[0];

wsc[1] = 0.4375 * wDummy[1];

wsc[2] = 0.3244;

wsc[3] = 0.2381;

Float_t dsc = 1.223;

// --- Crates boxes {Al, Cu, Fe, Cr, Ni} ---

Float_t acra[5] = {26.981539, 63.546, 55.845, 51.9961, 58.6934};

Float_t zcra[5] = {13., 29., 26., 24., 28.};

Float_t wcra[5] = {0.7, 0.2, 0.07, 0.018, 0.012};

Float_t dcra = 0.77;

// --- Polietilene CH2 ---

Float_t aPlastic[2] = {12.011, 1.00794};

Float_t zPlastic[2] = {6., 1.};

Float_t wPlastic[2] = {1., 2.};

// Float_t dPlastic = 0.92; // PDB value

Float_t dPlastic = 0.93; // (~+1.1%)

Int_t nwPlastic = -2;

Mixture(0, "Air$", aAir, zAir, dAir, 4, wAir);

Mixture(1, "Nomex$", anox, znox, dnox, nnox, wnox);

Mixture(2, "G10$", ag10, zg10, densg10, nlmatg10, wmatg10);

Mixture(3, "fibre glass$", afg, zfg, dfg, nfg, wfg);

Material(4, "Al $", 26.981539, 13., 2.7, -8.9, 999.);

Float_t factor = 0.4 / 1.5 * 2. / 3.;

Material(5, "Al honeycomb$", 26.981539, 13., 2.7 * factor, -8.9 / factor, 999.);

Mixture(6, "Freon$", afre, zfre, densfre, nfre, wfre);

Mixture(7, "Glass$", aq, zq, dq, nq, wq);

Mixture(8, "Water$", awa, zwa, dwa, nwa, wwa);

Mixture(9, "cables+tubes$", acbt, zcbt, decbt, 2, wcbt);

Material(10, "Cu $", 63.546, 29., 8.96, -1.43, 999.);

Mixture(11, "cable$", asc, zsc, dsc, 4, wsc);

Mixture(12, "Al+Cu+steel$", acra, zcra, dcra, 5, wcra);

Mixture(13, "plastic$", aPlastic, zPlastic, dPlastic, nwPlastic, wPlastic);

Float_t factorHoles = 1. / 36.5;

Material(14, "Al honey for holes$", 26.981539, 13., 2.7 * factorHoles, -8.9 / factorHoles, 999.);

Float_t epsil, stmin, deemax, stemax;

// STD data

// EPSIL = 0.1 ! Tracking precision,

// STEMAX = 0.1 ! Maximum displacement for multiple scattering

// DEEMAX = 0.1 ! Maximum fractional energy loss, DLS

// STMIN = 0.1

// TOF data

epsil = .001; // Tracking precision,

stemax = -1.; // Maximum displacement for multiple scattering

deemax = -.3; // Maximum fractional energy loss, DLS

stmin = -.8;

Medium(kAir, "Air$", 0, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kNomex, "Nomex$", 1, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kG10, "G10$", 2, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kFiberGlass, "fibre glass$", 3, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kAlFrame, "Al Frame$", 4, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kHoneycomb, "honeycomb$", 5, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kFre, "Fre$", 6, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kCuS, "Cu-S$", 10, 1, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kGlass, "Glass$", 7, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kWater, "Water$", 8, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kCable, "Cable$", 11, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kCableTubes, "Cables+Tubes$", 9, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kCopper, "Copper$", 10, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kPlastic, "Plastic$", 13, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kCrates, "Crates$", 12, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

Medium(kHoneyHoles, "honey_holes$", 14, 0, isxfld, sxmgmx, 10., stemax, deemax, epsil, stmin);

}

void Detector::MaterialMixer(Float_t* p, const Float_t* const a, const Float_t* const m, Int_t n) const

{

// a[] atomic weights vector (in)

// (atoms present in more compound appear separately)

// m[] number of corresponding atoms in the compound (in)

Float_t t = 0.;

for (Int_t i = 0; i < n; ++i) {

p[i] = a[i] * m[i];

t += p[i];

}

for (Int_t i = 0; i < n; ++i) {

p[i] = p[i] / t;

}

}

void Detector::ConstructGeometry()

{

CreateMaterials();

/*

xTof = 124.5;//fTOFGeometry->StripLength()+2.*(0.3+0.03); // cm, x-dimension of FTOA volume

yTof = fTOFGeometry->Rmax()-fTOFGeometry->Rmin(); // cm, y-dimension of FTOA volume

Float_t zTof = fTOFGeometry->ZlenA(); // cm, z-dimension of FTOA volume

*/

Float_t xTof = Geo::STRIPLENGTH + 2.5, yTof = Geo::RMAX - Geo::RMIN, zTof = Geo::ZLENA;

DefineGeometry(xTof, yTof, zTof);

LOG(info) << "Loaded TOF geometry";

}

void Detector::EndOfEvent() { Reset(); }

void Detector::DefineGeometry(Float_t xtof, Float_t ytof, Float_t zlenA)

{

//

// Definition of the Time Of Fligh Resistive Plate Chambers

//

Float_t xFLT, yFLT, zFLTA;

xFLT = xtof - 2. * Geo::MODULEWALLTHICKNESS;

yFLT = ytof * 0.5 - Geo::MODULEWALLTHICKNESS;

zFLTA = zlenA - 2. * Geo::MODULEWALLTHICKNESS;

createModules(xtof, ytof, zlenA, xFLT, yFLT, zFLTA);

makeStripsInModules(ytof, zlenA);

createModuleCovers(xtof, zlenA);

createBackZone(xtof, ytof, zlenA);

makeFrontEndElectronics(xtof);

makeFEACooling(xtof);

makeNinoMask(xtof);

makeSuperModuleCooling(xtof, ytof, zlenA);

makeSuperModuleServices(xtof, ytof, zlenA);

makeModulesInBTOFvolumes(ytof, zlenA);

makeCoversInBTOFvolumes();

makeBackInBTOFvolumes(ytof);

makeReadoutCrates(ytof);

}

void Detector::createModules(Float_t xtof, Float_t ytof, Float_t zlenA, Float_t xFLT, Float_t yFLT, Float_t zFLTA) const

{

//

// Create supermodule volume

// and wall volumes to separate 5 modules

//

Int_t idrotm[8];

for (Int_t ii = 0; ii < 8; ii++) {

idrotm[ii] = 0;

}

// Definition of the of fibre glass modules (FTOA, FTOB and FTOC)

Float_t par[3];

par[0] = xtof * 0.5;

par[1] = ytof * 0.25;

par[2] = zlenA * 0.5;

TVirtualMC::GetMC()->Gsvolu("FTOA", "BOX ", getMediumID(kFiberGlass), par, 3); // Fibre glass

if (mTOFHoles) {

par[0] = xtof * 0.5;

par[1] = ytof * 0.25;

par[2] = (zlenA * 0.5 - Geo::INTERCENTRMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gsvolu("FTOB", "BOX ", getMediumID(kFiberGlass), par, 3); // Fibre glass

TVirtualMC::GetMC()->Gsvolu("FTOC", "BOX ", getMediumID(kFiberGlass), par, 3); // Fibre glass

}

// Definition and positioning

// of the not sensitive volumes with Insensitive Freon (FLTA, FLTB and FLTC)

par[0] = xFLT * 0.5;

par[1] = yFLT * 0.5;

par[2] = zFLTA * 0.5;

TVirtualMC::GetMC()->Gsvolu("FLTA", "BOX ", getMediumID(kFre), par, 3); // Freon mix

Float_t xcoor, ycoor, zcoor;

xcoor = 0.;

ycoor = Geo::MODULEWALLTHICKNESS * 0.5;

zcoor = 0.;

TVirtualMC::GetMC()->Gspos("FLTA", 0, "FTOA", xcoor, ycoor, zcoor, 0, "ONLY");

if (mTOFHoles) {

par[2] = (zlenA * 0.5 - 2. * Geo::MODULEWALLTHICKNESS - Geo::INTERCENTRMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gsvolu("FLTB", "BOX ", getMediumID(kFre), par, 3); // Freon mix

TVirtualMC::GetMC()->Gsvolu("FLTC", "BOX ", getMediumID(kFre), par, 3); // Freon mix

// xcoor = 0.;

// ycoor = Geo::MODULEWALLTHICKNESS*0.5;

zcoor = Geo::MODULEWALLTHICKNESS;

TVirtualMC::GetMC()->Gspos("FLTB", 0, "FTOB", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FLTC", 0, "FTOC", xcoor, ycoor, -zcoor, 0, "ONLY");

}

// Definition and positioning

// of the fibre glass walls between central and intermediate modules (FWZ1 and FWZ2)

Float_t alpha, tgal, beta, tgbe, trpa[11];

// tgal = (yFLT - 2.*Geo::LENGTHINCEMODBORDER)/(Geo::INTERCENTRMODBORDER2 - Geo::INTERCENTRMODBORDER1);

tgal = (yFLT - Geo::LENGTHINCEMODBORDERU - Geo::LENGTHINCEMODBORDERD) /

(Geo::INTERCENTRMODBORDER2 - Geo::INTERCENTRMODBORDER1);

alpha = TMath::ATan(tgal);

beta = (TMath::Pi() * 0.5 - alpha) * 0.5;

tgbe = TMath::Tan(beta);

trpa[0] = xFLT * 0.5;

trpa[1] = 0.;

trpa[2] = 0.;

trpa[3] = 2. * Geo::MODULEWALLTHICKNESS;

// trpa[4] = (Geo::LENGTHINCEMODBORDER - 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

// trpa[5] = (Geo::LENGTHINCEMODBORDER + 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

trpa[4] = (Geo::LENGTHINCEMODBORDERD - 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[5] = (Geo::LENGTHINCEMODBORDERD + 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[6] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

trpa[7] = 2. * Geo::MODULEWALLTHICKNESS;

trpa[8] = (Geo::LENGTHINCEMODBORDERD - 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[9] = (Geo::LENGTHINCEMODBORDERD + 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

// trpa[8] = (Geo::LENGTHINCEMODBORDER - 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

// trpa[9] = (Geo::LENGTHINCEMODBORDER + 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

trpa[10] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

TVirtualMC::GetMC()->Gsvolu("FWZ1D", "TRAP", getMediumID(kFiberGlass), trpa, 11); // Fibre glass

Matrix(idrotm[0], 90., 90., 180., 0., 90., 180.);

Matrix(idrotm[1], 90., 90., 0., 0., 90., 0.);

// xcoor = 0.;

// ycoor = -(yFLT - Geo::LENGTHINCEMODBORDER)*0.5;

ycoor = -(yFLT - Geo::LENGTHINCEMODBORDERD) * 0.5;

zcoor = Geo::INTERCENTRMODBORDER1;

TVirtualMC::GetMC()->Gspos("FWZ1D", 1, "FLTA", xcoor, ycoor, zcoor, idrotm[0], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ1D", 2, "FLTA", xcoor, ycoor, -zcoor, idrotm[1], "ONLY");

Float_t y0B, ycoorB, zcoorB;

if (mTOFHoles) {

// y0B = Geo::LENGTHINCEMODBORDER - Geo::MODULEWALLTHICKNESS*tgbe;

y0B = Geo::LENGTHINCEMODBORDERD - Geo::MODULEWALLTHICKNESS * tgbe;

trpa[0] = xFLT * 0.5;

trpa[1] = 0.;

trpa[2] = 0.;

trpa[3] = Geo::MODULEWALLTHICKNESS;

trpa[4] = (y0B - Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[5] = (y0B + Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[6] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

trpa[7] = Geo::MODULEWALLTHICKNESS;

trpa[8] = (y0B - Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[9] = (y0B + Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[10] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

// xcoor = 0.;

ycoorB = ycoor - Geo::MODULEWALLTHICKNESS * 0.5 * tgbe;

zcoorB =

(zlenA * 0.5 - 2. * Geo::MODULEWALLTHICKNESS - Geo::INTERCENTRMODBORDER1) * 0.5 - 2. * Geo::MODULEWALLTHICKNESS;

TVirtualMC::GetMC()->Gsvolu("FWZAD", "TRAP", getMediumID(kFiberGlass), trpa, 11); // Fibre glass

TVirtualMC::GetMC()->Gspos("FWZAD", 1, "FLTB", xcoor, ycoorB, zcoorB, idrotm[1], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZAD", 2, "FLTC", xcoor, ycoorB, -zcoorB, idrotm[0], "ONLY");

}

tgal = (yFLT - Geo::LENGTHINCEMODBORDERU - Geo::LENGTHINCEMODBORDERD) /

(Geo::INTERCENTRMODBORDER2 - Geo::INTERCENTRMODBORDER1);

alpha = TMath::ATan(tgal);

beta = (TMath::Pi() * 0.5 - alpha) * 0.5;

tgbe = TMath::Tan(beta);

trpa[0] = xFLT * 0.5;

trpa[1] = 0.;

trpa[2] = 0.;

trpa[3] = 2. * Geo::MODULEWALLTHICKNESS;

// trpa[4] = (Geo::LENGTHINCEMODBORDER - 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

// trpa[5] = (Geo::LENGTHINCEMODBORDER + 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

trpa[4] = (Geo::LENGTHINCEMODBORDERU - 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[5] = (Geo::LENGTHINCEMODBORDERU + 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[6] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

trpa[7] = 2. * Geo::MODULEWALLTHICKNESS;

trpa[8] = (Geo::LENGTHINCEMODBORDERU - 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[9] = (Geo::LENGTHINCEMODBORDERU + 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

// trpa[8] = (Geo::LENGTHINCEMODBORDER - 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

// trpa[9] = (Geo::LENGTHINCEMODBORDER + 2.*Geo::MODULEWALLTHICKNESS*tgbe)*0.5;

trpa[10] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

TVirtualMC::GetMC()->Gsvolu("FWZ1U", "TRAP", getMediumID(kFiberGlass), trpa, 11); // Fibre glass

Matrix(idrotm[2], 90., 270., 0., 0., 90., 180.);

Matrix(idrotm[3], 90., 270., 180., 0., 90., 0.);

// xcoor = 0.;

// ycoor = (yFLT - Geo::LENGTHINCEMODBORDER)*0.5;

ycoor = (yFLT - Geo::LENGTHINCEMODBORDERU) * 0.5;

zcoor = Geo::INTERCENTRMODBORDER2;

TVirtualMC::GetMC()->Gspos("FWZ1U", 1, "FLTA", xcoor, ycoor, zcoor, idrotm[2], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ1U", 2, "FLTA", xcoor, ycoor, -zcoor, idrotm[3], "ONLY");

if (mTOFHoles) {

// y0B = Geo::LENGTHINCEMODBORDER + Geo::MODULEWALLTHICKNESS*tgbe;

y0B = Geo::LENGTHINCEMODBORDERU + Geo::MODULEWALLTHICKNESS * tgbe;

trpa[0] = xFLT * 0.5;

trpa[1] = 0.;

trpa[2] = 0.;

trpa[3] = Geo::MODULEWALLTHICKNESS;

trpa[4] = (y0B - Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[5] = (y0B + Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[6] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

trpa[7] = Geo::MODULEWALLTHICKNESS;

trpa[8] = (y0B - Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[9] = (y0B + Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[10] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

TVirtualMC::GetMC()->Gsvolu("FWZBU", "TRAP", getMediumID(kFiberGlass), trpa, 11); // Fibre glass

// xcoor = 0.;

ycoorB = ycoor - Geo::MODULEWALLTHICKNESS * 0.5 * tgbe;

zcoorB = (zlenA * 0.5 - 2. * Geo::MODULEWALLTHICKNESS - Geo::INTERCENTRMODBORDER1) * 0.5 -

(Geo::INTERCENTRMODBORDER2 - Geo::INTERCENTRMODBORDER1) - 2. * Geo::MODULEWALLTHICKNESS;

TVirtualMC::GetMC()->Gspos("FWZBU", 1, "FLTB", xcoor, ycoorB, zcoorB, idrotm[3], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZBU", 2, "FLTC", xcoor, ycoorB, -zcoorB, idrotm[2], "ONLY");

}

trpa[0] = 0.5 * (Geo::INTERCENTRMODBORDER2 - Geo::INTERCENTRMODBORDER1) / TMath::Cos(alpha);

trpa[1] = 2. * Geo::MODULEWALLTHICKNESS;

trpa[2] = xFLT * 0.5;

trpa[3] = -beta * TMath::RadToDeg();

trpa[4] = 0.;

trpa[5] = 0.;

TVirtualMC::GetMC()->Gsvolu("FWZ2", "PARA", getMediumID(kFiberGlass), trpa, 6); // Fibre glass

Matrix(idrotm[4], alpha * TMath::RadToDeg(), 90., 90. + alpha * TMath::RadToDeg(), 90., 90., 180.);

Matrix(idrotm[5], 180. - alpha * TMath::RadToDeg(), 90., 90. - alpha * TMath::RadToDeg(), 90., 90., 0.);

// xcoor = 0.;

// ycoor = 0.;

ycoor = (Geo::LENGTHINCEMODBORDERD - Geo::LENGTHINCEMODBORDERU) * 0.5;

zcoor = (Geo::INTERCENTRMODBORDER2 + Geo::INTERCENTRMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gspos("FWZ2", 1, "FLTA", xcoor, ycoor, zcoor, idrotm[4], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ2", 2, "FLTA", xcoor, ycoor, -zcoor, idrotm[5], "ONLY");

if (mTOFHoles) {

trpa[0] = 0.5 * (Geo::INTERCENTRMODBORDER2 - Geo::INTERCENTRMODBORDER1) / TMath::Cos(alpha);

trpa[1] = Geo::MODULEWALLTHICKNESS;

trpa[2] = xFLT * 0.5;

trpa[3] = -beta * TMath::RadToDeg();

trpa[4] = 0.;

trpa[5] = 0.;

TVirtualMC::GetMC()->Gsvolu("FWZC", "PARA", getMediumID(kFiberGlass), trpa, 6); // Fibre glass

// xcoor = 0.;

ycoorB = ycoor - Geo::MODULEWALLTHICKNESS * tgbe;

zcoorB = (zlenA * 0.5 - 2. * Geo::MODULEWALLTHICKNESS - Geo::INTERCENTRMODBORDER1) * 0.5 -

(Geo::INTERCENTRMODBORDER2 - Geo::INTERCENTRMODBORDER1) * 0.5 - 2. * Geo::MODULEWALLTHICKNESS;

TVirtualMC::GetMC()->Gspos("FWZC", 1, "FLTB", xcoor, ycoorB, zcoorB, idrotm[5], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZC", 2, "FLTC", xcoor, ycoorB, -zcoorB, idrotm[4], "ONLY");

}

// Definition and positioning

// of the fibre glass walls between intermediate and lateral modules (FWZ3 and FWZ4)

tgal = (yFLT - 2. * Geo::LENGTHEXINMODBORDER) / (Geo::EXTERINTERMODBORDER2 - Geo::EXTERINTERMODBORDER1);

alpha = TMath::ATan(tgal);

beta = (TMath::Pi() * 0.5 - alpha) * 0.5;

tgbe = TMath::Tan(beta);

trpa[0] = xFLT * 0.5;

trpa[1] = 0.;

trpa[2] = 0.;

trpa[3] = 2. * Geo::MODULEWALLTHICKNESS;

trpa[4] = (Geo::LENGTHEXINMODBORDER - 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[5] = (Geo::LENGTHEXINMODBORDER + 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[6] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

trpa[7] = 2. * Geo::MODULEWALLTHICKNESS;

trpa[8] = (Geo::LENGTHEXINMODBORDER - 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[9] = (Geo::LENGTHEXINMODBORDER + 2. * Geo::MODULEWALLTHICKNESS * tgbe) * 0.5;

trpa[10] =

TMath::ATan(tgbe * 0.5) * TMath::RadToDeg(); // TMath::ATan((trpa[5] - trpa[4])/(2.*trpa[3]))*TMath::RadToDeg();

TVirtualMC::GetMC()->Gsvolu("FWZ3", "TRAP", getMediumID(kFiberGlass), trpa, 11); // Fibre glass

// xcoor = 0.;

ycoor = (yFLT - Geo::LENGTHEXINMODBORDER) * 0.5;

zcoor = Geo::EXTERINTERMODBORDER1;

TVirtualMC::GetMC()->Gspos("FWZ3", 1, "FLTA", xcoor, ycoor, zcoor, idrotm[3], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ3", 2, "FLTA", xcoor, ycoor, -zcoor, idrotm[2], "ONLY");

if (mTOFHoles) {

// xcoor = 0.;

// ycoor = (yFLT - Geo::LENGTHEXINMODBORDER)*0.5;

zcoor =

-Geo::EXTERINTERMODBORDER1 + (zlenA * 0.5 + Geo::INTERCENTRMODBORDER1 - 2. * Geo::MODULEWALLTHICKNESS) * 0.5;

TVirtualMC::GetMC()->Gspos("FWZ3", 5, "FLTB", xcoor, ycoor, zcoor, idrotm[2], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ3", 6, "FLTC", xcoor, ycoor, -zcoor, idrotm[3], "ONLY");

}

// xcoor = 0.;

ycoor = -(yFLT - Geo::LENGTHEXINMODBORDER) * 0.5;

zcoor = Geo::EXTERINTERMODBORDER2;

TVirtualMC::GetMC()->Gspos("FWZ3", 3, "FLTA", xcoor, ycoor, zcoor, idrotm[1], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ3", 4, "FLTA", xcoor, ycoor, -zcoor, idrotm[0], "ONLY");

if (mTOFHoles) {

// xcoor = 0.;

// ycoor = -(yFLT - Geo::LENGTHEXINMODBORDER)*0.5;

zcoor =

-Geo::EXTERINTERMODBORDER2 + (zlenA * 0.5 + Geo::INTERCENTRMODBORDER1 - 2. * Geo::MODULEWALLTHICKNESS) * 0.5;

TVirtualMC::GetMC()->Gspos("FWZ3", 7, "FLTB", xcoor, ycoor, zcoor, idrotm[0], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ3", 8, "FLTC", xcoor, ycoor, -zcoor, idrotm[1], "ONLY");

}

trpa[0] = 0.5 * (Geo::EXTERINTERMODBORDER2 - Geo::EXTERINTERMODBORDER1) / TMath::Cos(alpha);

trpa[1] = 2. * Geo::MODULEWALLTHICKNESS;

trpa[2] = xFLT * 0.5;

trpa[3] = -beta * TMath::RadToDeg();

trpa[4] = 0.;

trpa[5] = 0.;

TVirtualMC::GetMC()->Gsvolu("FWZ4", "PARA", getMediumID(kFiberGlass), trpa, 6); // Fibre glass

Matrix(idrotm[6], alpha * TMath::RadToDeg(), 90., 90. + alpha * TMath::RadToDeg(), 90., 90., 180.);

Matrix(idrotm[7], 180. - alpha * TMath::RadToDeg(), 90., 90. - alpha * TMath::RadToDeg(), 90., 90., 0.);

// xcoor = 0.;

ycoor = 0.;

zcoor = (Geo::EXTERINTERMODBORDER2 + Geo::EXTERINTERMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gspos("FWZ4", 1, "FLTA", xcoor, ycoor, zcoor, idrotm[7], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ4", 2, "FLTA", xcoor, ycoor, -zcoor, idrotm[6], "ONLY");

if (mTOFHoles) {

// xcoor = 0.;

// ycoor = 0.;

zcoor = -(Geo::EXTERINTERMODBORDER2 + Geo::EXTERINTERMODBORDER1) * 0.5 +

(zlenA * 0.5 + Geo::INTERCENTRMODBORDER1 - 2. * Geo::MODULEWALLTHICKNESS) * 0.5;

TVirtualMC::GetMC()->Gspos("FWZ4", 3, "FLTB", xcoor, ycoor, zcoor, idrotm[6], "ONLY");

TVirtualMC::GetMC()->Gspos("FWZ4", 4, "FLTC", xcoor, ycoor, -zcoor, idrotm[7], "ONLY");

}

}

void Detector::makeStripsInModules(Float_t ytof, Float_t zlenA) const

{

//

// Define MRPC strip volume, called FSTR

// Insert FSTR volume in FLTA/B/C volumes

//

Float_t yFLT = ytof * 0.5 - Geo::MODULEWALLTHICKNESS;

///////////////// Detector itself //////////////////////

// new description for strip volume -double stack strip-

// -- all constants are expressed in cm

// height of different layers

constexpr Float_t HGLFY = Geo::HFILIY + 2. * Geo::HGLASSY; // height of GLASS Layer

constexpr Float_t LSENSMX = Geo::NPADX * Geo::XPAD; // length of Sensitive Layer

constexpr Float_t HSENSMY = Geo::HSENSMY; // height of Sensitive Layer

constexpr Float_t WSENSMZ = Geo::NPADZ * Geo::ZPAD; // width of Sensitive Layer

// height of the FSTR Volume (the strip volume)

constexpr Float_t HSTRIPY = 2. * Geo::HHONY + 2. * Geo::HPCBY + 4. * Geo::HRGLY + 2. * HGLFY + Geo::HCPCBY;

// width of the FSTR Volume (the strip volume)

constexpr Float_t WSTRIPZ = Geo::WCPCBZ;

// length of the FSTR Volume (the strip volume)

constexpr Float_t LSTRIPX = Geo::STRIPLENGTH;

// FSTR volume definition-filling this volume with non sensitive Gas Mixture

Float_t parfp[3] = {static_cast<Float_t>(LSTRIPX * 0.5), static_cast<Float_t>(HSTRIPY * 0.5),

static_cast<Float_t>(WSTRIPZ * 0.5)};

TVirtualMC::GetMC()->Gsvolu("FSTR", "BOX", getMediumID(kFre), parfp, 3); // Freon mix

Float_t posfp[3] = {0., 0., 0.};

// NOMEX (HONEYCOMB) Layer definition

// parfp[0] = LSTRIPX*0.5;

parfp[1] = Geo::HHONY * 0.5;

parfp[2] = Geo::WHONZ * 0.5;

TVirtualMC::GetMC()->Gsvolu("FHON", "BOX", getMediumID(kNomex), parfp, 3); // Nomex (Honeycomb)

// positioning 2 NOMEX Layers on FSTR volume

// posfp[0] = 0.;

posfp[1] = -HSTRIPY * 0.5 + parfp[1];

// posfp[2] = 0.;

TVirtualMC::GetMC()->Gspos("FHON", 1, "FSTR", 0., posfp[1], 0., 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FHON", 2, "FSTR", 0., -posfp[1], 0., 0, "ONLY");

// Lower PCB Layer definition

// parfp[0] = LSTRIPX*0.5;

parfp[1] = Geo::HPCBY * 0.5;

parfp[2] = Geo::WPCBZ1 * 0.5;

TVirtualMC::GetMC()->Gsvolu("FPC1", "BOX", getMediumID(kG10), parfp, 3); // G10

// Upper PCB Layer definition

// parfp[0] = LSTRIPX*0.5;

// parfp[1] = Geo::HPCBY*0.5;

parfp[2] = Geo::WPCBZ2 * 0.5;

TVirtualMC::GetMC()->Gsvolu("FPC2", "BOX", getMediumID(kG10), parfp, 3); // G10

// positioning 2 external PCB Layers in FSTR volume

// posfp[0] = 0.;

posfp[1] = -HSTRIPY * 0.5 + Geo::HHONY + parfp[1];

// posfp[2] = 0.;

TVirtualMC::GetMC()->Gspos("FPC1", 1, "FSTR", 0., -posfp[1], 0., 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FPC2", 1, "FSTR", 0., posfp[1], 0., 0, "ONLY");

// Central PCB layer definition

// parfp[0] = LSTRIPX*0.5;

parfp[1] = Geo::HCPCBY * 0.5;

parfp[2] = Geo::WCPCBZ * 0.5;

TVirtualMC::GetMC()->Gsvolu("FPCB", "BOX", getMediumID(kG10), parfp, 3); // G10

gGeoManager->GetVolume("FPCB")->VisibleDaughters(kFALSE);

// positioning the central PCB layer

TVirtualMC::GetMC()->Gspos("FPCB", 1, "FSTR", 0., 0., 0., 0, "ONLY");

// Sensitive volume definition

Float_t parfs[3] = {static_cast<Float_t>(LSENSMX * 0.5), static_cast<Float_t>(HSENSMY * 0.5),

static_cast<Float_t>(WSENSMZ * 0.5)};

TVirtualMC::GetMC()->Gsvolu("FSEN", "BOX", getMediumID(kCuS), parfs, 3); // Cu sensitive

// printf("check material\n");

// printf("ID used = %i\n",getMediumID(kCuS));

// printf("ID needed = %i\n",gGeoManager->GetMedium("TOF_Cu-S$")->GetId());

// getchar();

// dividing FSEN along z in Geo::NPADZ=2 and along x in Geo::NPADX=48

TVirtualMC::GetMC()->Gsdvn("FSEZ", "FSEN", Geo::NPADZ, 3);

TVirtualMC::GetMC()->Gsdvn("FPAD", "FSEZ", Geo::NPADX, 1);

// positioning sensitive layer inside FPCB

TVirtualMC::GetMC()->Gspos("FSEN", 1, "FPCB", 0., 0., 0., 0, "ONLY");

// RED GLASS Layer definition

// parfp[0] = LSTRIPX*0.5;

parfp[1] = Geo::HRGLY * 0.5;

parfp[2] = Geo::WRGLZ * 0.5;

TVirtualMC::GetMC()->Gsvolu("FRGL", "BOX", getMediumID(kGlass), parfp, 3); // red glass

// positioning 4 RED GLASS Layers in FSTR volume

// posfp[0] = 0.;

posfp[1] = -HSTRIPY * 0.5 + Geo::HHONY + Geo::HPCBY + parfp[1];

// posfp[2] = 0.;

TVirtualMC::GetMC()->Gspos("FRGL", 1, "FSTR", 0., posfp[1], 0., 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FRGL", 4, "FSTR", 0., -posfp[1], 0., 0, "ONLY");

// posfp[0] = 0.;

posfp[1] = (Geo::HCPCBY + Geo::HRGLY) * 0.5;

// posfp[2] = 0.;

TVirtualMC::GetMC()->Gspos("FRGL", 2, "FSTR", 0., -posfp[1], 0., 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FRGL", 3, "FSTR", 0., posfp[1], 0., 0, "ONLY");

// GLASS Layer definition

// parfp[0] = LSTRIPX*0.5;

parfp[1] = Geo::HGLASSY;

parfp[2] = Geo::WGLFZ * 0.5;

TVirtualMC::GetMC()->Gsvolu("FGLF", "BOX", getMediumID(kGlass), parfp, 3); // glass

// positioning 2 GLASS Layers in FSTR volume

// posfp[0] = 0.;

posfp[1] = (Geo::HCPCBY + HGLFY) * 0.5 + Geo::HRGLY;

// posfp[2] = 0.;

TVirtualMC::GetMC()->Gspos("FGLF", 1, "FSTR", 0., -posfp[1], 0., 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FGLF", 2, "FSTR", 0., posfp[1], 0., 0, "ONLY");

// Positioning the Strips (FSTR volumes) in the FLT volumes

Int_t maxStripNumbers[5] = {Geo::NSTRIPC, Geo::NSTRIPB, Geo::NSTRIPA, Geo::NSTRIPB, Geo::NSTRIPC};

Int_t idrotm[Geo::NSTRIPXSECTOR];

for (Int_t ii = 0; ii < Geo::NSTRIPXSECTOR; ii++) {

idrotm[ii] = 0;

}

Int_t totalStrip = 0;

Float_t xpos, zpos, ypos, ang;

for (Int_t iplate = 0; iplate < Geo::NPLATES; iplate++) {

if (iplate > 0) {

totalStrip += maxStripNumbers[iplate - 1];

}

for (Int_t istrip = 0; istrip < maxStripNumbers[iplate]; istrip++) {

ang = Geo::getAngles(iplate, istrip);

if (ang > 0.) {

Matrix(idrotm[istrip + totalStrip], 90., 0., 90. + ang, 90., ang, 90.);

} else if (ang == 0.) {

Matrix(idrotm[istrip + totalStrip], 90., 0., 90., 90., 0., 0.);

} else if (ang < 0.) {

Matrix(idrotm[istrip + totalStrip], 90., 0., 90. + ang, 90., -ang, 270.);

}

xpos = 0.;

ypos = Geo::getHeights(iplate, istrip) + yFLT * 0.5;

zpos = Geo::getDistances(iplate, istrip);

TVirtualMC::GetMC()->Gspos("FSTR", istrip + totalStrip + 1, "FLTA", xpos, ypos, -zpos,

idrotm[istrip + totalStrip], "ONLY");

if (mTOFHoles) {

if (istrip + totalStrip + 1 > 53) {

TVirtualMC::GetMC()->Gspos(

"FSTR", istrip + totalStrip + 1, "FLTC", xpos, ypos,

-zpos - (zlenA * 0.5 - 2. * Geo::MODULEWALLTHICKNESS + Geo::INTERCENTRMODBORDER1) * 0.5,

idrotm[istrip + totalStrip], "ONLY");

}

if (istrip + totalStrip + 1 < 39) {

TVirtualMC::GetMC()->Gspos(

"FSTR", istrip + totalStrip + 1, "FLTB", xpos, ypos,

-zpos + (zlenA * 0.5 - 2. * Geo::MODULEWALLTHICKNESS + Geo::INTERCENTRMODBORDER1) * 0.5,

idrotm[istrip + totalStrip], "ONLY");

}

}

}

}

}

void Detector::createModuleCovers(Float_t xtof, Float_t zlenA) const

{

//

// Create covers for module:

// per each module zone, defined according to

// fgkInterCentrModBorder2, fgkExterInterModBorder1 and zlenA+2 values,

// there is a frame of thickness 2cm in Al

// and the contained zones in honeycomb of Al.

// There is also an interface layer (1.6mm thichness)

// and plastic and Cu corresponding to the flat cables.

//

Float_t par[3];

par[0] = xtof * 0.5 + 2.;

par[1] = Geo::MODULECOVERTHICKNESS * 0.5;

par[2] = zlenA * 0.5 + 2.;

TVirtualMC::GetMC()->Gsvolu("FPEA", "BOX ", getMediumID(kAir), par, 3); // Air

if (mTOFHoles) {

TVirtualMC::GetMC()->Gsvolu("FPEB", "BOX ", getMediumID(kAir), par, 3); // Air

}

constexpr Float_t ALCOVERTHICKNESS = 1.5;

constexpr Float_t INTERFACECARDTHICKNESS = 0.16;

constexpr Float_t ALSKINTHICKNESS = 0.1;

constexpr Float_t PLASTICFLATCABLETHICKNESS = 0.25;

constexpr Float_t COPPERFLATCABLETHICKNESS = 0.01;

// par[0] = xtof*0.5 + 2.;

par[1] = ALCOVERTHICKNESS * 0.5;

// par[2] = zlenA*0.5 + 2.;

TVirtualMC::GetMC()->Gsvolu("FALT", "BOX ", getMediumID(kAlFrame), par, 3); // Al

if (mTOFHoles) {

TVirtualMC::GetMC()->Gsvolu("FALB", "BOX ", getMediumID(kAlFrame), par, 3); // Al

}

Float_t xcoor, ycoor, zcoor;

xcoor = 0.;

ycoor = 0.;

zcoor = 0.;

TVirtualMC::GetMC()->Gspos("FALT", 0, "FPEA", xcoor, ycoor, zcoor, 0, "ONLY");

if (mTOFHoles) {

TVirtualMC::GetMC()->Gspos("FALB", 0, "FPEB", xcoor, ycoor, zcoor, 0, "ONLY");

}

par[0] = xtof * 0.5;

// par[1] = ALCOVERTHICKNESS*0.5;

par[2] = Geo::INTERCENTRMODBORDER2 - 2.;

TVirtualMC::GetMC()->Gsvolu("FPE1", "BOX ", getMediumID(kHoneycomb), par, 3); // Al honeycomb

// xcoor = 0.;

// ycoor = 0.;

// zcoor = 0.;

TVirtualMC::GetMC()->Gspos("FPE1", 0, "FALT", xcoor, ycoor, zcoor, 0, "ONLY");

if (mTOFHoles) {

// par[0] = xtof*0.5;

par[1] = ALCOVERTHICKNESS * 0.5 - ALSKINTHICKNESS;

// par[2] = Geo::INTERCENTRMODBORDER2 - 2.;

TVirtualMC::GetMC()->Gsvolu("FPE4", "BOX ", getMediumID(kHoneyHoles), par, 3); // Al honeycomb for holes

// xcoor = 0.;

// ycoor = 0.;

// zcoor = 0.;

TVirtualMC::GetMC()->Gspos("FPE4", 0, "FALB", xcoor, ycoor, zcoor, 0, "ONLY");

}

// par[0] = xtof*0.5;

// par[1] = ALCOVERTHICKNESS*0.5;

par[2] = (Geo::EXTERINTERMODBORDER1 - Geo::INTERCENTRMODBORDER2) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FPE2", "BOX ", getMediumID(kHoneycomb), par, 3); // Al honeycomb

// xcoor = 0.;

// ycoor = 0.;

zcoor = (Geo::EXTERINTERMODBORDER1 + Geo::INTERCENTRMODBORDER2) * 0.5;

TVirtualMC::GetMC()->Gspos("FPE2", 1, "FALT", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FPE2", 2, "FALT", xcoor, ycoor, -zcoor, 0, "ONLY");

if (mTOFHoles) {

// xcoor = 0.;

// ycoor = 0.;

// zcoor = (Geo::EXTERINTERMODBORDER1 + Geo::INTERCENTRMODBORDER2)*0.5;

TVirtualMC::GetMC()->Gspos("FPE2", 1, "FALB", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FPE2", 2, "FALB", xcoor, ycoor, -zcoor, 0, "ONLY");

}

// par[0] = xtof*0.5;

// par[1] = ALCOVERTHICKNESS*0.5;

par[2] = (zlenA * 0.5 + 2. - Geo::EXTERINTERMODBORDER1) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FPE3", "BOX ", getMediumID(kHoneycomb), par, 3); // Al honeycomb

// xcoor = 0.;

// ycoor = 0.;

zcoor = (zlenA * 0.5 + 2. + Geo::EXTERINTERMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gspos("FPE3", 1, "FALT", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FPE3", 2, "FALT", xcoor, ycoor, -zcoor, 0, "ONLY");

if (mTOFHoles) {

// xcoor = 0.;

// ycoor = 0.;

zcoor = (zlenA * 0.5 + 2. + Geo::EXTERINTERMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gspos("FPE3", 1, "FALB", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FPE3", 2, "FALB", xcoor, ycoor, -zcoor, 0, "ONLY");

}

// volumes for Interface cards

par[0] = xtof * 0.5;

par[1] = INTERFACECARDTHICKNESS * 0.5;

par[2] = Geo::INTERCENTRMODBORDER2 - 2.;

TVirtualMC::GetMC()->Gsvolu("FIF1", "BOX ", getMediumID(kG10), par, 3); // G10

// xcoor = 0.;

ycoor = ALCOVERTHICKNESS * 0.5 + INTERFACECARDTHICKNESS * 0.5;

zcoor = 0.;

TVirtualMC::GetMC()->Gspos("FIF1", 0, "FPEA", xcoor, ycoor, zcoor, 0, "ONLY");

// par[0] = xtof*0.5;

// par[1] = INTERFACECARDTHICKNESS*0.5;

par[2] = (Geo::EXTERINTERMODBORDER1 - Geo::INTERCENTRMODBORDER2) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FIF2", "BOX ", getMediumID(kG10), par, 3); // G10

// xcoor = 0.;

// ycoor = ALCOVERTHICKNESS*0.5 + INTERFACECARDTHICKNESS*0.5;

zcoor = (Geo::EXTERINTERMODBORDER1 + Geo::INTERCENTRMODBORDER2) * 0.5;

TVirtualMC::GetMC()->Gspos("FIF2", 1, "FPEA", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FIF2", 2, "FPEA", xcoor, ycoor, -zcoor, 0, "ONLY");

if (mTOFHoles) {

TVirtualMC::GetMC()->Gspos("FIF2", 1, "FPEB", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FIF2", 2, "FPEB", xcoor, ycoor, -zcoor, 0, "ONLY");

}

// par[0] = xtof*0.5;

// par[1] = INTERFACECARDTHICKNESS*0.5;

par[2] = (zlenA * 0.5 + 2. - Geo::EXTERINTERMODBORDER1) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FIF3", "BOX ", getMediumID(kG10), par, 3); // G10

// xcoor = 0.;

// ycoor = ALCOVERTHICKNESS*0.5 + INTERFACECARDTHICKNESS*0.5;

zcoor = (zlenA * 0.5 + 2. + Geo::EXTERINTERMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gspos("FIF3", 1, "FPEA", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FIF3", 2, "FPEA", xcoor, ycoor, -zcoor, 0, "ONLY");

if (mTOFHoles) {

TVirtualMC::GetMC()->Gspos("FIF3", 1, "FPEB", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FIF3", 2, "FPEB", xcoor, ycoor, -zcoor, 0, "ONLY");

}

// volumes for flat cables

// plastic

par[0] = xtof * 0.5;

par[1] = PLASTICFLATCABLETHICKNESS * 0.5;

par[2] = Geo::INTERCENTRMODBORDER2 - 2.;

TVirtualMC::GetMC()->Gsvolu("FFC1", "BOX ", getMediumID(kPlastic), par, 3); // Plastic (CH2)

// xcoor = 0.;

ycoor = -ALCOVERTHICKNESS * 0.5 - PLASTICFLATCABLETHICKNESS * 0.5;

zcoor = 0.;

TVirtualMC::GetMC()->Gspos("FFC1", 0, "FPEA", xcoor, ycoor, zcoor, 0, "ONLY");

// par[0] = xtof*0.5;

// par[1] = PLASTICFLATCABLETHICKNESS*0.5;

par[2] = (Geo::EXTERINTERMODBORDER1 - Geo::INTERCENTRMODBORDER2) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FFC2", "BOX ", getMediumID(kPlastic), par, 3); // Plastic (CH2)

// xcoor = 0.;

// ycoor = -ALCOVERTHICKNESS*0.5 - PLASTICFLATCABLETHICKNESS*0.5;

zcoor = (Geo::EXTERINTERMODBORDER1 + Geo::INTERCENTRMODBORDER2) * 0.5;

TVirtualMC::GetMC()->Gspos("FFC2", 1, "FPEA", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FFC2", 2, "FPEA", xcoor, ycoor, -zcoor, 0, "ONLY");

if (mTOFHoles) {

TVirtualMC::GetMC()->Gspos("FFC2", 1, "FPEB", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FFC2", 2, "FPEB", xcoor, ycoor, -zcoor, 0, "ONLY");

}

// par[0] = xtof*0.5;

// par[1] = PLASTICFLATCABLETHICKNESS*0.5;

par[2] = (zlenA * 0.5 + 2. - Geo::EXTERINTERMODBORDER1) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FFC3", "BOX ", getMediumID(kPlastic), par, 3); // Plastic (CH2)

// xcoor = 0.;

// ycoor = -ALCOVERTHICKNESS*0.5 - PLASTICFLATCABLETHICKNESS*0.5;

zcoor = (zlenA * 0.5 + 2. + Geo::EXTERINTERMODBORDER1) * 0.5;

TVirtualMC::GetMC()->Gspos("FFC3", 1, "FPEA", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FFC3", 2, "FPEA", xcoor, ycoor, -zcoor, 0, "ONLY");

if (mTOFHoles) {

TVirtualMC::GetMC()->Gspos("FFC3", 1, "FPEB", xcoor, ycoor, zcoor, 0, "ONLY");

TVirtualMC::GetMC()->Gspos("FFC3", 2, "FPEB", xcoor, ycoor, -zcoor, 0, "ONLY");

}

// Cu

par[0] = xtof * 0.5;

par[1] = COPPERFLATCABLETHICKNESS * 0.5;

par[2] = Geo::INTERCENTRMODBORDER2 - 2.;

TVirtualMC::GetMC()->Gsvolu("FCC1", "BOX ", getMediumID(kCopper), par, 3); // Cu

TVirtualMC::GetMC()->Gspos("FCC1", 0, "FFC1", 0., 0., 0., 0, "ONLY");

// par[0] = xtof*0.5;

// par[1] = COPPERFLATCABLETHICKNESS*0.5;

par[2] = (Geo::EXTERINTERMODBORDER1 - Geo::INTERCENTRMODBORDER2) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FCC2", "BOX ", getMediumID(kCopper), par, 3); // Cu

TVirtualMC::GetMC()->Gspos("FCC2", 0, "FFC2", 0., 0., 0., 0, "ONLY");

// par[0] = xtof*0.5;

// par[1] = COPPERFLATCABLETHICKNESS*0.5;

par[2] = (zlenA * 0.5 + 2. - Geo::EXTERINTERMODBORDER1) * 0.5 - 2.;

TVirtualMC::GetMC()->Gsvolu("FCC3", "BOX ", getMediumID(kCopper), par, 3); // Cu

TVirtualMC::GetMC()->Gspos("FCC3", 0, "FFC3", 0., 0., 0., 0, "ONLY");

}

void Detector::createBackZone(Float_t xtof, Float_t ytof, Float_t zlenA) const

{