#include <exitHadronAbsorber.h>

Inheritance diagram for exitHadronAbsorber:

Collaboration diagram for exitHadronAbsorber:

[legend]

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

	exitHadronAbsorber ()

virtual	~exitHadronAbsorber ()

virtual void	Initialize ()

virtual Bool_t	ProcessHits (FairVolume *v=0)

virtual void	Register ()

virtual TClonesArray *	GetCollection (Int_t iColl) const

virtual void	Reset ()

void	ConstructGeometry ()

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 ()

vetoPoint *	AddHit (Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgcode, TVector3 Lpos, TVector3 Lmom)

void	SetEnergyCut (Float_t emax)

void	SetOnlyMuons ()

void	SetOpt4DP ()

void	SkipNeutrinos ()

void	SetZposition (Float_t x)

Private Attributes
Int_t	fTrackID

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	fzPos
	length

Bool_t	withNtuple
	zPos, optional

TNtuple *	fNtuple
	special option for Dark Photon physics studies

Float_t	EMax

Bool_t	fOnlyMuons
	max energy to transport

Bool_t	fSkipNeutrinos
	flag if only muons should be stored

TFile *	fout
	flag if neutrinos should be ignored

TClonesArray *	fElectrons

Int_t	index

TClonesArray *	fexitHadronAbsorberPointCollection

Detailed Description

Definition at line 16 of file exitHadronAbsorber.h.

Constructor & Destructor Documentation

◆ exitHadronAbsorber() [1/2]

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

Name : Detector Name Active: kTRUE for active detectors (ProcessHits() will be called) kFALSE for inactive detectors

◆ exitHadronAbsorber() [2/2]

exitHadronAbsorber::exitHadronAbsorber ( )

default constructor

Definition at line 40 of file exitHadronAbsorber.cxx.

  : FairDetector("exitHadronAbsorber", kTRUE, kVETO),
    fTrackID(-1),
    fVolumeID(-1),
    fPos(),
    fMom(),
    fTime(-1.),
    fLength(-1.),
    fOnlyMuons(kFALSE),
    fSkipNeutrinos(kFALSE),
    fzPos(3E8),
    withNtuple(kFALSE),
    fexitHadronAbsorberPointCollection(new TClonesArray("vetoPoint"))
{}

◆ ~exitHadronAbsorber()

exitHadronAbsorber::~exitHadronAbsorber ( )

virtual

destructor

Definition at line 55 of file exitHadronAbsorber.cxx.

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

Member Function Documentation

◆ AddHit()

vetoPoint * exitHadronAbsorber::AddHit	(	Int_t	trackID,
		Int_t	detID,
		TVector3	pos,
		TVector3	mom,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss,
		Int_t	pdgcode,
		TVector3	Lpos,
		TVector3	Lmom
	)

Definition at line 239 of file exitHadronAbsorber.cxx.

{
  TClonesArray& clref = *fexitHadronAbsorberPointCollection;
  Int_t size = clref.GetEntriesFast();
  return new(clref[size]) vetoPoint(trackID, detID, pos, mom,
         time, length, eLoss, pdgCode,Lpos,Lmom);
}

◆ BeginEvent()

virtual void exitHadronAbsorber::BeginEvent ( )

inlinevirtual

Definition at line 66 of file exitHadronAbsorber.h.

66{;}

◆ BeginPrimary()

virtual void exitHadronAbsorber::BeginPrimary ( )

inlinevirtual

Definition at line 63 of file exitHadronAbsorber.h.

63{;}

◆ ConstructGeometry()

void exitHadronAbsorber::ConstructGeometry ( )

Create the detector geometry

Definition at line 205 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==NULL)
     geoBuild->createMedium(ShipMedium);
   vac =gGeoManager->GetMedium("vacuums");
   TGeoVolume *top=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 (not nav->cd("/MuonShieldArea_1/CoatWall_1")) {
      nav->cd("/MuonShieldArea_1/MagnAbsorb2_MagRetL_1");
      distance = -1.;}
    TGeoBBox* tmp =  (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
   TGeoVolume *sensPlane = gGeoManager->MakeBox("sensPlane",vac,10.*m-1.*mm,10.*m-1.*mm,1.*mm);
   nav->cd("/MuonShieldArea_1/");
   nav->GetCurrentNode()->GetVolume()->AddNode(sensPlane, 1, new TGeoTranslation(0, 0, zLoc));
   AddSensitiveVolume(sensPlane);
}

◆ CopyClones()

virtual void exitHadronAbsorber::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 57 of file exitHadronAbsorber.h.

58 {;}

◆ EndOfEvent()

void exitHadronAbsorber::EndOfEvent ( )

virtual

Definition at line 132 of file exitHadronAbsorber.cxx.

{
 
  fexitHadronAbsorberPointCollection->Clear();
 
}

◆ FinishPrimary()

virtual void exitHadronAbsorber::FinishPrimary ( )

inlinevirtual

Definition at line 61 of file exitHadronAbsorber.h.

61{;}

◆ FinishRun()

void exitHadronAbsorber::FinishRun ( )

virtual

Definition at line 173 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 = "";key+=idhnu;
    TSeqCollection* fileList=gROOT->GetListOfFiles();
    ((TFile*)fileList->At(0))->cd();
    TH1D* Hidhnu = (TH1D*)fout->Get(key);
    Hidhnu->Write();
    key="";key+=idhnu+1000;
    TH2D* Hidhnu100 = (TH2D*)fout->Get(key);
    Hidhnu100->Write();
    key = "";key+=idhnu+2000;
    TH2D* Hidhnu200 = (TH2D*)fout->Get(key);
    Hidhnu200->Write();
   }
  }
  if(withNtuple){fNtuple->Write();}
}

◆ GetCollection()

TClonesArray * exitHadronAbsorber::GetCollection ( Int_t iColl ) const

virtual

Gets the produced collections

Definition at line 257 of file exitHadronAbsorber.cxx.

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

◆ Initialize()

void exitHadronAbsorber::Initialize ( )

virtual

Initialization of the detector is done here

Definition at line 88 of file exitHadronAbsorber.cxx.

{
  FairDetector::Initialize();
  TSeqCollection* fileList=gROOT->GetListOfFiles();
  fout = ((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 = "";key+=idhnu;
    TH1D* Hidhnu = new TH1D(key,title,400,0.,400.);
    title = name;title+="  log10-p vs log10-pt";
    key = "";key+=idhnu+1000;
    TH2D* Hidhnu100 = new TH2D(key,title,100,-0.3,1.7,100,-2.,0.5);
    title = name;title+="  log10-p vs log10-pt";
    key = "";key+=idhnu+2000;
    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");
  }
}

◆ PostTrack()

virtual void exitHadronAbsorber::PostTrack ( )

inlinevirtual

Definition at line 64 of file exitHadronAbsorber.h.

64{;}

◆ PreTrack()

void exitHadronAbsorber::PreTrack ( )

virtual

Definition at line 139 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; key+=idhnu;
         TH1D* h1 = (TH1D*)fout->Get(key);
         if (h1){h1->Fill(fMom.P(),wspill);}
         key="";key+=idhnu+1000;
         TH2D* h2 = (TH2D*)fout->Get(key);
         if (h2){h2->Fill(l10ptot,l10pt,wspill);}
         key="";key+=idhnu+2000;
         h2 = (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 or idabs==14 or idabs == 16 )){gMC->StopTrack();}
   }
}

◆ ProcessHits()

Bool_t exitHadronAbsorber::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 63 of file exitHadronAbsorber.cxx.

{
  if ( gMC->IsTrackEntering() ) {
    fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
    TParticle* p  = gMC->GetStack()->GetCurrentTrack();
    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(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 = (ShipStack*) gMC->GetStack();
      stack->AddPoint(kVETO);
      }
    }
  }
  gMC->StopTrack();
  return kTRUE;
}

◆ Register()

void exitHadronAbsorber::Register ( )

virtual

Registers the produced collections in FAIRRootManager.

Definition at line 250 of file exitHadronAbsorber.cxx.

{
 
  FairRootManager::Instance()->Register("vetoPoint", "veto",
                                        fexitHadronAbsorberPointCollection, kTRUE);
}

◆ Reset()

void exitHadronAbsorber::Reset ( )

virtual

has to be called after each event to reset the containers

Definition at line 262 of file exitHadronAbsorber.cxx.

{
  fexitHadronAbsorberPointCollection->Clear();
}

◆ SetEnergyCut()

void exitHadronAbsorber::SetEnergyCut ( Float_t emax )

inline

Definition at line 72 of file exitHadronAbsorber.h.

72{EMax=emax;}// min energy to be copied to Geant4

makeSNDGenieEvents.emax

emax

Definition makeSNDGenieEvents.py:140

◆ SetOnlyMuons()

void exitHadronAbsorber::SetOnlyMuons ( )

inline

Definition at line 73 of file exitHadronAbsorber.h.

73{fOnlyMuons=kTRUE;}

◆ SetOpt4DP()

void exitHadronAbsorber::SetOpt4DP ( )

inline

Definition at line 74 of file exitHadronAbsorber.h.

74{withNtuple=kTRUE;}

◆ SetSpecialPhysicsCuts()

virtual void exitHadronAbsorber::SetSpecialPhysicsCuts ( )

inlinevirtual

Definition at line 59 of file exitHadronAbsorber.h.

59{;}

◆ SetZposition()

void exitHadronAbsorber::SetZposition ( Float_t x )

inline

Definition at line 76 of file exitHadronAbsorber.h.

76{fzPos=x;}

geometry_config.x

x

Definition geometry_config.py:180

◆ SkipNeutrinos()

void exitHadronAbsorber::SkipNeutrinos ( )

inline

Definition at line 75 of file exitHadronAbsorber.h.

75{fSkipNeutrinos=kTRUE;}

Member Data Documentation

◆ EMax

Float_t exitHadronAbsorber::EMax

private

Definition at line 92 of file exitHadronAbsorber.h.

◆ fElectrons

TClonesArray* exitHadronAbsorber::fElectrons

private

Definition at line 96 of file exitHadronAbsorber.h.

◆ fexitHadronAbsorberPointCollection

TClonesArray* exitHadronAbsorber::fexitHadronAbsorberPointCollection

private

container for data points

Definition at line 99 of file exitHadronAbsorber.h.

◆ fLength

Double_t exitHadronAbsorber::fLength

private

time

Definition at line 88 of file exitHadronAbsorber.h.

◆ fMom

TLorentzVector exitHadronAbsorber::fMom

private

position at entrance

Definition at line 86 of file exitHadronAbsorber.h.

◆ fNtuple

TNtuple* exitHadronAbsorber::fNtuple

private

special option for Dark Photon physics studies

Definition at line 91 of file exitHadronAbsorber.h.

◆ fOnlyMuons

Bool_t exitHadronAbsorber::fOnlyMuons

private

max energy to transport

Definition at line 93 of file exitHadronAbsorber.h.

◆ fout

TFile* exitHadronAbsorber::fout

private

flag if neutrinos should be ignored

Definition at line 95 of file exitHadronAbsorber.h.

◆ fPos

TLorentzVector exitHadronAbsorber::fPos

private

volume id

Definition at line 85 of file exitHadronAbsorber.h.

◆ fSkipNeutrinos

Bool_t exitHadronAbsorber::fSkipNeutrinos

private

flag if only muons should be stored

Definition at line 94 of file exitHadronAbsorber.h.

◆ fTime

Double_t exitHadronAbsorber::fTime

private

momentum at entrance

Definition at line 87 of file exitHadronAbsorber.h.

◆ fTrackID

Int_t exitHadronAbsorber::fTrackID

private

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

Definition at line 83 of file exitHadronAbsorber.h.

◆ fVolumeID

Int_t exitHadronAbsorber::fVolumeID

private

track index

Definition at line 84 of file exitHadronAbsorber.h.

◆ fzPos

Double_t exitHadronAbsorber::fzPos

private

length

Definition at line 89 of file exitHadronAbsorber.h.

◆ index

Int_t exitHadronAbsorber::index

private

Definition at line 97 of file exitHadronAbsorber.h.

◆ withNtuple

Bool_t exitHadronAbsorber::withNtuple

private

zPos, optional

Definition at line 90 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

Detailed Description

Constructor & Destructor Documentation

◆ exitHadronAbsorber() [1/2]

◆ exitHadronAbsorber() [2/2]

◆ ~exitHadronAbsorber()

Member Function Documentation

◆ AddHit()

◆ BeginEvent()

◆ BeginPrimary()

◆ ConstructGeometry()

◆ CopyClones()

◆ EndOfEvent()

◆ FinishPrimary()

◆ FinishRun()

◆ GetCollection()

◆ Initialize()

◆ PostTrack()

◆ PreTrack()

◆ ProcessHits()

◆ Register()

◆ Reset()

◆ SetEnergyCut()

◆ SetOnlyMuons()

◆ SetOpt4DP()

◆ SetSpecialPhysicsCuts()

◆ SetZposition()

◆ SkipNeutrinos()

Member Data Documentation

◆ EMax

◆ fElectrons

◆ fexitHadronAbsorberPointCollection

◆ fLength

◆ fMom

◆ fNtuple

◆ fOnlyMuons

◆ fout

◆ fPos

◆ fSkipNeutrinos

◆ fTime

◆ fTrackID

◆ fVolumeID

◆ fzPos

◆ index

◆ withNtuple