CSCSensitiveDetectorCosmics Class Reference

#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
void	Initialize (G4HCofThisEvent *HCE) override final
	member functions
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
std::string	m_hitCollectionName
	member data
CSCSimHitCollection *	m_myCSCHitColl {nullptr}
AtlasG4EventUserInfo *	m_g4UserEventInfo {nullptr}
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 68 of file CSCSensitiveDetectorCosmics.h.

Constructor & Destructor Documentation

CSCSensitiveDetectorCosmics()

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

construction/destruction

Definition at line 22 of file CSCSensitiveDetectorCosmics.cxx.

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

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 34 of file CSCSensitiveDetectorCosmics.cxx.

{
  m_myCSCHitColl = nullptr;
  if (auto* eventInfo = AtlasG4EventUserInfo::GetEventUserInfo()) {
    m_myCSCHitColl = eventInfo->GetHitCollectionMap()->Find<CSCSimHitCollection>(m_hitCollectionName);
    m_g4UserEventInfo = eventInfo;
  }
  // 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

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 47 of file CSCSensitiveDetectorCosmics.cxx.

                                                                                   {
 
  if (!m_myCSCHitColl) {
    G4Exception("CSCSensitiveDetectorCosmics::ProcessHits", "CSCCosmicHitCollectionMissing", FatalException,
                "Hit collection not initialized; did SetupEvent run?");
    return false;
  }
 
  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(m_g4UserEventInfo ? m_g4UserEventInfo->GetEventStore() : nullptr), 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;
}

Member Data Documentation

m_currVertex

Amg::Vector3D CSCSensitiveDetectorCosmics::m_currVertex

private

Definition at line 83 of file CSCSensitiveDetectorCosmics.h.

m_g4UserEventInfo

AtlasG4EventUserInfo* CSCSensitiveDetectorCosmics::m_g4UserEventInfo {nullptr}

private

Definition at line 90 of file CSCSensitiveDetectorCosmics.h.

{nullptr};

m_globalTime

double CSCSensitiveDetectorCosmics::m_globalTime

private

Definition at line 85 of file CSCSensitiveDetectorCosmics.h.

m_globH

Amg::Vector3D CSCSensitiveDetectorCosmics::m_globH

private

Definition at line 84 of file CSCSensitiveDetectorCosmics.h.

m_hitCollectionName

std::string CSCSensitiveDetectorCosmics::m_hitCollectionName

private

member data

Definition at line 88 of file CSCSensitiveDetectorCosmics.h.

m_mom

Amg::Vector3D CSCSensitiveDetectorCosmics::m_mom

private

Definition at line 80 of file CSCSensitiveDetectorCosmics.h.

m_momMag

double CSCSensitiveDetectorCosmics::m_momMag

private

Definition at line 81 of file CSCSensitiveDetectorCosmics.h.

m_muonHelper

const CscHitIdHelper* CSCSensitiveDetectorCosmics::m_muonHelper

private

Definition at line 91 of file CSCSensitiveDetectorCosmics.h.

m_myCSCHitColl

CSCSimHitCollection* CSCSensitiveDetectorCosmics::m_myCSCHitColl {nullptr}

private

Definition at line 89 of file CSCSensitiveDetectorCosmics.h.

{nullptr};

m_vertex

Amg::Vector3D CSCSensitiveDetectorCosmics::m_vertex

private

Definition at line 82 of file CSCSensitiveDetectorCosmics.h.

---

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

• CSCSensitiveDetectorCosmics.h
• CSCSensitiveDetectorCosmics.cxx