MagneticSpectrometer Class Reference

#include <MagneticSpectrometer.h>

Inheritance diagram for MagneticSpectrometer:

Collaboration diagram for MagneticSpectrometer:

Public Member Functions
	MagneticSpectrometer (const char *name, const Double_t Zcenter, Bool_t Active, const char *Title="MagneticSpectrometer")
	MagneticSpectrometer ()
virtual	~MagneticSpectrometer ()
void	SetDesign (Int_t Design)
void	SetTotDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetFeDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetRpcDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetRpcStripDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetRpcGasDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetRpcElectrodeDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetRpcPETDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetNFeInArm (Int_t N)
void	SetNRpcInArm (Int_t N)
void	SetZDimensionArm (Double_t Z)
void	SetGapDownstream (Double_t Gap)
void	SetGapMiddle (Double_t Gap)
void	SetMagneticField (Double_t B)
void	SetReturnYokeDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetSmallerYokeDimensions (Double_t X, Double_t Y, Double_t Z)
void	SetCoilParameters (Double_t CoilH, Double_t CoilW, Int_t N, Double_t CoilG)
void	SetPillarDimensions (Double_t X, Double_t Y, Double_t Z)
void	ConstructGeometry ()
virtual void	Initialize ()
virtual Bool_t	ProcessHits (FairVolume *v=0)
virtual void	Register ()
virtual TClonesArray *	GetCollection (Int_t iColl) const
virtual void	Reset ()
ShipRpcPoint *	AddHit (Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgCode)
virtual void	CopyClones (TClonesArray *cl1, TClonesArray *cl2, Int_t offset)
virtual void	SetSpecialPhysicsCuts ()
virtual void	EndOfEvent ()
virtual void	FinishPrimary ()
virtual void	FinishRun ()
virtual void	BeginPrimary ()
virtual void	PostTrack ()
virtual void	PreTrack ()
virtual void	BeginEvent ()
void	DecodeVolumeID (Int_t detID, int &nARM, int &nRPC)

Protected Member Functions
	MagneticSpectrometer (const MagneticSpectrometer &)
MagneticSpectrometer &	operator= (const MagneticSpectrometer &)

Protected Attributes
Int_t	fDesign
Double_t	fZcenter
Double_t	fXtot
Double_t	fYtot
Double_t	fZtot
Int_t	fNFe
Int_t	fNRpc
Double_t	fXFe
Double_t	fXRpc
Double_t	fYFe
Double_t	fYRpc
Double_t	fZFe
Double_t	fZRpc
Double_t	fZArm
Double_t	fGapDown
Double_t	fGapMiddle
Double_t	fField
Double_t	fXRyoke
Double_t	fYRyoke
Double_t	fZRyoke
Double_t	fXRyoke_s
Double_t	fYRyoke_s
Double_t	fZRyoke_s
Double_t	fCoilGap
Double_t	fCoilH
Double_t	fCoilW
Int_t	fNCoil
Double_t	fXStrip
Double_t	fYStrip
Double_t	fZStrip
Double_t	fXPet
Double_t	fYPet
Double_t	fZPet
Double_t	fXEle
Double_t	fYEle
Double_t	fZEle
Double_t	fXGas
Double_t	fYGas
Double_t	fZGas
Double_t	fPillarX
Double_t	fPillarY
Double_t	fPillarZ

Private Member Functions
Int_t	InitMedium (const char *name)

Private Attributes
Int_t	fTrackID
Int_t	fPdgCode
	track index
Int_t	fVolumeID
	pdg code
TLorentzVector	fPos
	volume id
TLorentzVector	fMom
	position at entrance
Double32_t	fTime
	momentum at entrance
Double32_t	fLength
	time
Double32_t	fELoss
	length
TClonesArray *	fShipRpcPointCollection
	energy loss

Detailed Description

Definition at line 18 of file MagneticSpectrometer.h.

Constructor & Destructor Documentation

MagneticSpectrometer() [1/3]

MagneticSpectrometer::MagneticSpectrometer	(	const char *	name,
		const Double_t	Zcenter,
		Bool_t	Active,
		const char *	Title = "MagneticSpectrometer"
	)

Definition at line 59 of file MagneticSpectrometer.cxx.

  : FairDetector(name, Active, ktauRpc),
    fTrackID(-1),
    fPdgCode(),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fELoss(-1),
    fShipRpcPointCollection(new TClonesArray("ShipRpcPoint"))
{
  fZcenter = Zcenter;
}

MagneticSpectrometer() [2/3]

MagneticSpectrometer::MagneticSpectrometer

(

)

Definition at line 45 of file MagneticSpectrometer.cxx.

  : FairDetector("MagneticSpectrometer",kTRUE, ktauRpc),
    fTrackID(-1),
    fPdgCode(),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fELoss(-1),
    fShipRpcPointCollection(new TClonesArray("ShipRpcPoint"))
{
}

~MagneticSpectrometer()

MagneticSpectrometer::~MagneticSpectrometer ( )

virtual

Definition at line 190 of file MagneticSpectrometer.cxx.

{
  if (fShipRpcPointCollection) {
    fShipRpcPointCollection->Delete();
    delete fShipRpcPointCollection;
  }
}

MagneticSpectrometer() [3/3]

MagneticSpectrometer::MagneticSpectrometer ( const MagneticSpectrometer &  )

protected

Member Function Documentation

AddHit()

ShipRpcPoint * MagneticSpectrometer::AddHit	(	Int_t	trackID,
		Int_t	detID,
		TVector3	pos,
		TVector3	mom,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss,
		Int_t	pdgCode
	)

This method is an example of how to add your own point of type muonPoint to the clones array

Definition at line 597 of file MagneticSpectrometer.cxx.

{
  TClonesArray& clref = *fShipRpcPointCollection;
  Int_t size = clref.GetEntriesFast();
  //cout << "ShipRpctau hit called"<< pos.z()<<endl;
  //    return new(clref[size]) ShipRpcPoint(trackID, detID, pos, mom,time, length, eLoss, pdgCode,NArm, NRpc, NHpt);
  return new(clref[size]) ShipRpcPoint(trackID, detID, pos, mom,time, length, eLoss, pdgCode);
}

BeginEvent()

virtual void MagneticSpectrometer::BeginEvent ( )

inlinevirtual

Definition at line 85 of file MagneticSpectrometer.h.

{;}

BeginPrimary()

virtual void MagneticSpectrometer::BeginPrimary ( )

inlinevirtual

Definition at line 82 of file MagneticSpectrometer.h.

{;}

ConstructGeometry()

void MagneticSpectrometer::ConstructGeometry

(

)

Definition at line 224 of file MagneticSpectrometer.cxx.

{
  TGeoVolume *top = gGeoManager->GetTopVolume();
  TGeoVolumeAssembly *tTauNuDet = new TGeoVolumeAssembly("tTauNuDet");
  top->AddNode(tTauNuDet, 1, new TGeoTranslation(0, 0, 0));
 
  InitMedium("RPCgas");
  TGeoMedium *RPCmat =gGeoManager->GetMedium("RPCgas");
   
  InitMedium("vacuum");
  TGeoMedium *vacuum =gGeoManager->GetMedium("vacuum");
    
  InitMedium("Bakelite");
  TGeoMedium *bakelite =gGeoManager->GetMedium("Bakelite");  
 
  InitMedium("iron");
  TGeoMedium *Iron =gGeoManager->GetMedium("iron");
 
  InitMedium("steel");
  TGeoMedium *Steel =gGeoManager->GetMedium("steel");
    
  InitMedium("copper");
  TGeoMedium *Cu =gGeoManager->GetMedium("copper");
    
  InitMedium("Concrete");
  TGeoMedium *Conc =gGeoManager->GetMedium("Concrete");
    
  TGeoUniformMagField *retFieldU    = new TGeoUniformMagField(0.,0.,-fField); //magnetic field up return yoke
  TGeoUniformMagField *retFieldL   = new TGeoUniformMagField(0.,0.,fField); //magnetic field low return yoke
    
  Double_t d = 0;
 
  if(fDesign!=3)
    {
      TGeoBBox *MSBox = new TGeoBBox("MagneticSpectrometerBox", fXRyoke/2, fYtot/2, fZtot/2);
      TGeoVolume *volMSBox = new TGeoVolume("volMagneticSpectrometer", MSBox, vacuum);
      tTauNuDet->AddNode(volMSBox, 1, new TGeoTranslation(0,10*cm,fZcenter));
  
      
      TGeoBBox *UpYokeBox = new TGeoBBox("UpYokeBox", fXRyoke/2, fYRyoke/2, fZRyoke/2);
      TGeoVolume *volUpYoke = new TGeoVolume("volUpYoke",UpYokeBox,vacuum);
      volMSBox->AddNode(volUpYoke,1,new TGeoTranslation(0,fYtot/2 - fYRyoke/2,0));
      volUpYoke->SetField(retFieldU);
    
    
      TGeoBBox *FeYoke = new TGeoBBox("FeYoke",fXtot/2, fYRyoke/2, fZArm/2);
      TGeoVolume *volFeYoke = new TGeoVolume("volFeYoke",FeYoke,Iron);
      volFeYoke->SetLineColor(kGray+1);
  
      TGeoBBox *FeYoke1 = new TGeoBBox("FeYoke1",fXtot/2, fYRyoke_s/2, fZRyoke_s/2);
      TGeoVolume *volFeYoke1 = new TGeoVolume("volFeYoke1",FeYoke1,Iron);
      volFeYoke1->SetLineColor(kGray+1);
    
      TGeoBBox *CoilContainer = new TGeoBBox("CoilContainer",fXtot/2, fCoilH/2, 40*cm);
      TGeoVolume *volCoilContainer = new TGeoVolume("volCoilContainer",CoilContainer,vacuum);
    
      TGeoBBox *Coil = new TGeoBBox("Coil",fXtot/2, fCoilH/2, fCoilW/2);
      TGeoVolume *volCoil = new TGeoVolume("volCoil",Coil,Cu);
      volCoil->SetLineColor(kOrange -5);
      for(int i = 0; i < fNCoil; i++)
    {
      volCoilContainer->AddNode(volCoil, i, new TGeoTranslation(0,0, -40*cm + fCoilW/2 + i*(fCoilGap + fCoilW)));
    }
    
      //vertical coils
      TGeoBBox *CoilV = new TGeoBBox("CoilV",fCoilH/2, fYRyoke/2 , fCoilW/2);
      TGeoVolume *volCoilV = new TGeoVolume("volCoilV",CoilV,Cu);
      volCoilV->SetLineColor(kOrange -5);
      for(int i = 0; i < fNCoil; i++)
    {
      volUpYoke->AddNode(volCoilV, i, new TGeoTranslation(fXRyoke/2 - fCoilH/2,0, -40*cm + fCoilW/2 + i*(fCoilGap + fCoilW)));
    }
      for(int i = 0; i < fNCoil; i++)
    {
      volUpYoke->AddNode(volCoilV, i, new TGeoTranslation(-fXRyoke/2 + fCoilH/2,0, -40*cm + fCoilW/2 + i*(fCoilGap + fCoilW)));
    }
 
      //cout <<"fZArm: " << fZArm<< endl;
    
      volUpYoke->AddNode(volFeYoke,1, new TGeoTranslation(0,0,- (fZArm + fGapMiddle)/2));
      volUpYoke->AddNode(volFeYoke,2, new TGeoTranslation(0,0,(fZArm + fGapMiddle)/2));
      volUpYoke->AddNode(volFeYoke1,1,new TGeoTranslation(0,0,0));
      volUpYoke->AddNode(volCoilContainer,1,new TGeoTranslation(0,fYRyoke/2 - fCoilH/2,0)); //up
      volUpYoke->AddNode(volCoilContainer,2,new TGeoTranslation(0,-fYRyoke/2 + fCoilH/2,0)); //low
    
      TGeoBBox *LowYokeBox = new TGeoBBox("LowYokeBox", fXRyoke/2, fYRyoke/2, fZRyoke/2);
      TGeoVolume *volLowYoke = new TGeoVolume("volLowYoke",LowYokeBox,vacuum);
      volMSBox->AddNode(volLowYoke,1,new TGeoTranslation(0,-fYtot/2 + fYRyoke/2,0));
      volLowYoke->SetField(retFieldL);
   
      //vertical coils
      for(int i = 0; i < fNCoil; i++)
    {
      volLowYoke->AddNode(volCoilV, i, new TGeoTranslation(fXRyoke/2 - fCoilH/2,0, -40*cm + fCoilW/2 + i*(fCoilGap + fCoilW)));
    }
      for(int i = 0; i < fNCoil; i++)
    {
      volLowYoke->AddNode(volCoilV, i, new TGeoTranslation(-fXRyoke/2 + fCoilH/2,0, -40*cm + fCoilW/2 + i*(fCoilGap + fCoilW)));
    }
    
      volLowYoke->AddNode(volFeYoke,3, new TGeoTranslation(0,0,- (fZArm + fGapMiddle)/2));
      volLowYoke->AddNode(volFeYoke,4, new TGeoTranslation(0,0,(fZArm + fGapMiddle)/2));
      volLowYoke->AddNode(volFeYoke1,1,new TGeoTranslation(0,0,0));
      volLowYoke->AddNode(volCoilContainer,3,new TGeoTranslation(0,fYRyoke/2- fCoilH/2,0)); //up
      volLowYoke->AddNode(volCoilContainer,4,new TGeoTranslation(0,-fYRyoke/2 + fCoilH/2,0)); //low
 
      Int_t ArmNumber = 1;
      TGeoBBox *Arm1Box = new TGeoBBox("Arm1MSBox", fXFe/2, fYFe/2, fZArm/2);
      TGeoVolume *volArm1 = new TGeoVolume("volArm1MS", Arm1Box,vacuum);
      TGeoUniformMagField *magField1 = new TGeoUniformMagField(0.,-fField,0.); //magnetic field arm1
      volArm1->SetField(magField1);
      volMSBox ->AddNode(volArm1,ArmNumber,new TGeoTranslation(0,0,-(fGapMiddle+fZArm)/2));
    
      Int_t nr =  ArmNumber*1E4;
 
      TGeoBBox *IronLayer = new TGeoBBox("Iron",fXFe/2, fYFe/2, fZFe/2);
      TGeoVolume *volIron = new TGeoVolume("volIron",IronLayer,Iron);
      //volIron->SetField(magField1);
 
      for(Int_t i = 0; i < fNFe; i++)
    {
      volArm1->AddNode(volIron,nr + 100 + i, new TGeoTranslation(0, 0, -fZArm/2+i*(fZFe +fZRpc) +fZFe/2));
    }
    
      TGeoBBox *RpcContainer = new TGeoBBox("RpcContainer", fXRpc/2, fYRpc/2, fZRpc/2);
      TGeoVolume *volRpcContainer = new TGeoVolume("volRpcContainer",RpcContainer,vacuum);
  
      TGeoBBox *Strip = new TGeoBBox("Strip",fXStrip/2, fYStrip/2, fZStrip/2);
      TGeoVolume *volStrip = new TGeoVolume("volStrip",Strip,Cu);
      volStrip->SetLineColor(kRed);
      volRpcContainer->AddNode(volStrip,1,new TGeoTranslation (0,0,-3.25*mm));
      volRpcContainer->AddNode(volStrip,2,new TGeoTranslation (0,0,3.25*mm));
      TGeoBBox *PETinsulator = new TGeoBBox("PETinsulator", fXPet/2, fYPet/2, fZPet/2);
      TGeoVolume *volPETinsulator = new TGeoVolume("volPETinsulator", PETinsulator, bakelite);
      volPETinsulator->SetLineColor(kYellow);
      volRpcContainer->AddNode(volPETinsulator,1,new TGeoTranslation(0,0,-3.1*mm));
      volRpcContainer->AddNode(volPETinsulator,2,new TGeoTranslation(0,0, 3.1*mm));
      TGeoBBox *Electrode = new TGeoBBox("Electrode",fXEle/2, fYEle/2, fZEle/2);
      TGeoVolume *volElectrode = new TGeoVolume("volElectrode",Electrode,bakelite);
      volElectrode->SetLineColor(kGreen);
      volRpcContainer->AddNode(volElectrode,1,new TGeoTranslation(0,0,-2*mm));
      volRpcContainer->AddNode(volElectrode,2,new TGeoTranslation(0,0, 2*mm));
      TGeoBBox *RpcGas = new TGeoBBox("RpcGas", fXGas/2, fYGas/2, fZGas/2);
      TGeoVolume *volRpc = new TGeoVolume("volRpc",RpcGas,RPCmat);
      volRpc->SetLineColor(kCyan);
      volRpcContainer->AddNode(volRpc,1,new TGeoTranslation(0,0,0));
   
      AddSensitiveVolume(volRpc);
    
      for(Int_t i = 0; i < fNRpc; i++)
    {
      volArm1->AddNode(volRpcContainer,nr + i,new TGeoTranslation(0, -fYFe/2 + fYRpc/2, -fZArm/2+(i+1)*fZFe + i*fZRpc +fZRpc/2));
    }
    
      ArmNumber = 2;
      nr =  ArmNumber*1E4;
 
      TGeoBBox *Arm2Box = new TGeoBBox("Arm2MSBox",fXFe/2, fYFe/2, fZArm/2);
      TGeoVolume *volArm2 = new TGeoVolume("volArm2MS", Arm2Box,vacuum);
      TGeoUniformMagField *magField2 = new TGeoUniformMagField(0.,fField,0.); //magnetic field arm2
      volArm2->SetField(magField2);
      volMSBox ->AddNode(volArm2,1,new TGeoTranslation(0,0,(fGapMiddle+fZArm)/2));
            TGeoVolume *volIron2 = new TGeoVolume("volIron2",IronLayer,Iron);
 
     //different volumes for second arm
      
      TGeoVolume *volRpcContainer2 = new TGeoVolume("volRpcContainer2",RpcContainer,vacuum);
      TGeoVolume *volStrip2 = new TGeoVolume("volStrip2",Strip,Cu);
      
      volStrip2->SetLineColor(kRed);
      volRpcContainer2->AddNode(volStrip2,1,new TGeoTranslation (0,0,-3.25*mm));
      volRpcContainer2->AddNode(volStrip2,2,new TGeoTranslation (0,0,3.25*mm));
 
      TGeoVolume *volPETinsulator2 = new TGeoVolume("volPETinsulator2", PETinsulator, bakelite);
      volPETinsulator2->SetLineColor(kYellow);
      volRpcContainer2->AddNode(volPETinsulator2,1,new TGeoTranslation(0,0,-3.1*mm));
      volRpcContainer2->AddNode(volPETinsulator2,2,new TGeoTranslation(0,0, 3.1*mm));
 
      TGeoVolume *volElectrode2 = new TGeoVolume("volElectrode2",Electrode,bakelite);
      volElectrode2->SetLineColor(kGreen);
      volRpcContainer2->AddNode(volElectrode2,1,new TGeoTranslation(0,0,-2*mm));
      volRpcContainer2->AddNode(volElectrode2,2,new TGeoTranslation(0,0, 2*mm));
   
      TGeoVolume *volRpc2 = new TGeoVolume("volRpc2",RpcGas,RPCmat);
      volRpc2->SetLineColor(kCyan);
      volRpcContainer2->AddNode(volRpc2,1,new TGeoTranslation(0,0,0));
      AddSensitiveVolume(volRpc2);
   
      for(Int_t i = 0; i < fNFe; i++)
    {
      volArm2->AddNode(volIron2,nr + 100 + i,new TGeoTranslation(0, 0, -fZArm/2+i*(fZFe +fZRpc) +fZFe/2));
    }
    
      for(Int_t i = 0; i < fNRpc; i++)
    {
      volArm2->AddNode(volRpcContainer2, nr + i,new TGeoTranslation(0, -fYFe/2 + fYRpc/2, -fZArm/2+(i+1)*fZFe + i*fZRpc +fZRpc/2));
    }
    
      //10 cm of Concrete on which the whole Magnetic Spectrometer volume (HPT included) will be placed
      TGeoBBox *Base = new TGeoBBox("Base", fXtot/2, 10*cm/2, fZtot/2);
      TGeoVolume *volBase = new TGeoVolume("volBase",Base,Conc);
      volBase->SetLineColor(kYellow-3);
 
      tTauNuDet->AddNode(volBase,1, new TGeoTranslation(0,-fYtot/2 + 10*cm/2,fZcenter));
 
 
      TGeoBBox *Pillar1Box = new TGeoBBox(fPillarX/2,fPillarY/2, fPillarZ/2);
      TGeoVolume *Pillar1Vol = new TGeoVolume("Pillar1Vol",Pillar1Box,Steel);
      Pillar1Vol->SetLineColor(kGreen+3);
 
      tTauNuDet->AddNode(Pillar1Vol,1, new TGeoTranslation(-fXtot/2+fPillarX/2,-fYtot/2-fPillarY/2,fZcenter-fZArm/2 - fGapMiddle/2 +fPillarZ/2));
      tTauNuDet->AddNode(Pillar1Vol,2, new TGeoTranslation(fXtot/2-fPillarX/2,-fYtot/2-fPillarY/2,fZcenter-fZArm/2 - fGapMiddle/2 +fPillarZ/2));
      tTauNuDet->AddNode(Pillar1Vol,3, new TGeoTranslation(-fXtot/2+fPillarX/2,-fYtot/2-fPillarY/2,fZcenter+fZArm/2+fGapMiddle/2-fPillarZ/2));
      tTauNuDet->AddNode(Pillar1Vol,4, new TGeoTranslation(fXtot/2-fPillarX/2,-fYtot/2-fPillarY/2,fZcenter+fZArm/2+fGapMiddle/2-fPillarZ/2));
    }
  if(fDesign==3)
    {
      Int_t nr = 1E4;
      
      TGeoBBox *MSBox = new TGeoBBox("MagneticSpectrometerBox", fXtot/2, fYtot/2, fZtot/2);
      TGeoVolume *volMSBox = new TGeoVolume("volMagneticSpectrometer", MSBox, vacuum);
      tTauNuDet->AddNode(volMSBox, 1, new TGeoTranslation(0,0,fZcenter));
 
      TGeoBBox *IronLayer = new TGeoBBox("Iron",fXFe/2, fYFe/2, fZFe/2);
      TGeoVolume *volIron = new TGeoVolume("volIron",IronLayer,Iron);
      //volIron->SetField(magField1);
 
      for(Int_t i = 0; i < fNFe; i++)
    {
      volMSBox->AddNode(volIron,nr + 100 + i, new TGeoTranslation(0, 0,-fZtot/2+i*fZFe+fZFe/2+i*fZRpc));
    }
 
       TGeoBBox *RpcContainer = new TGeoBBox("RpcContainer", fXRpc/2, fYRpc/2, fZRpc/2);
      TGeoVolume *volRpcContainer = new TGeoVolume("volRpcContainer",RpcContainer,vacuum);
  
      TGeoBBox *Strip = new TGeoBBox("Strip",fXStrip/2, fYStrip/2, fZStrip/2);
      TGeoVolume *volStrip = new TGeoVolume("volStrip",Strip,Cu);
      volStrip->SetLineColor(kRed);
      volRpcContainer->AddNode(volStrip,1,new TGeoTranslation (0,0,-3.25*mm));
      volRpcContainer->AddNode(volStrip,2,new TGeoTranslation (0,0,3.25*mm));
      TGeoBBox *PETinsulator = new TGeoBBox("PETinsulator", fXPet/2, fYPet/2, fZPet/2);
      TGeoVolume *volPETinsulator = new TGeoVolume("volPETinsulator", PETinsulator, bakelite);
      volPETinsulator->SetLineColor(kYellow);
      volRpcContainer->AddNode(volPETinsulator,1,new TGeoTranslation(0,0,-3.1*mm));
      volRpcContainer->AddNode(volPETinsulator,2,new TGeoTranslation(0,0, 3.1*mm));
      TGeoBBox *Electrode = new TGeoBBox("Electrode",fXEle/2, fYEle/2, fZEle/2);
      TGeoVolume *volElectrode = new TGeoVolume("volElectrode",Electrode,bakelite);
      volElectrode->SetLineColor(kGreen);
      volRpcContainer->AddNode(volElectrode,1,new TGeoTranslation(0,0,-2*mm));
      volRpcContainer->AddNode(volElectrode,2,new TGeoTranslation(0,0, 2*mm));
      TGeoBBox *RpcGas = new TGeoBBox("RpcGas", fXGas/2, fYGas/2, fZGas/2);
      TGeoVolume *volRpc = new TGeoVolume("volRpc",RpcGas,RPCmat);
      volRpc->SetLineColor(kCyan);
      volRpcContainer->AddNode(volRpc,1,new TGeoTranslation(0,0,0));
   
      AddSensitiveVolume(volRpc);
    
      for(Int_t i = 0; i < fNRpc; i++)
    {
      volMSBox->AddNode(volRpcContainer,nr + i,new TGeoTranslation(0, 0, -fZtot/2+(i+1)*fZFe + i*fZRpc +fZRpc/2));
    }
      
      TGeoBBox *Pillar1Box = new TGeoBBox(fPillarX/2,fPillarY/2, fPillarZ/2);
      TGeoVolume *Pillar1Vol = new TGeoVolume("Pillar1Vol",Pillar1Box,Steel);
      Pillar1Vol->SetLineColor(kGreen+3);
 
      tTauNuDet->AddNode(Pillar1Vol,1, new TGeoTranslation(-fXtot/2+fPillarX/2,-fYtot/2-fPillarY/2,fZcenter-fZtot/2+fPillarZ/2));
      tTauNuDet->AddNode(Pillar1Vol,2, new TGeoTranslation(fXtot/2-fPillarX/2,-fYtot/2-fPillarY/2,fZcenter-fZtot/2 +fPillarZ/2));
      //      tTauNuDet->AddNode(Pillar1Vol,3, new TGeoTranslation(-fXtot/2+fPillarX/2,-fYtot/2-fPillarY/2,fZcenter+fZtot/2-fPillarZ/2));
      //tTauNuDet->AddNode(Pillar1Vol,4, new TGeoTranslation(fXtot/2-fPillarX/2,-fYtot/2-fPillarY/2,fZcenter+fZtot/2-fPillarZ/2));
    }
 
}

CopyClones()

virtual void MagneticSpectrometer::CopyClones ( TClonesArray *  cl1, TClonesArray *  cl2, Int_t  offset  )

inlinevirtual

The following methods can be implemented if you need to make any optional action in your detector during the transport.

Definition at line 76 of file MagneticSpectrometer.h.

                                            {;}

DecodeVolumeID()

void MagneticSpectrometer::DecodeVolumeID	(	Int_t	detID,
		int &	nARM,
		int &	nRPC
	)

Definition at line 579 of file MagneticSpectrometer.cxx.

{
  nARM =  detID/1E4;
  nRPC =  detID - nARM*1E4;
}

EndOfEvent()

void MagneticSpectrometer::EndOfEvent ( )

virtual

Definition at line 558 of file MagneticSpectrometer.cxx.

{
  fShipRpcPointCollection->Clear();
}

FinishPrimary()

virtual void MagneticSpectrometer::FinishPrimary ( )

inlinevirtual

Definition at line 80 of file MagneticSpectrometer.h.

{;}

FinishRun()

virtual void MagneticSpectrometer::FinishRun ( )

inlinevirtual

Definition at line 81 of file MagneticSpectrometer.h.

{;}

GetCollection()

TClonesArray * MagneticSpectrometer::GetCollection ( Int_t  iColl ) const

virtual

Gets the produced collections

Definition at line 585 of file MagneticSpectrometer.cxx.

{
  if (iColl == 0) { return fShipRpcPointCollection; }
  else { return NULL; }
}

Initialize()

void MagneticSpectrometer::Initialize ( )

virtual

Initialization of the detector is done here

Definition at line 198 of file MagneticSpectrometer.cxx.

{
  FairDetector::Initialize();
}

InitMedium()

Int_t MagneticSpectrometer::InitMedium ( const char *  name )

private

Definition at line 204 of file MagneticSpectrometer.cxx.

{
  static FairGeoLoader *geoLoad=FairGeoLoader::Instance();
  static FairGeoInterface *geoFace=geoLoad->getGeoInterface();
  static FairGeoMedia *media=geoFace->getMedia();
  static FairGeoBuilder *geoBuild=geoLoad->getGeoBuilder();
 
  FairGeoMedium *ShipMedium=media->getMedium(name);
 
  if (!ShipMedium)
    {
      Fatal("InitMedium","Material %s not defined in media file.", name);
      return -1111;
    }
  TGeoMedium* medium=gGeoManager->GetMedium(name);
  if (medium!=NULL)
    return ShipMedium->getMediumIndex();
  return geoBuild->createMedium(ShipMedium);
}

operator=()

MagneticSpectrometer & MagneticSpectrometer::operator= ( const MagneticSpectrometer &  )

protected

PostTrack()

virtual void MagneticSpectrometer::PostTrack ( )

inlinevirtual

Definition at line 83 of file MagneticSpectrometer.h.

{;}

PreTrack()

virtual void MagneticSpectrometer::PreTrack ( )

inlinevirtual

Definition at line 84 of file MagneticSpectrometer.h.

{;}

ProcessHits()

Bool_t MagneticSpectrometer::ProcessHits ( FairVolume *  v = 0 )

virtual

this method is called for each step during simulation (see FairMCApplication::Stepping())

This method is called from the MC stepping

Definition at line 499 of file MagneticSpectrometer.cxx.

{
  //Set parameters at entrance of volume. Reset ELoss.
  if ( gMC->IsTrackEntering() ) {
    fELoss  = 0.;
    fTime   = gMC->TrackTime() * 1.0e09;
    fLength = gMC->TrackLength();
    gMC->TrackPosition(fPos);
    gMC->TrackMomentum(fMom);
  }
  // Sum energy loss for all steps in the active volume
  fELoss += gMC->Edep();
    
  // Create muonPoint at exit of active volume
  if ( gMC->IsTrackExiting()    ||
       gMC->IsTrackStop()       ||
       gMC->IsTrackDisappeared()   ) {
    fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
    fVolumeID = vol->getMCid();
    if (fELoss == 0. ) { return kFALSE; }
    TParticle* p=gMC->GetStack()->GetCurrentTrack();
    Int_t pdgCode = p->GetPdgCode();
    Int_t detID=0;
    gMC->CurrentVolID(detID);
 
    // cout<< "detID = " << detID << endl;
    Int_t MaxLevel = gGeoManager->GetLevel();
    const Int_t MaxL = MaxLevel;
    //cout << "MaxLevel = " << MaxL << endl;
    //cout << gMC->CurrentVolPath()<< endl;
    Int_t NRpc =0;
    const char *name;
    name = gMC->CurrentVolName();
    //cout << name << endl;
    Int_t motherID = gGeoManager->GetMother(0)->GetNumber();
    const char *mumname = gMC->CurrentVolOffName(0);
    //cout<<mumname<<"   "<< motherID<<endl;
    detID = motherID;
    //cout<< "detID = " << detID << endl;
    //cout<<endl;
    fVolumeID = detID;
 
    TLorentzVector Pos; 
    gMC->TrackPosition(Pos); 
    Double_t xmean = (fPos.X()+Pos.X())/2. ;      
    Double_t ymean = (fPos.Y()+Pos.Y())/2. ;      
    Double_t zmean = (fPos.Z()+Pos.Z())/2. ;     
     
    AddHit(fTrackID, fVolumeID, TVector3(xmean, ymean,  zmean), TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength,fELoss, pdgCode);
        
    // Increment number of muon det points in TParticle
    ShipStack* stack = (ShipStack*) gMC->GetStack();
    stack->AddPoint(ktauRpc);
  }
    
  return kTRUE;
}

Register()

void MagneticSpectrometer::Register ( )

virtual

Registers the produced collections in FAIRRootManager.

This will create a branch in the output tree called ShipRpcPoint, setting the last parameter to kFALSE means: this collection will not be written to the file, it will exist only during the simulation.

Definition at line 564 of file MagneticSpectrometer.cxx.

{
    
  FairRootManager::Instance()->Register("ShipRpcPoint", "MagneticSpectrometer",
                    fShipRpcPointCollection, kTRUE);
}

Reset()

void MagneticSpectrometer::Reset ( )

virtual

has to be called after each event to reset the containers

Definition at line 591 of file MagneticSpectrometer.cxx.

{
  fShipRpcPointCollection->Clear();
}

SetCoilParameters()

void MagneticSpectrometer::SetCoilParameters	(	Double_t	CoilH,
		Double_t	CoilW,
		Int_t	N,
		Double_t	CoilG
	)

Definition at line 173 of file MagneticSpectrometer.cxx.

{
  fCoilH = CoilH;
  fCoilW = CoilW;
  fCoilGap = CoilG;
  fNCoil = N;
}

SetDesign()

void MagneticSpectrometer::SetDesign

(

Int_t

Design

)

Definition at line 74 of file MagneticSpectrometer.cxx.

{  
  fDesign = Design;
  cout <<" Mag Spectro Design "<< fDesign<<endl;
}

SetFeDimensions()

void MagneticSpectrometer::SetFeDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 87 of file MagneticSpectrometer.cxx.

{  
  fXFe = X;
  fYFe = Y;
  fZFe = Z;
}

SetGapDownstream()

void MagneticSpectrometer::SetGapDownstream

(

Double_t

Gap

)

Definition at line 148 of file MagneticSpectrometer.cxx.

{
  fGapDown = Gap;
}

SetGapMiddle()

void MagneticSpectrometer::SetGapMiddle

(

Double_t

Gap

)

Definition at line 153 of file MagneticSpectrometer.cxx.

{
  fGapMiddle = Gap;
}

SetMagneticField()

void MagneticSpectrometer::SetMagneticField

(

Double_t

B

)

Definition at line 143 of file MagneticSpectrometer.cxx.

{
  fField =B;
}

SetNFeInArm()

void MagneticSpectrometer::SetNFeInArm

(

Int_t

N

)

Definition at line 163 of file MagneticSpectrometer.cxx.

{
  fNFe = N;
}

SetNRpcInArm()

void MagneticSpectrometer::SetNRpcInArm

(

Int_t

N

)

Definition at line 168 of file MagneticSpectrometer.cxx.

{
  fNRpc = N;
}

SetPillarDimensions()

void MagneticSpectrometer::SetPillarDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 182 of file MagneticSpectrometer.cxx.

{
  fPillarX = X;
  fPillarY = Y;
  fPillarZ = Z;
}

SetReturnYokeDimensions()

void MagneticSpectrometer::SetReturnYokeDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 129 of file MagneticSpectrometer.cxx.

{  
  fXRyoke = X;
  fYRyoke = Y;
  fZRyoke = Z;
}

SetRpcDimensions()

void MagneticSpectrometer::SetRpcDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 94 of file MagneticSpectrometer.cxx.

{  
  fXRpc = X;
  fYRpc = Y;
  fZRpc = Z;
}

SetRpcElectrodeDimensions()

void MagneticSpectrometer::SetRpcElectrodeDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 115 of file MagneticSpectrometer.cxx.

{  
  fXEle = X;
  fYEle = Y;
  fZEle = Z;
}

SetRpcGasDimensions()

void MagneticSpectrometer::SetRpcGasDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 101 of file MagneticSpectrometer.cxx.

{  
  fXGas = X;
  fYGas = Y;
  fZGas = Z;
}

SetRpcPETDimensions()

void MagneticSpectrometer::SetRpcPETDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 108 of file MagneticSpectrometer.cxx.

{  
  fXPet = X;
  fYPet = Y;
  fZPet = Z;
}

SetRpcStripDimensions()

void MagneticSpectrometer::SetRpcStripDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 122 of file MagneticSpectrometer.cxx.

{  
  fXStrip = X;
  fYStrip = Y;
  fZStrip = Z;
}

SetSmallerYokeDimensions()

void MagneticSpectrometer::SetSmallerYokeDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 136 of file MagneticSpectrometer.cxx.

{  
  fXRyoke_s = X;
  fYRyoke_s = Y;
  fZRyoke_s = Z;
}

SetSpecialPhysicsCuts()

virtual void MagneticSpectrometer::SetSpecialPhysicsCuts ( )

inlinevirtual

Definition at line 78 of file MagneticSpectrometer.h.

{;}

SetTotDimensions()

void MagneticSpectrometer::SetTotDimensions	(	Double_t	X,
		Double_t	Y,
		Double_t	Z
	)

Definition at line 80 of file MagneticSpectrometer.cxx.

{  
  fXtot = X;
  fYtot = Y;
  fZtot = Z;
}

SetZDimensionArm()

void MagneticSpectrometer::SetZDimensionArm

(

Double_t

Z

)

Definition at line 158 of file MagneticSpectrometer.cxx.

{
  fZArm = Z;
}

Member Data Documentation

fCoilGap

Double_t MagneticSpectrometer::fCoilGap

protected

Definition at line 138 of file MagneticSpectrometer.h.

fCoilH

Double_t MagneticSpectrometer::fCoilH

protected

Definition at line 139 of file MagneticSpectrometer.h.

fCoilW

Double_t MagneticSpectrometer::fCoilW

protected

Definition at line 140 of file MagneticSpectrometer.h.

fDesign

Int_t MagneticSpectrometer::fDesign

protected

Definition at line 112 of file MagneticSpectrometer.h.

fELoss

Double32_t MagneticSpectrometer::fELoss

private

length

Definition at line 101 of file MagneticSpectrometer.h.

fField

Double_t MagneticSpectrometer::fField

protected

Definition at line 128 of file MagneticSpectrometer.h.

fGapDown

Double_t MagneticSpectrometer::fGapDown

protected

Definition at line 126 of file MagneticSpectrometer.h.

fGapMiddle

Double_t MagneticSpectrometer::fGapMiddle

protected

Definition at line 127 of file MagneticSpectrometer.h.

fLength

Double32_t MagneticSpectrometer::fLength

private

time

Definition at line 100 of file MagneticSpectrometer.h.

fMom

TLorentzVector MagneticSpectrometer::fMom

private

position at entrance

Definition at line 98 of file MagneticSpectrometer.h.

fNCoil

Int_t MagneticSpectrometer::fNCoil

protected

Definition at line 141 of file MagneticSpectrometer.h.

fNFe

Int_t MagneticSpectrometer::fNFe

protected

Definition at line 117 of file MagneticSpectrometer.h.

fNRpc

Int_t MagneticSpectrometer::fNRpc

protected

Definition at line 118 of file MagneticSpectrometer.h.

fPdgCode

Int_t MagneticSpectrometer::fPdgCode

private

track index

Definition at line 95 of file MagneticSpectrometer.h.

fPillarX

Double_t MagneticSpectrometer::fPillarX

protected

Definition at line 158 of file MagneticSpectrometer.h.

fPillarY

Double_t MagneticSpectrometer::fPillarY

protected

Definition at line 159 of file MagneticSpectrometer.h.

fPillarZ

Double_t MagneticSpectrometer::fPillarZ

protected

Definition at line 160 of file MagneticSpectrometer.h.

fPos

TLorentzVector MagneticSpectrometer::fPos

private

volume id

Definition at line 97 of file MagneticSpectrometer.h.

fShipRpcPointCollection

TClonesArray* MagneticSpectrometer::fShipRpcPointCollection

private

energy loss

container for data points

Definition at line 104 of file MagneticSpectrometer.h.

fTime

Double32_t MagneticSpectrometer::fTime

private

momentum at entrance

Definition at line 99 of file MagneticSpectrometer.h.

fTrackID

Int_t MagneticSpectrometer::fTrackID

private

Track information to be stored until the track leaves the active volume.

Definition at line 94 of file MagneticSpectrometer.h.

fVolumeID

Int_t MagneticSpectrometer::fVolumeID

private

pdg code

Definition at line 96 of file MagneticSpectrometer.h.

fXEle

Double_t MagneticSpectrometer::fXEle

protected

Definition at line 150 of file MagneticSpectrometer.h.

fXFe

Double_t MagneticSpectrometer::fXFe

protected

Definition at line 119 of file MagneticSpectrometer.h.

fXGas

Double_t MagneticSpectrometer::fXGas

protected

Definition at line 153 of file MagneticSpectrometer.h.

fXPet

Double_t MagneticSpectrometer::fXPet

protected

Definition at line 147 of file MagneticSpectrometer.h.

fXRpc

Double_t MagneticSpectrometer::fXRpc

protected

Definition at line 120 of file MagneticSpectrometer.h.

fXRyoke

Double_t MagneticSpectrometer::fXRyoke

protected

Definition at line 129 of file MagneticSpectrometer.h.

fXRyoke_s

Double_t MagneticSpectrometer::fXRyoke_s

protected

Definition at line 134 of file MagneticSpectrometer.h.

fXStrip

Double_t MagneticSpectrometer::fXStrip

protected

Definition at line 144 of file MagneticSpectrometer.h.

fXtot

Double_t MagneticSpectrometer::fXtot

protected

Definition at line 114 of file MagneticSpectrometer.h.

fYEle

Double_t MagneticSpectrometer::fYEle

protected

Definition at line 151 of file MagneticSpectrometer.h.

fYFe

Double_t MagneticSpectrometer::fYFe

protected

Definition at line 121 of file MagneticSpectrometer.h.

fYGas

Double_t MagneticSpectrometer::fYGas

protected

Definition at line 154 of file MagneticSpectrometer.h.

fYPet

Double_t MagneticSpectrometer::fYPet

protected

Definition at line 148 of file MagneticSpectrometer.h.

fYRpc

Double_t MagneticSpectrometer::fYRpc

protected

Definition at line 122 of file MagneticSpectrometer.h.

fYRyoke

Double_t MagneticSpectrometer::fYRyoke

protected

Definition at line 130 of file MagneticSpectrometer.h.

fYRyoke_s

Double_t MagneticSpectrometer::fYRyoke_s

protected

Definition at line 135 of file MagneticSpectrometer.h.

fYStrip

Double_t MagneticSpectrometer::fYStrip

protected

Definition at line 145 of file MagneticSpectrometer.h.

fYtot

Double_t MagneticSpectrometer::fYtot

protected

Definition at line 115 of file MagneticSpectrometer.h.

fZArm

Double_t MagneticSpectrometer::fZArm

protected

Definition at line 125 of file MagneticSpectrometer.h.

fZcenter

Double_t MagneticSpectrometer::fZcenter

protected

Definition at line 113 of file MagneticSpectrometer.h.

fZEle

Double_t MagneticSpectrometer::fZEle

protected

Definition at line 152 of file MagneticSpectrometer.h.

fZFe

Double_t MagneticSpectrometer::fZFe

protected

Definition at line 123 of file MagneticSpectrometer.h.

fZGas

Double_t MagneticSpectrometer::fZGas

protected

Definition at line 155 of file MagneticSpectrometer.h.

fZPet

Double_t MagneticSpectrometer::fZPet

protected

Definition at line 149 of file MagneticSpectrometer.h.

fZRpc

Double_t MagneticSpectrometer::fZRpc

protected

Definition at line 124 of file MagneticSpectrometer.h.

fZRyoke

Double_t MagneticSpectrometer::fZRyoke

protected

Definition at line 131 of file MagneticSpectrometer.h.

fZRyoke_s

Double_t MagneticSpectrometer::fZRyoke_s

protected

Definition at line 136 of file MagneticSpectrometer.h.

fZStrip

Double_t MagneticSpectrometer::fZStrip

protected

Definition at line 146 of file MagneticSpectrometer.h.

fZtot

Double_t MagneticSpectrometer::fZtot

protected

Definition at line 116 of file MagneticSpectrometer.h.

---

The documentation for this class was generated from the following files:

• nutaudet/MagneticSpectrometer.h
• nutaudet/MagneticSpectrometer.cxx