MuonGM::Mdt Class Reference

#include <Mdt.h>

Inheritance diagram for MuonGM::Mdt:

Collaboration diagram for MuonGM::Mdt:

Public Member Functions
	Mdt (const MYSQL &mysql, Component *s1, const std::string &s2)
GeoFullPhysVol *	build (StoredMaterialManager &matManager, const MYSQL &mysql)
GeoFullPhysVol *	build (StoredMaterialManager &matManager, const MYSQL &mysql, std::vector< Cutout * > &)
void	processCutouts (std::vector< Cutout * > &vcutdef)
virtual void	print () const override
void	setLogVolName (const std::string &str)

Public Attributes
double	width {0.}
double	length {0.}
double	thickness {0.}
double	longWidth {0.}
int	index {0}
double	tubelenStepSize {0.}
double	tubePitch {0.}
std::unique_ptr< MultiLayer >	layer {}
std::string	name {}
std::string	logVolName {}

Private Attributes
MdtComponent *	m_component {nullptr}

Detailed Description

Definition at line 23 of file Mdt.h.

Constructor & Destructor Documentation

Mdt()

MuonGM::Mdt::Mdt	(	const MYSQL &	mysql,
		Component *	s1,
		const std::string &	s2 )

Definition at line 32 of file Mdt.cxx.

                                                     : DetectorElement(ss->name) {
        logVolName = lVName;
        MdtComponent *s = static_cast<MdtComponent *>(ss);
        const MDT *thism = dynamic_cast<const MDT*>(mysql.GetTechnology(s->name));
 
        width = s->dx1;
        longWidth = s->dx2;
        thickness = s->GetThickness(mysql);
        length = s->dy;
        m_component = s;
        m_component->cutoutTubeXShift = 0.;
        layer = std::make_unique<MultiLayer>(mysql, s->name);
        layer->logVolName = lVName;
        layer->cutoutNsteps = 0;
        layer->width = width;
        layer->longWidth = longWidth;
        tubePitch = thism->pitch;
        layer->length = length;
        layer->nrOfTubes = (int)(layer->length / thism->pitch);
        if (longWidth > width)
            layer->nrOfSteps = int(length / s->tubelenStepSize);
        tubelenStepSize = s->tubelenStepSize;
        index = s->index;
    }

Member Function Documentation

build() [1/2]

GeoFullPhysVol * MuonGM::Mdt::build	(	StoredMaterialManager &	matManager,
		const MYSQL &	mysql )

Definition at line 60 of file Mdt.cxx.

                                                   {
        layer->cutoutNsteps = 1;
        return layer->build(matManager, mysql);
    }

build() [2/2]

GeoFullPhysVol * MuonGM::Mdt::build	(	StoredMaterialManager &	matManager,
		const MYSQL &	mysql,
		std::vector< Cutout * > &	vcutdef )

Definition at line 66 of file Mdt.cxx.

                                                             {
        MsgStream log(Athena::getMessageSvc(), "Mdt::build");
 
        int Ncuts = vcutdef.size();
        if (Ncuts > 0) {
 
        processCutouts(vcutdef);
        return layer->build(matManager, mysql);
 
        } else {
            return build(matManager, mysql);
        }
    }

print()

void MuonGM::Mdt::print ( ) const

overridevirtual

Implements MuonGM::DetectorElement.

Definition at line 82 of file Mdt.cxx.

                          {
        MsgStream log(Athena::getMessageSvc(), "MuonGM::Mdt");
        log << MSG::INFO << "Mdt " << name.c_str() << " :" << endmsg;
    }

processCutouts()

void MuonGM::Mdt::processCutouts

(

std::vector< Cutout * > &

vcutdef

)

Definition at line 87 of file Mdt.cxx.

                                                         {
      MsgStream log(Athena::getMessageSvc(), "Mdt::processCutouts");
      if (verbose_mdt) {
    log << MSG::VERBOSE << " mdt cutouts are on " << endmsg;
      }
      
      int Ncuts = vcutdef.size();
      bool cutoutsVsX {false}, cutoutsVsY{false};
      // first check whether there are several cutouts with different amdb x or y coordinates
      float lastX{FLT_MAX}, lastY{FLT_MAX};
      for (int i = 0; i < Ncuts; ++i) {
    log<<MSG::DEBUG<<__FILE__":"<<__LINE__
       <<" name: "<<name
       <<" logVolName: "<<logVolName
       <<" Ncuts: "<<Ncuts
       <<" lastX: "<<lastX
       <<" cutDef->dX "<<vcutdef[i]->dx
       <<" lastY: "<<lastY
       <<" dutDef->dY "<<vcutdef[i]->dy<<endmsg;
    
    if (lastX != FLT_MAX && lastX * vcutdef[i]->dx < 0)
      cutoutsVsX = true;
    lastX = vcutdef[i]->dx;
    if (lastY != FLT_MAX && lastY * vcutdef[i]->dy < 0)
      cutoutsVsY = true;
    lastY = vcutdef[i]->dy;
      }
      if (cutoutsVsX && cutoutsVsY) {               
    throw std::runtime_error(
                 Form("%s:%d \nMdt::build() - Found more than one cutout in amdb-x  direction and more than one cutout in amdb-y direction, currently not supported",
                      __FILE__, __LINE__));
      }
      
      if (!cutoutsVsX) {              // nominal case (used for BMS/BOG/BMG etc.)
    std::array<int, 5> cutoutNtubes{};        // Number of tubes in sub-multilayer[i]
    std::array<bool, 5> cutoutFullLength{};   // True if this region is outside the cutout
    std::array<double, 5> cutoutXtubes{};     // Location of tube center within sub-ml[i] along x-amdb
    std::array<double, 5> cutoutTubeLength{}; // Tube length
    std::array<double, 5> cutoutYmax{};
    
    for (int i = 0; i < 5; i++) {
      cutoutFullLength[i] = true;
      cutoutTubeLength[i] = width;
    }
    
    // Order cutouts by increasing dy
    for (int i = 0; i < Ncuts; i++) {
      for (int j = i + 1; j < Ncuts; j++) {
        if (vcutdef[j]->dy < vcutdef[i]->dy) {
          Cutout *c = vcutdef[i];
          vcutdef[i] = vcutdef[j];
          vcutdef[j] = c;
        }
      }
    }
    
    // Set up cut location code
    double top = length - tubePitch / 2.;
    int cutLocationCode[3] = {0, 0, 0};
    for (int i = 0; i < Ncuts; i++) {
      if (vcutdef[i]->dy <= 0)
        cutLocationCode[i] = -1;
      if (round(vcutdef[i]->dy + vcutdef[i]->lengthY, 2) >= round(top, 2))
        cutLocationCode[i] = 1;
    }
    
    // Calculate quantities needed by multilayer
    double twidth{0.}, xmin{0.}, xmax{0.};
    bool cutAtAngle = false;
    int Nsteps = 0;
    for (int i = 0; i < Ncuts; i++) {
      Cutout *c = vcutdef[i];
      if (c->dead1 > 1.)
        cutAtAngle = true;
      twidth = width + (longWidth - width) * c->dy / length;
      xmin = -twidth / 2. < c->dx - c->widthXs / 2. ? -twidth / 2. : c->dx + c->widthXs / 2.;
      xmax = twidth / 2. > c->dx + c->widthXs / 2. ? twidth / 2. : c->dx - c->widthXs / 2.;
      if (cutLocationCode[i] == -1) {
        cutoutYmax[Nsteps] = c->lengthY;
        cutoutTubeLength[Nsteps] = xmax - xmin;
        cutoutXtubes[Nsteps] = (xmax + xmin) / 2.;
        cutoutFullLength[Nsteps] = false;
        Nsteps++;
      } else if (cutLocationCode[i] == 1) {
        cutoutYmax[Nsteps] = c->dy;
        Nsteps++;
        cutoutTubeLength[Nsteps] = xmax - xmin;
        cutoutXtubes[Nsteps] = (xmax + xmin) / 2.;
        cutoutFullLength[Nsteps] = false;
      } else {
        cutoutYmax[Nsteps] = c->dy;
        Nsteps++;
        cutoutYmax[Nsteps] = c->dy + c->lengthY;
        cutoutTubeLength[Nsteps] = xmax - xmin;
        cutoutXtubes[Nsteps] = (xmax + xmin) / 2.;
        cutoutFullLength[Nsteps] = false;
        Nsteps++;
      }
    }
    cutoutYmax[Nsteps] = top;
    Nsteps++;
    
    double regionLength{0.}, low{0.};
    int fullLengthCounter{0}, tubeCounter{0};
    for (int i = 0; i < Nsteps; i++) {
      if (cutoutFullLength[i])
        fullLengthCounter++;
      
      regionLength = cutoutYmax[i] - low;
      cutoutNtubes[i] = int(regionLength / tubePitch);
      if ((regionLength / tubePitch - cutoutNtubes[i]) > 0.5)
        cutoutNtubes[i] += 1;
      
      if (fullLengthCounter > 1)
        cutoutNtubes[i]++;
      low = cutoutYmax[i];
      tubeCounter += cutoutNtubes[i];
    }
    if (tubeCounter > layer->nrOfTubes)
      --cutoutNtubes[Nsteps - 1];
    
    if (verbose_mdt) {
      for (int i = 0; i < Nsteps; i++) {
        log << MSG::VERBOSE << " cutoutYmax[" << i << "] = " << cutoutYmax[i] << " cutoutTubeLength[" << i << "] = " << cutoutTubeLength[i] << " cutoutXtubes[" << i
        << "] = " << cutoutXtubes[i] << " cutoutFullLength[" << i << "] = " << cutoutFullLength[i] << " cutoutNtubes[" << i << "] = " << cutoutNtubes[i]
        << endmsg;
      }
    }
    
    // encoding BMG chambers in Nsteps: Nsteps negative => BMG chamber
    // a/ 1st digit: eta index (1 to 3)
    // b/ 2nd digit: 1 == A side, 2 == C side
    // c/ 3rd digit: multilayer (1 or 2)
    // d/ last 2 digits: sector (12 or 14)
    
    if (logVolName.find("MDT10") != std::string::npos) {
      // multilayer 1 of BMG1A12, BMG2A12, BMG3A12, BMG1C14, BMG2C14, BMG3C14
      if (vcutdef[0]->icut == 1) { // these are BMG1A12 and BMG1C14
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -11112; // cut angle for A side positive BMG1A12
        else
          Nsteps = -12114; // cut angle for C side negative BMG1C14
      } else if (vcutdef[0]->icut == 2) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -21112; // cut angle for A side positive BMG2A12
        else
          Nsteps = -22114; // cut angle for C side negative BMG2C14
      } else if (vcutdef[0]->icut == 3) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -31112; // cut angle for A side positive BMG3A12
        else
          Nsteps = -32114; // cut angle for C side negative BMG3C14
      } else {
        log << MSG::ERROR << "massive error with MDT10 (BMG chambers)" << endmsg;
        std::abort();
      }
    } else if (logVolName.find("MDT11") != std::string::npos) {
      // multilayer 1 of BMG1A14, BMG2A14, BMG3A14, BMG1C12, BMG2C12, BMG3C12
      if (vcutdef[0]->icut == 1) { // these are BMG1A12 and BMG1C14
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -11114; // cut angle for A side positive BMG1A14
        else
          Nsteps = -12112; // cut angle for C side negative BMG1C12
      } else if (vcutdef[0]->icut == 2) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -21114; // cut angle for A side positive BMG2A14
        else
          Nsteps = -22112; // cut angle for C side negative BMG2C12
      } else if (vcutdef[0]->icut == 3) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -31114; // cut angle for A side positive BMG3A14
        else
          Nsteps = -32112; // cut angle for C side negative BMG3C12
      } else {
        log << MSG::ERROR << "massive error with MDT10 (BMG chambers)" << endmsg;
        std::abort();
      }
    } else if (logVolName.find("MDT12") != std::string::npos) {
      // multilayer 2 of BMG1A12, BMG2A12, BMG3A12, BMG1C14, BMG2C14, BMG3C14
      if (vcutdef[0]->icut == 1) { // these are BMG1A12 and BMG1C14
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -11212; // cut angle for A side positive BMG1A12
        else
          Nsteps = -12214; // cut angle for C side negative BMG1C14
      } else if (vcutdef[0]->icut == 2) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -21212; // cut angle for A side positive BMG2A12
        else
          Nsteps = -22214; // cut angle for C side negative BMG2C14
      } else if (vcutdef[0]->icut == 3) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -31212; // cut angle for A side positive BMG3A12
        else
          Nsteps = -32214; // cut angle for C side negative BMG3C14
      } else {
        log << MSG::ERROR << "massive error with MDT10 (BMG chambers)" << endmsg;
        std::abort();
      }
    } else if (logVolName.find("MDT13") != std::string::npos) {
      // multilayer 2 of BMG1A14, BMG2A14, BMG3A14, BMG1C12, BMG2C12, BMG3C12
      if (vcutdef[0]->icut == 1) { // these are BMG1A12 and BMG1C14
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -11214; // cut angle for A side positive BMG1A14
        else
          Nsteps = -12212; // cut angle for C side negative BMG1C12
      } else if (vcutdef[0]->icut == 2) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -21214; // cut angle for A side positive BMG2A14
        else
          Nsteps = -22212; // cut angle for C side negative BMG2C12
      } else if (vcutdef[0]->icut == 3) {
        if (vcutdef[0]->dead1 > 0.)
          Nsteps = -31214; // cut angle for A side positive BMG3A14
        else
          Nsteps = -32212; // cut angle for C side negative BMG3C12
      } else {
        log << MSG::ERROR << "massive error with MDT10 (BMG chambers)" << endmsg;
        std::abort();
      }
    }
    
    if (logVolName.find("MDT10") != std::string::npos || logVolName.find("MDT11") != std::string::npos || logVolName.find("MDT12") != std::string::npos ||
        logVolName.find("MDT13") != std::string::npos) {
      
      for (int i = 0; i < 5; i++) {
        cutoutNtubes[i] = 0;
        cutoutFullLength[i] = true;
        cutoutXtubes[i] = 0.;
        cutoutTubeLength[i] = width;
        cutoutYmax[i] = 0.;
      }
    }
    
    // Pass information to multilayer and MdtComponent
    layer->cutoutNsteps = Nsteps;
    m_component->cutoutTubeXShift = 0.;
    for (int i = 0; i < 5; i++) {
      layer->cutoutNtubes[i] = cutoutNtubes[i];
      layer->cutoutTubeLength[i] = cutoutTubeLength[i];
      layer->cutoutFullLength[i] = cutoutFullLength[i];
      layer->cutoutXtubes[i] = cutoutXtubes[i];
      layer->cutoutYmax[i] = cutoutYmax[i];
      if (!cutoutFullLength[i])
        m_component->cutoutTubeXShift = cutoutXtubes[i];
      // For now assume multiple cutouts have same width and take only the last value
    }
    
    layer->cutoutAtAngle = cutAtAngle;
      } else {
    // there are several cutouts along the amdb-x coordinate
    
    if (longWidth != width) {
      throw std::runtime_error(Form("File: %s, Line: %d\nMdt::build() - only support cutouts along amdb-x for rectangular chambers", __FILE__, __LINE__));
    }
    
    std::vector<std::pair<double, double>> nonCutoutXSteps;
    std::vector<std::pair<double, double>> nonCutoutYSteps;
    
    // Order cutouts by increasing dx
    std::sort(vcutdef.begin(),vcutdef.end(),[](const Cutout *a ,const Cutout* b ){
        return a->dx < b->dx;
      });
    // in amdb-coordinates
    double xminChamber = round(-width / 2, 2);
    double xmaxChamber = round(width / 2, 2);
    double yminChamber = 0;
    double ymaxChamber = round(length - tubePitch / 2, 2);
    
    double latestXMax = xminChamber;
    
    for (int i = 0; i < Ncuts; ++i) {
      Cutout *c = vcutdef[i];
      double lowerX = round(c->dx - c->widthXs / 2, 2);
      double xmin = std::max(lowerX, xminChamber);
      log<<MSG::DEBUG<<__FILE__<<":"<<__LINE__
         <<" name: "<<name
         <<" logVolName: "<<logVolName
         <<" xminChamber: "<<xminChamber<<" xmaxChamber: "<<xmaxChamber
         <<" ymaxChamber: "<<ymaxChamber
         <<" yminChamber: "<<yminChamber
         <<" latestXMax: "<<latestXMax
         <<" c->dx: "<<c->dx
         <<" c->widthXs/2 "<<c->widthXs / 2<<endmsg;
      if (xmin < latestXMax) {
        throw std::runtime_error(Form("File: %s, Line: %d\nMdt::build() - cannot have cutouts along amdb-x which overlap in amdb-x %f > %f",
                      __FILE__, __LINE__, latestXMax, xmin));
      }
      if (i == 0 && xmin > xminChamber) {
        // we start with a full slice without cutout
        nonCutoutXSteps.emplace_back(xminChamber, lowerX);
        nonCutoutYSteps.emplace_back(yminChamber, ymaxChamber);
      }
      double upperX = round(c->dx + c->widthXs / 2, 2);
      double xmax = (upperX >= xmaxChamber) ? xmaxChamber : upperX;
      
      double ymin = round(c->dy + c->lengthY, 2) < ymaxChamber ? c->dy + c->lengthY : 0;
      double ymax = ymaxChamber <= round(c->dy + c->lengthY, 2) ? c->dy : ymaxChamber;
      
      if (latestXMax < xmin) {
        // there is a full slice between latestXMax and xmin
        nonCutoutXSteps.emplace_back(latestXMax, xmin);
        nonCutoutYSteps.emplace_back(yminChamber, ymaxChamber);
      }
      
      nonCutoutXSteps.emplace_back(xmin, xmax);
      nonCutoutYSteps.emplace_back(ymin, ymax);
      
      if (i == Ncuts - 1 && xmax < xmaxChamber) {
        // we end with a full slice without cutout
        nonCutoutXSteps.emplace_back(xmax, xmaxChamber);
        nonCutoutYSteps.emplace_back(yminChamber, ymaxChamber);
      }
      
      latestXMax = xmax;
    }
    
    // Pass information to multilayer and MdtComponent
    m_component->cutoutTubeXShift = 0;
    layer->cutoutNsteps = nonCutoutXSteps.size();
    layer->m_nonCutoutXSteps = nonCutoutXSteps;
    layer->m_nonCutoutYSteps = nonCutoutYSteps;
    layer->cutoutAtAngle = false;
      }
    }

setLogVolName()

void MuonGM::DetectorElement::setLogVolName ( const std::string & str )

inlineinherited

Definition at line 22 of file DetectorElement.h.

{ logVolName = str; }

Member Data Documentation

index

int MuonGM::Mdt::index {0}

Definition at line 30 of file Mdt.h.

{0};

layer

std::unique_ptr<MultiLayer> MuonGM::Mdt::layer {}

Definition at line 34 of file Mdt.h.

{};

length

double MuonGM::Mdt::length {0.}

Definition at line 27 of file Mdt.h.

{0.};

logVolName

std::string MuonGM::DetectorElement::logVolName {}

inherited

Definition at line 18 of file DetectorElement.h.

{};

longWidth

double MuonGM::Mdt::longWidth {0.}

Definition at line 29 of file Mdt.h.

{0.}; // for trapezoidal layers

m_component

MdtComponent* MuonGM::Mdt::m_component {nullptr}

private

Definition at line 47 of file Mdt.h.

{nullptr};

name

std::string MuonGM::DetectorElement::name {}

inherited

Definition at line 17 of file DetectorElement.h.

{};

thickness

double MuonGM::Mdt::thickness {0.}

Definition at line 28 of file Mdt.h.

{0.};

tubelenStepSize

double MuonGM::Mdt::tubelenStepSize {0.}

Definition at line 31 of file Mdt.h.

{0.};

tubePitch

double MuonGM::Mdt::tubePitch {0.}

Definition at line 32 of file Mdt.h.

{0.};

width

double MuonGM::Mdt::width {0.}

Definition at line 26 of file Mdt.h.

{0.};

---

The documentation for this class was generated from the following files:

• Mdt.h
• Mdt.cxx