/* This file is part of MAUS: http://micewww.pp.rl.ac.uk:8080/projects/maus
 *
 * MAUS is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * MAUS is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with MAUS.  If not, see <http://www.gnu.org/licenses/>.
 *
 */
#include "Geant4/G4VPhysicalVolume.hh"
#include "Geant4/G4ThreeVector.hh"
#include "Geant4/G4Tubs.hh"

#include "src/common_cpp/DetModel/SciFi/DoubletFiberParam.hh"

// taking the parameters that define a fibre
DoubletFiberParam::DoubletFiberParam(G4double pCentralFiber,
                                     G4double pActiveRadius,
                                     G4double pOuterDiameter,
                                     G4double pInnerDiameter,
                                     G4double pFiberDiameter,
                                     G4double pFiberPitch ) {
//  sensitiveRadius = pSensitiveRadius;
  centralFiber = pCentralFiber;
  activeRadius = pActiveRadius;
  outerDiameter = pOuterDiameter;
  innerDiameter = pInnerDiameter;
  fiberDiameter = pFiberDiameter;
  fiberPitch = pFiberPitch;

  coreRotation = new G4RotationMatrix(CLHEP::HepRotationX(90.0*deg));

  // This is the spacing between fibres in z
//  zSpacing = sqrt(fiberDiameter*fiberDiameter*(1-fiberPitch*fiberPitch/4));
  zSpacing = sqrt(fiberDiameter*fiberDiameter - fiberPitch*fiberPitch/4);
}

void DoubletFiberParam::ComputeTransformation(const G4int copyNo, G4VPhysicalVolume* physVol)
                                              const {
  // These are the coordinates of the fiber placement
  // channel numbers run from right to left
//  G4double xPos = - copyNo*fiberDiameter*fiberPitch/2
//                  + (activeRadius-fiberDiameter/2);
  G4double xPos = - (copyNo-centralFiber+1)*(fiberPitch/2);
  G4double yPos = 0.0;
  G4double zPos = (copyNo%2) ? 0.5*zSpacing : -0.5*zSpacing;

  physVol->SetRotation(coreRotation);
  physVol->SetTranslation(G4ThreeVector(xPos, yPos, zPos));
}

void DoubletFiberParam::ComputeDimensions(G4Tubs& fiberElement, G4int copyNo,
                                                          const G4VPhysicalVolume* physVol) const {
//  G4double xPos = - copyNo * (fiberDiameter*fiberPitch/2)
//                  + (activeRadius-fiberDiameter/2);

  G4double xPos = - (copyNo-centralFiber+1)*(fiberPitch/2);
  G4double fiberHalfLen = 0.0;

//  if ( sensitiveRadius > fabs( xPos ) ) {
//    fiberHalfLen = sqrt(sensitiveRadius * sensitiveRadius - xPos * xPos);
//  }
  if (activeRadius > fabs(xPos)) {
    fiberHalfLen = sqrt(activeRadius * activeRadius - xPos * xPos);
  }

  fiberElement.SetInnerRadius(innerDiameter/2);
  fiberElement.SetOuterRadius(outerDiameter/2);
  fiberElement.SetZHalfLength(fiberHalfLen);
  fiberElement.SetStartPhiAngle(0.0*deg);
  fiberElement.SetDeltaPhiAngle(360.0*deg);
}