CCalEcal.cc

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CCalEcal.cc

//
// ********************************************************************
// * 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: CCalEcal.cc
// Description: CCalEcal Geometry factory class for crystal matrix
#include "CCalEcal.hh"

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

//#define debug

CCalEcal::~CCalEcal() {}

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

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

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

  //Let's read overall box dimensions and positions
  readName(is,genMat);
  is >> widBox >> lengBox >> xpos >> ypos >> zpos >> thetaX >> phiX
     >> thetaY >> phiY >> thetaZ >> phiZ >> jump;
#ifdef debug
  G4cout << tab << "General material: " << genMat << "\tBox dimensions "
       << widBox << ", " << lengBox << G4endl;
  G4cout << tab << "Positioned at (" << xpos << ", " << ypos << ", " << zpos
       << ") with rotation (" << thetaX << ", " << phiX << ", " << thetaY
       << ", " << phiY << ", " << thetaZ << ", " << phiZ << ")" << G4endl;
#endif

  //Then the layer positions
  G4int i=0;
  readName(is,layMat);
  is >> layNum >> layRadius >> layAngle >> lengFront;
  for (i=0; i<5; i++) 
    is >> layPar[i];
#ifdef debug
  G4cout << tab << "Layer material: " << layMat << " Number " << layNum
       << " Radius " << layRadius << " Angle " << layAngle/deg 
       << " front dist " << lengFront << " Parameters ";
  for (i=0; i<5; i++)
    G4cout << layPar[i] << " ";
  G4cout << G4endl;
#endif

  //Then the crystal positions
  readName(is,crystMat);
  is >> crystNum >> crystLength >> crystTol;
  for (i=0; i<5; i++) 
    is >> crystPar[i];
#ifdef debug
  G4cout << tab << "Crystal material: " << crystMat << " Number " << crystNum
       << " Length " << crystLength << " Tolerance " << crystTol
       << " Parameters ";
  for (i=0; i<5; i++)
    G4cout << crystPar[i] << " ";
  G4cout << G4endl;
#endif

  //Then the support material
  readName(is,suppMat);
  is >> dxSupp >> dySupp >> dzSupp >> distSupp >> jump;
#ifdef debug
  G4cout << tab << "Support material: " << suppMat << " Dimensions " << dxSupp
       << ", " << dySupp << ", " << dzSupp << " Distance " << distSupp << G4endl;
#endif
   
  // Close the file
  G4cout << " ==> Closing file " << File() << G4endl;
  is.close();

  return 1;

}

void CCalEcal::constructDaughters() {}