#include <ShipTAUMagneticSpectrometer.h>

Inheritance diagram for ShipTAUMagneticSpectrometer:

Collaboration diagram for ShipTAUMagneticSpectrometer:

[legend]

Public Member Functions
	ShipTAUMagneticSpectrometer (const char name, const Double_t zLS, const Double_t FeL, const Double_t AirL, const Double_t SpectroL, const Double_t GapV, const Double_t DGap, const Double_t MGap, const Double_t mf, const char Title="ShipTAUMagneticSpectrometer")

	ShipTAUMagneticSpectrometer ()

virtual	~ShipTAUMagneticSpectrometer ()

void	ConstructGeometry ()

void	Initialize ()

Public Attributes
ClassDef(ShipTAUMagneticSpectrometer, 1) protected Double_t	IronLenght

Double_t	AirLenght

Double_t	SpectrometerLenght

Double_t	GapFromVacuum

Double_t	DriftGap

Double_t	MiddleGap

Double_t	MagneticField

Detailed Description

Definition at line 10 of file ShipTAUMagneticSpectrometer.h.

Constructor & Destructor Documentation

◆ ShipTAUMagneticSpectrometer() [1/2]

ShipTAUMagneticSpectrometer::ShipTAUMagneticSpectrometer	(	const char *	name,
		const Double_t	zLS,
		const Double_t	FeL,
		const Double_t	AirL,
		const Double_t	SpectroL,
		const Double_t	GapV,
		const Double_t	DGap,
		const Double_t	MGap,
		const Double_t	mf,
		const char *	Title = `"ShipTAUMagneticSpectrometer"`
	)

Definition at line 34 of file ShipTAUMagneticSpectrometer.cxx.

  : FairModule(name ,Title)
{
    zLastSlab = zLS;
    IronLenght = FeL;
    AirLenght = AirL;
    SpectrometerLenght = SpectroL;
    GapFromVacuum = GapV;
    DriftGap = DGap;
    MiddleGap = MGap;
    MagneticField = mf;
}

◆ ShipTAUMagneticSpectrometer() [2/2]

ShipTAUMagneticSpectrometer::ShipTAUMagneticSpectrometer ( )

Definition at line 29 of file ShipTAUMagneticSpectrometer.cxx.

  : FairModule("ShipTAUMagneticSpectrometer", "")
{
}

◆ ~ShipTAUMagneticSpectrometer()

ShipTAUMagneticSpectrometer::~ShipTAUMagneticSpectrometer ( )

virtual

Definition at line 25 of file ShipTAUMagneticSpectrometer.cxx.

26{

27}

Member Function Documentation

◆ ConstructGeometry()

void ShipTAUMagneticSpectrometer::ConstructGeometry ( )

Definition at line 47 of file ShipTAUMagneticSpectrometer.cxx.

{
    Int_t NIronSlabs = 12;
    
    TGeoVolume *top=gGeoManager->GetTopVolume();
    TGeoMedium *Fe  = gGeoManager->GetMedium("iron");
    
    Double_t d = 0;
    
    TGeoUniformMagField *magField = new TGeoUniformMagField(0.,-MagneticField,0.);
    TGeoUniformMagField *RetField     = new TGeoUniformMagField(0.,MagneticField,0.);
    
    TGeoBBox *Layer = new TGeoBBox(225,400,IronLenght/2);
    TGeoVolume *volLayer = new TGeoVolume("volLayer",Layer,Fe);
    for(Int_t i = 0; i< NIronSlabs; i++)
    {
        d = zLastSlab - i*(IronLenght+AirLenght);
        top->AddNode(volLayer,i,new TGeoTranslation(0, 0, d));
    }
    volLayer->SetField(magField);
    
    cout <<"************************************" << endl;
    cout << " IronLenght+AirLenght = " << IronLenght+AirLenght << endl;
    cout << " NIronSlabs = " << NIronSlabs << endl;
    cout << " zLastSlab = "<< zLastSlab << endl;
    cout << d << endl;
    Double_t d1 = d- (MiddleGap + IronLenght); //z coord of the center of the last layer of the first spectrometer
    cout << d1 << endl;
    
    TGeoVolume *volLayer2 = new TGeoVolume("volLayer2",Layer,Fe);
    for(Int_t i = 0; i< NIronSlabs; i++)
    {
        Double_t d2 = d1-i*(IronLenght+AirLenght);
        top->AddNode(volLayer2,i,new TGeoTranslation(0, 0, d2));
    }
    volLayer2->SetField(RetField);
    
    cout <<"************************************" << endl;
    
    
}