#include <CSCSensitiveDetectorCosmics.h>

Inheritance diagram for CSCSensitiveDetectorCosmics:

Collaboration diagram for CSCSensitiveDetectorCosmics:

Public Member Functions
	CSCSensitiveDetectorCosmics (const std::string &name, const std::string &hitCollectionName)
	construction/destruction More...

	~CSCSensitiveDetectorCosmics ()

void	Initialize (G4HCofThisEvent *HCE) override final
	member functions More...

G4bool	ProcessHits (G4Step aStep, G4TouchableHistory ROhist) override final

Private Member Functions
	FRIEND_TEST (CSCSensitiveDetectorCosmicstest, Initialize)

	FRIEND_TEST (CSCSensitiveDetectorCosmicstest, ProcessHits)

Private Attributes
Amg::Vector3D	m_mom

double	m_momMag

Amg::Vector3D	m_vertex

Amg::Vector3D	m_currVertex

Amg::Vector3D	m_globH

double	m_globalTime

SG::WriteHandle< CSCSimHitCollection >	m_myCSCHitColl
	member data More...

const CscHitIdHelper *	m_muonHelper

Detailed Description

Author: ketev.nosp@m.i@bn.nosp@m.l.gov

Class methods and properties

The method CSCSensitiveDetectorCosmics::ProcessHits is executed by the G4 kernel each time a charged particle (or a geantino) crosses one of the CSC Gas Gaps.

Navigating with the toucahbleHistory method GetHistoryDepth() through the hierarchy of the volumes crossed by the particles, the CSCSensitiveDetectorCosmics determinates the correct set of geometry parameters to be folded in the Simulation Identifier associated to each hit. The CSC SimIDs are 32-bit unsigned integers, built using the MuonSimEvent/CscHitIdHelper class which inherits from the MuonHitIdHelper base class.

The modified version (Daniela Rebuzzi, 16/12/2005) calculates the tof for the cosmics, using as t0 the time to the closest approach of the track to (0,0,0). This implementation works also for slow particles, since the discrimination physics event/cosmics is done relying on the G4Track vertex (not on the tof).

We describe in the following, how each field of the identifier is retrieved.

1) stationName, stationEta, stationPhi: when a volume is found in the hierarchy, whose name contains the substring "station", the stationName is extracted from the volume's name; stationPhi and stationEta are calculated starting from the volume copy number, assigned by MuonGeoModel.

2) multilayer: when a volume is found in the hierarchy, whose name contains the substring "component", the multilayer is set to 1 or to 2, according to the component number (multilayer=1 if the component copy number is 1, multilayer=3 if the component copy number is 2).

3) wireLayer: when a volume is found in the hierarchy, whose name contains the substring "CscArCO2", the wire Layer is calculated, increasing with |z| Atlas global coordinate.

notes:

1) for each hit, the time of flight (the G4 globalTime), is recorded and associated to the hit.

2) the CSCHit object contains: the SimID, the globalTime, the energy deposited, the start and end point of the step in local coordinates (chamber reference system), the particle pdg and the track identifier.

Definition at line 66 of file CSCSensitiveDetectorCosmics.h.

Constructor & Destructor Documentation

◆ CSCSensitiveDetectorCosmics()

CSCSensitiveDetectorCosmics::CSCSensitiveDetectorCosmics	(	const std::string &	name,
		const std::string &	hitCollectionName
	)

construction/destruction

Definition at line 20 of file CSCSensitiveDetectorCosmics.cxx.

   : G4VSensitiveDetector( name )
   // START OF COSMICS-SPECIFIC CODE
   , m_momMag(0.)
   , m_globalTime(0.)
   // END OF COSMICS-SPECIFIC CODE
   , m_myCSCHitColl( hitCollectionName )
 {
   m_muonHelper = CscHitIdHelper::GetHelper();
 }

◆ ~CSCSensitiveDetectorCosmics()

CSCSensitiveDetectorCosmics::~CSCSensitiveDetectorCosmics ( )

inline

Definition at line 72 of file CSCSensitiveDetectorCosmics.h.

72 {}

Member Function Documentation

◆ FRIEND_TEST() [1/2]

CSCSensitiveDetectorCosmics::FRIEND_TEST	(	CSCSensitiveDetectorCosmicstest	,
		Initialize
	)

private

◆ FRIEND_TEST() [2/2]

CSCSensitiveDetectorCosmics::FRIEND_TEST	(	CSCSensitiveDetectorCosmicstest	,
		ProcessHits
	)

private

◆ Initialize()

void CSCSensitiveDetectorCosmics::Initialize ( G4HCofThisEvent * HCE )

finaloverride

member functions

Definition at line 32 of file CSCSensitiveDetectorCosmics.cxx.

 {
   if (!m_myCSCHitColl.isValid()) m_myCSCHitColl = std::make_unique<CSCSimHitCollection>();
   // START OF COSMICS-SPECIFIC CODE
   m_mom = Amg::Vector3D(0.,0.,0.);
   m_globH = Amg::Vector3D(0.,0.,0.);
   // END OF COSMICS-SPECIFIC CODE
 }

◆ ProcessHits()

G4bool CSCSensitiveDetectorCosmics::ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *	ROhist
	)

finaloverride

CSCs sensitive to charged particle only; plus energy threshold

attributes of the CSC identifier construction

hit information to be recorded

the track id

scan geometry tree to identify the gas layer

station: name, eta and phi (-> chamber!)

station name, station eta and station phi

station name

station eta and station phi

chamber layer

the wire layer number

get the particle ID

record global time

energy deposit

get the transformation

local position of HitStart and HitEnd

get the end of the step as return by G4 - could be at the beginning of the next volume

construct the hit identifier

insert hit in collection

Definition at line 41 of file CSCSensitiveDetectorCosmics.cxx.

                                                                                    {
  
   G4Track* currentTrack = aStep->GetTrack();
  
   auto trackDef = currentTrack->GetDefinition();
   if (trackDef->GetPDGCharge() == 0.0) {
     if (trackDef != G4Geantino::GeantinoDefinition()) return true;
     else if (trackDef == G4ChargedGeantino::ChargedGeantinoDefinition()) return true;
   }
  
   const G4TouchableHistory* touchHist = static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
   G4ThreeVector startPos=aStep->GetPreStepPoint()->GetPosition();
   G4ThreeVector endPos=aStep->GetPostStepPoint()->GetPosition();
   double kinEnergy = aStep->GetPreStepPoint()->GetKineticEnergy();
  
   std::string stationName="";
   int stationEta=0;
   int stationPhi=0;
   int multiLayer=0;
   int wireLayer=0;
  
   Amg::Vector3D HitStart = Amg::Vector3D(-1000,-1000,-1000);
   Amg::Vector3D HitEnd   = Amg::Vector3D(-1000,-1000,-1000);
   double globalTime   = -1;
   double energyDeposit= -1;
   G4int lundcode      = -1;
   //const int trackid(aStep->GetTrack()->GetTrackID());
  
   // START OF COSMICS-SPECIFIC CODE
   // global coordinates
   Amg::Vector3D globVrtx = Amg::Hep3VectorToEigen( aStep->GetPreStepPoint()->GetPosition() );
  
   // distance of the hit from (0,0,0) vertex - calculated from the PreStepPoint (approximation)
   const double dist = globVrtx.mag();
   const double lightspeed = 299.792458; /* in vacuo speed of light [mm/ns] */ //FIXME this should be taken from standard headers
   const double tOrigin = dist / lightspeed;
  
   m_currVertex = Amg::Hep3VectorToEigen( aStep->GetTrack()->GetVertexPosition() );
  
   // for cosmics: only primary muon tracks - track momentum when first entering the spectrometer (one muon per event)
   if ((m_currVertex != m_vertex) && (aStep->GetTrack()->GetTrackID() == 1)) {
     // after calculationg the momentum magnidude, normalize it
     m_mom = Amg::Hep3VectorToEigen( currentTrack->GetMomentum() );
     m_momMag = m_mom.mag();
     m_mom.normalize();
     // the direction of the primary mu is used to calculate the t0, the position to the t0, m_globH, is ONE for a track
     Amg::Vector3D globVrtxFix = globVrtx;
     const double AlphaGlobal = -1*(globVrtxFix[0]*m_mom[0] + globVrtxFix[1]*m_mom[1] + globVrtxFix[2]*m_mom[2])/(m_mom[0]*m_mom[0] + m_mom[1]*m_mom[1] + m_mom[2]*m_mom[2]);
     m_globH = globVrtxFix + AlphaGlobal*m_mom;
     //    G4cout << "COSMICS MAIN TRACK IN THES CSC!" << G4endl;
   }
   const double globalDist = sqrt((m_globH[0] - globVrtx[0])*(m_globH[0] - globVrtx[0]) +
                                  (m_globH[1] - globVrtx[1])*(m_globH[1] - globVrtx[1]) +
                                  (m_globH[2] - globVrtx[2])*(m_globH[2] - globVrtx[2]));
   double tof = globalDist / lightspeed;
   // END OF COSMICS SPECIFIC CODE
  
   bool isAssembly = false;
   for (int i = touchHist->GetHistoryDepth(); i>=0; --i) {
     std::string::size_type npos;
     std::string volName = touchHist->GetVolume(i)->GetName();
  
     if ((npos = volName.find("av_")) != std::string::npos &&
         (npos = volName.find("impr_")) != std::string::npos)  isAssembly = true;
  
     if ((npos = volName.find("station")) != std::string::npos && (!isAssembly)) {
  
       volName.resize(npos-2);
       int volCopyNo = touchHist->GetVolume(i)->GetCopyNo();
       stationName   = volName;
       //stationEta    = volCopyNo/100;
       // bug fix for 14.2.0 to cope with non-mirrored chambers as proposed by
       // S.Spagnolo (ported from CSCSensitiveDetector.cxx for rel 20.3.3)
       stationEta    = (volCopyNo%1000)/100;
       stationPhi    = abs(volCopyNo%100);
  
     } else if ((npos = volName.find("CSC")) != std::string::npos && isAssembly ) {
       // vol name for Assembly components are
       // av_WWW_impr_XXX_Muon::BMSxMDTxx_pv_ZZZ_NAME
       //         where WWW is ass. istance nr.
       //               XXX is comp. imprint nr.
       //               BMSxMDTxx is the name of the comp. log.Vol.
       //                  x station sub-type; xx technology subtype
       //               ZZZ is No of order of this daugther
       //               NAME is the comp. tag (geoIdentifierTag)
       //               for CSCs it is cl[1] or cl[2]
       // copy numbers for Ass.components are =
       //               CopyNoBase(= geoIdentifierTag of the assembly) + geoIdentifierTag of the component
       //               geoIdentifierTag of the component = Job
       //               geoIdentifierTag of the assembly = (sideC*10000 +
       //                             mirsign*1000 + abs(zi)*100 + fi+1)*100000;
       //                             mirsign*1000 + abs(zi)*100 + fi+1)*100000;
  
       std::string::size_type loc1,loc2;
       if ((loc1 = volName.find("Muon::")) != std::string::npos) {
         stationName = volName.substr(loc1+6,3); //type only
       }
  
       int volCopyNo = touchHist->GetVolume(i)->GetCopyNo();
       int copyNrBase = int(volCopyNo/100000);
       int sideC  = int(copyNrBase/10000);
       int zi     = int((copyNrBase%1000)/100);
       // int mirfl  = int((copyNrBase%10000)/1000);
       int fi     = int(copyNrBase%100);
       if (sideC == 1) zi = -zi;
       stationEta    = zi;
       stationPhi    = fi;
  
       // now get the geoIdentifierTag of the rpc components
       int gmID = 0;
       if ((loc1 = volName.find('[')) != std::string::npos) {
         if ((loc2 = volName.find(']', loc1+1)) != std::string::npos) {
           std::istringstream istrvar(volName.substr(loc1+1,loc2-loc1-1));
           istrvar>>gmID;
         }
       }
       multiLayer = gmID;
     } else if ((npos = volName.find("component")) != std::string::npos && (!isAssembly)) {
  
       multiLayer = touchHist->GetVolume(i)->GetCopyNo();
       if(multiLayer==3) multiLayer=2; //multilayer index
       // G4cout << "CSC:::::: multiLayer  "<<multiLayer << G4endl;
     } else if ((npos = volName.find("CscArCO2")) != std::string::npos) {
  
       wireLayer=touchHist->GetVolume(i)->GetCopyNo();
       wireLayer+=1;
       if(wireLayer==4) wireLayer=1;
       else if(wireLayer==3) wireLayer=2;
       else if(wireLayer==2) wireLayer=3;
       else if(wireLayer==1) wireLayer=4;
       // G4cout << "CSC:::::: wireLayer  "<<wireLayer << G4endl;
  
       G4String particle=aStep->GetTrack()->GetDefinition()->GetParticleName();
       // http://www.slac.stanford.edu/BFROOT/www/Computing/Environment/NewUser/htmlbug/node51.html #GEANT code
       if (particle=="gamma") lundcode=1; // FIXME was photon, but changed to gamma to match CSCSensitiveDetctor
       else if (particle=="e+")     lundcode=2;
       else if (particle=="e-")     lundcode=3;
       else if (particle=="mu+")    lundcode=5;
       else if (particle=="mu-")    lundcode=6;
       else if (particle=="pi+")    lundcode=8;
       else if (particle=="pi-")    lundcode=9;
       else if (particle=="kaon+")  lundcode=11;
       else if (particle=="kaon-")  lundcode=12;
       else if (particle=="proton") lundcode=14;
       else if (particle=="anti_proton") lundcode=15;
       else if (particle=="sigma+") lundcode=19;
       else if (particle=="sigma-") lundcode=21;
       else if (particle=="anti_sigma-") lundcode=27;
       else if (particle=="anti_sigma+") lundcode=29;
       else if (particle=="deuteron") lundcode=45;
       else if (particle=="geantino") lundcode=999;
       else lundcode=0;
  
       globalTime = aStep->GetPreStepPoint()->GetGlobalTime();
  
       energyDeposit = aStep->GetTotalEnergyDeposit();
  
       const G4AffineTransform transform = touchHist->GetHistory()->GetTopTransform();
  
       HitStart = Amg::Hep3VectorToEigen( transform.TransformPoint(startPos) );
  
       HitEnd = Amg::Hep3VectorToEigen( transform.TransformPoint(endPos) );
  
     }
   }
  
   // START OF COSMICS SPECIFIC CODE
   m_vertex = Amg::Hep3VectorToEigen( aStep->GetTrack()->GetVertexPosition() );
   // if the track vertex is far from (0,0,0), takes the tof, otherwise take the "usual" g4 globalTime
   (((m_vertex.mag() < 100) || ((fabs(globalTime - tOrigin)) < 0.1) ) ? (m_globalTime  = globalTime)
                                                                    : (m_globalTime = tof));
   // if m_globalTime  != globalTime and m_globalTime != tof in the output, this is due to multiple hits
   // before founding the good one (small approximation)
   // END OF COSMICS SPECIFIC CODE
  
   TrackHelper trHelp(aStep->GetTrack());
  
   HitID CSCid = m_muonHelper->BuildCscHitId(stationName, stationPhi,
                                           stationEta, multiLayer, wireLayer);
  
   m_myCSCHitColl->Emplace(CSCid, m_globalTime, energyDeposit,
                           HitStart, HitEnd, lundcode,
                           trHelp.GenerateParticleLink(), kinEnergy);
  
   // #ifndef CSCG4_DEBUG
   //
   // // printouts for cosmics
   //         G4cout << "------------------CSC------------------" << G4endl;
   //    G4cout << "Track "<<aStep->GetTrack()->GetTrackID() << G4endl;
   //    G4cout << "Track m_vertex "<<m_vertex[0]<<" " <<m_vertex[1]<<" " <<m_vertex[2] << G4endl;
 //      G4cout << "Global position of the hit " << globVrtx[0] <<" " << globVrtx[1] <<" " << globVrtx[2]  << G4endl;
 //         G4cout << "Distance from (0,0,0) and time " << dist << " " <<tOrigin  << G4endl;
 //         G4cout << "Momentum "<<m_momMag << G4endl;
 //      G4cout << "Momentum director cosines " <<m_mom[0]<<" " <<m_mom[1]<<" " <<m_mom[2] << G4endl;
 //      G4cout << "Eta and phi "<<m_mom.eta()<<" " <<m_mom.phi() << G4endl;
 //         G4cout << "Closest approach position and distance from (0,0,0) "
 //                << m_globH[0] <<" "<<m_globH[1]<<" "<<m_globH[2]<<" "
 //                << sqrt(m_globH[0]*m_globH[0] + m_globH[1]*m_globH[1] + m_globH[2]*m_globH[2])  << G4endl;
 //         G4cout << "Distance from t0 and tof " << globalDist <<" " <<tof  << G4endl;
 //         G4cout << "g4 globalTime " << globalTime  << G4endl;
 //         G4cout << "Time " << m_globalTime  << G4endl;
 //
 //      G4cout << "TUB "<<m_muonHelper->GetStationName(CSCid)
 //                  << " "<<m_muonHelper->GetFieldValue("PhiSector")
 //                  << " "<<m_muonHelper->GetFieldValue("ZSector")
 //                  << " "<<m_muonHelper->GetFieldValue("ChamberLayer")
 //                  << " "<<m_muonHelper->GetFieldValue("WireLayer") << G4endl;
 //
 //      G4cout << m_muonHelper->GetStationName(CSCid)<<" "<<newHit->print() << G4endl;
 //
 //  #endif
  
   return true;
 }