CDSimpleMC.h

Go to the documentation of this file.

// Program name: CDSimpleMC.h 
//                                                                     
// Package: CalDetTracker
//                                                                    
// Coded by Jeff Hartnell Mar/2004          
//                                                                    
// Purpose: To allow a simple MC simulation of stopping muons 
//                                                                    
// Contact: jeffrey.hartnell@physics.ox.ac.uk                         
//  Modified by ADM in April 2006 to use MinosMaterial  (in Conventions) rather than
//  StandardMaterial in MuELoss

#ifndef CDSIMPLEMC_H
#define CDSIMPLEMC_H

#include <map>
#include <string>
#include <vector>

#include "TMaterial.h"

#include "MuELoss/MuEnergyLoss.h"

class CDSimpleMC : public TObject
{
  
 public:  
  
  CDSimpleMC();
  ~CDSimpleMC();

  //methods that modify the object
  void RunMC();
  void CalcRatioScToFe();
  void SetThkFe(Double_t thkFe);
  void SetParticleMomentum(Double_t pMom) {fParticleMomentum=pMom;}
  
  //methods that do not modify the object
  Double_t GetParticleEnergy() const;
  Double_t GetThkFe() {return fThkFe;}
  Double_t GetTotalNumScPlanesHit() {return fTotalNumScPlanesHit;}
  Double_t GetTotalScEnLoss() {return fTotalScEnLoss;}
  void PrintDensities() const;
  void PrintThicknesses() const;

  //data members
  //must make these private sometime...
  std::vector<Double_t> fEn;
  std::vector<Double_t> fMom;//muon momentum
  std::vector<Double_t> fMomGeV;//muon momentum
  std::vector<Double_t> fdE_dxFe;
  std::vector<Double_t> fdE_dxAl;
  std::vector<Double_t> fdE_dxSc;

  std::vector<Double_t> fEnDeposited;
  std::vector<Double_t> fX;
  std::vector<Double_t> fXPlanes;

  std::vector<Double_t> fEnFeDeposited;
  std::vector<Double_t> fXFe;
  std::vector<Double_t> fXFePlanes;

  std::vector<Double_t> fEnAlDeposited;
  std::vector<Double_t> fXAl;
  std::vector<Double_t> fXAlPlanes;

  std::vector<Double_t> fEnScTiO2Deposited;
  std::vector<Double_t> fXScTiO2;
  std::vector<Double_t> fXScTiO2Planes;

  std::vector<Double_t> fEnNotScDeposited;
  std::vector<Double_t> fXNotSc;
  std::vector<Double_t> fXNotScPlanes;

  std::vector<Double_t> fEnScDeposited;
  std::vector<Double_t> fEnSc;//muon energy in the scintillator 
  std::vector<Double_t> fMomSc;//muon momentum in the scintillator 
  std::vector<Double_t> fMomScGeV;//muon momentum in the scintillator 
  std::vector<Double_t> fXSc;
  std::vector<Double_t> fXScPlanes;

  std::vector<Double_t> fEnDepositedPl;

  std::vector<Double_t> fEnRatio;
  std::vector<Double_t> fXRatio;
  std::vector<Double_t> fXPlanesRatio;

  //vectors to store the energy deposited in both Al and TiO2 layers
  std::vector<Double_t> fEnTiO2x2Deposited;
  std::vector<Double_t> fXTiO2x2;
  std::vector<Double_t> fEnAlx2Deposited;
  std::vector<Double_t> fXAlx2;

 private:

  //methods that do modify the object
  void CalcPlaneThk();

  //methods that do not modify the object
  Double_t dE_dxIterative(MuEnergyLoss& mat,
                          Double_t energyRemaining,
                          Double_t totalThk,
                          Double_t rho,
                          Double_t energyCutOff,
                          Int_t slices,
                          Double_t& fractionThrough) const;
  //data members
  Double_t fParticleMomentum;
  Double_t fTotalScEnLoss;
  Double_t fTotalNumScPlanesHit;
  std::string fS;//general purpose string

  //hold standard materials
  MinosMaterial::StdMaterial_t fStdMatFe;
  MinosMaterial::StdMaterial_t fStdMatAl;
  MinosMaterial::StdMaterial_t fStdMatSc;
  MinosMaterial::StdMaterial_t fStdMatScTiO2;

  //hold the process to use to calc dE_dx etc
  Process::Process_t fELossProcess;

  //all the densities
  Float_t fRoFe;
  Float_t fRoAl;
  Float_t fRoScTiO2;
  Float_t fRoSc;
  
  Float_t fMuMass;
  
  //variables to hold the thicknesses of materials
  Float_t fThkFe;
  Float_t fThkAl;
  Float_t fThkSc;
  Float_t fThkScTiO2;
  Float_t fThkAir;
  Float_t fThkPlane;//whole plane thickness
  
  //standard materials
  Material* fFe;
  Material* fAl;
  Material* fSc;
  
  //energyLoss objects
  MuEnergyLoss* fIonFe;
  MuEnergyLoss* fIonAl;
  MuEnergyLoss* fIonSc;
  MuEnergyLoss* fIonScTiO2;

  Int_t fSlices;//number of slices to use in traversing material
  ModelSelector fModelSelector;//use default model selections

  ClassDef(CDSimpleMC,0)
};

#endif // CDSIMPLEMC_H

---