source: trunk/examples/advanced/composite_calorimeter/src/CCalHcal.cc @ 1303

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
///////////////////////////////////////////////////////////////////////////////
// File: CCalHcal.cc
// Description: CCalHcal Geometry factory class for the hadron calorimeter
///////////////////////////////////////////////////////////////////////////////
#include "CCalHcal.hh"

#include <fstream>
#include "CCalutils.hh"

//#define debug

CCalHcal::CCalHcal(const G4String &name):
  CCalDetector(name),xposBox(0),typeLayerScnt(0),mothLayerScnt(0),
  xposLayerScnt(0),typeLayerAbs(0),mothLayerAbs(0),xposLayerAbs(0),
  dx_2Absorber(0),dy_2ScntLayer(0),dx_2ScntLayer(0),dx_2Wrapper(0),
  dx_2FrontPlastic(0),dx_2BackPlastic(0),dx_2Scintillator(0) {}

CCalHcal::~CCalHcal() {
  if (xposBox)
    delete[] xposBox;
  if (typeLayerScnt)
    delete[] typeLayerScnt;
  if (mothLayerScnt)
    delete[] mothLayerScnt;
  if (xposLayerScnt)
    delete[] xposLayerScnt;
  if (typeLayerAbs)
    delete[] typeLayerAbs;
  if (mothLayerAbs)
    delete[] mothLayerAbs;
  if (xposLayerAbs)
    delete[] xposLayerAbs;
  if (dx_2Absorber)
    delete[] dx_2Absorber;
  if (dy_2ScntLayer)
    delete[] dy_2ScntLayer;
  if (dx_2ScntLayer)
    delete[] dx_2ScntLayer;
  if (dx_2Wrapper)
    delete[] dx_2Wrapper;
  if (dx_2FrontPlastic)
    delete[] dx_2FrontPlastic;
  if (dx_2BackPlastic)
    delete[] dx_2BackPlastic;
  if (dx_2Scintillator)
    delete[] dx_2Scintillator;
}

int CCalHcal::readFile() {
  ///////////////////////////////////////////////////////////////
  //Let's open the file
  G4cout << " ==> Opening file " << File() << " to read elements..."
       << G4endl;

  std::ifstream is;
  bool ok = openGeomFile(is, pathName, File());
  if (!ok)
    return 0;

  // Find *DO HCal
  findDO(is, G4String("HCal"));

  // Calorimeter boundaries
  readName(is,genMaterial);
  is >> dy_2Cal >> dx_2Cal >> xposCal >> jump;
#ifdef debug
  G4cout << tab << "General material: " << genMaterial  << " Size " << dy_2Cal
       << ", " << dx_2Cal << " Position " << xposCal << G4endl;
#endif

  // Boxes
  readName(is,boxMaterial);
  is >> nBox >> dy_2Box >> dx_2Box >> wallThickBox;
  int i = 0;
  xposBox = new double[nBox];
  for (i=0; i<nBox; i++)
    is >> xposBox[i];
#ifdef debug
  G4cout << tab << "Box material: " << boxMaterial  << " Size " << dy_2Box
       << ", " << dx_2Box << " Wall Thickness " << wallThickBox << " number "
       << nBox << " position ";
  for (i=0; i<nBox; i++)
    G4cout << i << " " << xposBox[i] << " ";
  G4cout << G4endl;
#endif

  // Layers of scintillators
  G4String rubbish;
  readName(is,rubbish);
  is >> nLayerScnt;
  typeLayerScnt = new int[nLayerScnt];
  mothLayerScnt = new int[nLayerScnt];
  xposLayerScnt = new double[nLayerScnt];
  for (i=0; i<nLayerScnt; i++)
    is >> typeLayerScnt[i] >> mothLayerScnt[i] >> xposLayerScnt[i];
#ifdef debug
  G4cout << tab << nLayerScnt << " Layers of scintillators of type/mother box/"
       << "position" << G4endl;
  for (i=0; i<nLayerScnt; i++)
    G4cout << tab << i << " " << typeLayerScnt[i] << " " << mothLayerScnt[i] 
         << " " << xposLayerScnt[i] << G4endl;
#endif

  // Layers of absorbers
  readName(is,rubbish);
  is >> nLayerAbs;
  typeLayerAbs = new int[nLayerAbs];
  mothLayerAbs = new int[nLayerAbs];
  xposLayerAbs = new double[nLayerAbs];
  for (i=0; i<nLayerAbs; i++)
    is >> typeLayerAbs[i] >> mothLayerAbs[i] >> xposLayerAbs[i];
#ifdef debug
  G4cout << tab << nLayerAbs << " Layers of absorbers of type/mother box/"
       << "position" << G4endl;
  for (i=0; i<nLayerAbs; i++)
    G4cout << tab << i << " " << typeLayerAbs[i] << " " << mothLayerAbs[i] 
         << " " << xposLayerAbs[i] << G4endl;
#endif

  // Absorber parameters
  readName(is,absMaterial);
  is >> nAbsorber >> dy_2Absorber;
  dx_2Absorber = new double[nAbsorber];
  for (i=0; i<nAbsorber; i++)
    is >> dx_2Absorber[i];
#ifdef debug
  G4cout << "\tAbsorber mad of " << absMaterial << " with " << nAbsorber
       << " types and size " << dy_2Absorber;
  for (i=0;  i<nAbsorber; i++)
    G4cout << "  " << i << " " << dx_2Absorber[i];
  G4cout << G4endl;
#endif

  // Scintillator parameters
  readName(is,scntMaterial);
  readName(is,wrapMaterial);
  readName(is,plasMaterial);
  is >> nScintillator;
  dy_2ScntLayer    = new double[nScintillator];
  dx_2ScntLayer    = new double[nScintillator];
  dx_2Wrapper      = new double[nScintillator];
  dx_2FrontPlastic = new double[nScintillator];
  dx_2BackPlastic  = new double[nScintillator];
  dx_2Scintillator = new double[nScintillator];
  for (i=0; i<nScintillator; i++)
    is >> dy_2ScntLayer[i] >> dx_2ScntLayer[i] >> dx_2Wrapper[i]
       >> dx_2FrontPlastic[i] >> dx_2BackPlastic[i] >> dx_2Scintillator[i];
#ifdef debug
  G4cout << tab << nScintillator << " Scintillator layers made of " 
       << scntMaterial << " " << wrapMaterial << " and " << plasMaterial 
       << " of sizes " << G4endl;
  for (i=0; i<nScintillator; i++)
    G4cout << tab << i << " " << dy_2ScntLayer[i] << " " << dx_2ScntLayer[i] 
         << " " << dx_2Wrapper[i] << " " << dx_2FrontPlastic[i] << " " 
         << dx_2BackPlastic[i] << " " << dx_2Scintillator[i] << G4endl;
#endif
  
  ///////////////////////////////////////////////////////////////
  // Close the file
  G4cout << " ==> Closing file " << File() << G4endl;
  is.close();

  return 1;

}

void CCalHcal::constructDaughters() {
}