MagneticSpectrometer.h

Go to the documentation of this file.

#ifndef MAGNETICSPECTROMETER_H
#define MAGNETICSPECTROMETER_H
 
#include "FairModule.h"                 // for FairModule
#include "FairDetector.h"                  // for FairDetector
 
#include "Rtypes.h"                     // for ShipMuonShield::Class, Bool_t, etc
 
#include <string>                       // for string
 
#include "TVector3.h"
#include "TLorentzVector.h"
 
class ShipRpcPoint;
class FairVolume;
class TClonesArray;
 
class MagneticSpectrometer:public FairDetector
{
  public:
  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);
    
private:
    
    Int_t          fTrackID;           
    Int_t          fPdgCode;           
    Int_t          fVolumeID;          
    TLorentzVector fPos;               
    TLorentzVector fMom;               
    Double32_t     fTime;              
    Double32_t     fLength;            
    Double32_t     fELoss;             
    
    TClonesArray*  fShipRpcPointCollection;
    
    Int_t InitMedium(const char* name);
    
    
    
protected:
 
    Int_t fDesign;
    Double_t fZcenter; //z distance of the center of the spectrometer in cm from the center of the vacuum tube
    Double_t fXtot;
    Double_t fYtot;
    Double_t fZtot; //Dimension of the whole magnetic spectrometr (1st + 2nd arm + HPTs) alogn beam axis
    Int_t fNFe;
    Int_t fNRpc;
    Double_t fXFe;
    Double_t fXRpc;
    Double_t fYFe;
    Double_t fYRpc;
    Double_t fZFe; // Width of the Iron Slabs
    Double_t fZRpc; // Width of the Rpc planes
    Double_t fZArm; // Width of the Spectrometer Arms
    Double_t fGapDown; //distance between the end of the second arm of the spectrometer and the decay vessel
    Double_t fGapMiddle; // distance between the two arms of the spectrometer
    Double_t fField;
    Double_t fXRyoke;
    Double_t fYRyoke;
    Double_t fZRyoke;
    
    //Dimensions of the smaller part of the yoke
    Double_t fXRyoke_s;
    Double_t fYRyoke_s;
    Double_t fZRyoke_s;
 
    Double_t fCoilGap;
    Double_t fCoilH;
    Double_t fCoilW;
    Int_t fNCoil;
    
    //Dimension for detailed RPC simulation:
    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;
 
    //Dimension of the pillars
    Double_t fPillarX;
    Double_t fPillarY;
    Double_t fPillarZ;
 
 
    MagneticSpectrometer(const MagneticSpectrometer&);
    MagneticSpectrometer& operator=(const MagneticSpectrometer&);
    ClassDef(MagneticSpectrometer,5)
 
};
 
 
 
#endif //TAUMAGNETICSPECTROMETER_H