/data/refman/pasoursint/CMSSW_5_3_6/src/SimG4CMS/CherenkovAnalysis/interface/DreamSD.h

Go to the documentation of this file.

#ifndef SimG4CMS_DreamSD_h
#define SimG4CMS_DreamSD_h

#include "SimG4CMS/Calo/interface/CaloSD.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "G4String.hh"
#include "G4PhysicsOrderedFreeVector.hh"

#include <map>

const int MAXPHOTONS = 500; // Maximum number of photons we can store

class G4LogicalVolume;

class TTree;

class DreamSD : public CaloSD {

public:    

  DreamSD(G4String, const DDCompactView &, SensitiveDetectorCatalog &, 
          edm::ParameterSet const &, const SimTrackManager*);
  virtual ~DreamSD() {}
  virtual bool   ProcessHits(G4Step * step,G4TouchableHistory * tHistory);
  virtual uint32_t setDetUnitId(G4Step*);

protected:

  virtual G4bool getStepInfo(G4Step* aStep);
  virtual void   initRun();

private:    

  typedef std::pair<double,double> Doubles;
  typedef std::map<G4LogicalVolume*,Doubles> DimensionMap;

  void           initMap(G4String, const DDCompactView &);
  double         curve_LY(G4Step*, int); 
  const double   crystalLength(G4LogicalVolume*) const;
  const double   crystalWidth(G4LogicalVolume*) const;

  double         cherenkovDeposit_( G4Step* aStep );
  const double   getAverageNumberOfPhotons_(const double charge,
                                            const double beta,
                                            const G4Material* aMaterial,
                                            const G4MaterialPropertyVector* rIndex ) const;
  const double getPhotonEnergyDeposit_( const G4ParticleMomentum& p, const G4ThreeVector& x,
                                        const G4Step* aStep );
  bool setPbWO2MaterialProperties_( G4Material* aMaterial );

  bool         useBirk, doCherenkov_, readBothSide_;
  double       birk1, birk2, birk3;
  double       slopeLY;
  DimensionMap xtalLMap; // Store length and width

  int          side;

  std::auto_ptr<G4PhysicsOrderedFreeVector> chAngleIntegrals_;
  G4MaterialPropertiesTable*                materialPropertiesTable;
  // Histogramming
  TTree* ntuple_;
  int nphotons_;
  float px_[MAXPHOTONS],py_[MAXPHOTONS],pz_[MAXPHOTONS];
  float x_[MAXPHOTONS],y_[MAXPHOTONS],z_[MAXPHOTONS];

};

#endif // DreamSD_h