Ex03CalorimeterSD.cxx

Go to the documentation of this file.

//------------------------------------------------
// The Virtual Monte Carlo examples
// Copyright (C) 2014 - 2018 Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see geant4_vmc/LICENSE.
// Contact: root-vmc@cern.ch
//-------------------------------------------------
 
/// \file Ex03CalorimeterSD.cxx
/// \brief Implementation of the Ex03CalorimeterSD class
///
/// Geant4 ExampleN03 adapted to Virtual Monte Carlo \n
/// Id: ExN03CalorimeterSD.cc,v 1.6 2002/01/09 17:24:12 ranjard Exp \n
/// GEANT4 tag $Name:  $
///
/// \date 06/03/2002
/// \author I. Hrivnacova; IPN, Orsay
 
#include "Ex03CalorimeterSD.h"
#include "Ex03CalorHit.h"
#include "Ex03DetectorConstruction.h"
 
#include <Riostream.h>
#include <TLorentzVector.h>
#include <TMCRootManager.h>
#include <TTree.h>
#include <TVirtualMC.h>
 
/// \cond CLASSIMP
ClassImp(Ex03CalorimeterSD)
  /// \endcond
 
  using namespace std;
 
//_____________________________________________________________________________
Ex03CalorimeterSD::Ex03CalorimeterSD(
  const char* name, Ex03DetectorConstruction* detector)
  : TNamed(name, ""),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    fAbsorberVolId(0),
    fGapVolId(0),
    fAssemblyHierarchyPrinted(kFALSE),
    fVerboseLevel(1),
    fReflectedHit()
{
  /// Standard constructor.
  /// Create hits collection and an empty hit for each layer
  /// As the copy numbers may start from 0 or 1 (depending on
  /// geometry model, we create one more layer for this case.)
  /// \param name      The calorimeter hits collection name
  /// \param detector  The detector construction
 
  fCalCollection = new TClonesArray("Ex03CalorHit", 500);
  for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
    new ((*fCalCollection)[i]) Ex03CalorHit();
}
 
//_____________________________________________________________________________
Ex03CalorimeterSD::Ex03CalorimeterSD(
  const Ex03CalorimeterSD& origin, Ex03DetectorConstruction* detector)
  : TNamed(origin),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    fAbsorberVolId(origin.fAbsorberVolId),
    fGapVolId(origin.fGapVolId),
    fAssemblyHierarchyPrinted(kFALSE),
    fVerboseLevel(origin.fVerboseLevel),
    fReflectedHit()
{
  /// Copy constructor (for clonig on worker thread in MT mode).
  /// Create hits collection and an empty hit for each layer
  /// As the copy numbers may start from 0 or 1 (depending on
  /// geometry model, we create one more layer for this case.)
  /// \param origin    The source object (on master).
  /// \param detector  The detector construction
 
  fCalCollection = new TClonesArray("Ex03CalorHit", 500);
  for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
    new ((*fCalCollection)[i]) Ex03CalorHit();
}
 
//_____________________________________________________________________________
Ex03CalorimeterSD::Ex03CalorimeterSD()
  : TNamed(),
    fDetector(0),
    fCalCollection(0),
    fAbsorberVolId(0),
    fGapVolId(0),
    fAssemblyHierarchyPrinted(kFALSE),
    fVerboseLevel(1),
    fReflectedHit()
{
  /// Default constructor
}
 
//_____________________________________________________________________________
Ex03CalorimeterSD::~Ex03CalorimeterSD()
{
  /// Destructor
 
  if (fCalCollection) fCalCollection->Delete();
  delete fCalCollection;
}
 
//
// private methods
//
 
//_____________________________________________________________________________
Ex03CalorHit* Ex03CalorimeterSD::GetHit(Int_t i) const
{
  /// \return   The hit for the specified layer.
  /// \param i  The layer number
 
  return (Ex03CalorHit*)fCalCollection->At(i);
}
 
//_____________________________________________________________________________
void Ex03CalorimeterSD::ResetHits()
{
  /// Reset all hits in the hits collection.
 
  for (Int_t i = 0; i < fCalCollection->GetEntriesFast(); i++)
    GetHit(i)->Reset();
}
 
//
// public methods
//
 
//_____________________________________________________________________________
void Ex03CalorimeterSD::Initialize()
{
  /// Register hits collection in the Root manager;
  /// set sensitive volumes.
 
  if (TMCRootManager::Instance()) Register();
 
  // Keep the pointer to TVirtualMC object as a data member
  // to avoid a possible performance penalty due to a frequent retrieval
  // from the thread-local storage
  fMC = gMC;
 
  fAbsorberVolId = fMC->VolId("ABSO");
  fGapVolId = fMC->VolId("GAPX");
}
 
//_____________________________________________________________________________
Bool_t Ex03CalorimeterSD::ProcessHits()
{
  /// Account energy deposit and track lengths for each layer in its hit.
 
  Int_t copyNo;
  Int_t id = fMC->CurrentVolID(copyNo);
 
  if (id != fAbsorberVolId && id != fGapVolId) return false;
 
  if (fDetector->GetUseAssemblies() && !fAssemblyHierarchyPrinted) {
    fAssemblyHierarchyPrinted = kTRUE;
    cout << "Assembly volume hierarchy:" << endl;
    for (Int_t off = 0; off <= 6; ++off) {
      Int_t offCopyNo = -1;
      TString offName = fMC->CurrentVolOffName(off);
      Int_t offId = fMC->CurrentVolOffID(off, offCopyNo);
      cout << "  off=" << off << " name=" << offName << " id=" << offId
           << " copyNo=" << offCopyNo << endl;
    }
  }
 
  const Bool_t reflected = fDetector->GetUseReflection() && copyNo == 2;
  fMC->CurrentVolOffID(2, copyNo);
  // cout << "Got copyNo "<< copyNo << " " << fMC->CurrentVolPath() << endl;
 
  Double_t edep = fMC->Edep();
 
  Double_t step = 0.;
  if (fMC->TrackCharge() != 0.) step = fMC->TrackStep();
 
  if (!GetHit(copyNo)) {
    std::cerr << "No hit found for layer with copyNo = " << copyNo << endl;
    return false;
  }
 
  if (id == fAbsorberVolId) {
    GetHit(copyNo)->AddAbs(edep, step);
    if (reflected) fReflectedHit.AddAbs(edep, step);
  }
 
  if (id == fGapVolId) {
    GetHit(copyNo)->AddGap(edep, step);
    if (reflected) fReflectedHit.AddGap(edep, step);
  }
 
  return true;
}
 
//_____________________________________________________________________________
void Ex03CalorimeterSD::EndOfEvent()
{
  /// Print hits collection (if verbose) and reset hits afterwards.
 
  if (fVerboseLevel > 1) Print();
 
  if (fDetector->GetUseReflection()) {
    cout << "   Reflected placements (included in layer totals):" << endl;
    fReflectedHit.Print();
  }
  if (fDetector->GetRequireReflectedHits()) {
    if (fReflectedHit.GetEdepAbs() <= 0. || fReflectedHit.GetTrakAbs() <= 0. ||
        fReflectedHit.GetEdepGap() <= 0. || fReflectedHit.GetTrakGap() <= 0.) {
      Fatal("EndOfEvent", "Missing hits in reflected absorber or gap");
    }
    cout << "   Reflected hits test passed" << endl;
  }
  fReflectedHit.Reset();
 
  // Reset hits collection
  ResetHits();
}
 
//_____________________________________________________________________________
void Ex03CalorimeterSD::Register()
{
  /// Register the hits collection in Root manager.
 
  TMCRootManager::Instance()->Register("hits", "TClonesArray", &fCalCollection);
}
 
//_____________________________________________________________________________
void Ex03CalorimeterSD::Print(Option_t* /*option*/) const
{
  /// Print the hits collection.
 
  Int_t nofHits = fCalCollection->GetEntriesFast();
 
  cout << "\n-------->Hits Collection: in this event: " << endl;
 
  for (Int_t i = 0; i < nofHits; i++) (*fCalCollection)[i]->Print();
}
 
//_____________________________________________________________________________
void Ex03CalorimeterSD::PrintTotal() const
{
  /// Print the total values for all layers.
 
  Double_t totEAbs = 0.;
  Double_t totLAbs = 0.;
  Double_t totEGap = 0.;
  Double_t totLGap = 0.;
 
  Int_t nofHits = fCalCollection->GetEntriesFast();
  for (Int_t i = 0; i < nofHits; i++) {
    totEAbs += GetHit(i)->GetEdepAbs();
    totLAbs += GetHit(i)->GetTrakAbs();
    totEGap += GetHit(i)->GetEdepGap();
    totLGap += GetHit(i)->GetTrakGap();
  }
 
  cout << "   Absorber: total energy (MeV): " << setw(7) << totEAbs * 1.0e03
       << "       total track length (cm):  " << setw(7) << totLAbs << endl
       << "   Gap:      total energy (MeV): " << setw(7) << totEGap * 1.0e03
       << "       total track length (cm):  " << setw(7) << totLGap << endl;
}