#include <exitHadronAbsorber.h>

Inheritance diagram for exitHadronAbsorber:

[legend]

Collaboration diagram for exitHadronAbsorber:

[legend]

Public Member Functions
	exitHadronAbsorber (const char *Name, Bool_t Active)

	exitHadronAbsorber ()

void	Initialize () override

Bool_t	ProcessHits (FairVolume *v=0) override

void	Register () override

void	ConstructGeometry () override

void	FinishRun () override

void	PreTrack () override

void	SetEnergyCut (Float_t emax)

void	SetOnlyMuons ()

void	SetOpt4DP ()

void	SkipNeutrinos ()

void	SetZposition (Float_t x)

void	SetVetoPointName (TString name)

void	SetCylindricalPlane ()

void	SetUseCaveCoordinates ()

Public Member Functions inherited from SHiP::Detector< vetoPoint >
	Detector ()=default

	Detector (const char *Name, Bool_t Active, Int_t detID)

	Detector (const char *Name, Bool_t Active)

	~Detector () override=default

vetoPoint *	AddHit (Args &&... args)

void	ConstructGeometry () override=0

void	Initialize () override

void	Reset () override

void	EndOfEvent () override

void	Register () override

TClonesArray *	GetCollection (Int_t iColl) const override

void	UpdatePointTrackIndices (const std::map< Int_t, Int_t > &indexMap) override
	Update track indices in point collection after track filtering.

void	SetSpecialPhysicsCuts () override

void	FinishPrimary () override

void	FinishRun () override

void	BeginPrimary () override

void	PostTrack () override

void	PreTrack () override

void	BeginEvent () override

void	CopyClones (TClonesArray cl1, TClonesArray cl2, Int_t offset) override

Public Member Functions inherited from ISTLPointContainer
virtual void	UpdatePointTrackIndices (const std::map< Int_t, Int_t > &indexMap)=0
	Update track indices in point collection after track filtering.

virtual	~ISTLPointContainer ()=default

Private Attributes
Int_t	fUniqueID

Bool_t	fOnlyMuons

Bool_t	fSkipNeutrinos
	flag if only muons should be stored

TString	fVetoName
	flag if neutrinos should be ignored

Double_t	fzPos

Bool_t	withNtuple
	zPos, optional

TNtuple *	fNtuple
	special option for Dark Photon physics studies

Float_t	EMax

Bool_t	fCylindricalPlane
	max energy to transport

Bool_t	fUseCaveCoordinates
	cylindrical sensPlane flag

TFile *	fout
	set position from cave

TClonesArray *	fElectrons

Int_t	index

Additional Inherited Members
Protected Attributes inherited from SHiP::Detector< vetoPoint >
Int_t	fEventID

Int_t	fTrackID
	event index

Int_t	fVolumeID
	track index

TLorentzVector	fPos
	volume id

TLorentzVector	fMom
	position at entrance

Double_t	fTime
	momentum at entrance

Double_t	fLength
	time

Double_t	fELoss
	length

std::vector< vetoPoint > *	fDetPoints
	energy loss

TGeoVolume *	fDetector

Detailed Description

Definition at line 15 of file exitHadronAbsorber.h.

Constructor & Destructor Documentation

◆ exitHadronAbsorber() [1/2]

exitHadronAbsorber::exitHadronAbsorber	(	const char *	Name,
		Bool_t	Active
	)

Definition at line 37 of file exitHadronAbsorber.cxx.

    : Detector(Name, Active, kVETO),
      fOnlyMuons(kFALSE),
      fSkipNeutrinos(kFALSE),
      fVetoName("veto"),
      fzPos(3E8),
      withNtuple(kFALSE),
      fCylindricalPlane(kFALSE) {}

◆ exitHadronAbsorber() [2/2]

exitHadronAbsorber::exitHadronAbsorber ( )

Definition at line 46 of file exitHadronAbsorber.cxx.

    : Detector("exitHadronAbsorber", kTRUE, kVETO),
      fUniqueID(-1),
      fOnlyMuons(kFALSE),
      fSkipNeutrinos(kFALSE),
      fVetoName("veto"),
      fzPos(3E8),
      withNtuple(kFALSE),
      fCylindricalPlane(kFALSE),
      fUseCaveCoordinates(kFALSE) {}

Member Function Documentation

◆ ConstructGeometry()

void exitHadronAbsorber::ConstructGeometry ( )

overridevirtual

Create the detector geometry

Implements SHiP::Detector< vetoPoint >.

Definition at line 255 of file exitHadronAbsorber.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("vacuums");
  TGeoMedium* vac = gGeoManager->GetMedium("vacuums");
  if (vac == nullptr) geoBuild->createMedium(ShipMedium);
  vac = gGeoManager->GetMedium("vacuums");
  gGeoManager->GetTopVolume();
  TGeoNavigator* nav = gGeoManager->GetCurrentNavigator();
  Double_t zLoc;
  if (fzPos > 1E8) {
    // Add thin sensitive plane after hadron absorber
    Float_t distance = 1.;
    Double_t local[3] = {0, 0, 0};
    if (!nav->cd("/MuonShieldArea_1/AbsorberVol_1")) {
      nav->cd("/MuonShieldArea_1/MagnAbsorb_MagRetL_1");
      distance = -1.;
    }
    TGeoBBox* tmp =
        dynamic_cast<TGeoBBox*>(nav->GetCurrentNode()->GetVolume()->GetShape());
    local[2] = distance * tmp->GetDZ();
    Double_t global[3] = {0, 0, 0};
    nav->LocalToMaster(local, global);
    zLoc = global[2] + distance * 1. * cm;
  } else {
    zLoc = fzPos;
  }  // use external input
 
  TString myname(this->GetName());
  TString shapename_prefix(myname);  // Use a prefix for shapenames
 
  Double_t xLocPlane = 0.0;
  Double_t yLocPlane = 0.0;
  Double_t zLocPlane = zLoc;
 
  if (!fCylindricalPlane) {
    TGeoVolume* sensPlane =
        gGeoManager->MakeBox(shapename_prefix + fVetoName + "_box", vac,
                             3.56 * m - 1. * mm, 1.7 * m - 1. * mm, 1. * mm);
    std::cout << this->GetName()
              << ", ConstructGeometry(): Created Box with dimensions: "
                 "3.56*m-1.*mm,1.7*m-1.*mm,1.*mm"
              << std::endl;
 
    if (!fUseCaveCoordinates) {
      nav->cd("/MuonShieldArea_1/");
    } else {
      Double_t local[3] = {0, 0, zLoc};
      Double_t global[3];
      nav->cd("/target_vacuum_box_1/TargetArea_1");
      nav->LocalToMaster(local, global);
      nav->cd("/cave_1");
      yLocPlane = global[1];
      zLocPlane = global[2];
      xLocPlane = global[0];
    }
    sensPlane->SetLineColor(kBlue - 10);
    nav->GetCurrentNode()->GetVolume()->AddNode(
        sensPlane, 1, new TGeoTranslation(xLocPlane, yLocPlane, zLocPlane));
    AddSensitiveVolume(sensPlane);
  } else {  // add cylindrical sensPlane
    TGeoVolume* sensPlaneCyl =
        gGeoManager->MakeTube(shapename_prefix + fVetoName + "_tube", vac,
                              12.51 * cm, 12.52 * cm, zLoc / 2.);
    std::cout << this->GetName()
              << ", ConstructGeometry(): Created Tube with dimensions: "
                 "12.51*cm,12.52*cm,zLoc/2."
              << std::endl;
    nav->cd("/target_vacuum_box_1/TargetArea_1/HeVolume_1");
    nav->GetCurrentNode()->GetVolume()->AddNode(sensPlaneCyl, 1,
                                                new TGeoTranslation(0, 0, 0));
    AddSensitiveVolume(sensPlaneCyl);
  }
}

◆ FinishRun()

void exitHadronAbsorber::FinishRun ( )

override

Definition at line 206 of file exitHadronAbsorber.cxx.

                                   {
  for (Int_t idnu = 11; idnu < 23; idnu += 1) {
    // nu or anti-nu
    for (Int_t idadd = -1; idadd < 3; idadd += 2) {
      Int_t idw = idnu;
      if (idnu == 18) {
        idw = 22;
      }
      if (idnu == 19) {
        idw = 111;
      }
      if (idnu == 20) {
        idw = 221;
      }
      if (idnu == 21) {
        idw = 223;
      }
      if (idnu == 22) {
        idw = 331;
      }
      Int_t idhnu = 10000 + idw;
      if (idadd == -1) {
        if (idnu > 17) {
          continue;
        }
        idhnu += 10000;
        idw = -idnu;
      }
      TString key = fVetoName;
      key += idhnu;
      TSeqCollection* fileList = gROOT->GetListOfFiles();
      dynamic_cast<TFile*>(fileList->At(0))->cd();
      TH1D* Hidhnu = dynamic_cast<TH1D*>(fout->Get(key));
      Hidhnu->Write();
      key = fVetoName;
      key += idhnu + 1000;
      TH2D* Hidhnu100 = dynamic_cast<TH2D*>(fout->Get(key));
      Hidhnu100->Write();
      key = fVetoName;
      key += idhnu + 2000;
      TH2D* Hidhnu200 = dynamic_cast<TH2D*>(fout->Get(key));
      Hidhnu200->Write();
    }
  }
  if (withNtuple) {
    fNtuple->Write();
  }
}

◆ Initialize()

void exitHadronAbsorber::Initialize ( )

override

Definition at line 86 of file exitHadronAbsorber.cxx.

                                    {
  FairDetector::Initialize();
  TSeqCollection* fileList = gROOT->GetListOfFiles();
  fout = dynamic_cast<TFile*>(fileList->At(0));
  // book hists for Genie neutrino momentum distribution
  // add also leptons, and photon
  // add pi0 111 eta 221 eta' 331  omega 223 for DM production
  TDatabasePDG* PDG = TDatabasePDG::Instance();
  for (Int_t idnu = 11; idnu < 26; idnu += 1) {
    // nu or anti-nu
    for (Int_t idadd = -1; idadd < 3; idadd += 2) {
      Int_t idw = idnu;
      if (idnu == 18) {
        idw = 22;
      }
      if (idnu == 19) {
        idw = 111;
      }
      if (idnu == 20) {
        idw = 221;
      }
      if (idnu == 21) {
        idw = 223;
      }
      if (idnu == 22) {
        idw = 331;
      }
      if (idnu == 23) {
        idw = 211;
      }
      if (idnu == 24) {
        idw = 321;
      }
      if (idnu == 25) {
        idw = 2212;
      }
      Int_t idhnu = 10000 + idw;
      if (idadd == -1) {
        if (idnu > 17) {
          continue;
        }
        idhnu += 10000;
        idw = -idnu;
      }
      TString name = PDG->GetParticle(idw)->GetName();
      TString title = name;
      title += " momentum (GeV)";
      TString key = fVetoName;
      key += idhnu;
      [[maybe_unused]] TH1D* Hidhnu = new TH1D(key, title, 400, 0., 400.);
      title = name;
      title += "  log10-p vs log10-pt";
      key = fVetoName;
      key += idhnu + 1000;
      [[maybe_unused]] TH2D* Hidhnu100 =
          new TH2D(key, title, 100, -0.3, 1.7, 100, -2., 0.5);
      title = name;
      title += "  log10-p vs log10-pt";
      key = fVetoName;
      key += idhnu + 2000;
      [[maybe_unused]] TH2D* Hidhnu200 =
          new TH2D(key, title, 25, -0.3, 1.7, 100, -2., 0.5);
    }
  }
  if (withNtuple) {
    fNtuple = new TNtuple("4DP", "4DP", "id:px:py:pz:x:y:z");
  }
}

◆ PreTrack()

void exitHadronAbsorber::PreTrack ( )

override

Definition at line 155 of file exitHadronAbsorber.cxx.

                                  {
  gMC->TrackMomentum(fMom);
  if ((fMom.E() - fMom.M()) < EMax) {
    gMC->StopTrack();
    return;
  }
  TParticle* p = gMC->GetStack()->GetCurrentTrack();
  Int_t pdgCode = p->GetPdgCode();
  // record statistics for neutrinos, electrons and photons
  // add pi0 111 eta 221 eta' 331  omega 223
  Int_t idabs = TMath::Abs(pdgCode);
  if (idabs < 18 || idabs == 22 || idabs == 111 || idabs == 221 ||
      idabs == 223 || idabs == 331 || idabs == 211 || idabs == 321 ||
      idabs == 2212) {
    Double_t wspill = p->GetWeight();
    Int_t idhnu = idabs + 10000;
    if (pdgCode < 0) {
      idhnu += 10000;
    }
    Double_t l10ptot =
        TMath::Min(TMath::Max(TMath::Log10(fMom.P()), -0.3), 1.69999);
    Double_t l10pt =
        TMath::Min(TMath::Max(TMath::Log10(fMom.Pt()), -2.), 0.4999);
    TString key = fVetoName;
    key += idhnu;
    TH1D* h1 = dynamic_cast<TH1D*>(fout->Get(key));
    if (h1) {
      h1->Fill(fMom.P(), wspill);
    }
    key = fVetoName;
    key += idhnu + 1000;
    TH2D* h2 = dynamic_cast<TH2D*>(fout->Get(key));
    if (h2) {
      h2->Fill(l10ptot, l10pt, wspill);
    }
    key = fVetoName;
    key += idhnu + 2000;
    h2 = dynamic_cast<TH2D*>(fout->Get(key));
    if (h2) {
      h2->Fill(l10ptot, l10pt, wspill);
    }
    if (withNtuple) {
      fNtuple->Fill(pdgCode, fMom.Px(), fMom.Py(), fMom.Pz(), fPos.X(),
                    fPos.Y(), fPos.Z());
    }
    if (fSkipNeutrinos && (idabs == 12 || idabs == 14 || idabs == 16)) {
      gMC->StopTrack();
    }
  }
}

◆ ProcessHits()

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

override

This method is called from the MC stepping

Definition at line 57 of file exitHadronAbsorber.cxx.

                                                      {
  if (gMC->IsTrackEntering()) {
    fTrackID = gMC->GetStack()->GetCurrentTrackNumber();
    fEventID = gMC->CurrentEvent();
    TParticle* p = gMC->GetStack()->GetCurrentTrack();
    fUniqueID = p->GetUniqueID();
    Int_t pdgCode = p->GetPdgCode();
    gMC->TrackMomentum(fMom);
    if (!fOnlyMuons || TMath::Abs(pdgCode) == 13) {
      fTime = gMC->TrackTime() * 1.0e09;
      fLength = gMC->TrackLength();
      gMC->TrackPosition(fPos);
      if ((fMom.E() - fMom.M()) > EMax) {
        AddHit(fEventID, fTrackID, 111, TVector3(fPos.X(), fPos.Y(), fPos.Z()),
               TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength, 0,
               pdgCode, TVector3(p->Vx(), p->Vy(), p->Vz()),
               TVector3(p->Px(), p->Py(), p->Pz()));
        ShipStack* stack = dynamic_cast<ShipStack*>(gMC->GetStack());
        stack->AddPoint(kVETO);
      }
    }
  }
  if ((!fCylindricalPlane) && fzPos > 1E8) {
    gMC->StopTrack();
  }
  return kTRUE;
}

◆ Register()

void exitHadronAbsorber::Register ( )

override

Definition at line 333 of file exitHadronAbsorber.cxx.

                                  {
  fDetPoints = new std::vector<vetoPoint>();
  TString name = fVetoName + "Point";
  FairRootManager::Instance()->RegisterAny(name.Data(), fDetPoints, kTRUE);
}

◆ SetCylindricalPlane()

void exitHadronAbsorber::SetCylindricalPlane ( )

inline

Definition at line 37 of file exitHadronAbsorber.h.

37{ fCylindricalPlane = kTRUE; }

◆ SetEnergyCut()

void exitHadronAbsorber::SetEnergyCut ( Float_t emax )

inline

Definition at line 31 of file exitHadronAbsorber.h.

31{ EMax = emax; }

◆ SetOnlyMuons()

void exitHadronAbsorber::SetOnlyMuons ( )

inline

Definition at line 32 of file exitHadronAbsorber.h.

32{ fOnlyMuons = kTRUE; }

◆ SetOpt4DP()

void exitHadronAbsorber::SetOpt4DP ( )

inline

Definition at line 33 of file exitHadronAbsorber.h.

33{ withNtuple = kTRUE; }

◆ SetUseCaveCoordinates()

void exitHadronAbsorber::SetUseCaveCoordinates ( )

inline

Definition at line 38 of file exitHadronAbsorber.h.

38{ fUseCaveCoordinates = kTRUE; }

◆ SetVetoPointName()

void exitHadronAbsorber::SetVetoPointName ( TString name )

inline

Definition at line 36 of file exitHadronAbsorber.h.

36{ fVetoName = name; }

◆ SetZposition()

void exitHadronAbsorber::SetZposition ( Float_t x )

inline

Definition at line 35 of file exitHadronAbsorber.h.

35{ fzPos = x; }

run_fixedTarget.x

ROOT x

Definition: run_fixedTarget.py:488

◆ SkipNeutrinos()

void exitHadronAbsorber::SkipNeutrinos ( )

inline

Definition at line 34 of file exitHadronAbsorber.h.

34{ fSkipNeutrinos = kTRUE; }

Member Data Documentation

◆ EMax

Float_t exitHadronAbsorber::EMax

private

Definition at line 48 of file exitHadronAbsorber.h.

◆ fCylindricalPlane

Bool_t exitHadronAbsorber::fCylindricalPlane

private

max energy to transport

Definition at line 49 of file exitHadronAbsorber.h.

◆ fElectrons

TClonesArray* exitHadronAbsorber::fElectrons

private

Definition at line 53 of file exitHadronAbsorber.h.

◆ fNtuple

TNtuple* exitHadronAbsorber::fNtuple

private

special option for Dark Photon physics studies

Definition at line 47 of file exitHadronAbsorber.h.

◆ fOnlyMuons

Bool_t exitHadronAbsorber::fOnlyMuons

private

Definition at line 42 of file exitHadronAbsorber.h.

◆ fout

TFile* exitHadronAbsorber::fout

private

set position from cave

Definition at line 52 of file exitHadronAbsorber.h.

◆ fSkipNeutrinos

Bool_t exitHadronAbsorber::fSkipNeutrinos

private

flag if only muons should be stored

Definition at line 43 of file exitHadronAbsorber.h.

◆ fUniqueID

Int_t exitHadronAbsorber::fUniqueID

private

Definition at line 41 of file exitHadronAbsorber.h.

◆ fUseCaveCoordinates

Bool_t exitHadronAbsorber::fUseCaveCoordinates

private

cylindrical sensPlane flag

Definition at line 50 of file exitHadronAbsorber.h.

◆ fVetoName

TString exitHadronAbsorber::fVetoName

private

flag if neutrinos should be ignored

Definition at line 44 of file exitHadronAbsorber.h.

◆ fzPos

Double_t exitHadronAbsorber::fzPos

private

Definition at line 45 of file exitHadronAbsorber.h.

◆ index

Int_t exitHadronAbsorber::index

private

Definition at line 54 of file exitHadronAbsorber.h.

◆ withNtuple

Bool_t exitHadronAbsorber::withNtuple

private

zPos, optional

Definition at line 46 of file exitHadronAbsorber.h.

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

muonShieldOptimization/exitHadronAbsorber.h
muonShieldOptimization/exitHadronAbsorber.cxx

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ exitHadronAbsorber() [1/2]

◆ exitHadronAbsorber() [2/2]

Member Function Documentation

◆ ConstructGeometry()

◆ FinishRun()

◆ Initialize()

◆ PreTrack()

◆ ProcessHits()

◆ Register()

◆ SetCylindricalPlane()

◆ SetEnergyCut()

◆ SetOnlyMuons()

◆ SetOpt4DP()

◆ SetUseCaveCoordinates()

◆ SetVetoPointName()

◆ SetZposition()

◆ SkipNeutrinos()

Member Data Documentation

◆ EMax

◆ fCylindricalPlane

◆ fElectrons

◆ fNtuple

◆ fOnlyMuons

◆ fout

◆ fSkipNeutrinos

◆ fUniqueID

◆ fUseCaveCoordinates

◆ fVetoName

◆ fzPos

◆ index

◆ withNtuple