veto.h

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
// SPDX-FileCopyrightText: Copyright CERN for the benefit of the SHiP
// Collaboration
 
#ifndef VETO_VETO_H_
#define VETO_VETO_H_
 
#include <map>
#include <vector>
 
#include "Detector.h"
#include "TGeoVolume.h"
#include "TLorentzVector.h"
#include "TVector3.h"
 
class vetoPoint;
class FairVolume;
class TClonesArray;
 
class veto : public SHiP::Detector<vetoPoint> {
 public:
  veto();
 
  Bool_t ProcessHits(FairVolume* v = 0) override;
 
  void SetFastMuon() { fFastMuon = true; }  // kill all tracks except of muons
  void SetFollowMuon() {
    fFollowMuon = true;
  }  // make muon shield active to follow muons
  void SetVesselDimensions(Float_t x1, Float_t x2, Float_t y1, Float_t y2,
                           Float_t vesselstart) {
    VetoStartInnerX = x1;
    VetoEndInnerX = x2;
    VetoStartInnerY = y1;
    VetoEndInnerY = y2;
    zStartDecayVol = vesselstart;
  }
 
  void ConstructGeometry() override;
 
  void SetVesselStructure(Float_t a, Float_t b, Float_t c, TString d, Float_t l,
                          TString e, TString f, TString v, Float_t r) {
    f_InnerSupportThickness = a;
    f_VetoThickness = b;
    f_OuterSupportThickness = c;
    supportMedIn_name = d;
    f_LidThickness = l;
    vetoMed_name = e;
    supportMedOut_name = f;
    decayVolumeMed_name = v;
    f_RibThickness = r;
  }
 
  void PreTrack() override;
 
  inline void SetUseSupport(Int_t use = 1) { fUseSupport = use; }
  inline Int_t GetUseSupport() const { return fUseSupport; }
 
  inline void SetLiquidVeto(Int_t liquid = 1) { fLiquidVeto = liquid; }
  inline Int_t GetLiquidVeto() const { return fLiquidVeto; }
 
 private:
  Bool_t fFastMuon, fFollowMuon;
 
  Float_t f_InnerSupportThickness;
  Float_t f_OuterSupportThickness;
 
  Float_t f_LidThickness;
  Float_t f_VetoThickness;
  Float_t f_RibThickness;
 
  TString vetoMed_name;
  TString supportMedIn_name;
  TString supportMedOut_name;
  TString decayVolumeMed_name;
 
  TGeoMedium* vetoMed;
  TGeoMedium* supportMedIn;
  TGeoMedium* supportMedOut;
  TGeoMedium* decayVolumeMed;
 
  Float_t VetoStartInnerX;
  Float_t VetoStartInnerY;
  Float_t VetoEndInnerX;
  Float_t VetoEndInnerY;
 
  Float_t zStartDecayVol;
 
  Int_t fUseSupport;
  Int_t fLiquidVeto;
 
  veto(const veto&) = delete;
  veto& operator=(const veto&) = delete;
  TGeoVolume* GeoTrapezoid(TString xname, Double_t wz, Double_t wX_start,
                           Double_t wX_end, Double_t wY_start, Double_t wY_end,
                           Int_t color, TGeoMedium* material, Bool_t sens);
 
  TGeoVolume* GeoTrapezoidHollow(TString xname, Double_t thick, Double_t wz,
                                 Double_t wX_start, Double_t wX_end,
                                 Double_t wY_start, Double_t wY_end,
                                 Int_t color, TGeoMedium* material,
                                 Bool_t sens);
 
  void AddBlock(TGeoVolumeAssembly* tInnerWall,
                TGeoVolumeAssembly* tDecayVacuum,
                TGeoVolumeAssembly* tOuterWall, TGeoVolumeAssembly* tLongitRib,
                TGeoVolumeAssembly* tVerticalRib, TGeoVolumeAssembly* ttLiSc,
                int& liScCounter, int blockNr, int nx, int ny, double z1,
                double z2, double Zshift, double dist, double wallThick,
                double liscThick1, double liscThick2, double ribThick);
  TGeoVolumeAssembly* GeoCornerRib(TString xname, double ribThick, double lt1,
                                   double lt2, double dz, double slopeX,
                                   double slopeY, Int_t color,
                                   TGeoMedium* material, Bool_t sens);
 
  int makeId(double z, double x, double y);
  int liscId(TString ShapeTypeName, int blockNr, int Zlayer, int number,
             int position);
  double wx(double z);
  double wy(double z);
  TGeoVolume* GeoSideObj(TString xname, double dz, double a1, double b1,
                         double a2, double b2, double dA, double dB,
                         Int_t color, TGeoMedium* material, Bool_t sens);
  TGeoVolume* GeoCornerLiSc1(TString xname, double dz, bool isClockwise,
                             double a1, double a2, double b1, double b2,
                             double dA, double dB, Int_t color,
                             TGeoMedium* material, Bool_t sens);
  TGeoVolume* GeoCornerLiSc2(TString xname, double dz, bool isClockwise,
                             double a1, double a2, double b1, double b2,
                             double dA, double dB, Int_t color,
                             TGeoMedium* material, Bool_t sens);
 
  TGeoVolume* MakeSegments();
 
  ClassDefOverride(veto, 3)
};
 
#endif  // VETO_VETO_H_