CSCSensitiveDetector Class Reference

#include <CSCSensitiveDetector.h>

Inheritance diagram for CSCSensitiveDetector:

Collaboration diagram for CSCSensitiveDetector:

Public Member Functions
	CSCSensitiveDetector (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 (CSCSensitiveDetectortest, Initialize)
	FRIEND_TEST (CSCSensitiveDetectortest, ProcessHits)

Private Attributes
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 CSCSensitiveDetector::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 CSCSensitiveDetector 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.

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 63 of file CSCSensitiveDetector.h.

Constructor & Destructor Documentation

CSCSensitiveDetector()

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

construction/destruction

Definition at line 20 of file CSCSensitiveDetector.cxx.

  : G4VSensitiveDetector( name )
  , m_hitCollectionName( hitCollectionName )
{
  m_muonHelper = CscHitIdHelper::GetHelper();
}

Member Function Documentation

FRIEND_TEST() [1/2]

CSCSensitiveDetector::FRIEND_TEST ( CSCSensitiveDetectortest , Initialize  )

private

FRIEND_TEST() [2/2]

CSCSensitiveDetector::FRIEND_TEST ( CSCSensitiveDetectortest , ProcessHits  )

private

Initialize()

void CSCSensitiveDetector::Initialize ( G4HCofThisEvent * HCE )

finaloverride

member functions

Definition at line 29 of file CSCSensitiveDetector.cxx.

{
  m_myCSCHitColl = nullptr;
  if (auto* eventInfo = AtlasG4EventUserInfo::GetEventUserInfo()) {
    m_myCSCHitColl = eventInfo->GetHitCollectionMap()->Find<CSCSimHitCollection>(m_hitCollectionName);
    m_g4UserEventInfo = eventInfo;
  }
}

ProcessHits()

G4bool CSCSensitiveDetector::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

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

the track 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 38 of file CSCSensitiveDetector.cxx.

                                                                            {
 
  if (!m_myCSCHitColl) {
    G4Exception("CSCSensitiveDetector::ProcessHits", "CSCHitCollectionMissing", 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=1;
  int stationPhi=1;
  int multiLayer=1;
  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;
  // int trackid         = -1;

 
  bool isAssembly = false;
  for (int i = touchHist->GetHistoryDepth(); i>=0; --i) {
    std::string::size_type npos;
    std::string volName = touchHist->GetVolume(i)->GetName();
 
    //G4cout << "Csc SD: swimming through the tree: level "<<i-touchHist->GetHistoryDepth() << G4endl;
    //G4cout << "Csc SD: name "<<volName << G4endl;

    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
      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
    } 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;
      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=999;
 
      if (lundcode < 0)
        G4cout << "WP CSCSensitiveDetector.cxx have exceptional particlue out of selection::  "
               << particle<< G4endl;

      // trackid=aStep->GetTrack()->GetTrackID();

      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) );
 
    }
  }
 
  TrackHelper trHelp(aStep->GetTrack());

  HitID CSCid = m_muonHelper->BuildCscHitId(stationName, stationPhi,
                                          stationEta, multiLayer, wireLayer);
  m_myCSCHitColl->Emplace(CSCid, globalTime, energyDeposit,
                          HitStart, HitEnd, lundcode,
                          trHelp.GenerateParticleLink(m_g4UserEventInfo ? m_g4UserEventInfo->GetEventStore() : nullptr),kinEnergy);
 
  return true;
}

Member Data Documentation

m_g4UserEventInfo

AtlasG4EventUserInfo* CSCSensitiveDetector::m_g4UserEventInfo {nullptr}

private

Definition at line 79 of file CSCSensitiveDetector.h.

{nullptr};

m_hitCollectionName

std::string CSCSensitiveDetector::m_hitCollectionName

private

member data

Definition at line 77 of file CSCSensitiveDetector.h.

m_muonHelper

const CscHitIdHelper* CSCSensitiveDetector::m_muonHelper

private

Definition at line 80 of file CSCSensitiveDetector.h.

m_myCSCHitColl

CSCSimHitCollection* CSCSensitiveDetector::m_myCSCHitColl {nullptr}

private

Definition at line 78 of file CSCSensitiveDetector.h.

{nullptr};

---

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

• CSCSensitiveDetector.h
• CSCSensitiveDetector.cxx