MollerGlobalMagnetField.cc

Go to the documentation of this file.

//=============================================================================
// 
//   ---------------------------
//  | Doxygen File Information |
//  ---------------------------
// 
/**
 
   \file MollerGlobalMagnetField.cc

   $Revision: 1.3 $ 
   $Date: 2006/03/01 17:03:16 $

   \authors Klaus Hans Grimm Juliette M. Mammei

*/
//=============================================================================

//=============================================================================
//   -----------------------
//  | CVS File Information |
//  -----------------------
// 
//  Last Update:      $Author: mammei $
//  Update Date:      $Date: 2010/08/17 09:48:30 $
//  CVS/RCS Revision: $Revision: 1.3 $
//  Status:           $State: Exp $
// 
// ===================================
//  CVS Revision Log at end of file !!
// ===================================
//
//============================================================================

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
              
#include "MollerGlobalMagnetField.hh"
#include "MollerGlobalMagnetFieldMessenger.hh"

#include "G4UImanager.hh"
#include "G4UIcommand.hh"

#include "globals.hh"

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

MollerGlobalMagnetField::MollerGlobalMagnetField() 
{ 

  G4cout << G4endl << "###### Calling MollerGlobalMagnetField::MollerGlobalMagnetField() " << G4endl << G4endl;

  fieldMessenger = new MollerGlobalMagnetFieldMessenger(this); 

  G4String fPath = ".";
  if (getenv("MOLLERGEANTDIR"))
    fPath = getenv("MOLLERGEANTDIR");
  else
    G4cout << "You do not have MOLLERGEANTDIR defined in your environment!" << G4endl;

  
  fMagneticField_MainMagnet = NULL;
  fMagneticField_MiniMagnet = NULL;

}


void MollerGlobalMagnetField::ReadMagneticField()
{

  G4String fPath = ".";
  if (getenv("MOLLERGEANTDIR"))
    fPath = getenv("MOLLERGEANTDIR");
  else
    G4cout << "You do not have MOLLERGEANTDIR defined in your environment!" << G4endl;

  G4String upstreamFieldMap = (fPath+"/map_directory/"+upstreamFilename+".txt");
  G4String hybridFieldMap = (fPath+"/map_directory/"+hybridFilename+".txt");

  G4cout<<"Field maps:"<<G4endl;
  G4cout<<upstreamFieldMap<<G4endl;
  G4cout<<hybridFieldMap<<G4endl;
  
  BFieldScalingFactor_MainMagnet = 1.0; // default is 1.0 here
  fMagneticField_MainMagnet = new MollerMagnetFieldMap();

  if(fMagneticField_MainMagnet){

    fMagneticField_MainMagnet->SetFieldMap_RMin(0.015*m);      
    fMagneticField_MainMagnet->SetFieldMap_RMax(0.375*m); 
    fMagneticField_MainMagnet->SetFieldMap_RStepsize(0.002*m);
  
    fMagneticField_MainMagnet->SetFieldMap_ZMin(9.5*m); 
    fMagneticField_MainMagnet->SetFieldMap_ZMax(17.0*m);
    fMagneticField_MainMagnet->SetFieldMap_ZStepsize(0.1*m);
  
    fMagneticField_MainMagnet->SetFieldMap_PhiMin(-25*degree);
    fMagneticField_MainMagnet->SetFieldMap_PhiMax(25*degree);
    fMagneticField_MainMagnet->SetFieldMap_PhiStepsize(1*degree);
  
    fMagneticField_MainMagnet->InitializeGrid();
  
    fMagneticField_MainMagnet->ReadFieldMap(hybridFieldMap);
    //fMagneticField_MainMagnet->PrintFieldCheck("magfieldcheck_back.txt");
  }

  //---------------------------------------------------------------
  //    configuration values of the upstream field map
  //---------------------------------------------------------------
        
  BFieldScalingFactor_MiniMagnet =  1.25; // default is 1.25 here
  fMagneticField_MiniMagnet = new MollerMagnetFieldMap();

  if(fMagneticField_MiniMagnet){
  
    fMagneticField_MiniMagnet->SetFieldMap_RMin(0.015*m);      
    fMagneticField_MiniMagnet->SetFieldMap_RMax(0.375*m); 
    fMagneticField_MiniMagnet->SetFieldMap_RStepsize(0.002*m);
  
    fMagneticField_MiniMagnet->SetFieldMap_ZMin(6.0*m); 
    fMagneticField_MiniMagnet->SetFieldMap_ZMax(8.0*m);
    fMagneticField_MiniMagnet->SetFieldMap_ZStepsize(0.1*m);
  
    fMagneticField_MiniMagnet->SetFieldMap_PhiMin(-25*degree);
    fMagneticField_MiniMagnet->SetFieldMap_PhiMax(25*degree);
    fMagneticField_MiniMagnet->SetFieldMap_PhiStepsize(1*degree);
  
    fMagneticField_MiniMagnet->InitializeGrid();

    fMagneticField_MiniMagnet->ReadFieldMap(upstreamFieldMap);
  
    //fMagneticField_MiniMagnet->PrintGridCheck("maggridcheck_front.txt");
  }

  G4cout << G4endl << "###### Leaving MollerGlobalMagnetField::MollerGlobalMagnetField " << G4endl << G4endl;
}   
void MollerGlobalMagnetField::SetHybridMapFile(const G4String& nam)
{
  hybridFilename = nam;
}

void MollerGlobalMagnetField::SetUpstreamMapFile(const G4String& nam)
{
  upstreamFilename = nam;
}
    
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
MollerGlobalMagnetField::~MollerGlobalMagnetField()
{  
  if (fMagneticField_MainMagnet) delete fMagneticField_MainMagnet;  
  if (fMagneticField_MiniMagnet) delete fMagneticField_MiniMagnet;
    
  delete fieldMessenger;
} 
    
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
void MollerGlobalMagnetField::GetFieldValue(const G4double Point[4], G4double *Bfield ) const
{

  G4double myLocalPointInMainMagnet[4];
  G4double myLocalPointInMiniMagnet[4];
  
  G4double myLocalBfieldInMainMagnet[3];
  G4double myLocalBfieldInMiniMagnet[3];
  
  //---------------------------------------------------------------
  // translation from global Point[4] into local magnet coordinates
  //---------------------------------------------------------------
  //
  // No translation needed because for now, map is created with actual z positions
  // If add command to move fields, would do it here.

  myLocalPointInMainMagnet[0] = Point[0];             // x-pos
  myLocalPointInMainMagnet[1] = Point[1];             // y-pos
  myLocalPointInMainMagnet[2] = Point[2];             // z-pos
  myLocalPointInMainMagnet[3] = Point[3];             // time

  myLocalPointInMiniMagnet[0] = Point[0];             // x-pos
  myLocalPointInMiniMagnet[1] = Point[1];             // y-pos
  myLocalPointInMiniMagnet[2] = Point[2];             // z-pos
  myLocalPointInMiniMagnet[3] = Point[3];             // time


  //     G4cout << " X pos [cm] = " << Point[0]/cm << " ,Y pos [cm] = " << Point[1]/cm << " ,Z pos [cm] = " << Point[2]/cm << G4endl;


  if(fMagneticField_MiniMagnet && fMagneticField_MainMagnet){
    fMagneticField_MiniMagnet->GetFieldValue( myLocalPointInMiniMagnet, myLocalBfieldInMiniMagnet );
    fMagneticField_MainMagnet->GetFieldValue( myLocalPointInMainMagnet, myLocalBfieldInMainMagnet );
    
    Bfield[0] = myLocalBfieldInMainMagnet[0]*BFieldScalingFactor_MainMagnet + myLocalBfieldInMiniMagnet[0]*BFieldScalingFactor_MiniMagnet;
    Bfield[1] = myLocalBfieldInMainMagnet[1]*BFieldScalingFactor_MainMagnet + myLocalBfieldInMiniMagnet[1]*BFieldScalingFactor_MiniMagnet;
    Bfield[2] = myLocalBfieldInMainMagnet[2]*BFieldScalingFactor_MainMagnet + myLocalBfieldInMiniMagnet[2]*BFieldScalingFactor_MiniMagnet;
    //     printf("%f %f %f %f %f %f %f %f %f %f %f %f\n",Point[0],Point[1],Point[2],
    //     myLocalBfieldInMainMagnet[0]*100,myLocalBfieldInMiniMagnet[0]*100,Bfield[0]*100,
    //     myLocalBfieldInMainMagnet[1]*100,myLocalBfieldInMiniMagnet[1]*100,Bfield[1]*100,
    //     myLocalBfieldInMainMagnet[2]*100,myLocalBfieldInMiniMagnet[2]*100,Bfield[2]*100);
  }
  else if(fMagneticField_MainMagnet && !fMagneticField_MiniMagnet){
    fMagneticField_MainMagnet->GetFieldValue( myLocalPointInMainMagnet, myLocalBfieldInMainMagnet );
    
    Bfield[0] = myLocalBfieldInMainMagnet[0]*BFieldScalingFactor_MainMagnet;
    Bfield[1] = myLocalBfieldInMainMagnet[1]*BFieldScalingFactor_MainMagnet;
    Bfield[2] = myLocalBfieldInMainMagnet[2]*BFieldScalingFactor_MainMagnet;
  }
  else if(!fMagneticField_MainMagnet && fMagneticField_MiniMagnet){
    fMagneticField_MiniMagnet->GetFieldValue( myLocalPointInMiniMagnet, myLocalBfieldInMiniMagnet );

    Bfield[0] =  myLocalBfieldInMiniMagnet[0]*BFieldScalingFactor_MiniMagnet;
    Bfield[1] =  myLocalBfieldInMiniMagnet[1]*BFieldScalingFactor_MiniMagnet;
    Bfield[2] =  myLocalBfieldInMiniMagnet[2]*BFieldScalingFactor_MiniMagnet;    
  }
  else{
    Bfield[0] = 0.0*kilogauss;
    Bfield[1] = 0.0*kilogauss;
    Bfield[2] = 0.0*kilogauss;
  }
  
  //   G4cout << "------------------------------------------------------------------------" << G4endl;
  //   G4cout << " X pos [cm] = " << Point[0]/cm << "  , with X BField [kG] = " << Bfield[0] << G4endl;
  //   G4cout << " Y pos [cm] = " << Point[1]/cm << "  , with Y BField [kG] = " << Bfield[1] << G4endl;
  //   G4cout << " Z pos [cm] = " << Point[2]/cm << "  , with Z BField [kG] = " << Bfield[2] << G4endl;

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

//=======================================================
//   -----------------------
//  | CVS File Information |
//  -----------------------
// 
//      $Revisions$  
//      $Log: MollerGlobalMagnetField.cc,v $
//      Revision 1.3  2006/03/01 17:03:16  grimm
//      Changed BFieldScalingFactor_MiniMagnet =  -0.15434; // according to Juliette Mammei
//      Before the scaling factor was ~0.07
//
//      Revision 1.2  2005/12/27 19:10:21  grimm
//      - Redesign of Doxygen header containing CVS info like revision and date
//      - Added CVS revision log at the end of file
//
//