CocoaToDDLMgr Class Reference

#include <CocoaToDDLMgr.h>

Public Member Functions
ALIint	buildRotationNumber (OpticalObject *opto)
	CocoaToDDLMgr ()
void	lv (OpticalObject *opto)
void	ma (CocoaMaterialElementary *ma)
ALIbool	materialIsRepeated (CocoaMaterialElementary *ma)
void	measurementsAsSpecPars ()
void	newPartPost (std::string name, std::string extension)
void	newPartPre (std::string name)
void	newSectPost (std::string name)
void	newSectPost_lv (std::string name)
void	newSectPost_ma (std::string name)
void	newSectPost_pv (std::string name)
void	newSectPost_ro (std::string name)
void	newSectPost_so (std::string name)
void	newSectPost_specPar (std::string name)
void	newSectPre (std::string name, std::string type)
void	newSectPre_lv (std::string name)
void	newSectPre_ma (std::string name)
void	newSectPre_pv (std::string name)
void	newSectPre_ro (std::string name)
void	newSectPre_so (std::string name)
void	newSectPre_specPar (std::string name)
void	pv (OpticalObject *opto)
void	ro (const CLHEP::HepRotation &ro, int n)
std::string	scrubString (const std::string &s)
void	so (OpticalObject *opto)
void	specPar (OpticalObject *opto)
void	writeDDDFile (ALIstring filename)
void	writeHeader (ALIstring filename)
void	writeLogicalVolumes ()
void	writeMaterials ()
void	writePhysicalVolumes ()
void	writeRotations ()
void	writeSolids ()
void	writeSpecPars ()
void	writeSpecParsCocoa ()
	~CocoaToDDLMgr ()

Static Public Member Functions
static CocoaToDDLMgr *	getInstance ()

Private Attributes
ALIFileOut	file_
std::string	filename_
std::vector< CocoaMaterialElementary * >	theMaterialList
std::vector< CLHEP::HepRotation >	theRotationList

Static Private Attributes
static CocoaToDDLMgr *	instance = nullptr

Detailed Description

Definition at line 25 of file CocoaToDDLMgr.h.

Constructor & Destructor Documentation

CocoaToDDLMgr::CocoaToDDLMgr ( )

inline

Definition at line 28 of file CocoaToDDLMgr.h.

Referenced by getInstance().

{};

CocoaToDDLMgr::~CocoaToDDLMgr ( )

inline

Definition at line 29 of file CocoaToDDLMgr.h.

References buildRotationNumber(), edmStreamStallGrapher::extension, corrVsCorr::filename, getInstance(), lv(), ma(), materialIsRepeated(), measurementsAsSpecPars(), dqmiodumpmetadata::n, Skims_PA_cff::name, newPartPost(), newPartPre(), newSectPost(), newSectPost_lv(), newSectPost_ma(), newSectPost_pv(), newSectPost_ro(), newSectPost_so(), newSectPost_specPar(), newSectPre(), newSectPre_lv(), newSectPre_ma(), newSectPre_pv(), newSectPre_ro(), newSectPre_so(), newSectPre_specPar(), pv(), ro(), alignCSCRings::s, scrubString(), so(), specPar(), AlCaHLTBitMon_QueryRunRegistry::string, writeDDDFile(), writeHeader(), writeLogicalVolumes(), writeMaterials(), writePhysicalVolumes(), writeRotations(), writeSolids(), writeSpecPars(), and writeSpecParsCocoa().

{};

Member Function Documentation

ALIint CocoaToDDLMgr::buildRotationNumber

(

OpticalObject *

opto

)

Definition at line 741 of file CocoaToDDLMgr.cc.

References OpticalObject::rmLocal(), and theRotationList.

Referenced by pv(), and ~CocoaToDDLMgr().

                                                             {
  ALIint rotnum = -1;

  if (opto->rmLocal().isIdentity())
    return rotnum;

  std::vector<CLHEP::HepRotation>::const_iterator ite;

  int nc = 0;
  for (ite = theRotationList.begin(); ite != theRotationList.end(); ++ite) {
    if ((*ite) == opto->rmLocal()) {
      rotnum = nc;
      break;
    }
    nc++;
  }

  if (rotnum == -1) {
    theRotationList.push_back(opto->rmLocal());
    rotnum = theRotationList.size() - 1;
  }

  return rotnum;
}

CocoaToDDLMgr * CocoaToDDLMgr::getInstance ( )

static

Definition at line 23 of file CocoaToDDLMgr.cc.

References CocoaToDDLMgr(), and instance.

Referenced by ~CocoaToDDLMgr().

                                          {
  if (!instance) {
    instance = new CocoaToDDLMgr;
  }
  return instance;
}

void CocoaToDDLMgr::lv

(

OpticalObject *

opto

)

Definition at line 425 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, file_, OpticalObject::getMaterial(), CocoaMaterialElementary::getName(), Skims_PA_cff::name, OpticalObject::name(), OpticalObject::shortName(), AlCaHLTBitMon_QueryRunRegistry::string, and OpticalObject::type().

Referenced by writeLogicalVolumes(), and ~CocoaToDDLMgr().

                                          {
  std::string name = opto->shortName();
  std::string rSolid = opto->shortName();
  std::string sensitive = "unspecified";

  if (opto->type() == "system") {
    file_ << " <LogicalPart name=\"" << name << "\" category=\"" << sensitive << "\">" << std::endl
          << "  <rSolid name=\"" << rSolid << "\"/>" << std::endl
          << "  <rMaterial name=\"Hydrogen\""
          << "/>" << std::endl
          << " </LogicalPart>" << std::endl;
    return;
  }

#ifdef gdebug_v
  cout << "xml:lv " << opto->name() << std::endl;
#endif
  file_ << " <LogicalPart name=\"" << name << "\" category=\"" << sensitive << "\">" << std::endl
        << "  <rSolid name=\"" << rSolid << "\"/>" << std::endl
        << "  <rMaterial name=\"" << opto->getMaterial()->getName() << "\""
        << "/>" << std::endl
        << " </LogicalPart>" << std::endl;
}

void CocoaToDDLMgr::ma

(

CocoaMaterialElementary *

ma

)

g*cm3;

Definition at line 198 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, fastSimProducer_cff::density, file_, CocoaMaterialElementary::getA(), CocoaMaterialElementary::getDensity(), CocoaMaterialElementary::getName(), CocoaMaterialElementary::getSymbol(), CocoaMaterialElementary::getZ(), theMaterialList, and UC.

Referenced by materialIsRepeated(), writeMaterials(), and ~CocoaToDDLMgr().

                                                  {
  theMaterialList.push_back(ma);

#ifdef gdebug
  cout << "  ma:" << ma->getName() << std::endl;
#endif

  ALIfloat density = ma->getDensity();  

  // start tag
  file_ << "  <ElementaryMaterial";
  ALIstring cSymbol = ma->getSymbol();

  // name attribute
  file_ << " name=\"" << ma->getName() << "\"";

  // put out common attributes
  //  file_ << " density=\"" << UnitConverter(density,"Volumic Mass") << "\"";
  file_ << " density=\"" << UC(density, "Volumic Mass") << "\"";
  file_ << " symbol=\"" << ma->getSymbol() << "\"";

  // finish last two attributes and end material element
  file_ << " atomicWeight=\"" << (ma->getA()) << "*g/mole\""
        << " atomicNumber=\"" << ma->getZ() << "\""
        << "/>" << std::endl;
}

ALIbool CocoaToDDLMgr::materialIsRepeated

(

CocoaMaterialElementary *

ma

)

Definition at line 698 of file CocoaToDDLMgr.cc.

References ma(), and theMaterialList.

Referenced by writeMaterials(), and ~CocoaToDDLMgr().

                                                                     {
  ALIbool isRepeated = false;
  std::vector<CocoaMaterialElementary*>::const_iterator ite;

  for (ite = theMaterialList.begin(); ite != theMaterialList.end(); ++ite) {
    if (*(*ite) == *ma) {
      isRepeated = true;
      break;
    }
  }

  return isRepeated;
}

void CocoaToDDLMgr::measurementsAsSpecPars

(

)

Definition at line 617 of file CocoaToDDLMgr.cc.

References Measurement::dim(), file_, cuy::ii, Model::MeasurementList(), OpticalObject::name(), Measurement::name(), Measurement::OptONameList(), Measurement::sigma(), AlCaHLTBitMon_QueryRunRegistry::string, Measurement::type(), Measurement::value(), Measurement::valueIsSimulated(), and Measurement::valueType().

Referenced by writeSpecPars(), and ~CocoaToDDLMgr().

                                           {
  std::vector<Measurement*> measlist = Model::MeasurementList();
  std::vector<Measurement*>::iterator mite;
  std::vector<ALIstring>::iterator site;
  std::multimap<OpticalObject*, Measurement*> optoMeasMap;
  for (mite = measlist.begin(); mite != measlist.end(); ++mite) {
    auto& optolist = (*mite)->OptOList();
    OpticalObject* opto = optolist[optolist.size() - 1];
    optoMeasMap.insert(std::multimap<OpticalObject*, Measurement*>::value_type(opto, *mite));
  }

  typedef std::multimap<OpticalObject*, Measurement*>::const_iterator itemom;
  itemom omite;
  std::pair<itemom, itemom> omitep;
  itemom omite2, omite3;

  for (omite = optoMeasMap.begin(); omite != optoMeasMap.end(); ++omite) {
    omitep = optoMeasMap.equal_range((*omite).first);
    if (omite != optoMeasMap.begin() && (*omite).first == (*omite3).first)
      continue;  // check that it is not the same OptO than previous one
    omite3 = omite;
    for (omite2 = omitep.first; omite2 != omitep.second; ++omite2) {
      OpticalObject* opto = (*(omite2)).first;
      Measurement* meas = (*(omite2)).second;
      std::vector<ALIstring> namelist = meas->OptONameList();
      if (omite2 == omitep.first) {
        file_ << " <SpecPar name=\"" << meas->name() << "_MEASUREMENT\">" << std::endl;
        file_ << "   <PartSelector path=\"/" << opto->name() << "\"/> " << std::endl;
      }

      file_ << "   <Parameter name=\"" << std::string("meas_name") << "\" value=\"" << meas->name()
            << "\"  eval=\"false\" /> " << std::endl;
      file_ << "   <Parameter name=\"" << std::string("meas_type") << "\" value=\"" << meas->type()
            << "\"  eval=\"false\" /> " << std::endl;
      for (site = namelist.begin(); site != namelist.end(); ++site) {
        file_ << "   <Parameter name=\"" << std::string("meas_object_name_") + meas->name() << "\" value=\"" << (*site)
              << "\"  eval=\"false\" /> " << std::endl;
      }
      for (ALIuint ii = 0; ii < meas->dim(); ii++) {
        file_ << "   <Parameter name=\"" << std::string("meas_value_name_") + meas->name() << "\" value=\""
              << meas->valueType(ii) << "\"  eval=\"false\" /> " << std::endl;
        file_ << "   <Parameter name=\"" << std::string("meas_value_") + meas->name() << "\" value=\""
              << meas->value(ii) << "\"  eval=\"true\" /> " << std::endl;
        file_ << "   <Parameter name=\"" << std::string("meas_sigma_") + meas->name() << "\" value=\""
              << meas->sigma(ii) << "\"  eval=\"true\" /> " << std::endl;
        file_ << "   <Parameter name=\"" << std::string("meas_is_simulated_value_") + meas->name() << "\" value=\""
              << meas->valueIsSimulated(ii) << "\"  eval=\"true\" /> " << std::endl;
      }
    }
    file_ << " </SpecPar>" << std::endl;
  }
}

void CocoaToDDLMgr::newPartPost	(	std::string	name,
		std::string	extension
	)

Definition at line 183 of file CocoaToDDLMgr.cc.

References file_.

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

                                                                   {
  file_ << std::endl << "</DDDefinition>" << std::endl;
  file_.close();
}

void CocoaToDDLMgr::newPartPre

(

std::string

name

)

Definition at line 165 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, file_, filename_, and Skims_PA_cff::name.

Referenced by writeHeader(), and ~CocoaToDDLMgr().

                                             {
  filename_ = name;
  file_.open(filename_.c_str());
  file_.precision(8);
  file_ << "<?xml version=\"1.0\"?>" << std::endl;

  // all files get schema references and namespaces.
  file_ << "<DDDefinition xmlns=\"http://www.cern.ch/cms/DDL\""
        << " xmlns:xsi=\"http://www.w3.org/2001/XMLSchema-instance\""
        << " xsi:schemaLocation=\"http://www.cern.ch/cms/DDL ../../DDLSchema/DDLSchema.xsd\">" << std::endl
        << std::endl;

#ifdef gdebug
  cout << "part-pre:" << name << std::endl;
#endif
}

void CocoaToDDLMgr::newSectPost

(

std::string

name

)

Definition at line 695 of file CocoaToDDLMgr.cc.

References file_.

Referenced by newSectPost_lv(), newSectPost_ma(), newSectPost_pv(), newSectPost_ro(), newSectPost_so(), newSectPost_specPar(), and ~CocoaToDDLMgr().

{ file_ << "</" << name << ">" << std::endl << std::endl; }

void CocoaToDDLMgr::newSectPost_lv

(

std::string

name

)

Definition at line 450 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, and newSectPost().

Referenced by writeLogicalVolumes(), and ~CocoaToDDLMgr().

                                                 {
#ifdef gdebug
  cout << " sect-lv-post:" << name << '-' << std::endl;
#endif
  newSectPost("LogicalPartSection");
}

void CocoaToDDLMgr::newSectPost_ma

(

std::string

name

)

Definition at line 226 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, and newSectPost().

Referenced by writeMaterials(), and ~CocoaToDDLMgr().

                                                 {
#ifdef gdebug
  cout << " sect-mat-post:" << name << '-' << std::endl;
#endif
  newSectPost("MaterialSection");
}

void CocoaToDDLMgr::newSectPost_pv

(

std::string

name

)

Definition at line 506 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, and newSectPost().

Referenced by writePhysicalVolumes(), and ~CocoaToDDLMgr().

                                                 {
#ifdef gdebug
  cout << " sect-pv-post:" << name << '-' << std::endl;
#endif
  newSectPost("PosPartSection");
}

void CocoaToDDLMgr::newSectPost_ro

(

std::string

name

)

Definition at line 544 of file CocoaToDDLMgr.cc.

References newSectPost().

Referenced by writeRotations(), and ~CocoaToDDLMgr().

{ newSectPost("RotationSection"); }

void CocoaToDDLMgr::newSectPost_so

(

std::string

name

)

Definition at line 408 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, and newSectPost().

Referenced by writeSolids(), and ~CocoaToDDLMgr().

                                                 {
#ifdef gdebug
  cout << " sect-so-post:" << name << '-' << std::endl;
#endif
  newSectPost("SolidSection");
}

void CocoaToDDLMgr::newSectPost_specPar

(

std::string

name

)

Definition at line 685 of file CocoaToDDLMgr.cc.

References newSectPost().

Referenced by writeSpecPars(), and ~CocoaToDDLMgr().

{ newSectPost("SpecParSection"); }

void CocoaToDDLMgr::newSectPre	(	std::string	name,
		std::string	type
	)

Definition at line 690 of file CocoaToDDLMgr.cc.

References file_.

Referenced by newSectPre_lv(), newSectPre_ma(), newSectPre_pv(), newSectPre_ro(), newSectPre_so(), and ~CocoaToDDLMgr().

                                                             {
  file_ << "<" << type << " label=\"" << name << "\">" << std::endl;
}

void CocoaToDDLMgr::newSectPre_lv

(

std::string

name

)

Definition at line 417 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, filename_, newSectPre(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by writeLogicalVolumes(), and ~CocoaToDDLMgr().

                                                {
#ifdef gdebug
  cout << " sect-lv-pre:" << name << '-' << std::endl;
#endif
  newSectPre(filename_, std::string("LogicalPartSection"));
}

void CocoaToDDLMgr::newSectPre_ma

(

std::string

name

)

Definition at line 190 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, filename_, newSectPre(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by writeMaterials(), and ~CocoaToDDLMgr().

                                                {
#ifdef gdebug
  std::cout << " sect-mat-pre:" << name << '-' << std::endl;
#endif
  newSectPre(filename_, std::string("MaterialSection"));
}

void CocoaToDDLMgr::newSectPre_pv

(

std::string

name

)

Definition at line 459 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, filename_, newSectPre(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by writePhysicalVolumes(), and ~CocoaToDDLMgr().

                                                {
#ifdef gdebug
  cout << " sect-pv-pre:" << name << '-' << std::endl;
#endif
  newSectPre(filename_, std::string("PosPartSection"));
}

void CocoaToDDLMgr::newSectPre_ro

(

std::string

name

)

Definition at line 515 of file CocoaToDDLMgr.cc.

References filename_, newSectPre(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by writeRotations(), and ~CocoaToDDLMgr().

{ newSectPre(filename_, std::string("RotationSection")); }

void CocoaToDDLMgr::newSectPre_so

(

std::string

name

)

Definition at line 235 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, filename_, newSectPre(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by writeSolids(), and ~CocoaToDDLMgr().

                                                {
#ifdef gdebug
  cout << " sect-so-pre:" << name << '-' << std::endl;
#endif
  newSectPre(filename_, std::string("SolidSection"));
}

void CocoaToDDLMgr::newSectPre_specPar

(

std::string

name

)

Definition at line 548 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, file_, and filename_.

Referenced by writeSpecPars(), and ~CocoaToDDLMgr().

                                                     {
#ifdef gdebug
  cout << " sect-lv-pre:" << name << '-' << std::endl;
#endif
  //-   newSectPre(filename_,std::string("SpecParSection"));
  file_ << "<SpecParSection label=\"" << filename_ << "\" eval=\"true\">" << std::endl;
}

void CocoaToDDLMgr::pv

(

OpticalObject *

opto

)

Definition at line 467 of file CocoaToDDLMgr.cc.

References buildRotationNumber(), OpticalObject::centreLocal(), gather_cfg::cout, file_, OpticalObject::name(), OpticalObject::parent(), OpticalObject::shortName(), OrderedSet::t, and UC.

Referenced by writePhysicalVolumes(), and ~CocoaToDDLMgr().

                                          {
#ifdef gdebug_v
  cout << "  pv:" << opto->name() << ':' << opto->parent()->name() << std::endl;
#endif

  //   file_ << " <PosPart copyNumber=\"" << pv->GetCopyNo() << "\">" << std::endl;
  file_ << " <PosPart copyNumber=\""
        << "1"
        << "\">" << std::endl;
  file_ << "   <rParent name=\"";

  //t   if (file!=filename_) file_ << file << ":";

  file_ << opto->parent()->shortName() << "\"/>" << std::endl;

  file_ << "   <rChild name=\"";
  //t  if (file_d != filename_)  file_<<  file_d << ":";
  file_ << opto->shortName();
  file_ << "\"/>" << std::endl;

  int rotNumber = buildRotationNumber(opto);
  //CocoaDDLRotation* rot = buildRotationNotRepeated( opto );

  if (rotNumber != -1)
    file_ << "  <rRotation name=\"R" << rotNumber << "\"/>" << std::endl;

  CLHEP::Hep3Vector t = opto->centreLocal();
  if (t != CLHEP::Hep3Vector()) {  //if (0,0,0) write nothing
    const CLHEP::Hep3Vector t = opto->centreLocal();

    file_ << "  <Translation x=\"" << UC(t[0], "Length") << "\""
          << " y=\"" << UC(t[1], "Length") << "\""
          << " z=\"" << UC(t[2], "Length") << "\" />" << std::endl;
  }

  file_ << " </PosPart>" << std::endl;
}

void CocoaToDDLMgr::ro	(	const CLHEP::HepRotation &	ro,
		int	n
	)

Definition at line 519 of file CocoaToDDLMgr.cc.

References file_, identity(), GlobalPosition_Frontier_DevDB_cff::tag, and UC.

Referenced by writeRotations(), and ~CocoaToDDLMgr().

                                                        {
  CLHEP::HepRotation roinv = inverseOf(ro);
  //-    G4ThreeVector v(1.,1.,1.);
  //-     G4ThreeVector a;
  //-   a = (*ro)*v;
  bool identity = false;
  ALIstring tag = " <Rotation name=\"R";
  identity = roinv.isIdentity();

  //---- DDD convention is to use the inverse matrix, COCOA is the direct one!!!
  if (!identity) {
    file_ << tag << n << "\"";
    file_ << " phiX=\"" << UC(roinv.phiX(), "Angle") << "\""
          << " thetaX=\"" << UC(roinv.thetaX(), "Angle") << "\""
          << " phiY=\"" << UC(roinv.phiY(), "Angle") << "\""
          << " thetaY=\"" << UC(roinv.thetaY(), "Angle") << "\""
          << " phiZ=\"" << UC(roinv.phiZ(), "Angle") << "\""
          << " thetaZ=\"" << UC(roinv.thetaZ(), "Angle")
          << "\""
          //<< " angleUnit=\"degree\"/>"
          << " />" << std::endl;
  }
}

std::string CocoaToDDLMgr::scrubString

(

const std::string &

s

)

Definition at line 713 of file CocoaToDDLMgr.cc.

References runTheMatrix::ret, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~CocoaToDDLMgr().

                                                       {
  std::string::const_iterator ampat;
  static const std::string amp = "_";  //"&";
  std::string ret = "";
  for (ampat = s.begin(); ampat != s.end(); ++ampat) {
    if (*ampat == '&')
      ret = ret + amp;
    else if (*ampat == '/')
      ret = ret + ";";
    else if (*ampat == ':')
      ret = ret + '_';
    else
      ret = ret + *ampat;
  }
  // this works when used alone.  when in this file it fails.  i don't know why.
  //for (ampat = s.begin(); ampat != s.end(); ampat++)
  //  {
  //    if (*ampat == '&'){
  //    s.insert(ampat+1, amp.begin(), amp.end());
  //    }
  //  }
  //  replace(s.begin(), s.end(), '/', ';');
  //return s;
  //cout << "AMP: " << ret << endl;
  return ret;
}

void CocoaToDDLMgr::so

(

OpticalObject *

opto

)

Definition at line 243 of file CocoaToDDLMgr.cc.

References beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, cppFunctionSkipper::exception, file_, CocoaSolidShapeTubs::getDeltaPhiAngle(), CocoaSolidShapeTubs::getInnerRadius(), CocoaSolidShapeTubs::getOuterRadius(), CocoaSolidShapeTubs::getStartPhiAngle(), CocoaSolidShape::getType(), CocoaSolidShapeBox::getXHalfLength(), CocoaSolidShapeBox::getYHalfLength(), CocoaSolidShapeBox::getZHalfLength(), CocoaSolidShapeTubs::getZHalfLength(), Skims_PA_cff::name, ALIFileOut::name(), AlCaHLTBitMon_QueryRunRegistry::string, and UC.

Referenced by writeSolids(), and ~CocoaToDDLMgr().

                                          {
  std::cout << " CocoaToDDLMgr::so( " << opto;
  std::cout << " " << opto->shortName() << std::endl;

  std::string name = opto->shortName();

  if (opto->type() == "system") {
    //    file_ << " <Box name=\"" << name << "\"";
    file_ << " <Box name=\"" << opto->name() << "\"";
    file_ << " dx=\"10.*m"
          << "\" dy=\"10.*m"
          << "\" dz=\"10.*m"
          << "\"/>" << std::endl;
    return;
  }

  CocoaSolidShape* so = opto->getSolidShape();

  std::cout << " CocoaToDDLMgr::so( so " << so << std::endl;
  std::string solidType = so->getType();

  if (solidType == "Box") {
    file_ << " <" << solidType << " name=\"" << name << "\"";
    CocoaSolidShapeBox* sb = dynamic_cast<CocoaSolidShapeBox*>(so);
    file_ << " dx=\"" << UC(sb->getXHalfLength(), "Length") << "\" dy=\"" << UC(sb->getYHalfLength(), "Length")
          << "\" dz=\"" << UC(sb->getZHalfLength(), "Length") << "\"/>" << std::endl;
  } else if (solidType == "Tubs") {
    CocoaSolidShapeTubs* tu = dynamic_cast<CocoaSolidShapeTubs*>(so);
    file_ << " <" << solidType << " name=\"" << name << "\""
          << " rMin=\"" << UC(tu->getInnerRadius(), "Length") << "\""
          << " rMax=\"" << UC(tu->getOuterRadius(), "Length") << "\""
          << " dz=\"" << UC(tu->getZHalfLength(), "Length") << "\""
          << " startPhi=\"" << UC(tu->getStartPhiAngle(), "Angle") << "\""
          << " deltaPhi=\"" << UC(tu->getDeltaPhiAngle(), "Angle") << "\""
          << "/>" << std::endl;
  }
  /*  else if (solidType == "Cons")
    {
      G4Cons * cn = dynamic_cast < G4Cons * > (so);
      file_ << " <" << solidType 
            << " name=\""     << name                                              << "\""
        << " dz=\""       << UC(cn->getZHalfLength(),"Length")       << "\""
        << " rMin1=\""    << UC(cn->getInnerRadiusMinusZ(),"Length") << "\"" 
        << " rMax1=\""    << UC(cn->getOuterRadiusMinusZ(),"Length") << "\""
        << " rMin2=\""    << UC(cn->getInnerRadiusPlusZ(),"Length")  << "\""        
        << " rMax2=\""    << UC(cn->getOuterRadiusPlusZ(),"Length")  << "\""
            << " startPhi=\"" << UC(cn->getStartPhiAngle(),"Angle")     << "\""
        << " deltaPhi=\"" << UC(cn->getDeltaPhiAngle(),"Angle")    << "\""
      //<< " lengthUnit=\"mm\" angleUnit=\"degree\"/>" 
        << " />" << std::endl;
    }
  else if (solidType == "Polycone")
    {
      G4Polycone * pc = dynamic_cast < G4Polycone * > (so);
      file_ << " <Polycone name=\"" << name<< "\"";
      bool isOpen = pc->IsOpen();
      G4int numRZCorner = (dynamic_cast < G4Polycone * > (so))->getNumRZCorner();

      file_ << " startPhi=\"" << UC(pc->getStartPhi(),"Angle") << "\""
        //<< " deltaPhi=\"" << UC(fabs((pc->getEndPhi()/deg - pc->getStartPhi()/deg))*deg,"Angle")   << "\"" 
        //<< " deltaPhi=\"" << UC(pc->getEndPhi(),"Angle")   << "\"" 
        << " deltaPhi=\"" << UC(pc->original_parameters->Opening_angle,"Angle")   << "\""
        //<< " angleUnit=\"degree\">" 
        << " >" << std::endl;

      G4PolyconeSideRZ rz;

      //liendl: FIXME put a switch which decides whether RZ or Rmin,Rmax,Z types should
      //              by generated ....
      //outPolyZSections(rz, (dynamic_cast < G4Polycone * > (so)), numRZCorner);
      G4double * zVal;
      G4double * rmin;
      G4double * rmax;
      G4int zPlanes;
      zPlanes = pc->original_parameters->Num_z_planes;
      zVal = pc->original_parameters->Z_values;
      rmin = pc->original_parameters->Rmin;
      rmax = pc->original_parameters->Rmax;
      outPolySections(zPlanes, zVal, rmin, rmax);
      file_ << " </Polycone> " << std::endl;
    }
  else if (solidType == "Polyhedra")
    {
      //      bool isOpen = (dynamic_cast < G4Polyhedra * > (so))->IsOpen();
      G4Polyhedra * ph = (dynamic_cast < G4Polyhedra * > (so));
      G4int numRZCorner = ph->getNumRZCorner();

      file_ << " <Polyhedra name=\"" << name<< "\""
            << " numSide=\"" << ph->getNumSide() << "\""
            << " startPhi=\"" << UC(ph->getStartPhi(),"Angle") << "\""
        //<< " deltaPhi=\""   <<  UC(fabs((ph->getEndPhi()/deg - ph->getStartPhi()/deg))*deg,"Angle")   << "\""
        << " deltaPhi=\"" << UC(ph->original_parameters->Opening_angle,"Angle")   << "\"" 
        << " >" << std::endl;

      G4PolyhedraSideRZ rz;
      //liendl: FIXME put a switch which decides whether RZ or Rmin,Rmax,Z types should
      //              by generated ....
      // outPolyZSections(rz, (dynamic_cast < G4Polyhedra * > (so)), numRZCorner);
      G4double * zVal;
      G4double * rmin;
      G4double * rmax;
      // convertRad of ctor of G4Polyhedra(..) ....
      G4double strangeG3G4Factor = cos(0.5*ph->original_parameters->Opening_angle/G4double(ph->getNumSide()));
      G4int zPlanes;
      zPlanes = ph->original_parameters->Num_z_planes;
      zVal = ph->original_parameters->Z_values;
      rmin = ph->original_parameters->Rmin;
      rmax = ph->original_parameters->Rmax;
      for (int i=0; i<zPlanes;++i) {
       *(rmin+i) = *(rmin+i) * strangeG3G4Factor;
       *(rmax+i) = *(rmax+i) * strangeG3G4Factor;
      }
      outPolySections(zPlanes, zVal, rmin, rmax);

      file_ << " </Polyhedra>" << std::endl;
    }
  else if (solidType == "Trapezoid")
    {
      // DDL fields
      // ALP1, ALP2, Bl1, Bl2, Dz, H1, H2, Phi, Thet, TL1, TL2, lengthUnit, angleUnit    
      // Phi and Theta are !NOT!optional.
       G4Trap * trp = dynamic_cast < G4Trap * > (so);
       G4ThreeVector symAxis(trp->getSymAxis());
       double theta, phi;
       theta = symAxis.theta();
       phi = symAxis.phi();

      file_ << " <" << solidType 
        << " name=\"" << name                                                       << "\""
        << " dz=\""   << UC(trp->getZHalfLength(),"Length")      << "\""
        << " alp1=\"" << UC(atan(trp->getTanAlpha1()/rad),"Angle") << "\"" 
        << " bl1=\""  << UC(trp->getXHalfLength1(),"Length")     << "\""
        << " tl1=\""  << UC(trp->getXHalfLength2(),"Length")     << "\""        
        << " h1=\""   << UC(trp->getYHalfLength1(),"Length")    << "\""
        << " alp2=\"" << UC(atan(trp->getTanAlpha2()/rad),"Angle") << "\""
        << " bl2=\""  << UC(trp->getXHalfLength3(),"Length")     << "\""
        << " tl2=\""  << UC(trp->getXHalfLength4(),"Length")    << "\""     
        << " h2=\""   << UC(trp->getYHalfLength2(),"Length")     << "\""
        << " phi=\""  << UC(phi,"Angle") << "\""
        << " theta=\"" << UC(theta,"Angle") << "\"" 
        << " />" << std::endl ;
    }
  else if (solidType == "Trd1")
    {
      G4Trd * tr = dynamic_cast < G4Trd * > (so);
      file_ << " <" << solidType  
        << " name=\"" << name                                                  << "\""
            << " dz=\""   << UC(tr->getZHalfLength(),"Length")  << "\""
        << " dy1=\""  << UC(tr->getYHalfLength1(),"Length") << "\""
            << " dy2=\""  << UC(tr->getYHalfLength2(),"Length")<< "\"" 
        << " dx1=\""  << UC(tr->getXHalfLength1(),"Length") << "\""
        << " dx2=\""  << UC(tr->getXHalfLength2(),"Length") << "\""
        //<< " lengthUnit=\"mm\"/>" 
        << " />" << std::endl;
    }
  */
  else {
    std::cerr << " <!-- NOT HANDLED: " << solidType << " name=\"" << name << "\""
              << ">" << std::endl
              << " </" << solidType << "> -->" << std::endl;
    std::exception();
  }
}

void CocoaToDDLMgr::specPar

(

OpticalObject *

opto

)

Definition at line 557 of file CocoaToDDLMgr.cc.

References OpticalObject::CoordinateEntryList(), OpticalObject::ExtraEntryList(), file_, OpticalObject::getCmsswID(), cuy::ii, Entry::name(), OpticalObject::name(), Entry::quality(), Entry::sigma(), AlCaHLTBitMon_QueryRunRegistry::string, Entry::type(), OpticalObject::type(), UC, and Entry::value().

Referenced by writeSpecPars(), and ~CocoaToDDLMgr().

                                               {
  file_ << " <SpecPar name=\"" << opto->name() << "_PARAMS\">" << std::endl;
  file_ << "   <PartSelector path=\"/" << opto->name() << "\"/> " << std::endl;
  file_ << "   <Parameter name=\"cocoa_type\""
        << " value=\"" << opto->type() << "\"   eval=\"false\" /> " << std::endl;
  file_ << "   <Parameter name=\"cmssw_ID\""
        << " value=\"" << opto->getCmsswID() << "\" /> " << std::endl;

  const std::vector<Entry*>& coord = opto->CoordinateEntryList();
  for (int ii = 3; ii < 6; ii++) {
    Entry* ent = coord[ii];
    file_ << "   <Parameter name=\"" << ent->name() + std::string("_value") << "\" value=\"";
    file_ << UC(ent->value(), "Angle");
    file_ << "\" /> " << std::endl;
  }
  for (int ii = 0; ii < 6; ii++) {
    Entry* ent = coord[ii];
    file_ << "   <Parameter name=\"" << ent->name() + std::string("_sigma") << "\" value=\"";
    if (ii < 3) {
      file_ << UC(ent->sigma(), "Length");
    } else {
      file_ << UC(ent->sigma(), "Angle");
    }
    file_ << "\" /> " << std::endl;
    file_ << "   <Parameter name=\"" << ent->name() + std::string("_quality") << "\" value=\"" << ent->quality()
          << "\" /> " << std::endl;
  }

  const std::vector<Entry*>& extraEnt = opto->ExtraEntryList();
  for (ALIuint ii = 0; ii < extraEnt.size(); ii++) {
    Entry* ent = extraEnt[ii];
    file_ << "   <Parameter name=\"extra_entry\" value=\"" << ent->name() << "\"  eval=\"false\" /> " << std::endl;
    file_ << "   <Parameter name=\"dimType\" value=\"" << ent->type() << "\"  eval=\"false\" /> " << std::endl;
    file_ << "   <Parameter name=\"value\" value=\"";
    if (ent->type() == "nodim") {
      file_ << ent->value();
    } else if (ent->type() == "length") {
      file_ << UC(ent->value(), "Length");
    } else if (ent->type() == "angle") {
      file_ << UC(ent->value(), "Angle");
    }
    file_ << "\"  eval=\"true\" /> " << std::endl;

    file_ << "   <Parameter name=\"sigma\" value=\"";
    if (ent->type() == "nodim") {
      file_ << ent->sigma();
    } else if (ent->type() == "length") {
      file_ << UC(ent->sigma(), "Length");
    } else if (ent->type() == "angle") {
      file_ << UC(ent->sigma(), "Angle");
    }
    file_ << "\"  eval=\"true\" /> " << std::endl;

    file_ << "   <Parameter name=\"quality\" value=\"" << ent->quality() << "\"  eval=\"true\" /> " << std::endl;
  }

  file_ << " </SpecPar>" << std::endl;
}

void CocoaToDDLMgr::writeDDDFile

(

ALIstring

filename

)

Definition at line 31 of file CocoaToDDLMgr.cc.

References newPartPost(), writeHeader(), writeLogicalVolumes(), writeMaterials(), writePhysicalVolumes(), writeRotations(), writeSolids(), and writeSpecPars().

Referenced by ~CocoaToDDLMgr().

                                                   {
  //---- Write header
  writeHeader(filename);

  //---- Write materials
  writeMaterials();

  //---- Write solids
  writeSolids();

  //---- Write logical volumes
  writeLogicalVolumes();

  //---- Write physical volumes
  writePhysicalVolumes();

  //---- Write rotations
  writeRotations();

  //---- Write SpecPar's
  writeSpecPars();

  newPartPost(filename, "");
}

void CocoaToDDLMgr::writeHeader

(

ALIstring

filename

)

Definition at line 57 of file CocoaToDDLMgr.cc.

References newPartPre().

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

{ newPartPre(filename); }

void CocoaToDDLMgr::writeLogicalVolumes

(

)

Definition at line 99 of file CocoaToDDLMgr.cc.

References lv(), newSectPost_lv(), newSectPre_lv(), and Model::OptOList().

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

                                        {
  newSectPre_lv("");

  auto& optolist = Model::OptOList();
  for (auto ite = optolist.begin(); ite != optolist.end(); ite++) {
    bool alreadyWritten = false;
    for (auto ite2 = optolist.begin(); ite2 != ite; ite2++) {
      if ((*ite)->shortName() == (*ite2)->shortName()) {
        alreadyWritten = true;
      }
    }
    if (!alreadyWritten)
      lv(*ite);
  }

  newSectPost_lv("");
}

void CocoaToDDLMgr::writeMaterials

(

)

Definition at line 60 of file CocoaToDDLMgr.cc.

References ma(), materialIsRepeated(), newSectPost_ma(), newSectPre_ma(), and Model::OptOList().

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

                                   {
  newSectPre_ma("");
  auto& optolist = Model::OptOList();
  for (auto ite = optolist.begin(); ite != optolist.end(); ite++) {
    if ((*ite)->type() == "system")
      continue;
    CocoaMaterialElementary* mat = (*ite)->getMaterial();
    //-    std::cout << " mat of opto " << (*ite)->name() << " = " << mat->getName() << std::endl;
    if (mat) {
      if (!materialIsRepeated(mat))
        ma(mat);
    }
  }

  newSectPost_ma("");
}

void CocoaToDDLMgr::writePhysicalVolumes

(

)

Definition at line 118 of file CocoaToDDLMgr.cc.

References newSectPost_pv(), newSectPre_pv(), Model::OptOList(), and pv().

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

                                         {
  newSectPre_pv("");

  auto& optolist = Model::OptOList();
  for (auto ite = optolist.begin(); ite != optolist.end(); ite++) {
    if ((*ite)->type() == "system")
      continue;
    pv(*ite);
  }

  newSectPost_pv("");
}

void CocoaToDDLMgr::writeRotations

(

)

Definition at line 132 of file CocoaToDDLMgr.cc.

References newSectPost_ro(), newSectPre_ro(), ro(), and theRotationList.

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

                                   {
  newSectPre_ro("");
  std::vector<CLHEP::HepRotation>::const_iterator ite;
  int nc = 0;
  for (ite = theRotationList.begin(); ite != theRotationList.end(); ++ite) {
    //-  std::cout << nc << " rot size " <<  theRotationList.size() << std::endl;
    ro(*ite, nc);
    nc++;
  }
  newSectPost_ro("");
}

void CocoaToDDLMgr::writeSolids

(

)

Definition at line 78 of file CocoaToDDLMgr.cc.

References gather_cfg::cout, newSectPost_so(), newSectPre_so(), Model::OptOList(), and so().

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

                                {
  newSectPre_so("");

  auto& optolist = Model::OptOList();
  for (auto ite = optolist.begin(); ite != optolist.end(); ite++) {
    bool alreadyWritten = false;
    for (auto ite2 = optolist.begin(); ite2 != ite; ite2++) {
      if ((*ite)->shortName() == (*ite2)->shortName()) {
        alreadyWritten = true;
      }
    }
    std::cout << " CocoaToDDLMgr::writeSolids() " << alreadyWritten << *ite;
    std::cout << (*ite)->name() << std::endl;
    if (!alreadyWritten)
      so(*ite);
  }

  newSectPost_so("");
}

void CocoaToDDLMgr::writeSpecPars

(

)

Definition at line 145 of file CocoaToDDLMgr.cc.

References measurementsAsSpecPars(), newSectPost_specPar(), newSectPre_specPar(), Model::OptOList(), specPar(), and writeSpecParsCocoa().

Referenced by writeDDDFile(), and ~CocoaToDDLMgr().

                                  {
  newSectPre_specPar("");

  auto& optolist = Model::OptOList();
  for (auto ite = optolist.begin(); ite != optolist.end(); ite++) {
    if ((*ite)->type() == "system")
      continue;
    specPar(*ite);
  }

  writeSpecParsCocoa();

  //---- Write Measurements's
  measurementsAsSpecPars();

  newSectPost_specPar("");
}

void CocoaToDDLMgr::writeSpecParsCocoa

(

)

Definition at line 671 of file CocoaToDDLMgr.cc.

References file_, ALIFileOut::name(), and Model::OptOList().

Referenced by writeSpecPars(), and ~CocoaToDDLMgr().

                                       {
  file_ << "<!--    Define volumes as COCOA objects --> " << std::endl << "  <SpecPar name=\"COCOA\"> " << std::endl;

  auto& optolist = Model::OptOList();
  for (auto ite = optolist.begin(); ite != optolist.end(); ite++) {
    if ((*ite)->type() == "system")
      continue;
    file_ << "    <PartSelector path=\"/" << (*ite)->name() << "\"/> " << std::endl;
  }

  file_ << "   <String name=\"COCOA\" value=\"COCOA\"/> " << std::endl << "  </SpecPar> " << std::endl;
}

Member Data Documentation

ALIFileOut CocoaToDDLMgr::file_

private

Definition at line 73 of file CocoaToDDLMgr.h.

Referenced by lv(), ma(), measurementsAsSpecPars(), newPartPost(), newPartPre(), newSectPost(), newSectPre(), newSectPre_specPar(), pv(), ro(), so(), specPar(), and writeSpecParsCocoa().

std::string CocoaToDDLMgr::filename_

private

Definition at line 74 of file CocoaToDDLMgr.h.

Referenced by newPartPre(), newSectPre_lv(), newSectPre_ma(), newSectPre_pv(), newSectPre_ro(), newSectPre_so(), and newSectPre_specPar().

CocoaToDDLMgr * CocoaToDDLMgr::instance = nullptr

staticprivate

Definition at line 71 of file CocoaToDDLMgr.h.

Referenced by getInstance(), and production_tasks.Task::getname().

std::vector<CocoaMaterialElementary*> CocoaToDDLMgr::theMaterialList

private

Definition at line 76 of file CocoaToDDLMgr.h.

Referenced by ma(), and materialIsRepeated().

std::vector<CLHEP::HepRotation> CocoaToDDLMgr::theRotationList

private

Definition at line 77 of file CocoaToDDLMgr.h.

Referenced by buildRotationNumber(), and writeRotations().