/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/geometry/Detector.cpp

Go to the documentation of this file.

/* @version $Revision: 1703 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2012-05-03 10:11:12 +0200 (Thu, 03 May 2012) $ */

#include "geometry/Detector.hh"
#include "geometry/Chamber.hh"
#include "geometry/GassiplexChamber1.hh"
#include "geometry/GassiplexChamber4.hh"
#include "geometry/GassiplexStripChamber1.hh"
#include "geometry/Hardroc1Chamber.hh"
#include "geometry/Hardroc2Chamber1.hh"
#include "geometry/Hardroc2Chamber2.hh"
#include "geometry/Hardroc2Chamber6.hh"
#include "geometry/Hardroc2IPNLChamber6.hh"
#include "geometry/MicrorocChamber1.hh"
#include "geometry/MicrorocChamber6.hh"
#include "geometry/MicrorocTestChamber.hh"
#include "geometry/MetreCarreChamber.hh"
#include "geometry/DiracChamber1.hh"
#include "geometry/Board.hh"
#include "geometry/Hardroc1Board.hh"
#include "geometry/Hardroc2Board.hh"
#include "geometry/MicrorocBoard.hh"
#include "geometry/DiracBoard.hh"
#include "geometry/Dif.hh"
#include "geometry/Chip.hh"
#include "geometry/GassiplexChip.hh"
#include "geometry/Hardroc1Chip.hh"
#include "geometry/Hardroc2Chip.hh"
#include "geometry/DiracChip.hh"
#include "tools/SteerDesc.hh"

#include "tools/MicroException.hh"
#include "tools/Log.hh"
#include <iostream>

using namespace std;

std::ostream& operator << (std::ostream& out, const Detector& x) {
  return(x.operator <<(out));
}

//-------------------------------------------//
Detector::Detector()
:monDev(0)
{
  //FILE_LOG(logDEBUG) << "-----Detector constructor" << endl;
}

//--------------- Destructeur -------------//
Detector::~Detector()
{
  //FILE_LOG(logDEBUG)<< "----Detector destructor" << endl;
  for (ChamberMap_t::const_iterator it = chambers.begin(); it != chambers.end(); ++it) {
    delete (it->second);
  }

}

//-------------------------------------------//
std::ostream& Detector::operator << (std::ostream& out) const {
    out << description << endl ;
  for (ChamberMap_t::const_iterator it = chambers.begin(); it != chambers.end(); ++it) {
    out << *(it->second);
  }
  return(out);
}

//-------------------------------------------//
bool Detector::addChamber(Chamber* aChamber)
{
  if (!aChamber)
    return false;

  chambers.insert(make_pair(aChamber->getId(), aChamber));
  return true;
}

//-------------------------------------------//
Chamber& Detector::getChamberById(const i32 _id) const
{
  ChamberMap_t::const_iterator result =  chambers.find(_id);
  if (result != chambers.end())
  {
    return *(result->second);
  }
  throw MicroException("No Chamber correspond for this detector and this id.");
}

//-------------------------------------------//
Dif& Detector::getDifById(const i32 _id) const
{
  const ChamberMap_t &chambers = getChambers();
  for (ChamberMap_t::const_iterator it = chambers.begin(); it != chambers.end(); ++it) {
    Chamber& chamber = *(it->second);
    DifMap_t difs = chamber.getDifs();
    DifMap_t::iterator difIt = difs.find(_id);
    if (difIt != difs.end())
      return(*difIt->second);
  }
  throw MicroException("No Dif correspond for this detector and this id.");
}

//-------------------------------------------//
unsigned int Detector::getNumberOfGaziplex96(void) const
{
  unsigned int result = 0;

  ChamberMap_t::const_iterator iter = chambers.begin();

  for ( iter = chambers.begin(); iter != chambers.end() ; iter++ )
  {
//cout << "chamber[" << iter->second->getId() << "]nbGassiplex[" <<  iter->second->getNumberOfBoards() << "]" << endl;
    // only count GassiplexBoards
    /*
      if ( iter->second->getType().compare(GASSIPLEXCHAMBER4)==0
                || iter->second->getType().compare(GASSIPLEXCHAMBER1)==0
                || iter->second->getType().compare(GASSIPLEXSTRIPCHAMBER1)==0)
    */
    if ( iter->second->getType().find("GASSIPLEX") != string::npos )
      {
          result = result + iter->second->getNumberOfBoards();
      }
  }
  return result;
}

//-------------------------------------------//
int Detector::build(const SteerDesc& sd) {


  int nbChamber = 0;

  monDev = sd.monitoringDif; 

  for (int index = 0 ; index < sd.chambers.size(); index++)
  {
    int brow = 0;
    const chamber_info_t& chamberInfo = sd.chambers[index];


    // default unit is micrometer 
    ui32  unit = 1;                               
    if ( chamberInfo.unit.compare("m")==0 )       { unit = 1000000 ;}  // 1 m =  1 000 000  um 
    else if ( chamberInfo.unit.compare("dm")==0 ) { unit = 100000 ;}   // 1 dm = 100 000 um 
    else if ( chamberInfo.unit.compare("cm")==0 ) { unit = 10000 ;}    // 1 cm = 10 000 um 
    else if ( chamberInfo.unit.compare("mm")==0 ) { unit = 1000 ;}     // 1 mm = 1000 um 
    else if ( chamberInfo.unit.compare("um")==0 ) { unit = 1 ;}   
  
    float xPos = chamberInfo.xPos * unit;
    float yPos = chamberInfo.yPos * unit;
    float zPos = chamberInfo.zPos * unit;
    bool xRotation = chamberInfo.xRotation;
    bool yRotation = chamberInfo.yRotation;
    i16 zRotation = chamberInfo.zRotation;
    bool stripHorizontal = chamberInfo.stripHorizontal;

    //Get dif for this chamber
    std::vector<DifInfo_t>difs;  // Steer desc infos:  Chamber Id and  Id
    std::vector<BoardInfo_t>boards;  // Steer desc infos:  Chamber Id dif Id and id
    std::vector<ChipInfo_t>chips;  // Steer desc infos:  Chamber Id dif Id Board Id and id

    // fill Dif, Boards and chip for current chamber
    sd.fillDifByChamberId(difs,chamberInfo.id);
    sd.fillBoardByChamberId(boards,chamberInfo.id);
    sd.fillChipByChamberId(chips,chamberInfo.id);

    FILE_LOG(logDEBUG1)<< "Detector.build : chamber " << chamberInfo.id << " ,chamberInfo.xPos: " << xPos
                      << "um ,chamberInfo.yPos: " << yPos
                      << "um ,chamberInfo.zPos: " << zPos 
                      << "um ,chamberInfo.stripHorizontal: " << chamberInfo.stripHorizontal << endl;
    if ( !chamberInfo.type.compare(GASSIPLEXCHAMBER1))
    {
      GassiplexChamber1 *chamber = new GassiplexChamber1(xPos, yPos, zPos,  xRotation, yRotation, zRotation, 1, 1, chamberInfo.id);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }

    else if ( !chamberInfo.type.compare(GASSIPLEXCHAMBER4))
    {
      GassiplexChamber4 *chamber = new GassiplexChamber4(xPos, yPos, zPos,   xRotation,yRotation, zRotation,2, 2, chamberInfo.id);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }

    }

    else if ( !chamberInfo.type.compare(GASSIPLEXSTRIPCHAMBER1))
    {
      GassiplexStripChamber1 *chamber = new GassiplexStripChamber1(xPos, yPos, zPos,   xRotation,yRotation, zRotation,1, 1, chamberInfo.id,chamberInfo.stripHorizontal);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }

    }


    else if ( !chamberInfo.type.compare(HARDROC1CHAMBER))
    /*
    ASU Test Beam 2009 et calibration: 1 Chambre (Hardroc1Chamber) composee de 1 PCB ( Hardroc1Board)
    composee de 4 Chip HR1 (Hardroc1Chip ) composee de 64 Channel ( Hardroc1Channel )
    */
    {
      Hardroc1Chamber *chamber = new Hardroc1Chamber(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }

    /*
    * ASU "Calibration 2009" : 1 Chambre ( Hardroc2Chamber1 ) composee de 1 PCB ( Hardroc2Board24 )
    * composee de 24 Chip HR2 ( * Hardroc2Chip ) composee  de 64 Channel ( Hardroc2Channel )
    */
    else if (  !chamberInfo.type.compare(HARDROC2CHAMBER1))
    {
      Hardroc2Chamber1 *chamber =  new Hardroc2Chamber1(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }
    /*
    * ASU  : 1 Chambre ( Hardroc2Chamber2 ) composee de 2 PCB ( Hardroc2Board24 )
    * composee de 24 Chip HR2 ( * Hardroc2Chip ) composee  de 64 Channel ( Hardroc2Channel )
    */

    /*
    * ASU  : 1 Chambre ( Hardroc2Chamber6 ) composee de 6 PCB ( Hardroc2Board24 )
    * composee de 24 Chip HR2 ( * Hardroc2Chip ) composee  de 64 Channel ( Hardroc2Channel )
    */
    else if (  !chamberInfo.type.compare(HARDROC2CHAMBER2))
    {
      Hardroc2Chamber2 *chamber = new Hardroc2Chamber2(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }

    /*
    * 1 Chambre ( Hardroc2Chamber6 ) composee de 6 PCB ( Hardroc2Board24 )
    * composee de 24 Chip HR2 ( * Hardroc2Chip ) composee  de 64 Channel ( Hardroc2Channel )
    */
    else if (  !chamberInfo.type.compare(HARDROC2CHAMBER6))
    {
      Hardroc2Chamber6 *chamber = new Hardroc2Chamber6(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }

    /*
    * ASU "IPNL " : 1 Chambre ( Hardroc2Chamber1 ) composee de 1 PCB ( Hardroc2Board24 )
    * composee de 24 Chip HR2 IPNL ( * Hardroc2IPNLChip ) composee  de 64 Channel ( Hardroc2Channel )
    */
    else if (  !chamberInfo.type.compare(HARDROC2IPNLCHAMBER6))
    {
      Hardroc2IPNLChamber6 *chamber = new Hardroc2IPNLChamber6(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }

        /*
    // Mettre Carre HR@ pour test beam june 2010
    * composee de 5 ASU de 24 Chip HR2 ( * Hardroc2Chip ) composee  de 64 Channel ( Hardroc2Channel )
    */
    else if (  !chamberInfo.type.compare(METRECARRECHAMBER))
    {
      MetreCarreChamber *chamber = new MetreCarreChamber(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }


     /*
    // Calibration  MICROROC@ pour test beam  2011
    * composee de 1 ASU de 24 Chip Microroc ( * MicrorocChip ) composee  de 64 Channel ( MicrorocChannel )
    */
    else if (  !chamberInfo.type.compare(MICROROCCHAMBER1))
    {
      MicrorocChamber1 *chamber = new MicrorocChamber1(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }
     /*
    // Mettre Carre MICROROC@ pour test beam  2011
    * composee de 6 ASU de 24 Chip Microroc ( * MicrorocChip ) composee  de 64 Channel ( MicrorocChannel )
    */
    else if (  !chamberInfo.type.compare(MICROROCCHAMBER6))
    {
      MicrorocChamber6 *chamber = new MicrorocChamber6(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }

     /*
    // Carte de test MICROROC
    * composee de 1 ASU de 14 Chip Microroc ( * MicrorocChip ) composee  de 64 Channel ( MicrorocChannel )
    */
    else if (  !chamberInfo.type.compare(MICROROCTESTCHAMBER))
    {
      MicrorocTestChamber *chamber = new MicrorocTestChamber(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    }
   else if ( !chamberInfo.type.compare(DIRACCHAMBER1))
    {
      DiracChamber1 *chamber = new DiracChamber1(xPos, yPos, zPos,    xRotation,yRotation, zRotation, chamberInfo.id,difs,boards,chips);
      if ( chamber->isValid()) {
        addChamber(chamber);
        nbChamber++;
      }
    } // endif DIRAC1
  }

  //FILE_LOG(logDEBUG) << *this;
  return nbChamber;
} // Detector.build

---