CalorimetryManager Class Reference

#include <CalorimetryManager.h>

Public Member Functions
	CalorimetryManager ()
	CalorimetryManager (FSimEvent *aSimEvent, const edm::ParameterSet &fastCalo, const edm::ParameterSet &MuonECALPars, const edm::ParameterSet &MuonHCALPars, const edm::ParameterSet &fastGflash)
CaloGeometryHelper *	getCalorimeter () const
FastHFShowerLibrary *	getHFShowerLibrary () const
void	loadFromEcalBarrel (edm::PCaloHitContainer &c) const
void	loadFromEcalEndcap (edm::PCaloHitContainer &c) const
void	loadFromHcal (edm::PCaloHitContainer &c) const
void	loadFromPreshower (edm::PCaloHitContainer &c) const
void	loadMuonSimTracks (edm::SimTrackContainer &m) const
void	reconstruct (RandomEngineAndDistribution const *)
	~CalorimetryManager ()

Private Member Functions
void	clean ()
void	EMShowerSimulation (const FSimTrack &myTrack, RandomEngineAndDistribution const *)
void	HDShowerSimulation (const FSimTrack &myTrack, RandomEngineAndDistribution const *)
	Hadronic Shower Simulation. More...
void	MuonMipSimulation (const FSimTrack &myTrack, RandomEngineAndDistribution const *)
void	readParameters (const edm::ParameterSet &fastCalo)
void	reconstructHCAL (const FSimTrack &myTrack, RandomEngineAndDistribution const *)
void	respCorr (double)
void	updateECAL (const std::map< CaloHitID, float > &hitMap, int onEcal, int trackID=0, float corr=1.0)
void	updateHCAL (const std::map< CaloHitID, float > &hitMap, int trackID=0, float corr=1.0)
void	updatePreshower (const std::map< CaloHitID, float > &hitMap, int trackID=0, float corr=1.0)

Private Attributes
GammaFunctionGenerator *	aGammaGenerator
const LandauFluctuationGenerator *	aLandauGenerator
double	aTerm
bool	bFixedLength_
double	bTerm
double	crackPadSurvivalProbability_
bool	debug_
std::vector< std::pair< CaloHitID, float > >	EBMapping_
std::unique_ptr< KKCorrectionFactors >	ecalCorrection
bool	EcalDigitizer_
double	ecorr
std::vector< std::pair< CaloHitID, float > >	EEMapping_
std::vector< std::pair< CaloHitID, float > >	ESMapping_
std::vector< unsigned int >	evtsToDebug_
int	gridSize_
bool	HcalDigitizer_
double	hcorr
int	hdGridSize_
int	hdSimMethod_
std::vector< std::pair< CaloHitID, float > >	HMapping_
int	ietaShiftHB_
int	ietaShiftHE_
int	ietaShiftHF_
int	ietaShiftHO_
bool	initialized_
std::vector< double >	k_e
std::vector< double >	k_h
std::vector< double >	mipValues_
std::vector< FSimTrack >	muonSimTracks
CaloGeometryHelper *	myCalorimeter_
RawParticle	myElec
	A few pointers to save time. More...
HCALResponse *	myHDResponse_
Histos *	myHistos
HSParameters *	myHSParameters_
RawParticle	myPart
RawParticle	myPosi
FSimEvent *	mySimEvent
int	optionHDSim_
std::vector< double >	p_knots
double	pulledPadSurvivalProbability_
double	radiusFactorEB_
double	radiusFactorEE_
std::vector< double >	radiusPreshowerCorrections_
double	RCFactor_
std::vector< double >	rsp
double	RTFactor_
std::vector< double >	samplingHBHE_
std::vector< double >	samplingHF_
std::vector< double >	samplingHO_
bool	simulatePreshower_
double	spotFraction_
GflashAntiProtonShowerProfile *	theAntiProtonProfile
std::vector< double >	theCoreIntervals_
FastHFShowerLibrary *	theHFShowerLibrary
MaterialEffects *	theMuonEcalEffects
MaterialEffects *	theMuonHcalEffects
GflashPiKShowerProfile *	thePiKProfile
GflashHadronShowerProfile *	theProfile
GflashProtonShowerProfile *	theProtonProfile
std::vector< double >	theTailIntervals_
std::vector< double >	timeShiftHB_
std::vector< double >	timeShiftHE_
std::vector< double >	timeShiftHF_
std::vector< double >	timeShiftHO_
bool	unfoldedMode_
bool	useCorrectionSL
bool	useShowerLibrary

Static Private Attributes
static std::vector< std::pair< int, float > >	myZero_

Detailed Description

Definition at line 46 of file CalorimetryManager.h.

Constructor & Destructor Documentation

CalorimetryManager::CalorimetryManager

(

)

Definition at line 71 of file CalorimetryManager.cc.

                                       : 
  myCalorimeter_(nullptr),
  //  myHistos(0),
  initialized_(false)
{;}

CalorimetryManager::CalorimetryManager	(	FSimEvent *	aSimEvent,
		const edm::ParameterSet &	fastCalo,
		const edm::ParameterSet &	MuonECALPars,
		const edm::ParameterSet &	MuonHCALPars,
		const edm::ParameterSet &	fastGflash
	)

Definition at line 77 of file CalorimetryManager.cc.

References aGammaGenerator, aLandauGenerator, ecalCorrection, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), TrajectoryFactories_cff::MaterialEffects, myCalorimeter_, myHDResponse_, myHSParameters_, readParameters(), theAntiProtonProfile, theHFShowerLibrary, theMuonEcalEffects, theMuonHcalEffects, thePiKProfile, theProfile, and theProtonProfile.

  : 
  mySimEvent(aSimEvent), 
  initialized_(false),
  theMuonEcalEffects(nullptr), theMuonHcalEffects(nullptr), bFixedLength_(false)
{
  
  aLandauGenerator = new LandauFluctuationGenerator;
  aGammaGenerator = new GammaFunctionGenerator;
  
  //Gflash
  theProfile = new GflashHadronShowerProfile(parGflash);
  thePiKProfile = new GflashPiKShowerProfile(parGflash);
  theProtonProfile = new GflashProtonShowerProfile(parGflash);
  theAntiProtonProfile = new GflashAntiProtonShowerProfile(parGflash);
  
  // FastHFShowerLibrary
  theHFShowerLibrary = new FastHFShowerLibrary(fastCalo);
  
  readParameters(fastCalo);
  
  //  myHistos = 0; 
  
  //   myHistos = Histos::instance();
  //   myHistos->book("h10",140,-3.5,3.5,100,-0.5,99.5);
  //   myHistos->book("h20",150,0,150.,100,-0.5,99.5);
  //   myHistos->book("h100",140,-3.5,3.5,100,0,0.1);
  //   myHistos->book("h110",140,-3.5,3.5,100,0,10.);
  //   myHistos->book("h120",200,-5.,5.,100,0,0.5);
  
  //   myHistos->book("h200",300,0,3.,100,0.,35.);
  //   myHistos->book("h210",720,-M_PI,M_PI,100,0,35.);
  //   myHistos->book("h212",720,-M_PI,M_PI,100,0,35.);
  
  //   myHistos->bookByNumber("h30",0,7,300,-3.,3.,100,0.,35.);
  //   myHistos->book("h310",75,-3.,3.,"");
  //   myHistos->book("h400",100,-10.,10.,100,0.,35.);
  //   myHistos->book("h410",720,-M_PI,M_PI);
  
  myCalorimeter_ = 
    new CaloGeometryHelper(fastCalo);
  myHDResponse_ = 
    new HCALResponse(fastCalo.getParameter<edm::ParameterSet>("HCALResponse"));
  myHSParameters_ = 
    new HSParameters(fastCalo.getParameter<edm::ParameterSet>("HSParameters"));
  
  // Material Effects for Muons in ECAL (only EnergyLoss implemented so far)
  
  if ( fastMuECAL.getParameter<bool>("PairProduction") || 
       fastMuECAL.getParameter<bool>("Bremsstrahlung") ||
       fastMuECAL.getParameter<bool>("MuonBremsstrahlung") ||
       fastMuECAL.getParameter<bool>("EnergyLoss") || 
       fastMuECAL.getParameter<bool>("MultipleScattering") )
    theMuonEcalEffects = new MaterialEffects(fastMuECAL);
  
  // Material Effects for Muons in HCAL (only EnergyLoss implemented so far)
  
  if ( fastMuHCAL.getParameter<bool>("PairProduction") || 
       fastMuHCAL.getParameter<bool>("Bremsstrahlung") ||
       fastMuHCAL.getParameter<bool>("MuonBremsstrahlung") ||
       fastMuHCAL.getParameter<bool>("EnergyLoss") || 
       fastMuHCAL.getParameter<bool>("MultipleScattering") )
    theMuonHcalEffects = new MaterialEffects(fastMuHCAL);

  if( fastCalo.exists("ECALResponseScaling") ) {
    ecalCorrection = std::unique_ptr<KKCorrectionFactors>( new KKCorrectionFactors( fastCalo.getParameter<edm::ParameterSet>("ECALResponseScaling") ) );
  }

}

CalorimetryManager::~CalorimetryManager

(

)

Definition at line 160 of file CalorimetryManager.cc.

References myCalorimeter_, myHDResponse_, theHFShowerLibrary, theMuonEcalEffects, theMuonHcalEffects, and theProfile.

{
  if(myCalorimeter_) delete myCalorimeter_;
  if(myHDResponse_) delete myHDResponse_;
  
  if ( theMuonEcalEffects ) delete theMuonEcalEffects;
  if ( theMuonHcalEffects ) delete theMuonHcalEffects;
  
  if ( theProfile ) delete theProfile;
  
  if ( theHFShowerLibrary ) delete theHFShowerLibrary;
}

Member Function Documentation

void CalorimetryManager::clean ( )

private

Definition at line 151 of file CalorimetryManager.cc.

References EBMapping_, EEMapping_, ESMapping_, HMapping_, and muonSimTracks.

Referenced by reconstruct().

{
  EBMapping_.clear();
  EEMapping_.clear();
  HMapping_.clear();
  ESMapping_.clear();
  muonSimTracks.clear();
}

void CalorimetryManager::EMShowerSimulation ( const FSimTrack &  myTrack, RandomEngineAndDistribution const *  random  )

private

Definition at line 252 of file CalorimetryManager.cc.

References funct::abs(), aGammaGenerator, aLandauGenerator, aTerm, bFixedLength_, bTerm, EMShower::compute(), crackPadSurvivalProbability_, debug_, particleFlowClusterECALTimeSelected_cfi::depth, GetRecoTauVFromDQM_MC_cff::dir2, MillePedeFileConverter_cfg::e, ecalCorrection, FSimTrack::ecalEntrance(), Calorimeter::ecalProperties(), RawParticle::eta(), RandomEngineAndDistribution::flatShoot(), CaloGeometryHelper::getClosestCell(), PreshowerHitMaker::getHits(), HcalHitMaker::getHits(), EcalHitMaker::getHits(), EMShower::getMaximumOfShower(), gridSize_, FSimTrack::hcalEntrance(), Calorimeter::hcalProperties(), FSimTrack::id(), FSimTrack::layer1Entrance(), Calorimeter::layer1Properties(), FSimTrack::layer2Entrance(), Calorimeter::layer2Properties(), cmsBatch::log, particleFlowClusterECALTimeSelected_cfi::maxEnergy, mipValues_, myCalorimeter_, myElec, myPart, myPosi, FSimTrack::onEcal(), FSimTrack::onHcal(), FSimTrack::onLayer1(), FSimTrack::onLayer2(), pulledPadSurvivalProbability_, radiusFactorEB_, radiusFactorEE_, ECALProperties::radLenIncm(), RCFactor_, RTFactor_, Scenarios_cff::scale, EcalHitMaker::setCrackPadSurvivalProbability(), EMShower::setGrid(), EMShower::setHcal(), EMShower::setPreshower(), EcalHitMaker::setPreshowerPresent(), EcalHitMaker::setPulledPadSurvivalProbability(), EcalHitMaker::setRadiusFactor(), EcalHitMaker::setTrackParameters(), RawParticle::setVertex(), simulatePreshower_, findQualityFiles::size, theCoreIntervals_, theTailIntervals_, CoreSimTrack::type(), updateECAL(), updateHCAL(), updatePreshower(), and RawParticle::vertex().

Referenced by reconstruct().

                                                                                       {
  std::vector<const RawParticle*> thePart;
  double X0depth;
  
  if (debug_) {
    LogInfo("FastCalorimetry") << " EMShowerSimulation "  <<myTrack << std::endl;      
  }
  
  //  std::cout << " Simulating " << myTrack << std::endl;
  
  // The Particle at ECAL entrance
  //  std::cout << " Before ecalEntrance " << std::endl;
  myPart = myTrack.ecalEntrance(); 
  
  // protection against infinite loop.  
  if ( myTrack.type() == 22 && myPart.e()<0.055) return; 
  
  
  // Barrel or Endcap ?
  int onEcal = myTrack.onEcal();
  int onHcal = myTrack.onHcal();
  int onLayer1 = myTrack.onLayer1();
  int onLayer2 = myTrack.onLayer2();
  
  // The entrance in ECAL
  XYZPoint ecalentrance = myPart.vertex().Vect();
  
  //  std::cout << " Ecal entrance " << ecalentrance << std::endl;
  
  // The preshower
  PreshowerHitMaker * myPreshower = nullptr;
  if(simulatePreshower_ && (onLayer1 || onLayer2))
    {
      XYZPoint layer1entrance,layer2entrance;
      XYZVector dir1,dir2;
      if(onLayer1) 
    {
      layer1entrance = XYZPoint(myTrack.layer1Entrance().vertex().Vect());
      dir1 = XYZVector(myTrack.layer1Entrance().Vect().Unit());
    }
      if(onLayer2) 
    {
      layer2entrance = XYZPoint(myTrack.layer2Entrance().vertex().Vect());
      dir2 = XYZVector(myTrack.layer2Entrance().Vect().Unit());
    }
      //      std::cout << " Layer1entrance " << layer1entrance << std::endl;
      //      std::cout << " Layer2entrance " << layer2entrance << std::endl;
      myPreshower = new PreshowerHitMaker(myCalorimeter_,
                      layer1entrance,
                      dir1,
                      layer2entrance,
                      dir2,
                      aLandauGenerator,
                                          random);
      myPreshower->setMipEnergy(mipValues_[0],mipValues_[1]);
    }
  
  // The ECAL Properties
  EMECALShowerParametrization 
    showerparam(myCalorimeter_->ecalProperties(onEcal), 
        myCalorimeter_->hcalProperties(onHcal), 
        myCalorimeter_->layer1Properties(onLayer1), 
        myCalorimeter_->layer2Properties(onLayer2),
        theCoreIntervals_,
        theTailIntervals_,
        RCFactor_,
        RTFactor_);
  
  // Photons : create an e+e- pair
  if ( myTrack.type() == 22 ) {
    
    // Depth for the first e+e- pair creation (in X0)
    X0depth = -log(random->flatShoot()) * (9./7.);
    
    // Initialization
    double eMass = 0.000510998902; 
    double xe=0;
    double xm=eMass/myPart.e();
    double weight = 0.;
    
    // Generate electron energy between emass and eGamma-emass
    do {
      xe = random->flatShoot()*(1.-2.*xm) + xm;
      weight = 1. - 4./3.*xe*(1.-xe);
    } while ( weight < random->flatShoot() );
    
    // Protection agains infinite loop in Famos Shower
    if ( myPart.e()*xe < 0.055 || myPart.e()*(1.-xe) < 0.055 ) {
      
      if ( myPart.e() > 0.055 ) thePart.push_back(&myPart);
      
    } else {      
      
      myElec = (myPart) * xe;
      myPosi = (myPart) * (1.-xe);
      myElec.setVertex(myPart.vertex());
      myPosi.setVertex(myPart.vertex());
      thePart.push_back(&myElec);
      thePart.push_back(&myPosi);
    }
    // Electrons
  } else {  
    
    X0depth = 0.;
    if ( myPart.e() > 0.055 ) thePart.push_back(&myPart);
    
  } 
  
  // After the different protections, this shouldn't happen. 
  if(thePart.empty()) 
    { 
      if(myPreshower==nullptr) return; 
      delete myPreshower; 
      return; 
    } 
  
  // find the most energetic particle
  double maxEnergy=-1.;
  for(unsigned ip=0;ip < thePart.size();++ip)
    if(thePart[ip]->e() > maxEnergy) maxEnergy = thePart[ip]->e();
  
  // Initialize the Grid in ECAL
  int size = gridSize_;
  if(maxEnergy>100) size=11;
  //  if ( maxEnergy < threshold5x5 ) size = 5;
  //  if ( maxEnergy < threshold3x3 ) size = 3;
  
  
  EMShower theShower(random,aGammaGenerator,&showerparam,&thePart,nullptr,nullptr,bFixedLength_);
  
  
  double maxShower = theShower.getMaximumOfShower();
  if (maxShower > 20.) maxShower = 2.; // simple pivot-searching protection 
  
  double depth((X0depth + maxShower) * 
           myCalorimeter_->ecalProperties(onEcal)->radLenIncm());
  XYZPoint meanShower = ecalentrance + myPart.Vect().Unit()*depth;
  
  //  if(onEcal!=1) return ; 
  
  // The closest crystal
  DetId pivot(myCalorimeter_->getClosestCell(meanShower, true, onEcal==1));
  
  if(pivot.subdetId() == 0) {   // further protection against avbsence of pivot
    edm::LogWarning("CalorimetryManager") <<  "Pivot for egamma  e = "  << myTrack.hcalEntrance().e() << " is not found at depth " << depth << " and meanShower coordinates = " << meanShower << std::endl; 
    if(myPreshower) delete myPreshower;
    return;
  }
  
  EcalHitMaker myGrid(myCalorimeter_,ecalentrance,pivot,onEcal,size,0,random);
  //                                             ^^^^
  //                                         for EM showers
  myGrid.setPulledPadSurvivalProbability(pulledPadSurvivalProbability_);
  myGrid.setCrackPadSurvivalProbability(crackPadSurvivalProbability_);
  
  //maximumdepth dependence of the radiusfactorbehindpreshower
  //First tuning: Shilpi Jain (Mar-Apr 2010); changed after tuning - Feb-July - Shilpi Jain
  /* **************
     myGrid.setRadiusFactor(radiusFactor_);
     if(onLayer1 || onLayer2)
     {
     float b               = radiusPreshowerCorrections_[0];
     float a               = radiusFactor_*( 1.+radiusPreshowerCorrections_[1]*radiusPreshowerCorrections_[0] );
     float maxdepth        = X0depth+theShower.getMaximumOfShower();
     float newRadiusFactor = radiusFactor_;
     if(myPart.e()<=250.)
     {
     newRadiusFactor = a/(1.+b*maxdepth); 
     }
     myGrid.setRadiusFactor(newRadiusFactor);
     }
     else // otherwise use the normal radius factor
     {
     myGrid.setRadiusFactor(radiusFactor_);
     }
     ************** */
  if(myTrack.onEcal() == 2) // if on EE  
    {
      if( (onLayer1 || onLayer2) && myPart.e()<=250.)
    {
      double maxdepth        = X0depth+theShower.getMaximumOfShower();
      double newRadiusFactor = radiusFactorEE_ * aTerm/(1.+bTerm*maxdepth);
      myGrid.setRadiusFactor(newRadiusFactor);
    }
      else // otherwise use the normal radius factor
    {
      myGrid.setRadiusFactor(radiusFactorEE_);
    }
    }//if(myTrack.onEcal() == 2)
  else                      // else if on EB
    {
      myGrid.setRadiusFactor(radiusFactorEB_);
    }
  //(end of) changed after tuning - Feb-July - Shilpi Jain
  
  myGrid.setPreshowerPresent(simulatePreshower_);
  
  // The shower simulation
  myGrid.setTrackParameters(myPart.Vect().Unit(),X0depth,myTrack);
  
  //  std::cout << " PS ECAL GAP HCAL X0 " << myGrid.ps1TotalX0()+myGrid.ps2TotalX0() << " " << myGrid.ecalTotalX0();
  //  std::cout << " " << myGrid.ecalHcalGapTotalX0() << " " << myGrid.hcalTotalX0() << std::endl;
  //  std::cout << " PS ECAL GAP HCAL L0 " << myGrid.ps1TotalL0()+myGrid.ps2TotalL0() << " " << myGrid.ecalTotalL0();
  //  std::cout << " " << myGrid.ecalHcalGapTotalL0() << " " << myGrid.hcalTotalL0() << std::endl;
  //  std::cout << "ECAL-HCAL " << myTrack.momentum().eta() << " " <<  myGrid.ecalHcalGapTotalL0() << std::endl;
  //
  //  std::cout << " Grid created " << std::endl;
  if(myPreshower) theShower.setPreshower(myPreshower);
  
  HcalHitMaker myHcalHitMaker(myGrid,(unsigned)0); 
  
  theShower.setGrid(&myGrid);
  theShower.setHcal(&myHcalHitMaker);
  theShower.compute();
  //myHistos->fill("h502", myPart->eta(),myGrid.totalX0());

  // calculate the total simulated energy for this particle
  float simE = 0;
  for( const auto& mapIterator : myGrid.getHits() ) {
    simE += mapIterator.second;
  }

  auto scale = ecalCorrection ? ecalCorrection->getScale( myTrack.ecalEntrance().e(),
     std::abs( myTrack.ecalEntrance().eta() ), simE ) : 1.;

  // Save the hits !
  updateECAL( myGrid.getHits(), onEcal,myTrack.id(), scale );

  // Now fill the HCAL hits
  updateHCAL(myHcalHitMaker.getHits(),myTrack.id());
  
  // delete the preshower
  if(myPreshower!=nullptr) {
    updatePreshower(myPreshower->getHits(),myTrack.id());
    delete myPreshower;
    //  std::cout << " Deleting myPreshower " << std::endl;
  }
  
}

CaloGeometryHelper* CalorimetryManager::getCalorimeter ( ) const

inline

Definition at line 61 of file CalorimetryManager.h.

Referenced by FamosManager::setupGeometryAndField().

{return myCalorimeter_;}

FastHFShowerLibrary* CalorimetryManager::getHFShowerLibrary ( ) const

inline

Definition at line 64 of file CalorimetryManager.h.

References EnergyCorrector::c, clean, corr, and funct::m.

Referenced by FamosManager::setupGeometryAndField().

{return theHFShowerLibrary;}

void CalorimetryManager::HDShowerSimulation ( const FSimTrack &  myTrack, RandomEngineAndDistribution const *  random  )

private

Hadronic Shower Simulation.

Definition at line 557 of file CalorimetryManager.cc.

References funct::abs(), HcalHitMaker::addHit(), FSimTrack::charge(), ALCARECOTkAlJpsiMuMu_cff::charge, HFShower::compute(), HDShower::compute(), HDRShower::computeShower(), HCALResponse::correctHF(), debug_, FSimTrack::ecalEntrance(), Calorimeter::ecalProperties(), ecorr, objects.autophobj::float, Gflash::getCalorimeterNumber(), CaloGeometryHelper::getClosestCell(), GflashHadronShowerProfile::getGflashHitList(), GflashHadronShowerProfile::getGflashShowino(), GflashTrajectory::getGflashTrajectoryPoint(), HCALResponse::getHCALEnergyResponse(), Calorimeter::getHcalGeometry(), GflashShowino::getHelix(), HcalHitMaker::getHits(), FastHFShowerLibrary::getHitsMap(), HDShower::getmip(), HCALResponse::getMIPfraction(), GflashShowino::getPathLengthAtShower(), GflashShowino::getPathLengthOnEcal(), GeV, GflashHadronShowerProfile::hadronicParameterization(), FSimTrack::hcalEntrance(), Calorimeter::hcalProperties(), hcorr, hdGridSize_, hdSimMethod_, FSimTrack::id(), GflashHadronShowerProfile::initialize(), Gflash::intLength, Gflash::kENCA, Gflash::kESPM, Gflash::kHB, Gflash::kHE, Gflash::kNULL, GflashHadronShowerProfile::loadParameters(), SiStripPI::max, FSimTrack::momentum(), myCalorimeter_, myHDResponse_, myHSParameters_, myPart, DetId::null(), FSimTrack::onEcal(), FSimTrack::onHcal(), FSimTrack::onVFcal(), optionHDSim_, objects.autophobj::particleType, random, DetId::rawId(), FastHFShowerLibrary::recoHFShowerLibrary(), respCorr(), HCALResponse::responseHCAL(), HcalHitMaker::setDepth(), FastHFShowerLibrary::SetRandom(), HcalHitMaker::setSpotEnergy(), EcalHitMaker::setTrackParameters(), mps_update::status, theAntiProtonProfile, theHFShowerLibrary, thePiKProfile, theProfile, theProtonProfile, CoreSimTrack::type(), updateECAL(), updateHCAL(), useShowerLibrary, RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by reconstruct().

                                                                                                              {//,
  // const edm::ParameterSet& fastCalo){
  
  theHFShowerLibrary->SetRandom(random);

  //  TimeMe t(" FASTEnergyReconstructor::HDShower");
  const XYZTLorentzVector& moment = myTrack.momentum();
  
  if(debug_)
    LogInfo("FastCalorimetry") 
      << "CalorimetryManager::HDShowerSimulation - track param."
      << std::endl
      << "  eta = " << moment.eta() << std::endl
      << "  phi = " << moment.phi() << std::endl
      << "   et = " << moment.Et()  << std::endl
      << "   e  = " << myTrack.hcalEntrance().e() << std::endl;
  
  if (debug_) {
    LogInfo("FastCalorimetry") << " HDShowerSimulation "  << myTrack << std::endl;      
  }
  
  
  int hit;
  //  int pid = abs(myTrack.type());
  
  XYZTLorentzVector trackPosition;
  if ( myTrack.onEcal() ) {
    trackPosition=myTrack.ecalEntrance().vertex();
    hit = myTrack.onEcal()-1;                               //
    myPart = myTrack.ecalEntrance();
  } else if ( myTrack.onVFcal()) {
    trackPosition=myTrack.vfcalEntrance().vertex();
    hit = 2;
    myPart = myTrack.vfcalEntrance();
  }
  else
    {
      LogInfo("FastCalorimetry") << " The particle is not in the acceptance " << std::endl;
      return;
    }
  
  // int onHCAL = hit + 1; - specially for myCalorimeter->hcalProperties(onHCAL)
  // (below) to get VFcal properties ...
  int onHCAL = hit + 1;
  int onECAL = myTrack.onEcal();
  
  double pathEta   = trackPosition.eta();
  double pathPhi   = trackPosition.phi();   
  //  double pathTheta = trackPosition.theta();
  
  double eint  = moment.e();
  double eGen  = myTrack.hcalEntrance().e();
  
  double emeas = 0.;
  double pmip= myHDResponse_->getMIPfraction(eGen, pathEta);
  //  std::cout << " CalorimetryManager onHcal " <<  pmip  << std::endl;
  
  
  //double emeas = -0.000001; 
  
  //===========================================================================
  if(eGen > 0.) {  
    
    // ECAL and HCAL properties to get
    HDShowerParametrization 
      theHDShowerparam(myCalorimeter_->ecalProperties(onECAL),
               myCalorimeter_->hcalProperties(onHCAL),
               myHSParameters_);
    
    //Making ECAL Grid (and segments calculation)
    XYZPoint caloentrance;
    XYZVector direction;
    if(myTrack.onEcal()) 
      { 
    caloentrance = myTrack.ecalEntrance().vertex().Vect();
    direction = myTrack.ecalEntrance().Vect().Unit();
      }
    else if(myTrack.onHcal())
      {
    caloentrance = myTrack.hcalEntrance().vertex().Vect();
    direction = myTrack.hcalEntrance().Vect().Unit();
      }
    else
      {
    caloentrance = myTrack.vfcalEntrance().vertex().Vect();
    direction = myTrack.vfcalEntrance().Vect().Unit();
      }
    
    if(debug_)
      LogInfo("FastCalorimetry") 
    << "CalorimetryManager::HDShowerSimulation - on-calo 1 "
    << std::endl
    << "  onEcal    = " <<  myTrack.onEcal()  << std::endl
    << "  onHcal    = " <<  myTrack.onHcal()  << std::endl
    << "  onVFcal   = " <<  myTrack.onVFcal() << std::endl
    << "  position  = " << caloentrance << std::endl;
    
    
    DetId pivot;
    if(myTrack.onEcal())
      {
    pivot=myCalorimeter_->getClosestCell(caloentrance,
                         true, myTrack.onEcal()==1);
      }
    else if(myTrack.onHcal())
      {
    //  std::cout << " CalorimetryManager onHcal " <<  myTrack.onHcal() << " caloentrance" << caloentrance  << std::endl;
    pivot=myCalorimeter_->getClosestCell(caloentrance,                       
                         false, false);
      }
    
    EcalHitMaker myGrid(myCalorimeter_,caloentrance,pivot,
            pivot.null()? 0 : myTrack.onEcal(),hdGridSize_,1,
            random);
    // 1=HAD shower
    
    myGrid.setTrackParameters(direction,0,myTrack);
    // Build the FAMOS HCAL 
    HcalHitMaker myHcalHitMaker(myGrid,(unsigned)1); 
    
    // Shower simulation
    bool status = false;
    int  mip = 2;
    // Use HFShower for HF
    if ( !myTrack.onEcal() && !myTrack.onHcal() ) {
      //      std::cout << "CalorimetryManager::HDShowerSimulation(): track entrance = "
      //        << myTrack.vfcalEntrance().vertex().X() << " "
      //        << myTrack.vfcalEntrance().vertex().Y() << " "
      //        << myTrack.vfcalEntrance().vertex().Z() << " "
      //        << " , Energy (Gen/Scale) = " << eGen << " " << e << std::endl;
      
      // Warning : We give here the particle energy with the response
      //           but without the resolution/gaussian smearing
      //           For HF, the resolution is due to the PE statistic
      
      if(useShowerLibrary) {
    theHFShowerLibrary->recoHFShowerLibrary(myTrack);  
    status = true;
      } else {
        HFShower theShower(random,
                   &theHDShowerparam,
               &myGrid,
               &myHcalHitMaker,
               onECAL,
               eGen);
    //           eGen);
    //           e); // PV Warning : temporarly set the energy to the generated E
    
        status = theShower.compute();
      }
    } else { 
      if(hdSimMethod_ == 0) {
    HDShower theShower(random,
               &theHDShowerparam,
               &myGrid,
               &myHcalHitMaker,
               onECAL,
               eGen,
                           pmip);
    status = theShower.compute();
        mip    = theShower.getmip();
      }
      else if (hdSimMethod_ == 1) {
    HDRShower theShower(random,
                &theHDShowerparam,
                &myGrid,
                &myHcalHitMaker,
                onECAL,
                eGen);
    status = theShower.computeShower();
        mip = 2;
      }
      else if (hdSimMethod_ == 2 ) {
    //        std::cout << "Using GflashHadronShowerProfile hdSimMethod_ == 2" << std::endl;
    
        //dynamically loading a corresponding profile by the particle type
        int particleType = myTrack.type();
        theProfile = thePiKProfile;
        if(particleType == -2212) theProfile = theAntiProtonProfile;
        else if(particleType == 2212) theProfile = theProtonProfile;
    
        //input variables for GflashHadronShowerProfile
        int showerType = 99 + myTrack.onEcal();
        double globalTime = 150.0; // a temporary reference hit time in nanosecond
        float charge = (float)(myTrack.charge());
        Gflash3Vector gfpos(trackPosition.X(),trackPosition.Y(),trackPosition.Z());
        Gflash3Vector gfmom(moment.X(),moment.Y(),moment.Z());
    
        theProfile->initialize(showerType,eGen,globalTime,charge,gfpos,gfmom);
        theProfile->loadParameters();
        theProfile->hadronicParameterization();
    
        //make hits
    std::vector<GflashHit>& gflashHitList = theProfile->getGflashHitList();
    std::vector<GflashHit>::const_iterator spotIter    = gflashHitList.begin();
    std::vector<GflashHit>::const_iterator spotIterEnd = gflashHitList.end();
    
    Gflash::CalorimeterNumber whichCalor = Gflash::kNULL;
    
        for( ; spotIter != spotIterEnd; spotIter++){
      
          double pathLength = theProfile->getGflashShowino()->getPathLengthAtShower()
            + (30*100/eGen)*(spotIter->getTime() - globalTime);
      
          double currentDepth = std::max(0.0,pathLength - theProfile->getGflashShowino()->getPathLengthOnEcal());
      
          //find the the showino position at the currentDepth
          GflashTrajectoryPoint trajectoryPoint;
          theProfile->getGflashShowino()->getHelix()->getGflashTrajectoryPoint(trajectoryPoint,pathLength);
          Gflash3Vector positionAtCurrentDepth = trajectoryPoint.getPosition();
          //find radial distrance
          Gflash3Vector lateralDisplacement = positionAtCurrentDepth - spotIter->getPosition()/CLHEP::cm;
          double rShower = lateralDisplacement.r();
          double azimuthalAngle = lateralDisplacement.phi();
      
          whichCalor = Gflash::getCalorimeterNumber(positionAtCurrentDepth);
      
          if(whichCalor==Gflash::kESPM || whichCalor==Gflash::kENCA) {
            bool statusPad = myGrid.getPads(currentDepth,true);
            if(!statusPad) continue;
            myGrid.setSpotEnergy(1.2*spotIter->getEnergy()/CLHEP::GeV);
            myGrid.addHit(rShower/Gflash::intLength[Gflash::kESPM],azimuthalAngle,0);
          }
          else if(whichCalor==Gflash::kHB || whichCalor==Gflash::kHE) {
            bool setHDdepth = myHcalHitMaker.setDepth(currentDepth,true);
            if(!setHDdepth) continue;
            myHcalHitMaker.setSpotEnergy(1.4*spotIter->getEnergy()/CLHEP::GeV);
            myHcalHitMaker.addHit(rShower/Gflash::intLength[Gflash::kHB],azimuthalAngle,0);
          }
        }
        status = true;
      }
      else {
    edm::LogInfo("FastSimulationCalorimetry") << " SimMethod " << hdSimMethod_ <<" is NOT available ";
      }
    }
    
    
    if(status) {
      
      // Here to switch between simple formulae and parameterized response
      if(optionHDSim_ == 1) {
        emeas = myHDResponse_->getHCALEnergyResponse(eGen, hit, random);
      }
      else { // optionHDsim == 2
        emeas = myHDResponse_->responseHCAL(mip, eGen, pathEta, 1, random); // 1=hadron
      }
      
      double correction = emeas / eGen;
      
      // RespCorrP factors (ECAL and HCAL separately) calculation
      respCorr(eint);     
      
      if(debug_)
    LogInfo("FastCalorimetry") 
      << "CalorimetryManager::HDShowerSimulation - on-calo 2" << std::endl
      << "   eta  = " << pathEta << std::endl
      << "   phi  = " << pathPhi << std::endl
      << "  Egen  = " << eGen << std::endl
      << " Emeas  = " << emeas << std::endl
      << "  corr  = " << correction << std::endl
      << "   mip  = " << mip << std::endl;  
      
      if(myTrack.onEcal() > 0) {
    // Save ECAL hits
    updateECAL(myGrid.getHits(),onECAL,myTrack.id(),correction*ecorr);
      }
      
      // Save HCAL hits
      if(myTrack.onVFcal() && useShowerLibrary) {
        myHDResponse_->correctHF(eGen,abs(myTrack.type()));
        updateHCAL(theHFShowerLibrary->getHitsMap(),myTrack.id());
      } 
      else 
        updateHCAL(myHcalHitMaker.getHits(),myTrack.id(),correction*hcorr);
      
    }      
    else {  // shower simulation failed  
      //      std::cout << " Shower simulation failed " << trackPosition.Vect() << std::endl;
      //      std::cout << " The FSimTrack " << myTrack << std::endl;
      //      std::cout << " HF entrance on VFcal" << myTrack.onVFcal() << std::endl;
      //      std::cout << " trackPosition.eta() " << trackPosition.eta() << std::endl;
      if(myTrack.onHcal() || myTrack.onVFcal())
    {
      DetId cell = myCalorimeter_->getClosestCell(trackPosition.Vect(),false,false);
      double tof = (((HcalGeometry*)(myCalorimeter_->getHcalGeometry()))->getPosition(cell).mag())/29.98;//speed of light
      CaloHitID current_id(cell.rawId(),tof,myTrack.id());
      std::map<CaloHitID,float> hitMap;
      hitMap[current_id] = emeas;
      updateHCAL(hitMap,myTrack.id());
      if(debug_)
        LogInfo("FastCalorimetry") << " HCAL simple cell "   
                       << cell.rawId() << " added    E = " 
                       << emeas << std::endl;  
    }
    }
    
  } // e > 0. ...
  
  if(debug_)
    LogInfo("FastCalorimetry") << std::endl << " FASTEnergyReconstructor::HDShowerSimulation  finished "
                   << std::endl;
}

void CalorimetryManager::loadFromEcalBarrel

(

edm::PCaloHitContainer &

c

)

const

Definition at line 1364 of file CalorimetryManager.cc.

References EBMapping_, plotBeamSpotDB::first, mps_fire::i, edm::second(), and EBDetId::unhashIndex().

Referenced by FamosProducer::produce().

{ 
  c.reserve(c.size()+EBMapping_.size());
  for(unsigned i=0; i<EBMapping_.size(); i++) {
    c.push_back(PCaloHit(EBDetId::unhashIndex(EBMapping_[i].first.unitID()),EBMapping_[i].second,EBMapping_[i].first.timeSlice(),EBMapping_[i].first.trackID()));
  }
}

void CalorimetryManager::loadFromEcalEndcap

(

edm::PCaloHitContainer &

c

)

const

Definition at line 1372 of file CalorimetryManager.cc.

References EEMapping_, plotBeamSpotDB::first, mps_fire::i, edm::second(), and EEDetId::unhashIndex().

Referenced by FamosProducer::produce().

{
  c.reserve(c.size()+EEMapping_.size());
  for(unsigned i=0; i<EEMapping_.size(); i++) {
    c.push_back(PCaloHit(EEDetId::unhashIndex(EEMapping_[i].first.unitID()),EEMapping_[i].second,EEMapping_[i].first.timeSlice(),EEMapping_[i].first.trackID()));
  }
}

void CalorimetryManager::loadFromHcal

(

edm::PCaloHitContainer &

c

)

const

Definition at line 1380 of file CalorimetryManager.cc.

References plotBeamSpotDB::first, HMapping_, mps_fire::i, and edm::second().

Referenced by FamosProducer::produce().

{
  c.reserve(c.size()+HMapping_.size());
  for(unsigned i=0; i<HMapping_.size(); i++) {
    c.push_back(PCaloHit(DetId(HMapping_[i].first.unitID()),HMapping_[i].second,HMapping_[i].first.timeSlice(),HMapping_[i].first.trackID()));
  }
}

void CalorimetryManager::loadFromPreshower

(

edm::PCaloHitContainer &

c

)

const

Definition at line 1389 of file CalorimetryManager.cc.

References ESMapping_, plotBeamSpotDB::first, and mps_fire::i.

Referenced by FamosProducer::produce().

{
  c.reserve(c.size()+ESMapping_.size());
  for(unsigned i=0; i<ESMapping_.size(); i++) {
    c.push_back(PCaloHit(ESMapping_[i].first.unitID(),ESMapping_[i].second,ESMapping_[i].first.timeSlice(),ESMapping_[i].first.trackID()));
  }
}

void CalorimetryManager::loadMuonSimTracks

(

edm::SimTrackContainer &

m

)

const

Definition at line 1399 of file CalorimetryManager.cc.

References funct::abs(), mps_fire::i, muonSimTracks, and findQualityFiles::size.

Referenced by FamosProducer::produce().

{
  unsigned size=muons.size();
  for(unsigned i=0; i<size;++i)
    {
      int id=muons[i].trackId();
      if(abs(muons[i].type())!=13) continue;
      // identify the corresponding muon in the local collection
      
      std::vector<FSimTrack>::const_iterator itcheck=find_if(muonSimTracks.begin(),muonSimTracks.end(),FSimTrackEqual(id));
      if(itcheck!=muonSimTracks.end())
    {
      muons[i].setTkPosition(itcheck->trackerSurfacePosition());
      muons[i].setTkMomentum(itcheck->trackerSurfaceMomentum());
      //      std::cout << " Found the SimTrack " << std::endl;
      //      std::cout << *itcheck << std::endl;
      //      std::cout << "SimTrack Id "<< id << " " << muons[i] << " " << std::endl;
    }
      //      else
      //    {
      //      std::cout << " Calorimetery Manager : this should really not happen " << std::endl;
      //      std::cout << " Was looking for " << id << " " << muons[i] << std::endl;
      //      for(unsigned i=0;i<muonSimTracks.size();++i)
      //        std::cout << muonSimTracks[i] << std::endl;
      //    }
    }
  
}

void CalorimetryManager::MuonMipSimulation ( const FSimTrack &  myTrack, RandomEngineAndDistribution const *  random  )

private

// /// TEMPORARY First attempt to include HCAL (with straight-line extrapolation):

Definition at line 862 of file CalorimetryManager.cc.

References HcalHitMaker::addHit(), ALCARECOTkAlJpsiMuMu_cff::charge, FSimTrack::charge(), debug_, EnergyLossSimulator::deltaMom(), FSimTrack::ecalEntrance(), MaterialEffects::energyLossSimulator(), objects.autophobj::float, CaloGeometryHelper::getClosestCell(), HcalHitMaker::getHits(), CaloSegment::HCAL, FSimTrack::hcalEntrance(), hdGridSize_, FSimTrack::id(), FSimTrack::momentum(), muonSimTracks, myCalorimeter_, myPart, DetId::null(), FSimTrack::onEcal(), FSimTrack::onHcal(), FSimTrack::onVFcal(), CaloSegment::PbWO4, random, reconstructHCAL(), RawParticle::setID(), HcalHitMaker::setSpotEnergy(), SimTrack::setTkMomentum(), SimTrack::setTkPosition(), EcalHitMaker::setTrackParameters(), theMuonEcalEffects, theMuonHcalEffects, updateECAL(), updateHCAL(), MaterialEffectsSimulator::updateState(), RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by reconstruct().

{
  //  TimeMe t(" FASTEnergyReconstructor::HDShower");
  XYZTLorentzVector moment = myTrack.momentum();
  
  // Backward compatibility behaviour
  if(!theMuonHcalEffects) 
    {
      if(myTrack.onHcal() || myTrack.onVFcal() ) 
    reconstructHCAL(myTrack, random);
      
      return;
    }
  
  if(debug_)
    LogInfo("FastCalorimetry") << "CalorimetryManager::MuonMipSimulation - track param."
                   << std::endl
                   << "  eta = " << moment.eta() << std::endl
                   << "  phi = " << moment.phi() << std::endl
                   << "   et = " << moment.Et()  << std::endl;
  
  //  int hit;
  //  int pid = abs(myTrack.type());
  
  XYZTLorentzVector trackPosition;
  if ( myTrack.onEcal() ) {
    trackPosition=myTrack.ecalEntrance().vertex();
    //    hit = myTrack.onEcal()-1;                               //
    myPart = myTrack.ecalEntrance();
  } else if ( myTrack.onVFcal()) {
    trackPosition=myTrack.vfcalEntrance().vertex();
    //    hit = 2;
    myPart = myTrack.vfcalEntrance();
  }
  else
    {
      LogInfo("FastCalorimetry") << " The particle is not in the acceptance " << std::endl;
      return;
    }

  // int onHCAL = hit + 1; - specially for myCalorimeter->hcalProperties(onHCAL)
  // (below) to get VFcal properties ...
  // not needed ? 
  //  int onHCAL = hit + 1; 
  int onECAL = myTrack.onEcal();
  
  //  double pathEta   = trackPosition.eta();
  //  double pathPhi   = trackPosition.phi();   
  //  double pathTheta = trackPosition.theta();
  
  //===========================================================================
  
  // ECAL and HCAL properties to get
  
  //Making ECAL Grid (and segments calculation)
  XYZPoint caloentrance;
  XYZVector direction;
  if(myTrack.onEcal()) 
    {   
      caloentrance = myTrack.ecalEntrance().vertex().Vect();
      direction = myTrack.ecalEntrance().Vect().Unit();
    }
  else if(myTrack.onHcal())
    {
      caloentrance = myTrack.hcalEntrance().vertex().Vect();
      direction = myTrack.hcalEntrance().Vect().Unit();
    }
  else
    {
      caloentrance = myTrack.vfcalEntrance().vertex().Vect();
      direction = myTrack.vfcalEntrance().Vect().Unit();
    }
  
  DetId pivot;
  if(myTrack.onEcal())
    {
      pivot=myCalorimeter_->getClosestCell(caloentrance,
                       true, myTrack.onEcal()==1);
    }
  else if(myTrack.onHcal())
    {
      //    std::cout << " CalorimetryManager onHcal " <<  myTrack.onHcal() << " caloentrance" << caloentrance  << std::endl;
      pivot=myCalorimeter_->getClosestCell(caloentrance,                         
                       false, false);
    }
  
  EcalHitMaker myGrid(myCalorimeter_,caloentrance,pivot,
              pivot.null()? 0 : myTrack.onEcal(),hdGridSize_,0,
              random);
  // 0 =EM shower -> Unit = X0
  
  myGrid.setTrackParameters(direction,0,myTrack);
  
  // Now get the path in the Preshower, ECAL and HCAL along a straight line extrapolation 
  // but only those in the ECAL are used 
  
  const std::vector<CaloSegment>& segments=myGrid.getSegments();
  unsigned nsegments=segments.size();
  
  int ifirstHcal=-1;
  int ilastEcal=-1;
  
  EnergyLossSimulator* energyLossECAL = (theMuonEcalEffects) ?
    theMuonEcalEffects->energyLossSimulator() : nullptr;
  //  // Muon brem in ECAL
  //  MuonBremsstrahlungSimulator* muonBremECAL = 0;
  //  if (theMuonEcalEffects) muonBremECAL = theMuonEcalEffects->muonBremsstrahlungSimulator();
  
  for(unsigned iseg=0;iseg<nsegments&&ifirstHcal<0;++iseg)
    {
      
      // in the ECAL, there are two types of segments: PbWO4 and GAP
      float segmentSizeinX0=segments[iseg].X0length();
      
      // Martijn - insert your computations here
      float energy=0.0;
      if (segmentSizeinX0>0.001 && segments[iseg].material()==CaloSegment::PbWO4 ) {
    // The energy loss simulator
    float charge = (float)(myTrack.charge());
    ParticlePropagator theMuon(moment,trackPosition,charge,nullptr);
    theMuon.setID(-(int)charge*13);
    if ( energyLossECAL ) { 
      energyLossECAL->updateState(theMuon, segmentSizeinX0, random);
      energy = energyLossECAL->deltaMom().E();
      moment -= energyLossECAL->deltaMom();
    }
      } 
      // that's all for ECAL, Florian
      // Save the hit only if it is a crystal
      if(segments[iseg].material()==CaloSegment::PbWO4)
    {
      myGrid.getPads(segments[iseg].sX0Entrance()+segmentSizeinX0*0.5);
      myGrid.setSpotEnergy(energy);
      myGrid.addHit(0.,0.);
      ilastEcal=iseg;
    }
      // Check for end of loop:
      if(segments[iseg].material()==CaloSegment::HCAL)
    {
      ifirstHcal=iseg;
    }
    }
  
  
  // Build the FAMOS HCAL 
  HcalHitMaker myHcalHitMaker(myGrid,(unsigned)2);     
  // float mipenergy=0.1;
  // Create the helix with the stepping helix propagator
  // to add a hit, just do
  // myHcalHitMaker.setSpotEnergy(mipenergy);
  // math::XYZVector hcalEntrance;
  // if(ifirstHcal>=0) hcalEntrance=segments[ifirstHcal].entrance();
  // myHcalHitMaker.addHit(hcalEntrance);
  int ilastHcal=-1;
  float mipenergy=0.0;
  
  EnergyLossSimulator* energyLossHCAL =  (theMuonHcalEffects) ?
    theMuonHcalEffects->energyLossSimulator() : nullptr;
  //  // Muon Brem effect
  //  MuonBremsstrahlungSimulator* muonBremHCAL = 0;
  //  if (theMuonHcalEffects) muonBremHCAL = theMuonHcalEffects->muonBremsstrahlungSimulator(); 
  
  if(ifirstHcal>0 && energyLossHCAL){
    for(unsigned iseg=ifirstHcal;iseg<nsegments;++iseg)
      {
    float segmentSizeinX0=segments[iseg].X0length();
    if(segments[iseg].material()==CaloSegment::HCAL) {
      ilastHcal=iseg;
      if (segmentSizeinX0>0.001) {
        // The energy loss simulator
        float charge = (float)(myTrack.charge());
        ParticlePropagator theMuon(moment,trackPosition,charge,nullptr);
        theMuon.setID(-(int)charge*13);
        energyLossHCAL->updateState(theMuon, segmentSizeinX0, random);
        mipenergy = energyLossHCAL->deltaMom().E();
        moment -= energyLossHCAL->deltaMom();
        myHcalHitMaker.setSpotEnergy(mipenergy);
        myHcalHitMaker.addHit(segments[iseg].entrance());
      }
    } 
      }
  }
  
  //
  
  
  
  // Copy the muon SimTrack (Only for Energy loss)
  FSimTrack muonTrack(myTrack);
  if(energyLossHCAL && ilastHcal>=0) {
    math::XYZVector hcalExit=segments[ilastHcal].exit();
    muonTrack.setTkPosition(hcalExit);
    muonTrack.setTkMomentum(moment);
  } else if(energyLossECAL && ilastEcal>=0) {
    math::XYZVector ecalExit=segments[ilastEcal].exit();
    muonTrack.setTkPosition(ecalExit);
    muonTrack.setTkMomentum(moment);
  } // else just leave tracker surface position and momentum...  
  
  muonSimTracks.push_back(muonTrack);
  
  
  // no need to change below this line
  std::map<CaloHitID,float>::const_iterator mapitr;
  std::map<CaloHitID,float>::const_iterator endmapitr;
  if(myTrack.onEcal() > 0) {
    // Save ECAL hits
    updateECAL(myGrid.getHits(),onECAL,myTrack.id());
  }
  
  // Save HCAL hits
  updateHCAL(myHcalHitMaker.getHits(),myTrack.id());
  
  if(debug_)
    LogInfo("FastCalorimetry") << std::endl << " FASTEnergyReconstructor::MipShowerSimulation  finished "
                   << std::endl;
}

void CalorimetryManager::readParameters ( const edm::ParameterSet &  fastCalo )

private

Definition at line 1088 of file CalorimetryManager.cc.

References aTerm, bFixedLength_, bTerm, crackPadSurvivalProbability_, debug_, EcalDigitizer_, evtsToDebug_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gridSize_, HcalDigitizer_, hdGridSize_, hdSimMethod_, mps_fire::i, ietaShiftHB_, ietaShiftHE_, ietaShiftHF_, ietaShiftHO_, k_e, k_h, mipValues_, optionHDSim_, p_knots, pulledPadSurvivalProbability_, radiusFactorEB_, radiusFactorEE_, radiusPreshowerCorrections_, RCFactor_, rsp, RTFactor_, samplingHBHE_, samplingHF_, samplingHO_, simulatePreshower_, spotFraction_, theCoreIntervals_, theTailIntervals_, timeShiftHB_, timeShiftHE_, timeShiftHF_, timeShiftHO_, useCorrectionSL, and useShowerLibrary.

Referenced by CalorimetryManager().

                                                                       {
  
  edm::ParameterSet ECALparameters = fastCalo.getParameter<edm::ParameterSet>("ECAL");
  
  evtsToDebug_ = fastCalo.getUntrackedParameter<std::vector<unsigned int> >("EvtsToDebug",std::vector<unsigned>());
  debug_ = fastCalo.getUntrackedParameter<bool>("Debug");
  
  bFixedLength_ = ECALparameters.getParameter<bool>("bFixedLength");
  //   std::cout << "bFixedLength_ = " << bFixedLength_ << std::endl;
  
  gridSize_ = ECALparameters.getParameter<int>("GridSize");
  spotFraction_ = ECALparameters.getParameter<double>("SpotFraction");
  pulledPadSurvivalProbability_ = ECALparameters.getParameter<double>("FrontLeakageProbability");
  crackPadSurvivalProbability_ = ECALparameters.getParameter<double>("GapLossProbability");
  theCoreIntervals_ = ECALparameters.getParameter<std::vector<double> >("CoreIntervals");
  theTailIntervals_ = ECALparameters.getParameter<std::vector<double> >("TailIntervals");
  
  RCFactor_ = ECALparameters.getParameter<double>("RCFactor");
  RTFactor_ = ECALparameters.getParameter<double>("RTFactor");
  //changed after tuning - Feb-July - Shilpi Jain
  //  radiusFactor_ = ECALparameters.getParameter<double>("RadiusFactor");
  radiusFactorEE_ = ECALparameters.getParameter<double>("RadiusFactorEE");
  radiusFactorEB_ = ECALparameters.getParameter<double>("RadiusFactorEB");
  //(end of) changed after tuning - Feb-July - Shilpi Jain
  radiusPreshowerCorrections_ = ECALparameters.getParameter<std::vector<double> >("RadiusPreshowerCorrections");
  aTerm = 1.+radiusPreshowerCorrections_[1]*radiusPreshowerCorrections_[0];
  bTerm = radiusPreshowerCorrections_[0];
  mipValues_ = ECALparameters.getParameter<std::vector<double> >("MipsinGeV");
  simulatePreshower_ = ECALparameters.getParameter<bool>("SimulatePreshower");
  
  if(gridSize_ <1) gridSize_= 7;
  if(pulledPadSurvivalProbability_ <0. || pulledPadSurvivalProbability_>1 ) pulledPadSurvivalProbability_= 1.;
  if(crackPadSurvivalProbability_ <0. || crackPadSurvivalProbability_>1 ) crackPadSurvivalProbability_= 0.9;
  
  LogInfo("FastCalorimetry") << " Fast ECAL simulation parameters " << std::endl;
  LogInfo("FastCalorimetry") << " =============================== " << std::endl;
  if(simulatePreshower_)
    LogInfo("FastCalorimetry") << " The preshower is present " << std::endl;
  else
    LogInfo("FastCalorimetry") << " The preshower is NOT present " << std::endl;
  LogInfo("FastCalorimetry") << " Grid Size : " << gridSize_  << std::endl; 
  if(spotFraction_>0.) 
    LogInfo("FastCalorimetry") << " Spot Fraction : " << spotFraction_ << std::endl;
  else
    {
      LogInfo("FastCalorimetry") << " Core of the shower " << std::endl;
      for(unsigned ir=0; ir < theCoreIntervals_.size()/2;++ir)
    {
      LogInfo("FastCalorimetry") << " r < " << theCoreIntervals_[ir*2] << " R_M : " << theCoreIntervals_[ir*2+1] << "        ";
    }
      LogInfo("FastCalorimetry") << std::endl;
      
      LogInfo("FastCalorimetry") << " Tail of the shower " << std::endl;
      for(unsigned ir=0; ir < theTailIntervals_.size()/2;++ir)
    {
      LogInfo("FastCalorimetry") << " r < " << theTailIntervals_[ir*2] << " R_M : " << theTailIntervals_[ir*2+1] << "        ";
    }
      //changed after tuning - Feb-July - Shilpi Jain
      //      LogInfo("FastCalorimetry") << "Radius correction factor " << radiusFactor_ << std::endl;
      LogInfo("FastCalorimetry") << "Radius correction factors:  EB & EE " << radiusFactorEB_ << " : "<< radiusFactorEE_ << std::endl;
      //(end of) changed after tuning - Feb-July - Shilpi Jain
      LogInfo("FastCalorimetry") << std::endl;
      if(mipValues_.size()>2)   {
    LogInfo("FastCalorimetry") << "Improper number of parameters for the preshower ; using 95keV" << std::endl;
    mipValues_.clear();
    mipValues_.resize(2,0.000095);
      }
    }
  
  LogInfo("FastCalorimetry") << " FrontLeakageProbability : " << pulledPadSurvivalProbability_ << std::endl;
  LogInfo("FastCalorimetry") << " GapLossProbability : " << crackPadSurvivalProbability_ << std::endl;
  
  
  // RespCorrP: p (momentum), ECAL and HCAL corrections = f(p)
  edm::ParameterSet CalorimeterParam = fastCalo.getParameter<edm::ParameterSet>("CalorimeterProperties");
  
  rsp = CalorimeterParam.getParameter<std::vector<double> >("RespCorrP");
  LogInfo("FastCalorimetry") << " RespCorrP (rsp) size " << rsp.size() << std::endl;
  
  if( rsp.size()%3 !=0 )  {
    LogInfo("FastCalorimetry") 
      << " RespCorrP size is wrong -> no corrections applied !!!" 
      << std::endl;
    
    p_knots.push_back(14000.);
    k_e.push_back    (1.);
    k_h.push_back    (1.);
  }
  else {
    for(unsigned i = 0; i < rsp.size(); i += 3) { 
      LogInfo("FastCalorimetry") << "i = " << i/3 << "   p = " << rsp [i] 
                 << "   k_e(p) = " << rsp[i+1] 
                 << "   k_e(p) = " << rsp[i+2] << std::endl; 
      
      p_knots.push_back(rsp[i]);
      k_e.push_back    (rsp[i+1]);
      k_h.push_back    (rsp[i+2]); 
    }
  }  
  
  
  //FR
  edm::ParameterSet HCALparameters = fastCalo.getParameter<edm::ParameterSet>("HCAL");
  optionHDSim_ = HCALparameters.getParameter<int>("SimOption");
  hdGridSize_  = HCALparameters.getParameter<int>("GridSize");
  hdSimMethod_ = HCALparameters.getParameter<int>("SimMethod");
  //RF
  
  EcalDigitizer_    = ECALparameters.getUntrackedParameter<bool>("Digitizer",false);
  HcalDigitizer_    = HCALparameters.getUntrackedParameter<bool>("Digitizer",false);
  samplingHBHE_ = HCALparameters.getParameter< std::vector<double> >("samplingHBHE");
  samplingHF_   = HCALparameters.getParameter< std::vector<double> >("samplingHF");
  samplingHO_   = HCALparameters.getParameter< std::vector<double> >("samplingHO");
  ietaShiftHB_  = HCALparameters.getParameter< int >("ietaShiftHB");
  ietaShiftHE_  = HCALparameters.getParameter< int >("ietaShiftHE");
  ietaShiftHF_  = HCALparameters.getParameter< int >("ietaShiftHF");
  ietaShiftHO_  = HCALparameters.getParameter< int >("ietaShiftHO");
  timeShiftHB_  = HCALparameters.getParameter< std::vector<double> >("timeShiftHB");
  timeShiftHE_  = HCALparameters.getParameter< std::vector<double> >("timeShiftHE");
  timeShiftHF_  = HCALparameters.getParameter< std::vector<double> >("timeShiftHF");
  timeShiftHO_  = HCALparameters.getParameter< std::vector<double> >("timeShiftHO");
  
  // FastHFShowerLibrary
  edm::ParameterSet m_HS = fastCalo.getParameter<edm::ParameterSet>("HFShowerLibrary");
  useShowerLibrary       = m_HS.getUntrackedParameter<bool>("useShowerLibrary",false);
  useCorrectionSL        = m_HS.getUntrackedParameter<bool>("useCorrectionSL",false);
}

void CalorimetryManager::reconstruct

(

RandomEngineAndDistribution const *

random

)

Definition at line 173 of file CalorimetryManager.cc.

References funct::abs(), clean(), HCALResponse::correctHF(), gather_cfg::cout, debug_, EMShowerSimulation(), edm::EventID::event(), evtsToDebug_, spr::find(), FastHFShowerLibrary::getHitsMap(), FSimTrack::hcalEntrance(), HDShowerSimulation(), FSimEvent::id(), FSimTrack::id(), initialized_, createfilelist::int, MuonMipSimulation(), myCalorimeter_, myHDResponse_, mySimEvent, FSimTrack::noEndVertex(), FSimEvent::nTracks(), FSimTrack::onEcal(), FSimTrack::onHcal(), FSimTrack::onVFcal(), optionHDSim_, sysUtil::pid, CaloGeometryHelper::preshowerPresent(), FBaseSimEvent::print(), FastHFShowerLibrary::recoHFShowerLibrary(), reconstructHCAL(), FastHFShowerLibrary::SetRandom(), simulatePreshower_, theHFShowerLibrary, FBaseSimEvent::track(), CoreSimTrack::type(), updateHCAL(), and useShowerLibrary.

Referenced by FamosManager::reconstruct().

{
  theHFShowerLibrary->SetRandom(random);
  if(!evtsToDebug_.empty())
    {
      std::vector<unsigned int>::const_iterator itcheck=find(evtsToDebug_.begin(),evtsToDebug_.end(),mySimEvent->id().event());
      debug_=(itcheck!=evtsToDebug_.end());
      if(debug_)
    mySimEvent->print();
    }
  // Clear the content of the calorimeters 
  if(!initialized_)
    {
      
      // Check if the preshower is really available
      if(simulatePreshower_ && !myCalorimeter_->preshowerPresent())
    {
      std::cout << " WARNING " << std::endl;
      std::cout << " The preshower simulation has been turned on; but no preshower geometry is available " << std::endl;
      std::cout << " Disabling the preshower simulation " << std::endl;
      simulatePreshower_ = false;
    }
      
      initialized_=true;
    }
  clean();

  
  
  LogInfo("FastCalorimetry") << "Reconstructing " << (int) mySimEvent->nTracks() << " tracks." << std::endl;
  for( int fsimi=0; fsimi < (int) mySimEvent->nTracks() ; ++fsimi) {
    
    FSimTrack& myTrack = mySimEvent->track(fsimi);
    
    int pid = abs(myTrack.type());
    
    if (debug_) {
      LogInfo("FastCalorimetry") << " ===> pid = "  << pid << std::endl;      
    }
    
    
    // Check that the particle hasn't decayed
    if(myTrack.noEndVertex()) {
      // Simulate energy smearing for photon and electrons
      if ( pid == 11 || pid == 22 ) {
    
    if ( myTrack.onEcal() ) 
      EMShowerSimulation(myTrack, random);
    else if ( myTrack.onVFcal() ) {
          if(useShowerLibrary) {
            theHFShowerLibrary->recoHFShowerLibrary(myTrack);  
            myHDResponse_->correctHF(myTrack.hcalEntrance().e(),abs(myTrack.type()));
            updateHCAL(theHFShowerLibrary->getHitsMap(),myTrack.id());
          } 
          else reconstructHCAL(myTrack, random);
    }   
      } // electron or photon
      else if (pid==13)
    {
          MuonMipSimulation(myTrack, random);
    }
      // Simulate energy smearing for hadrons (i.e., everything 
      // but muons... and SUSY particles that deserve a special 
      // treatment.
      else if ( pid < 1000000 ) {
    if ( myTrack.onHcal() || myTrack.onVFcal() ) {    
      if(optionHDSim_ == 0 )  reconstructHCAL(myTrack, random);
      else HDShowerSimulation(myTrack, random);
    }
      } // pid < 1000000 
    } // myTrack.noEndVertex()
  } // particle loop
  //  LogInfo("FastCalorimetry") << " Number of  hits (barrel)" << EBMapping_.size() << std::endl;
  //  LogInfo("FastCalorimetry") << " Number of  hits (Hcal)" << HMapping_.size() << std::endl;
  //  std::cout << " Nombre de hit (endcap)" << EEMapping_.size() << std::endl;
  
} // reconstruct

void CalorimetryManager::reconstructHCAL ( const FSimTrack &  myTrack, RandomEngineAndDistribution const *  random  )

private

Definition at line 493 of file CalorimetryManager.cc.

References funct::abs(), debug_, CaloGeometryHelper::getClosestCell(), HCALResponse::getHCALEnergyResponse(), Calorimeter::getHcalGeometry(), FSimTrack::hcalEntrance(), FSimTrack::id(), myCalorimeter_, myHDResponse_, FSimTrack::onHcal(), sysUtil::pid, DetId::rawId(), HCALResponse::responseHCAL(), CoreSimTrack::type(), updateHCAL(), RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by MuonMipSimulation(), and reconstruct().

{
  int hit;
  int pid = abs(myTrack.type());
  if (debug_) {
    LogInfo("FastCalorimetry") << " reconstructHCAL "  << myTrack << std::endl;      
  }
  
  //  FSimTrack myTrack = mySimEvent.track(fsimi);
  
  //  int pid=abs(myTrack.type());
  //  std::cout << "reconstructHCAL " << std::endl;
  
  XYZTLorentzVector trackPosition;
  if (myTrack.onHcal()) {
    trackPosition=myTrack.hcalEntrance().vertex();
    hit = myTrack.onHcal()-1;
  } else {
    trackPosition=myTrack.vfcalEntrance().vertex();
    hit = 2;
  }
  
  double pathEta   = trackPosition.eta();
  double pathPhi   = trackPosition.phi();   
  //  double pathTheta = trackPosition.theta();
  
  double EGen  = myTrack.hcalEntrance().e();
  double emeas = 0.;
  //double emeas = -0.0001;
  
  if(pid == 13) { 
    //    std::cout << " We should not be here " << std::endl;
    emeas = myHDResponse_->responseHCAL(0, EGen, pathEta, 2, random); // 2=muon
    if(debug_)
      LogInfo("FastCalorimetry") << "CalorimetryManager::reconstructHCAL - MUON !!!" << std::endl;
  }
  else if( pid == 22 || pid == 11) {
    emeas = myHDResponse_->responseHCAL(0, EGen, pathEta, 0, random); // last par. = 0 = e/gamma
    //  cout <<  "CalorimetryManager::reconstructHCAL - e/gamma !!!" << std::endl;
    if(debug_)
      LogInfo("FastCalorimetry") << "CalorimetryManager::reconstructHCAL - e/gamma !!!" << std::endl;
  }
  else {
    emeas = myHDResponse_->getHCALEnergyResponse(EGen, hit, random);
  }
  
  if(debug_)
    LogInfo("FastCalorimetry") << "CalorimetryManager::reconstructHCAL - on-calo "   
                   << "  eta = " << pathEta 
                   << "  phi = " << pathPhi 
                   << "  Egen = " << EGen 
                   << "  Emeas = " << emeas << std::endl;
  
  if(emeas > 0.) {  
    DetId cell = myCalorimeter_->getClosestCell(trackPosition.Vect(),false,false);
    double tof = (((HcalGeometry*)(myCalorimeter_->getHcalGeometry()))->getPosition(cell).mag())/29.98;//speed of light
    CaloHitID current_id(cell.rawId(),tof,myTrack.id());
    std::map<CaloHitID,float> hitMap;
    hitMap[current_id] = emeas;
    updateHCAL(hitMap,myTrack.id());
  }
}

void CalorimetryManager::respCorr ( double  p )

private

Definition at line 1216 of file CalorimetryManager.cc.

References debug_, ecorr, hcorr, mps_fire::i, k_e, k_h, p_knots, globals_cff::x1, and globals_cff::x2.

Referenced by HDShowerSimulation().

                                          {
  
  int sizeP = p_knots.size();
  
  if(sizeP <= 1) {
    ecorr = 1.;
    hcorr = 1.;
  }
  else {
    int ip = -1;    
    for (int i = 0; i < sizeP; i++) { 
      if (p < p_knots[i]) { ip = i; break;}
    }
    if (ip == 0) {
      ecorr = k_e[0];
      hcorr = k_h[0];
    }
    else {
      if(ip == -1) {
    ecorr = k_e[sizeP-1];
    hcorr = k_h[sizeP-1];
      } 
      else {
    double x1 =  p_knots[ip-1];
    double x2 =  p_knots[ip];
    double y1 =  k_e[ip-1];
    double y2 =  k_e[ip];
    
    if(x1 == x2) {
      //        std::cout << " equal p_knots values!!! " << std::endl;
    }   
    
    ecorr = (y1 + (y2 - y1) * (p - x1)/(x2 - x1));
    
    y1 =  k_h[ip-1];
    y2 =  k_h[ip];
    hcorr = (y1 + (y2 - y1) * (p - x1)/(x2 - x1)); 
    
      }
    }
  }
  
  if(debug_)
    LogInfo("FastCalorimetry") << " p, ecorr, hcorr = " << p << " "  
                   << ecorr << "  " << hcorr << std::endl;
  
}

void CalorimetryManager::updateECAL ( const std::map< CaloHitID, float > &  hitMap, int  onEcal, int  trackID = 0, float  corr = 1.0  )

private

Definition at line 1264 of file CalorimetryManager.cc.

References corr, EBMapping_, and EEMapping_.

Referenced by EMShowerSimulation(), HDShowerSimulation(), and MuonMipSimulation().

{
  std::map<CaloHitID,float>::const_iterator mapitr;
  std::map<CaloHitID,float>::const_iterator endmapitr=hitMap.end();
  if(onEcal==1) {
    EBMapping_.reserve(EBMapping_.size()+hitMap.size());
    endmapitr=hitMap.end();
    for(mapitr=hitMap.begin();mapitr!=endmapitr;++mapitr) {
      //correct energy
      float energy = mapitr->second;
      energy *= corr;
      
      //make finalized CaloHitID
      CaloHitID current_id(mapitr->first.unitID(),mapitr->first.timeSlice(),trackID);
      
      EBMapping_.push_back(std::pair<CaloHitID,float>(current_id,energy));
    }
  }
  else if(onEcal==2) {
    EEMapping_.reserve(EEMapping_.size()+hitMap.size());
    endmapitr=hitMap.end();
    for(mapitr=hitMap.begin();mapitr!=endmapitr;++mapitr) {
      //correct energy
      float energy = mapitr->second;
      energy *= corr;
      
      //make finalized CaloHitID
      CaloHitID current_id(mapitr->first.unitID(),mapitr->first.timeSlice(),trackID);
      
      EEMapping_.push_back(std::pair<CaloHitID,float>(current_id,energy));
    }
  }
  
}

void CalorimetryManager::updateHCAL ( const std::map< CaloHitID, float > &  hitMap, int  trackID = 0, float  corr = 1.0  )

private

Definition at line 1299 of file CalorimetryManager.cc.

References corr, HcalDetId::depth(), HCALResponse::getCorrHFem(), HCALResponse::getCorrHFhad(), HcalBarrel, HcalDigitizer_, HcalEndcap, HcalForward, HcalOuter, HMapping_, HcalDetId::ietaAbs(), ietaShiftHB_, ietaShiftHE_, ietaShiftHF_, ietaShiftHO_, myHDResponse_, samplingHBHE_, samplingHF_, samplingHO_, DetId::subdetId(), ntuplemaker::time, timeShiftHB_, timeShiftHE_, timeShiftHF_, timeShiftHO_, useCorrectionSL, and useShowerLibrary.

Referenced by EMShowerSimulation(), HDShowerSimulation(), MuonMipSimulation(), reconstruct(), and reconstructHCAL().

{
  std::vector<double> hfcorrEm = myHDResponse_->getCorrHFem();
  std::vector<double> hfcorrHad = myHDResponse_->getCorrHFhad();
  std::map<CaloHitID,float>::const_iterator mapitr;
  std::map<CaloHitID,float>::const_iterator endmapitr=hitMap.end();
  HMapping_.reserve(HMapping_.size()+hitMap.size());
  for(mapitr=hitMap.begin(); mapitr!=endmapitr; ++mapitr) {
    //correct energy
    float energy = mapitr->second;
    energy *= corr;
    
    float time = mapitr->first.timeSlice();
    //put energy into uncalibrated state for digitizer && correct timing
    if(HcalDigitizer_){
      HcalDetId hdetid = HcalDetId(mapitr->first.unitID());
      if (hdetid.subdetId()== HcalBarrel){
    energy /= samplingHBHE_[hdetid.ietaAbs()-1]; //re-convert to GeV
    time = timeShiftHB_[hdetid.ietaAbs()-ietaShiftHB_];
      }
      else if (hdetid.subdetId()== HcalEndcap){
    energy /= samplingHBHE_[hdetid.ietaAbs()-1]; //re-convert to GeV
    time = timeShiftHE_[hdetid.ietaAbs()-ietaShiftHE_];
      }
      else if (hdetid.subdetId()== HcalForward){
        if(useShowerLibrary) {
      if(useCorrectionSL) {
             if(hdetid.depth()== 1) energy *= hfcorrEm[hdetid.ietaAbs()-ietaShiftHF_];
             if(hdetid.depth()== 2) energy *= hfcorrHad[hdetid.ietaAbs()-ietaShiftHF_];
          }
    } else {
      if(hdetid.depth()== 1) energy *= samplingHF_[0];
      if(hdetid.depth()== 2) energy *= samplingHF_[1];
          time = timeShiftHF_[hdetid.ietaAbs()-ietaShiftHF_];
    } 
      }
      else if (hdetid.subdetId()== HcalOuter){
    energy /= samplingHO_[hdetid.ietaAbs()-1];
    time = timeShiftHO_[hdetid.ietaAbs()-ietaShiftHO_];
      }
    }   
    
    //make finalized CaloHitID
    CaloHitID current_id(mapitr->first.unitID(),time,trackID);
    HMapping_.push_back(std::pair<CaloHitID,float>(current_id,energy));
  }
}

void CalorimetryManager::updatePreshower ( const std::map< CaloHitID, float > &  hitMap, int  trackID = 0, float  corr = 1.0  )

private

Definition at line 1347 of file CalorimetryManager.cc.

References corr, and ESMapping_.

Referenced by EMShowerSimulation().

{
  std::map<CaloHitID,float>::const_iterator mapitr;
  std::map<CaloHitID,float>::const_iterator endmapitr=hitMap.end();
  ESMapping_.reserve(ESMapping_.size()+hitMap.size());
  for(mapitr=hitMap.begin();mapitr!=endmapitr;++mapitr) {
    //correct energy
    float energy = mapitr->second;
    energy *= corr;
    
    //make finalized CaloHitID
    CaloHitID current_id(mapitr->first.unitID(),mapitr->first.timeSlice(),trackID);
    
    ESMapping_.push_back(std::pair<CaloHitID,float>(current_id,energy));
  }
}

Member Data Documentation

GammaFunctionGenerator* CalorimetryManager::aGammaGenerator

private

Definition at line 154 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), and EMShowerSimulation().

const LandauFluctuationGenerator* CalorimetryManager::aLandauGenerator

private

Definition at line 153 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), and EMShowerSimulation().

double CalorimetryManager::aTerm

private

Definition at line 143 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

bool CalorimetryManager::bFixedLength_

private

Definition at line 179 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

double CalorimetryManager::bTerm

private

Definition at line 143 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

double CalorimetryManager::crackPadSurvivalProbability_

private

Definition at line 138 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

bool CalorimetryManager::debug_

private

Definition at line 114 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), HDShowerSimulation(), MuonMipSimulation(), readParameters(), reconstruct(), reconstructHCAL(), and respCorr().

std::vector<std::pair<CaloHitID,float> > CalorimetryManager::EBMapping_

private

Definition at line 109 of file CalorimetryManager.h.

Referenced by clean(), loadFromEcalBarrel(), and updateECAL().

std::unique_ptr<KKCorrectionFactors> CalorimetryManager::ecalCorrection

private

Definition at line 192 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), and EMShowerSimulation().

bool CalorimetryManager::EcalDigitizer_

private

Definition at line 120 of file CalorimetryManager.h.

Referenced by readParameters().

double CalorimetryManager::ecorr

private

Definition at line 164 of file CalorimetryManager.h.

Referenced by HDShowerSimulation(), and respCorr().

std::vector<std::pair<CaloHitID,float> > CalorimetryManager::EEMapping_

private

Definition at line 110 of file CalorimetryManager.h.

Referenced by clean(), loadFromEcalEndcap(), and updateECAL().

std::vector<std::pair<CaloHitID,float> > CalorimetryManager::ESMapping_

private

Definition at line 112 of file CalorimetryManager.h.

Referenced by clean(), loadFromPreshower(), and updatePreshower().

std::vector<unsigned int> CalorimetryManager::evtsToDebug_

private

Definition at line 115 of file CalorimetryManager.h.

Referenced by readParameters(), and reconstruct().

int CalorimetryManager::gridSize_

private

Definition at line 145 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

bool CalorimetryManager::HcalDigitizer_

private

Definition at line 121 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

double CalorimetryManager::hcorr

private

Definition at line 165 of file CalorimetryManager.h.

Referenced by HDShowerSimulation(), and respCorr().

int CalorimetryManager::hdGridSize_

private

Definition at line 149 of file CalorimetryManager.h.

Referenced by HDShowerSimulation(), MuonMipSimulation(), and readParameters().

int CalorimetryManager::hdSimMethod_

private

Definition at line 149 of file CalorimetryManager.h.

Referenced by HDShowerSimulation(), and readParameters().

std::vector<std::pair<CaloHitID,float> > CalorimetryManager::HMapping_

private

Definition at line 111 of file CalorimetryManager.h.

Referenced by clean(), loadFromHcal(), and updateHCAL().

int CalorimetryManager::ietaShiftHB_

private

Definition at line 125 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

int CalorimetryManager::ietaShiftHE_

private

Definition at line 125 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

int CalorimetryManager::ietaShiftHF_

private

Definition at line 125 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

int CalorimetryManager::ietaShiftHO_

private

Definition at line 125 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

bool CalorimetryManager::initialized_

private

Definition at line 168 of file CalorimetryManager.h.

Referenced by reconstruct().

std::vector<double> CalorimetryManager::k_e

private

Definition at line 162 of file CalorimetryManager.h.

Referenced by readParameters(), and respCorr().

std::vector<double> CalorimetryManager::k_h

private

Definition at line 163 of file CalorimetryManager.h.

Referenced by readParameters(), and respCorr().

std::vector<double> CalorimetryManager::mipValues_

private

Definition at line 144 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

std::vector<FSimTrack> CalorimetryManager::muonSimTracks

private

Definition at line 170 of file CalorimetryManager.h.

Referenced by clean(), loadMuonSimTracks(), and MuonMipSimulation().

CaloGeometryHelper* CalorimetryManager::myCalorimeter_

private

Definition at line 102 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), EMShowerSimulation(), HDShowerSimulation(), MuonMipSimulation(), reconstruct(), reconstructHCAL(), and ~CalorimetryManager().

RawParticle CalorimetryManager::myElec

private

A few pointers to save time.

Definition at line 132 of file CalorimetryManager.h.

Referenced by EMShowerSimulation().

HCALResponse* CalorimetryManager::myHDResponse_

private

Definition at line 106 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), HDShowerSimulation(), reconstruct(), reconstructHCAL(), updateHCAL(), and ~CalorimetryManager().

Histos* CalorimetryManager::myHistos

private

Definition at line 104 of file CalorimetryManager.h.

HSParameters* CalorimetryManager::myHSParameters_

private

Definition at line 107 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), and HDShowerSimulation().

RawParticle CalorimetryManager::myPart

private

Definition at line 134 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), HDShowerSimulation(), and MuonMipSimulation().

RawParticle CalorimetryManager::myPosi

private

Definition at line 133 of file CalorimetryManager.h.

Referenced by EMShowerSimulation().

FSimEvent* CalorimetryManager::mySimEvent

private

Definition at line 101 of file CalorimetryManager.h.

Referenced by reconstruct().

std::vector< std::pair< int, float > > CalorimetryManager::myZero_

staticprivate

Initial value:

= std::vector<std::pair<int,float> >
  (1,std::pair<int,float>(0,0.))

Definition at line 156 of file CalorimetryManager.h.

int CalorimetryManager::optionHDSim_

private

Definition at line 149 of file CalorimetryManager.h.

Referenced by HDShowerSimulation(), readParameters(), and reconstruct().

std::vector<double> CalorimetryManager::p_knots

private

Definition at line 161 of file CalorimetryManager.h.

Referenced by readParameters(), and respCorr().

double CalorimetryManager::pulledPadSurvivalProbability_

private

Definition at line 137 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

double CalorimetryManager::radiusFactorEB_

private

Definition at line 141 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

double CalorimetryManager::radiusFactorEE_

private

Definition at line 141 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

std::vector<double> CalorimetryManager::radiusPreshowerCorrections_

private

Definition at line 142 of file CalorimetryManager.h.

Referenced by readParameters().

double CalorimetryManager::RCFactor_

private

Definition at line 147 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

std::vector<double> CalorimetryManager::rsp

private

Definition at line 160 of file CalorimetryManager.h.

Referenced by readParameters().

double CalorimetryManager::RTFactor_

private

Definition at line 147 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

std::vector<double> CalorimetryManager::samplingHBHE_

private

Definition at line 122 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

std::vector<double> CalorimetryManager::samplingHF_

private

Definition at line 123 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

std::vector<double> CalorimetryManager::samplingHO_

private

Definition at line 124 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

bool CalorimetryManager::simulatePreshower_

private

Definition at line 150 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), readParameters(), and reconstruct().

double CalorimetryManager::spotFraction_

private

Definition at line 139 of file CalorimetryManager.h.

Referenced by readParameters().

GflashAntiProtonShowerProfile* CalorimetryManager::theAntiProtonProfile

private

Definition at line 185 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), and HDShowerSimulation().

std::vector<double> CalorimetryManager::theCoreIntervals_

private

Definition at line 146 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

FastHFShowerLibrary* CalorimetryManager::theHFShowerLibrary

private

Definition at line 190 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), HDShowerSimulation(), reconstruct(), and ~CalorimetryManager().

MaterialEffects* CalorimetryManager::theMuonEcalEffects

private

Definition at line 171 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), MuonMipSimulation(), and ~CalorimetryManager().

MaterialEffects* CalorimetryManager::theMuonHcalEffects

private

Definition at line 172 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), MuonMipSimulation(), and ~CalorimetryManager().

GflashPiKShowerProfile* CalorimetryManager::thePiKProfile

private

Definition at line 183 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), and HDShowerSimulation().

GflashHadronShowerProfile* CalorimetryManager::theProfile

private

Definition at line 182 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), HDShowerSimulation(), and ~CalorimetryManager().

GflashProtonShowerProfile* CalorimetryManager::theProtonProfile

private

Definition at line 184 of file CalorimetryManager.h.

Referenced by CalorimetryManager(), and HDShowerSimulation().

std::vector<double> CalorimetryManager::theTailIntervals_

private

Definition at line 146 of file CalorimetryManager.h.

Referenced by EMShowerSimulation(), and readParameters().

std::vector<double> CalorimetryManager::timeShiftHB_

private

Definition at line 126 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

std::vector<double> CalorimetryManager::timeShiftHE_

private

Definition at line 127 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

std::vector<double> CalorimetryManager::timeShiftHF_

private

Definition at line 128 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

std::vector<double> CalorimetryManager::timeShiftHO_

private

Definition at line 129 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

bool CalorimetryManager::unfoldedMode_

private

Definition at line 117 of file CalorimetryManager.h.

bool CalorimetryManager::useCorrectionSL

private

Definition at line 189 of file CalorimetryManager.h.

Referenced by readParameters(), and updateHCAL().

bool CalorimetryManager::useShowerLibrary

private

Definition at line 188 of file CalorimetryManager.h.

Referenced by HDShowerSimulation(), readParameters(), reconstruct(), and updateHCAL().