PPSDiamondOrganization.cc

Go to the documentation of this file.

// ////////////////////////////////////
// Author
// Seyed Mohsen Etesami  setesami@cern.ch
// //////////////////////////////////////////

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "SimG4CMS/PPS/interface/PPSDiamondOrganization.h"
#include "DataFormats/CTPPSDetId/interface/CTPPSDiamondDetId.h"
#include "G4VPhysicalVolume.hh"
#include "G4VTouchable.hh"
#include "G4Step.hh"

#include <iostream>

//******************************************************************** Constructor and destructor

PPSDiamondOrganization ::PPSDiamondOrganization()
    : theArm_(-1), theStation_(-1), theRoman_pot_(-1), thePlane_(-1), theDetector_(-1) {}

uint32_t PPSDiamondOrganization::unitID(const G4Step* aStep) {
  G4VPhysicalVolume* physVol;
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();

  for (int ii = 0; ii < touch->GetHistoryDepth(); ii++) {
    physVol = touch->GetVolume(ii);

    if (physVol->GetName() == "CTPPS_Diamond_Segment" || physVol->GetName() == "CTPPS_UFSD_Segment") {
      theDetector_ = physVol->GetCopyNo() % 100;
      thePlane_ = physVol->GetCopyNo() / 100;
      LogDebug("PPSSimDiamond") << "\n---------------------CTPPS_Diamond_Segment-------------------------------------"
                                   "------------------------------";
      LogDebug("PPSSimDiamond") << "\t\t\t\t\tDetector name " << physVol->GetName()
                                << " copynumber= " << physVol->GetCopyNo();
      LogDebug("PPSSimDiamond") << "\t\t\t\t\tdetector= " << theDetector_ << " plane= " << thePlane_ << " ii = " << ii;
    }

    else if (physVol->GetName() == "Primary_Vacuum") {
      int cpy_no = physVol->GetCopyNo();
      theArm_ = (cpy_no / 100) % 10;
      theStation_ = (cpy_no / 10) % 10;
      theRoman_pot_ = cpy_no % 10;
    }
    LogDebug("PPSSimDiamond") << "Diamond"
                              << "physVol =" << physVol->GetName() << ", level=" << ii
                              << ", physVol->GetCopyNo()=" << physVol->GetCopyNo() << std::endl;
  }
  return CTPPSDiamondDetId(theArm_, theStation_, theRoman_pot_, thePlane_, theDetector_).rawId();
}