SiPixelDigitizerAlgorithm Class Reference

#include <SiPixelDigitizerAlgorithm.h>

Classes
class	Amplitude
class	CalParameters
class	EnergyDepositUnit
class	SignalPoint

Public Member Functions
void	fillDeadModules (const edm::EventSetup &es)
void	fillLorentzAngle (const edm::EventSetup &es)
void	fillMapandGeom (const edm::EventSetup &es)
void	init (const edm::EventSetup &es)
std::vector< PixelDigiSimLink >	make_link ()
edm::DetSet< PixelDigi > ::collection_type	run (const std::vector< PSimHit > &input, PixelGeomDetUnit *pixdet, GlobalVector)
	SiPixelDigitizerAlgorithm (const edm::ParameterSet &conf, CLHEP::HepRandomEngine &)
	~SiPixelDigitizerAlgorithm ()

Private Types
typedef GloballyPositioned < double >	Frame
typedef std::vector < edm::ParameterSet >	Parameters
typedef signal_map_type::iterator	signal_map_iterator
typedef std::map< int, Amplitude, std::less< int > >	signal_map_type
typedef std::map< unsigned int, std::vector< float > , std::less< unsigned int > >	simlink_map

Private Member Functions
void	add_noise ()
std::vector< PixelDigi >	digitize (PixelGeomDetUnit *det)
void	drift (const PSimHit &hit)
LocalVector	DriftDirection ()
void	fluctuateEloss (int particleId, float momentum, float eloss, float length, int NumberOfSegments, float elossVector[])
void	induce_signal (const PSimHit &hit)
void	make_digis ()
float	missCalibrate (int col, int row, float amp) const
void	module_killing_conf ()
void	module_killing_DB ()
void	pixel_inefficiency ()
void	pixel_inefficiency_db ()
void	primary_ionization (const PSimHit &hit)

Private Attributes
GlobalVector	_bfield
std::vector< SignalPoint >	_collection_points
const PixelGeomDetUnit *	_detp
std::vector< EnergyDepositUnit >	_ionization_points
std::vector< PSimHit >	_PixelHits
signal_map_type	_signal
bool	addChargeVCALSmearing
bool	addNoise
bool	addNoisyPixels
bool	addPixelInefficiency
bool	addThresholdSmearing
bool	alpha2Order
float	BPix_p0
float	BPix_p1
float	BPix_p2
float	BPix_p3
std::map< int, CalParameters, std::less< int > >	calmap
float	ClusterWidth
edm::ParameterSet	conf_
Parameters	DeadModules
uint32_t	detID
float	Dist300
bool	doMissCalibrate
double	electronsPerVCAL
double	electronsPerVCAL_Offset
CLHEP::RandFlat *	flatDistribution_
SiG4UniversalFluctuation *	fluctuate
bool	fluctuateCharge
float	FPix_p0
float	FPix_p1
float	FPix_p2
float	FPix_p3
CLHEP::RandGaussQ *	gaussDistribution_
CLHEP::RandGaussQ *	gaussDistributionVCALNoise_
edm::ESHandle< TrackerGeometry >	geom_
float	GeVperElectron
std::vector< PixelDigi >	internal_coll
std::vector< PixelDigiSimLink >	link_coll
edm::ESHandle < SiPixelFedCablingMap >	map_
float	moduleThickness
int	NumberOfSegments
int	numColumns
int	numRows
PixelIndices *	pIndexConverter
float	PixelChipEff
float	PixelChipEfficiency
float	PixelColEff
float	PixelColEfficiency
float	PixelEff
float	PixelEfficiency
bool	pixelInefficiency
CLHEP::HepRandomEngine &	rndEngine
float	Sigma0
simlink_map	simi
edm::ESHandle< SiPixelQuality >	SiPixelBadModule_
edm::ESHandle < SiPixelLorentzAngle >	SiPixelLorentzAngle_
CLHEP::RandGaussQ *	smearedChargeDistribution_
CLHEP::RandGaussQ *	smearedThreshold_BPix_
CLHEP::RandGaussQ *	smearedThreshold_FPix_
float	tanLorentzAnglePerTesla_BPix
float	tanLorentzAnglePerTesla_FPix
int	theAdcFullScale
int	theColsInChip
float	theElectronPerADC
float	theGainSmearing
CLHEP::RandGaussQ *	theGaussianDistribution
float	theNoiseInElectrons
GaussianTailNoiseGenerator *	theNoiser
float	theOffsetSmearing
float	thePixelChipEfficiency [6]
float	thePixelColEfficiency [6]
float	thePixelEfficiency [6]
int	thePixelLuminosity
float	thePixelThreshold
float	thePixelThresholdInE
float	theReadoutNoise
int	theRowsInChip
SiPixelGainCalibrationOfflineSimService *	theSiPixelGainCalibrationService_
float	theSmearedChargeRMS
float	theThresholdInE_BPix
float	theThresholdInE_FPix
double	theThresholdSmearing_BPix
double	theThresholdSmearing_FPix
float	theTofLowerCut
float	theTofUpperCut
double	tMax
const PixelTopology *	topol
bool	use_deadmodule_DB_
bool	use_ineff_from_db_
bool	use_LorentzAngle_DB_
bool	use_module_killing_

Detailed Description

Definition at line 59 of file SiPixelDigitizerAlgorithm.h.

Member Typedef Documentation

typedef GloballyPositioned<double> SiPixelDigitizerAlgorithm::Frame

private

Definition at line 231 of file SiPixelDigitizerAlgorithm.h.

typedef std::vector<edm::ParameterSet> SiPixelDigitizerAlgorithm::Parameters

private

Definition at line 90 of file SiPixelDigitizerAlgorithm.h.

typedef signal_map_type::iterator SiPixelDigitizerAlgorithm::signal_map_iterator

private

Definition at line 228 of file SiPixelDigitizerAlgorithm.h.

typedef std::map< int, Amplitude, std::less<int> > SiPixelDigitizerAlgorithm::signal_map_type

private

Definition at line 227 of file SiPixelDigitizerAlgorithm.h.

typedef std::map<unsigned int, std::vector<float>,std::less<unsigned int> > SiPixelDigitizerAlgorithm::simlink_map

private

Definition at line 230 of file SiPixelDigitizerAlgorithm.h.

Constructor & Destructor Documentation

SiPixelDigitizerAlgorithm::SiPixelDigitizerAlgorithm	(	const edm::ParameterSet &	conf,
		CLHEP::HepRandomEngine &	eng
	)

Definition at line 140 of file SiPixelDigitizerAlgorithm.cc.

References addChargeVCALSmearing, addNoise, addNoisyPixels, addThresholdSmearing, alpha2Order, BPix_p0, BPix_p1, BPix_p2, BPix_p3, calmap, benchmark_cfg::cerr, PixelIndices::channelToPixelROC(), ClusterWidth, conf_, gather_cfg::cout, Dist300, doMissCalibrate, electronsPerVCAL, electronsPerVCAL_Offset, lut2db_cfg::filename, flatDistribution_, fluctuate, fluctuateCharge, FPix_p0, FPix_p1, FPix_p2, FPix_p3, gaussDistribution_, gaussDistributionVCALNoise_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), GeVperElectron, i, recoMuon::in, geometryCSVtoXML::line, LogDebug, mergeVDriftHistosByStation::name, NumberOfSegments, L1TEmulatorMonitor_cff::p, SiPixelDigitizerAlgorithm::CalParameters::p0, SiPixelDigitizerAlgorithm::CalParameters::p1, SiPixelDigitizerAlgorithm::CalParameters::p2, SiPixelDigitizerAlgorithm::CalParameters::p3, fitWZ::par0, pixelInefficiency, PixelIndices::pixelToChannelROC(), rndEngine, Sigma0, smearedThreshold_BPix_, smearedThreshold_FPix_, theAdcFullScale, theElectronPerADC, theGainSmearing, theNoiseInElectrons, theNoiser, theOffsetSmearing, thePixelChipEfficiency, thePixelColEfficiency, thePixelEfficiency, thePixelLuminosity, theReadoutNoise, theThresholdInE_BPix, theThresholdInE_FPix, theThresholdSmearing_BPix, theThresholdSmearing_FPix, theTofLowerCut, theTofUpperCut, and tMax.

                                                                                                           :
  conf_(conf) , fluctuate(0), theNoiser(0), pIndexConverter(0),
  use_ineff_from_db_(conf_.getParameter<bool>("useDB")),
  use_module_killing_(conf_.getParameter<bool>("killModules")), // boolean to kill or not modules
  use_deadmodule_DB_(conf_.getParameter<bool>("DeadModules_DB")), // boolean to access dead modules from DB
  use_LorentzAngle_DB_(conf_.getParameter<bool>("LorentzAngle_DB")), // boolean to access Lorentz angle from DB 
  theSiPixelGainCalibrationService_(0),rndEngine(eng)
{
  using std::cout;
  using std::endl;

  // Common pixel parameters
  // These are parameters which are not likely to be changed
  NumberOfSegments = 20; // Default number of track segment divisions
  ClusterWidth = 3.;     // Charge integration spread on the collection plane
  GeVperElectron = 3.61E-09; // 1 electron =3.61eV, 1keV=277e, mod 9/06 d.k.
  Sigma0 = 0.00037;           // Charge diffusion constant 7->3.7 
  Dist300 = 0.0300;          //   normalized to 300micron Silicon

  alpha2Order = conf_.getParameter<bool>("Alpha2Order");   // switch on/off of E.B effect   

  // get external parameters:
  // ADC calibration 1adc count = 135e.
  // Corresponds to 2adc/kev, 270[e/kev]/135[e/adc]=2[adc/kev]
  // Be carefull, this parameter is also used in SiPixelDet.cc to 
  // calculate the noise in adc counts from noise in electrons.
  // Both defaults should be the same.
  theElectronPerADC=conf_.getParameter<double>("ElectronPerAdc");

  // ADC saturation value, 255=8bit adc.
  //theAdcFullScale=conf_.getUntrackedParameter<int>("AdcFullScale",255);
  theAdcFullScale=conf_.getParameter<int>("AdcFullScale");

  // Pixel threshold in units of noise:
  // thePixelThreshold=conf_.getParameter<double>("ThresholdInNoiseUnits");
  // Pixel threshold in electron units.
  theThresholdInE_FPix=conf_.getParameter<double>("ThresholdInElectrons_FPix");
  theThresholdInE_BPix=conf_.getParameter<double>("ThresholdInElectrons_BPix");

  // Add threshold gaussian smearing:
  addThresholdSmearing = conf_.getParameter<bool>("AddThresholdSmearing");
  theThresholdSmearing_FPix = conf_.getParameter<double>("ThresholdSmearing_FPix");
  theThresholdSmearing_BPix = conf_.getParameter<double>("ThresholdSmearing_BPix");

  // signal response new parameterization: split Fpix and BPix 
  FPix_p0 = conf_.getParameter<double>("FPix_SignalResponse_p0");
  FPix_p1 = conf_.getParameter<double>("FPix_SignalResponse_p1");
  FPix_p2 = conf_.getParameter<double>("FPix_SignalResponse_p2");
  FPix_p3 = conf_.getParameter<double>("FPix_SignalResponse_p3");

  BPix_p0 = conf_.getParameter<double>("BPix_SignalResponse_p0");
  BPix_p1 = conf_.getParameter<double>("BPix_SignalResponse_p1");
  BPix_p2 = conf_.getParameter<double>("BPix_SignalResponse_p2");
  BPix_p3 = conf_.getParameter<double>("BPix_SignalResponse_p3");

  // electrons to VCAL conversion needed in misscalibrate()
  electronsPerVCAL = conf_.getParameter<double>("ElectronsPerVcal");
  electronsPerVCAL_Offset = conf_.getParameter<double>("ElectronsPerVcal_Offset");

  // Add noise   
  addNoise=conf_.getParameter<bool>("AddNoise");

  // Smear the pixel charge with a gaussian which RMS is a function of the
  // pixel charge (Danek's study)
  addChargeVCALSmearing=conf_.getParameter<bool>("ChargeVCALSmearing");

  // Add noisy pixels 
  addNoisyPixels=conf_.getParameter<bool>("AddNoisyPixels");
  // Noise in electrons:
  // Pixel cell noise, relevant for generating noisy pixels 
  theNoiseInElectrons=conf_.getParameter<double>("NoiseInElectrons");
  // Fill readout noise, including all readout chain, relevant for smearing
  //theReadoutNoise=conf_.getUntrackedParameter<double>("ReadoutNoiseInElec",500.);
  theReadoutNoise=conf_.getParameter<double>("ReadoutNoiseInElec");

  //theTofCut 12.5, cut in particle TOD +/- 12.5ns
  //theTofCut=conf_.getUntrackedParameter<double>("TofCut",12.5);
  theTofLowerCut=conf_.getParameter<double>("TofLowerCut");
  theTofUpperCut=conf_.getParameter<double>("TofUpperCut");
  
  // Fluctuate charge in track subsegments
  fluctuateCharge=conf_.getUntrackedParameter<bool>("FluctuateCharge",true);

  // delta cutoff in MeV, has to be same as in OSCAR=0.030/cmsim=1.0 MeV
  //tMax = 0.030; // In MeV.  
  //tMax =conf_.getUntrackedParameter<double>("DeltaProductionCut",0.030);  
  tMax =conf_.getParameter<double>("DeltaProductionCut");  
 
  // Control the pixel inefficiency
  thePixelLuminosity=conf_.getParameter<int>("AddPixelInefficiency");

  // Get the constants for the miss-calibration studies
  doMissCalibrate=conf_.getParameter<bool>("MissCalibrate"); // Enable miss-calibration
  theGainSmearing=conf_.getParameter<double>("GainSmearing"); // sigma of the gain smearing
  theOffsetSmearing=conf_.getParameter<double>("OffsetSmearing"); //sigma of the offset smearing

  //pixel inefficiency
  // the first 3 settings [0],[1],[2] are for the barrel pixels
  // the next  3 settings [3],[4],[5] are for the endcaps (undecided how)  

  if (thePixelLuminosity==-1) {  // No inefficiency, all 100% efficient
    pixelInefficiency=false;
    for (int i=0; i<6;i++) {
      thePixelEfficiency[i]     = 1.;  // pixels = 100%
      thePixelColEfficiency[i]  = 1.;  // columns = 100%
      thePixelChipEfficiency[i] = 1.; // chips = 100%
    }
    
    // include only the static (non rate depedent) efficiency 
    // Usefull for very low rates (luminosity)
  } else if (thePixelLuminosity==0) { // static effciency
    pixelInefficiency=true;
    // Default efficiencies 
    for (int i=0; i<6;i++) {
      if(i<3) {  // For the barrel
    // Assume 1% inefficiency for single pixels, 
    // this is given by faulty bump-bonding and seus.  
    thePixelEfficiency[i]     = 1.-0.001;  // pixels = 99.9%
    // For columns make 0.1% default.
    thePixelColEfficiency[i]  = 1.-0.001;  // columns = 99.9%
    // A flat 0.1% inefficiency due to lost rocs
    thePixelChipEfficiency[i] = 1.-0.001; // chips = 99.9%
      } else { // For the endcaps
    // Assume 1% inefficiency for single pixels, 
    // this is given by faulty bump-bonding and seus.  
    thePixelEfficiency[i]     = 1.-0.001;  // pixels = 99.9%
    // For columns make 0.1% default.
    thePixelColEfficiency[i]  = 1.-0.001;  // columns = 99.9%
    // A flat 0.1% inefficiency due to lost rocs
    thePixelChipEfficiency[i] = 1.-0.001; // chips = 99.9%
      }
    }
    
    // Include also luminosity rate dependent inefficieny
  } else if (thePixelLuminosity>0) { // Include effciency
    pixelInefficiency=true;
    // Default efficiencies 
    for (int i=0; i<6;i++) {
      if(i<3) { // For the barrel
    // Assume 1% inefficiency for single pixels, 
    // this is given by faulty bump-bonding and seus.  
    thePixelEfficiency[i]     = 1.-0.01;  // pixels = 99%
    // For columns make 1% default.
    thePixelColEfficiency[i]  = 1.-0.01;  // columns = 99%
    // A flat 0.25% inefficiency due to lost data packets from TBM
    thePixelChipEfficiency[i] = 1.-0.0025; // chips = 99.75%
      } else { // For the endcaps
    // Assume 1% inefficiency for single pixels, 
    // this is given by faulty bump-bonding and seus.  
    thePixelEfficiency[i]     = 1.-0.01;  // pixels = 99%
    // For columns make 1% default.
    thePixelColEfficiency[i]  = 1.-0.01;  // columns = 99%
    // A flat 0.25% inefficiency due to lost data packets from TBM
    thePixelChipEfficiency[i] = 1.-0.0025; // chips = 99.75%
      }
    }
   
    // Special cases ( High-lumi for 4cm layer) where the readout losses are higher
    if(thePixelLuminosity==10) { // For high luminosity, bar layer 1
      thePixelColEfficiency[0] = 1.-0.034; // 3.4% for r=4 only
      thePixelEfficiency[0]    = 1.-0.015; // 1.5% for r=4
    }
    
  } // end the pixel inefficiency part

  // Init the random number services  
  if(addNoise || thePixelLuminosity || fluctuateCharge || addThresholdSmearing ) {
    gaussDistribution_ = new CLHEP::RandGaussQ(rndEngine, 0., theReadoutNoise);
    gaussDistributionVCALNoise_ = new CLHEP::RandGaussQ(rndEngine, 0., 1.);
    flatDistribution_ = new CLHEP::RandFlat(rndEngine, 0., 1.);
    
    if(addNoise) { 
      theNoiser = new GaussianTailNoiseGenerator(rndEngine);
    }
    
    if(fluctuateCharge) {
      fluctuate = new SiG4UniversalFluctuation(rndEngine);
    }

    // Threshold smearing with gaussian distribution:
    if(addThresholdSmearing) {
      smearedThreshold_FPix_ = new CLHEP::RandGaussQ(rndEngine, theThresholdInE_FPix , theThresholdSmearing_FPix);
      smearedThreshold_BPix_ = new CLHEP::RandGaussQ(rndEngine, theThresholdInE_BPix , theThresholdSmearing_BPix);
    }
    } //end Init the random number services


  // Prepare for the analog amplitude miss-calibration
  if(doMissCalibrate) {
    LogDebug ("PixelDigitizer ") 
      << " miss-calibrate the pixel amplitude "; 
    
    const bool ReadCalParameters = false;
    if(ReadCalParameters) {   // Read the calibration files from file
      // read the calibration constants from a file (testing only)
      ifstream in_file;  // data file pointer
      char filename[80] = "phCalibrationFit_C0.dat";
      
      in_file.open(filename, ios::in ); // in C++
      if (in_file.bad()) {
    cout << " File not found " << endl;
    return; // signal error
      }
      cout << " file opened : " << filename << endl;
      
      char line[500];
      for (int i = 0; i < 3; i++) {
    in_file.getline(line, 500,'\n');
    cout<<line<<endl;
      }
      
      cout << " test map" << endl;
      
      float par0,par1,par2,par3;
      int colid,rowid;
      string name;
      // Read MC tracks
      for(int i=0;i<(52*80);i++)  { // loop over tracks    
    in_file >> par0 >> par1 >> par2 >> par3 >> name >> colid
        >> rowid;
    if (in_file.bad()) { // check for errors
      cerr << "Cannot read data file" << endl;
      return;
    }
    if ( in_file.eof() != 0 ) {
      cerr << in_file.eof() << " " << in_file.gcount() << " "
           << in_file.fail() << " " << in_file.good() << " end of file "
           << endl;
      return;
    }
    
    //cout << " line " << i << " " <<par0<<" "<<par1<<" "<<par2<<" "<<par3<<" "
    //   <<colid<<" "<<rowid<<endl;
    
    CalParameters onePix;
    onePix.p0=par0;
    onePix.p1=par1;
    onePix.p2=par2;
    onePix.p3=par3;
    
    // Convert ROC pixel index to channel 
    int chan = PixelIndices::pixelToChannelROC(rowid,colid);
    calmap.insert(pair<int,CalParameters>(chan,onePix));
    
    // Testing the index conversion, can be skipped 
    pair<int,int> p = PixelIndices::channelToPixelROC(chan);
    if(rowid!=p.first) cout<<" wrong channel row "<<rowid<<" "<<p.first<<endl;
    if(colid!=p.second) cout<<" wrong channel col "<<colid<<" "<<p.second<<endl;
          
      } // pixel loop in a ROC
 
      cout << " map size  " << calmap.size() <<" max "<<calmap.max_size() << " "
       <<calmap.empty()<< endl;
      
//     cout << " map size  " << calmap.size()  << endl;
//     map<int,CalParameters,less<int> >::iterator ix,it;
//     map<int,CalParameters,less<int> >::const_iterator ip;
//     for (ix = calmap.begin(); ix != calmap.end(); ++ix) {
//       int i = (*ix).first;
//       pair<int,int> p = channelToPixelROC(i);
//       it  = calmap.find(i);
//       CalParameters y  = (*it).second;
//       CalParameters z = (*ix).second;
//       cout << i <<" "<<p.first<<" "<<p.second<<" "<<y.p0<<" "<<z.p0<<" "<<calmap[i].p0<<endl; 
      
//       //int dummy=0;
//       //cin>>dummy;
//     }

    } // end if readparameters
  } // end if missCalibration 

  LogInfo ("PixelDigitizer ") <<"SiPixelDigitizerAlgorithm constructed"
                  <<"Configuration parameters:" 
                  << "Threshold/Gain = "  
                  << "threshold in electron FPix = " 
                  << theThresholdInE_FPix
                  << "threshold in electron BPix = " 
                  << theThresholdInE_BPix 
                  <<" " << theElectronPerADC << " " << theAdcFullScale 
                  << " The delta cut-off is set to " << tMax
                  << " pix-inefficiency "<<thePixelLuminosity;

}

SiPixelDigitizerAlgorithm::~SiPixelDigitizerAlgorithm

(

)

Definition at line 425 of file SiPixelDigitizerAlgorithm.cc.

References addNoise, addThresholdSmearing, flatDistribution_, fluctuate, fluctuateCharge, gaussDistribution_, gaussDistributionVCALNoise_, LogDebug, smearedThreshold_BPix_, smearedThreshold_FPix_, theNoiser, and theSiPixelGainCalibrationService_.

                                                      {

  LogDebug ("PixelDigitizer")<<"SiPixelDigitizerAlgorithm deleted";
  
  // Destructor
  delete gaussDistribution_;
  delete gaussDistributionVCALNoise_;
  delete flatDistribution_;

  delete theSiPixelGainCalibrationService_;

  // Threshold gaussian smearing:
  if(addThresholdSmearing) {
    delete smearedThreshold_FPix_;
    delete smearedThreshold_BPix_;
  }

  if(addNoise) delete theNoiser;
  if(fluctuateCharge) delete fluctuate;

}

Member Function Documentation

void SiPixelDigitizerAlgorithm::add_noise ( )

private

Definition at line 1144 of file SiPixelDigitizerAlgorithm.cc.

References _signal, addChargeVCALSmearing, addNoisyPixels, gaussDistribution_, gaussDistributionVCALNoise_, GaussianTailNoiseGenerator::generate(), i, LogDebug, numColumns, numRows, PixelDigi::pixelToChannel(), theNoiseInElectrons, theNoiser, thePixelThreshold, theSmearedChargeRMS, theThresholdInE_BPix, and theThresholdInE_FPix.

Referenced by digitize().

                                          {

#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") << " enter add_noise " << theNoiseInElectrons;
#endif
 
  // First add noise to hit pixels
  
  for ( signal_map_iterator i = _signal.begin(); i != _signal.end(); i++) {
    
         if(addChargeVCALSmearing) 
      {
    if((*i).second < 3000)
      {
        theSmearedChargeRMS = 543.6 - (*i).second * 0.093;
      } else if((*i).second < 6000){
        theSmearedChargeRMS = 307.6 - (*i).second * 0.01;
      } else{
        theSmearedChargeRMS = -432.4 +(*i).second * 0.123;
    }

    // Noise from Vcal smearing:
    float noise_ChargeVCALSmearing = theSmearedChargeRMS * gaussDistributionVCALNoise_->fire() ;
    // Noise from full readout:
    float noise  = gaussDistribution_->fire() ;

        if(((*i).second + Amplitude(noise+noise_ChargeVCALSmearing,0,-1.)) < 0. ) {
          (*i).second.set(0);}
        else{
    (*i).second +=Amplitude(noise+noise_ChargeVCALSmearing,0,-1.);
        }

      } // End if addChargeVCalSmearing
     else
     {
    // Noise: ONLY full READOUT Noise.
    // Use here the FULL readout noise, including TBM,ALT,AOH,OPT-REC.

    float noise  = gaussDistribution_->fire() ;
        if(((*i).second + Amplitude(noise,0,-1.)) < 0. ) {
          (*i).second.set(0);}
        else{
    (*i).second +=Amplitude(noise,0,-1.);
        }
     } // end if only Noise from full readout
    
  }

  if(!addNoisyPixels)  // Option to skip noise in non-hit pixels
    return;

  // Add noise on non-hit pixels
  // Use here the pixel noise 
  int numberOfPixels = (numRows * numColumns);
  map<int,float, less<int> > otherPixels;
  map<int,float, less<int> >::iterator mapI;

  theNoiser->generate(numberOfPixels, 
                      thePixelThreshold, //thr. in un. of nois
              theNoiseInElectrons, // noise in elec. 
                      otherPixels );

#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") 
    <<  " Add noisy pixels " << numRows << " " 
    << numColumns << " " << theNoiseInElectrons << " " 
    << theThresholdInE_FPix << theThresholdInE_BPix <<" "<< numberOfPixels<<" " 
    << otherPixels.size() ;
#endif
  
  // Add noisy pixels
  for (mapI = otherPixels.begin(); mapI!= otherPixels.end(); mapI++) {
    int iy = ((*mapI).first) / numRows;
    int ix = ((*mapI).first) - (iy*numRows);

    // Keep for a while for testing.
    if( iy < 0 || iy > (numColumns-1) ) 
      LogWarning ("Pixel Geometry") << " error in iy " << iy ;
    if( ix < 0 || ix > (numRows-1) )
      LogWarning ("Pixel Geometry")  << " error in ix " << ix ;

    int chan = PixelDigi::pixelToChannel(ix, iy);

#ifdef TP_DEBUG
    LogDebug ("Pixel Digitizer")
      <<" Storing noise = " << (*mapI).first << " " << (*mapI).second 
      << " " << ix << " " << iy << " " << chan ;
#endif
 
    if(_signal[chan] == 0){
      //      float noise = float( (*mapI).second );
      int noise=int( (*mapI).second );
      _signal[chan] = Amplitude (noise, 0,-1.);
    }
  }

}

vector< PixelDigi > SiPixelDigitizerAlgorithm::digitize ( PixelGeomDetUnit *  det )

private

Definition at line 482 of file SiPixelDigitizerAlgorithm.cc.

References _collection_points, _PixelHits, _signal, add_noise(), addNoise, addNoisyPixels, addThresholdSmearing, BoundSurface::bounds(), detID, doMissCalibrate, drift(), induce_signal(), internal_coll, LogDebug, make_digis(), module_killing_conf(), module_killing_DB(), moduleThickness, PixelTopology::ncolumns(), PixelTopology::nrows(), numColumns, numRows, pIndexConverter, pixel_inefficiency(), pixel_inefficiency_db(), PixelSubdetector::PixelBarrel, pixelInefficiency, primary_ionization(), smearedThreshold_BPix_, smearedThreshold_FPix_, GeomDet::specificSurface(), PixelGeomDetUnit::specificTopology(), DetId::subdetId(), theNoiseInElectrons, thePixelThreshold, thePixelThresholdInE, theThresholdInE_BPix, theThresholdInE_FPix, theTofLowerCut, theTofUpperCut, Bounds::thickness(), topol, use_deadmodule_DB_, use_ineff_from_db_, and use_module_killing_.

Referenced by run().

                                                                          {
  
  if( _PixelHits.size() > 0 || addNoisyPixels) {
    
    topol=&det->specificTopology(); // cache topology
    numColumns = topol->ncolumns();  // det module number of cols&rows
    numRows = topol->nrows();
    
    // full detector thickness
    moduleThickness = det->specificSurface().bounds().thickness(); 
    
    // The index converter is only needed when inefficiencies or misscalibration
    // are simulated.
    if((pixelInefficiency>0) || doMissCalibrate ) {  // Init pixel indices
      pIndexConverter = new PixelIndices(numColumns,numRows);
    }
    
    // Noise already defined in electrons
    //thePixelThresholdInE = thePixelThreshold * theNoiseInElectrons ; 
    // Find the threshold in noise units, needed for the noiser.

  unsigned int Sub_detid=DetId(detID).subdetId();

  if(theNoiseInElectrons>0.){
    if(Sub_detid == PixelSubdetector::PixelBarrel){ // Barrel modules
      if(addThresholdSmearing) { 
    thePixelThresholdInE = smearedThreshold_BPix_->fire(); // gaussian smearing 
      } else {
    thePixelThresholdInE = theThresholdInE_BPix; // no smearing
      }
      
      thePixelThreshold = thePixelThresholdInE/theNoiseInElectrons; 
      
    } else { // Forward disks modules 
      if(addThresholdSmearing) {
    thePixelThresholdInE = smearedThreshold_FPix_->fire(); // gaussian smearing
      } else {
    thePixelThresholdInE = theThresholdInE_FPix; // no smearing
      }
      
      thePixelThreshold = thePixelThresholdInE/theNoiseInElectrons;
      
    }
  } else {
    thePixelThreshold = 0.;
  }


#ifdef TP_DEBUG
    LogDebug ("PixelDigitizer") 
      << " PixelDigitizer "  
      << numColumns << " " << numRows << " " << moduleThickness;
#endif

    // produce SignalPoint's for all SimHit's in detector
    // Loop over hits
 
    vector<PSimHit>::const_iterator ssbegin; 
    for (ssbegin= _PixelHits.begin();ssbegin !=_PixelHits.end(); ++ssbegin) {
      
#ifdef TP_DEBUG
      LogDebug ("Pixel Digitizer") 
    << (*ssbegin).particleType() << " " << (*ssbegin).pabs() << " " 
    << (*ssbegin).energyLoss() << " " << (*ssbegin).tof() << " " 
    << (*ssbegin).trackId() << " " << (*ssbegin).processType() << " " 
    << (*ssbegin).detUnitId()  
    << (*ssbegin).entryPoint() << " " << (*ssbegin).exitPoint() ; 
#endif      
      
      _collection_points.clear();  // Clear the container
      // fill _collection_points for this SimHit, indpendent of topology
      // Check the TOF cut
      //if (std::abs( (*ssbegin).tof() )<theTofCut){ // old cut
      if ( ((*ssbegin).tof() >= theTofLowerCut) && ((*ssbegin).tof() <= theTofUpperCut) ) {
    primary_ionization(*ssbegin); // fills _ionization_points   
    drift(*ssbegin);  // transforms _ionization_points to _collection_points    
    // compute induced signal on readout elements and add to _signal
    induce_signal(*ssbegin); // *ihit needed only for SimHit<-->Digi link
      } //  end if 
    } // end for 
  
    if(addNoise) add_noise();  // generate noise
    // Do only if needed 

    if((pixelInefficiency>0) && (_signal.size()>0)) 
      pixel_inefficiency(); // Kill some pixels

    if(use_ineff_from_db_ && (_signal.size()>0))
      pixel_inefficiency_db();

    delete pIndexConverter;

    if(use_module_killing_ && use_deadmodule_DB_) // remove dead modules using DB
      module_killing_DB();

    if(use_module_killing_ && !use_deadmodule_DB_) // remove dead modules using the list in cfg file
      module_killing_conf();

  }

  make_digis();

  return internal_coll;
}

void SiPixelDigitizerAlgorithm::drift ( const PSimHit &  hit )

private

Definition at line 717 of file SiPixelDigitizerAlgorithm.cc.

References _collection_points, _ionization_points, alpha2Order, Dist300, DriftDirection(), i, LogDebug, moduleThickness, FitTarget::Sigma, Sigma0, mathSSE::sqrt(), PSimHit::tof(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by digitize().

                                                       {


#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") << " enter drift " ;
#endif

  _collection_points.resize( _ionization_points.size()); // set size
  
  LocalVector driftDir=DriftDirection();  // get the charge drift direction
  if(driftDir.z() ==0.) {
    LogWarning("Magnetic field") << " pxlx: drift in z is zero ";
    return;
  }  

  // tangent of Lorentz angle
  //float TanLorenzAngleX = driftDir.x()/driftDir.z(); 
  //float TanLorenzAngleY = 0.; // force to 0, driftDir.y()/driftDir.z();

  float TanLorenzAngleX, TanLorenzAngleY,dir_z, CosLorenzAngleX,
    CosLorenzAngleY;
  if ( alpha2Order) {
    
    TanLorenzAngleX = driftDir.x(); // tangen of Lorentz angle
    TanLorenzAngleY = driftDir.y();
    dir_z = driftDir.z(); // The z drift direction
    CosLorenzAngleX = 1./sqrt(1.+TanLorenzAngleX*TanLorenzAngleX); //cosine
    CosLorenzAngleY = 1./sqrt(1.+TanLorenzAngleY*TanLorenzAngleY); //cosine;
    
  } else{
    
    TanLorenzAngleX = driftDir.x();
    TanLorenzAngleY = 0.; // force to 0, driftDir.y()/driftDir.z();
    dir_z = driftDir.z(); // The z drift direction
    CosLorenzAngleX = 1./sqrt(1.+TanLorenzAngleX*TanLorenzAngleX); //cosine to estimate the path length
    CosLorenzAngleY = 1.;
  }
  
  
#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") 
    << " Lorentz Tan " << TanLorenzAngleX << " " << TanLorenzAngleY <<" "
    << CosLorenzAngleX << " " << CosLorenzAngleY << " "
    << moduleThickness*TanLorenzAngleX << " " << driftDir;
#endif  
  
  float Sigma_x = 1.;  // Charge spread 
  float Sigma_y = 1.;
  float DriftDistance; // Distance between charge generation and collection 
  float DriftLength;   // Actual Drift Lentgh
  float Sigma;
  
  
  for (unsigned int i = 0; i != _ionization_points.size(); i++) {
    
    float SegX, SegY, SegZ; // position
    SegX = _ionization_points[i].x();
    SegY = _ionization_points[i].y();
    SegZ = _ionization_points[i].z();
    
    // Distance from the collection plane
    //DriftDistance = (moduleThickness/2. + SegZ); // Drift to -z 
    // Include explixitely the E drift direction (for CMS dir_z=-1)
    DriftDistance = moduleThickness/2. - (dir_z * SegZ); // Drift to -z 
    
    //if( DriftDistance <= 0.) 
    //cout<<" <=0 "<<DriftDistance<<" "<<i<<" "<<SegZ<<" "<<dir_z<<" "
    //  <<SegX<<" "<<SegY<<" "<<(moduleThickness/2)<<" "
    //  <<_ionization_points[i].energy()<<" "
    //  <<hit.particleType()<<" "<<hit.pabs()<<" "<<hit.energyLoss()<<" "
    //  <<hit.entryPoint()<<" "<<hit.exitPoint()
    //  <<endl;

    if( DriftDistance < 0.) {
      DriftDistance = 0.;
    } else if ( DriftDistance > moduleThickness )
      DriftDistance = moduleThickness;
    
    // Assume full depletion now, partial depletion will come later.
    float XDriftDueToMagField = DriftDistance * TanLorenzAngleX;
    float YDriftDueToMagField = DriftDistance * TanLorenzAngleY;
    
    // Shift cloud center
    float CloudCenterX = SegX + XDriftDueToMagField;
    float CloudCenterY = SegY + YDriftDueToMagField;

    // Calculate how long is the charge drift path
    DriftLength = sqrt( DriftDistance*DriftDistance + 
                        XDriftDueToMagField*XDriftDueToMagField +
                        YDriftDueToMagField*YDriftDueToMagField );

    // What is the charge diffusion after this path
    Sigma = sqrt(DriftLength/Dist300) * Sigma0;

    // Project the diffusion sigma on the collection plane
    Sigma_x = Sigma / CosLorenzAngleX ;
    Sigma_y = Sigma / CosLorenzAngleY ;

    SignalPoint sp( CloudCenterX, CloudCenterY,
     Sigma_x, Sigma_y, hit.tof(), _ionization_points[i].energy() );

    // Load the Charge distribution parameters
    _collection_points[i] = (sp);

  } // loop over ionization points, i.
 
} // end drift

LocalVector SiPixelDigitizerAlgorithm::DriftDirection ( )

private

Definition at line 1445 of file SiPixelDigitizerAlgorithm.cc.

References _bfield, _detp, alpha2Order, detID, LogDebug, PixelSubdetector::PixelBarrel, GloballyPositioned< T >::position(), GloballyPositioned< T >::rotation(), SiPixelLorentzAngle_, DetId::subdetId(), GeomDet::surface(), tanLorentzAnglePerTesla_BPix, tanLorentzAnglePerTesla_FPix, and use_LorentzAngle_DB_.

Referenced by drift().

                                                     {
  Frame detFrame(_detp->surface().position(),_detp->surface().rotation());
  LocalVector Bfield=detFrame.toLocal(_bfield);
  
  float alpha2_FPix;
  float alpha2_BPix;
  float alpha2;

  //float dir_x = -tanLorentzAnglePerTesla * Bfield.y();
  //float dir_y = +tanLorentzAnglePerTesla * Bfield.x();
  //float dir_z = -1.; // E field always in z direction, so electrons go to -z
  // The dir_z has to be +/- 1. !
  // LocalVector theDriftDirection = LocalVector(dir_x,dir_y,dir_z);

  float dir_x = 0.0;
  float dir_y = 0.0;
  float dir_z = 0.0;
  float scale = 0.0;
  
  unsigned int Sub_detid=DetId(detID).subdetId();

  // Read Lorentz angle from cfg file:**************************************************************
  
  if(!use_LorentzAngle_DB_){ 
    
    if ( alpha2Order) {
      alpha2_FPix = tanLorentzAnglePerTesla_FPix*tanLorentzAnglePerTesla_FPix;
      alpha2_BPix = tanLorentzAnglePerTesla_BPix*tanLorentzAnglePerTesla_BPix;
    }else {
      alpha2_FPix = 0.0;
      alpha2_BPix = 0.0;
    }
    
    if (Sub_detid == PixelSubdetector::PixelBarrel){// barrel layers
      dir_x = -( tanLorentzAnglePerTesla_BPix * Bfield.y() + alpha2_BPix* Bfield.z()* Bfield.x() );
      dir_y = +( tanLorentzAnglePerTesla_BPix * Bfield.x() - alpha2_BPix* Bfield.z()* Bfield.y() );
      dir_z = -(1 + alpha2_BPix* Bfield.z()*Bfield.z() );
      scale = (1 + alpha2_BPix* Bfield.z()*Bfield.z() );
      
    } else {// forward disks
      dir_x = -( tanLorentzAnglePerTesla_FPix * Bfield.y() + alpha2_FPix* Bfield.z()* Bfield.x() );
      dir_y = +( tanLorentzAnglePerTesla_FPix * Bfield.x() - alpha2_FPix* Bfield.z()* Bfield.y() );
      dir_z = -(1 + alpha2_FPix* Bfield.z()*Bfield.z() );
      scale = (1 + alpha2_FPix* Bfield.z()*Bfield.z() );
    }
    } // end: Read LA from cfg file.
  
  //Read Lorentz angle from DB:********************************************************************
  if(use_LorentzAngle_DB_){ 
    std::map<unsigned int,float> detid_la= SiPixelLorentzAngle_->getLorentzAngles();
    std::map<unsigned int,float>::const_iterator it;
    
    
    for (it=detid_la.begin();it!=detid_la.end();it++)
    {
      if (detID==it->first) {
    if (alpha2Order) {
      alpha2 = it->second * it->second;
    }  
    else {
      alpha2 = 0.0;
    } 
    //  std::cout << "detID is: " << it->first << "The LA per tesla is: " << it->second << std::endl;
    dir_x = -( it->second * Bfield.y() + alpha2 * Bfield.z()* Bfield.x() );
    dir_y = +( it->second * Bfield.x() - alpha2 * Bfield.z()* Bfield.y() );
    dir_z = -(1 + alpha2 * Bfield.z()*Bfield.z() );
    scale = (1 + alpha2 * Bfield.z()*Bfield.z() );
      }
    } 
  }// end: Read LA from DataBase.
  
  LocalVector theDriftDirection = LocalVector(dir_x/scale, dir_y/scale, dir_z/scale );  
  
#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") << " The drift direction in local coordinate is "   
                   << theDriftDirection ;
#endif
  
  return theDriftDirection;
}

void SiPixelDigitizerAlgorithm::fillDeadModules

(

const edm::EventSetup &

es

)

Definition at line 109 of file SiPixelDigitizerAlgorithm.cc.

References edm::EventSetup::get().

                                                                      {
  if(!use_deadmodule_DB_){
    DeadModules = conf_.getParameter<Parameters>("DeadModules"); // get dead module from cfg file
  }
  else{  // Get dead module from DB record 
    // ESHandle was defined in the header file   edm::ESHandle<SiPixelQuality> SiPixelBadModule_;
    es.get<SiPixelQualityRcd>().get(SiPixelBadModule_);
  }
}

void SiPixelDigitizerAlgorithm::fillLorentzAngle

(

const edm::EventSetup &

es

)

Definition at line 119 of file SiPixelDigitizerAlgorithm.cc.

References edm::EventSetup::get().

                                                                       {
  if(!use_LorentzAngle_DB_){
    // Get the Lorentz angle from the cfg file:
    tanLorentzAnglePerTesla_FPix=conf_.getParameter<double>("TanLorentzAnglePerTesla_FPix");
    tanLorentzAnglePerTesla_BPix=conf_.getParameter<double>("TanLorentzAnglePerTesla_BPix");
  } 
  else {
    // Get Lorentz angle from DB record 
    // ESHandle was defined in the header file edm::ESHandle<SiPixelLorentzAngle> SiPixelLorentzAngle_;
    es.get<SiPixelLorentzAngleSimRcd>().get(SiPixelLorentzAngle_);
  }
}

void SiPixelDigitizerAlgorithm::fillMapandGeom

(

const edm::EventSetup &

es

)

Definition at line 133 of file SiPixelDigitizerAlgorithm.cc.

References edm::EventSetup::get().

                                                                     {
   es.get<SiPixelFedCablingMapRcd>().get(map_);
   es.get<TrackerDigiGeometryRecord>().get(geom_);
}

void SiPixelDigitizerAlgorithm::fluctuateEloss ( int  particleId, float  momentum, float  eloss, float  length, int  NumberOfSegments, float  elossVector[]  )

private

Definition at line 662 of file SiPixelDigitizerAlgorithm.cc.

References abs, fluctuate, i, SiG4UniversalFluctuation::SampleFluctuations(), and tMax.

Referenced by primary_ionization().

                                                            {

  // Get dedx for this track
  float dedx;
  if( length > 0.) dedx = eloss/length;
  else dedx = eloss;

  double particleMass = 139.6; // Mass in MeV, Assume pion
  pid = abs(pid);
  if(pid!=211) {       // Mass in MeV
    if(pid==11)        particleMass = 0.511;        
    else if(pid==13)   particleMass = 105.7;
    else if(pid==321)  particleMass = 493.7;
    else if(pid==2212) particleMass = 938.3;
  }
  // What is the track segment length.
  float segmentLength = length/NumberOfSegs;

  // Generate charge fluctuations.
  float de=0.;
  float sum=0.;
  double segmentEloss = (1000.*eloss)/NumberOfSegs; //eloss in MeV
  for (int i=0;i<NumberOfSegs;i++) {
    //       material,*,   momentum,energy,*, *,  mass
    //myglandz_(14.,segmentLength,2.,2.,dedx,de,0.14);
    // The G4 routine needs momentum in MeV, mass in Mev, delta-cut in MeV,
    // track segment length in mm, segment eloss in MeV 
    // Returns fluctuated eloss in MeV
    double deltaCutoff = tMax; // the cutoff is sometimes redefined inside, so fix it.
    de = fluctuate->SampleFluctuations(double(particleMomentum*1000.),
                      particleMass, deltaCutoff, 
                      double(segmentLength*10.),
                      segmentEloss )/1000.; //convert to GeV
   
    elossVector[i]=de;
    sum +=de;
  }
  
  if(sum>0.) {  // If fluctuations give eloss>0.
    // Rescale to the same total eloss
    float ratio = eloss/sum;
   
    for (int ii=0;ii<NumberOfSegs;ii++) elossVector[ii]= ratio*elossVector[ii];
  } else {  // If fluctuations gives 0 eloss
    float averageEloss = eloss/NumberOfSegs;
    for (int ii=0;ii<NumberOfSegs;ii++) elossVector[ii]= averageEloss; 
  }
  return;
}

void SiPixelDigitizerAlgorithm::induce_signal ( const PSimHit &  hit )

private

Definition at line 827 of file SiPixelDigitizerAlgorithm.cc.

References _collection_points, _signal, Topology::channel(), PixelDigi::channelToPixel(), ClusterWidth, i, Topology::localPosition(), LogDebug, Topology::measurementPosition(), numColumns, numRows, PixelTopology::pitch(), PixelDigi::pixelToChannel(), query::result, ntuplemaker::status, topol, x, PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), detailsBasic3DVector::y, PV2DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by digitize().

                                                                 {

  // X  - Rows, Left-Right, 160, (1.6cm)   for barrel
  // Y  - Columns, Down-Up, 416, (6.4cm)

#ifdef TP_DEBUG
    LogDebug ("Pixel Digitizer") 
      << " enter induce_signal, " 
      << topol->pitch().first << " " << topol->pitch().second; //OK
#endif

   // local map to store pixels hit by 1 Hit.      
   typedef map< int, float, less<int> > hit_map_type;
   hit_map_type hit_signal;

   // map to store pixel integrals in the x and in the y directions
   map<int, float, less<int> > x,y; 
   
   // Assign signals to readout channels and store sorted by channel number
   
   // Iterate over collection points on the collection plane
   for ( vector<SignalPoint>::const_iterator i=_collection_points.begin();
     i != _collection_points.end(); i++) {
     
     float CloudCenterX = i->position().x(); // Charge position in x
     float CloudCenterY = i->position().y(); //                 in y
     float SigmaX = i->sigma_x();            // Charge spread in x
     float SigmaY = i->sigma_y();            //               in y
     float Charge = i->amplitude();          // Charge amplitude


     //if(SigmaX==0 || SigmaY==0) {
     //cout<<SigmaX<<" "<<SigmaY
     //   << " cloud " << i->position().x() << " " << i->position().y() << " " 
     //   << i->sigma_x() << " " << i->sigma_y() << " " << i->amplitude()<<endl;
     //}

#ifdef TP_DEBUG
       LogDebug ("Pixel Digitizer") 
     << " cloud " << i->position().x() << " " << i->position().y() << " " 
     << i->sigma_x() << " " << i->sigma_y() << " " << i->amplitude();
#endif      
     
     // Find the maximum cloud spread in 2D plane , assume 3*sigma
     float CloudRight = CloudCenterX + ClusterWidth*SigmaX;
     float CloudLeft  = CloudCenterX - ClusterWidth*SigmaX;
     float CloudUp    = CloudCenterY + ClusterWidth*SigmaY;
     float CloudDown  = CloudCenterY - ClusterWidth*SigmaY;
     
     // Define 2D cloud limit points
     LocalPoint PointRightUp  = LocalPoint(CloudRight,CloudUp);
     LocalPoint PointLeftDown = LocalPoint(CloudLeft,CloudDown);
     
     // This points can be located outside the sensor area.
     // The conversion to measurement point does not check for that
     // so the returned pixel index might be wrong (outside range).
     // We rely on the limits check below to fix this.
     // But remember whatever we do here THE CHARGE OUTSIDE THE ACTIVE
     // PIXEL ARE IS LOST, it should not be collected.

     // Convert the 2D points to pixel indices
     MeasurementPoint mp = topol->measurementPosition(PointRightUp ); //OK
     
     int IPixRightUpX = int( floor( mp.x()));
     int IPixRightUpY = int( floor( mp.y()));
     
#ifdef TP_DEBUG
     LogDebug ("Pixel Digitizer") << " right-up " << PointRightUp << " " 
                  << mp.x() << " " << mp.y() << " "
                  << IPixRightUpX << " " << IPixRightUpY ;
#endif
 
     mp = topol->measurementPosition(PointLeftDown ); //OK
    
     int IPixLeftDownX = int( floor( mp.x()));
     int IPixLeftDownY = int( floor( mp.y()));
     
#ifdef TP_DEBUG
     LogDebug ("Pixel Digitizer") << " left-down " << PointLeftDown << " " 
                  << mp.x() << " " << mp.y() << " "
                  << IPixLeftDownX << " " << IPixLeftDownY ;
#endif

     // Check detector limits to correct for pixels outside range.
     IPixRightUpX = numRows>IPixRightUpX ? IPixRightUpX : numRows-1 ;
     IPixRightUpY = numColumns>IPixRightUpY ? IPixRightUpY : numColumns-1 ;
     IPixLeftDownX = 0<IPixLeftDownX ? IPixLeftDownX : 0 ;
     IPixLeftDownY = 0<IPixLeftDownY ? IPixLeftDownY : 0 ;
     
     x.clear(); // clear temporary integration array
     y.clear();

     // First integrate cahrge strips in x
     int ix; // TT for compatibility
     for (ix=IPixLeftDownX; ix<=IPixRightUpX; ix++) {  // loop over x index
       float xUB, xLB, UpperBound, LowerBound;
      
       // Why is set to 0 if ix=0, does it meen that we accept charge 
       // outside the sensor? CHeck How it was done in ORCA? 
       //if (ix == 0) LowerBound = 0.;
       if (ix == 0 || SigmaX==0. )  // skip for surface segemnts 
     LowerBound = 0.;
       else {
     mp = MeasurementPoint( float(ix), 0.0);
     xLB = topol->localPosition(mp).x();
     gsl_sf_result result;
     int status = gsl_sf_erf_Q_e( (xLB-CloudCenterX)/SigmaX, &result);
     if (status != 0)  
       LogWarning ("Integration")<<"could not compute gaussian probability";
     LowerBound = 1-result.val;
       }
     
       if (ix == numRows-1 || SigmaX==0. ) 
     UpperBound = 1.;
       else {
     mp = MeasurementPoint( float(ix+1), 0.0);
     xUB = topol->localPosition(mp).x();
     gsl_sf_result result;
     int status = gsl_sf_erf_Q_e( (xUB-CloudCenterX)/SigmaX, &result);
     if (status != 0)  
       LogWarning ("Integration")<<"could not compute gaussian probability";
     UpperBound = 1. - result.val;
       }
       
       float   TotalIntegrationRange = UpperBound - LowerBound; // get strip
       x[ix] = TotalIntegrationRange; // save strip integral 
       //if(SigmaX==0 || SigmaY==0) 
       //cout<<TotalIntegrationRange<<" "<<ix<<endl;

     }

    // Now integarte strips in y
    int iy; // TT for compatibility
    for (iy=IPixLeftDownY; iy<=IPixRightUpY; iy++) { //loope over y ind  
      float yUB, yLB, UpperBound, LowerBound;

      if (iy == 0 || SigmaY==0.) 
    LowerBound = 0.;
      else {
        mp = MeasurementPoint( 0.0, float(iy) );
        yLB = topol->localPosition(mp).y();
    gsl_sf_result result;
    int status = gsl_sf_erf_Q_e( (yLB-CloudCenterY)/SigmaY, &result);
    if (status != 0)  
      LogWarning ("Integration")<<"could not compute gaussian probability";
    LowerBound = 1. - result.val;
      }

      if (iy == numColumns-1 || SigmaY==0. ) 
    UpperBound = 1.;
      else {
        mp = MeasurementPoint( 0.0, float(iy+1) );
        yUB = topol->localPosition(mp).y();
    gsl_sf_result result;
    int status = gsl_sf_erf_Q_e( (yUB-CloudCenterY)/SigmaY, &result);
    if (status != 0)  
      LogWarning ("Integration")<<"could not compute gaussian probability";
    UpperBound = 1. - result.val;
      }

      float   TotalIntegrationRange = UpperBound - LowerBound;
      y[iy] = TotalIntegrationRange; // save strip integral
      //if(SigmaX==0 || SigmaY==0) 
      //cout<<TotalIntegrationRange<<" "<<iy<<endl;
    }       

    // Get the 2D charge integrals by folding x and y strips
    int chan;
    for (ix=IPixLeftDownX; ix<=IPixRightUpX; ix++) {  // loop over x index
      for (iy=IPixLeftDownY; iy<=IPixRightUpY; iy++) { //loope over y ind  
    
        float ChargeFraction = Charge*x[ix]*y[iy];

        if( ChargeFraction > 0. ) {
      chan = PixelDigi::pixelToChannel( ix, iy);  // Get index 
          // Load the amplitude                      
          hit_signal[chan] += ChargeFraction;
    } // endif

    
    mp = MeasurementPoint( float(ix), float(iy) );
    LocalPoint lp = topol->localPosition(mp);
    chan = topol->channel(lp);

#ifdef TP_DEBUG
    LogDebug ("Pixel Digitizer") 
      << " pixel " << ix << " " << iy << " - "<<" "
      << chan << " " << ChargeFraction<<" "
      << mp.x() << " " << mp.y() <<" "
      << lp.x() << " " << lp.y() << " "  // givex edge position
      << chan; // edge belongs to previous ?
#endif

      } // endfor iy
    } //endfor ix
   

    // Test conversions (THIS IS FOR TESTING ONLY) comment-out.
    //     mp = topol->measurementPosition( i->position() ); //OK
    //     LocalPoint lp = topol->localPosition(mp);     //OK
    //     pair<float,float> p = topol->pixel( i->position() );  //OK
    //     chan = PixelDigi::pixelToChannel( int(p.first), int(p.second));
    //     pair<int,int> ip = PixelDigi::channelToPixel(chan);
    //     MeasurementPoint mp1 = MeasurementPoint( float(ip.first),
    //                       float(ip.second) );
    //     LogDebug ("Pixel Digitizer") << " Test "<< mp.x() << " " << mp.y() 
    //               << " "<< lp.x() << " " << lp.y() << " "<<" "
    //               <<p.first <<" "<<p.second<<" "<<chan<< " "
    //               <<" " << ip.first << " " << ip.second << " "
    //               << mp1.x() << " " << mp1.y() << " " //OK
    //               << topol->localPosition(mp1).x() << " "  //OK
    //               << topol->localPosition(mp1).y() << " "
    //               << topol->channel( i->position() ); //OK
    
    
  } // loop over charge distributions
   
  // Fill the global map with all hit pixels from this event
   
  for ( hit_map_type::const_iterator im = hit_signal.begin();
    im != hit_signal.end(); im++) {
    _signal[(*im).first] += Amplitude( (*im).second, &hit, (*im).second);
    
    int chan =  (*im).first; 
    pair<int,int> ip = PixelDigi::channelToPixel(chan);

#ifdef TP_DEBUG
    LogDebug ("Pixel Digitizer") 
      << " pixel " << ip.first << " " << ip.second << " "
      << _signal[(*im).first];    
#endif
  }

} // end induce_signal

void SiPixelDigitizerAlgorithm::init

(

const edm::EventSetup &

es

)

Definition at line 96 of file SiPixelDigitizerAlgorithm.cc.

Referenced by cms::SiPixelDigitizer::produce().

                                                           {
  if(use_ineff_from_db_){// load gain calibration service fromdb...
    theSiPixelGainCalibrationService_= new SiPixelGainCalibrationOfflineSimService(conf_);
    theSiPixelGainCalibrationService_->setESObjects( es );
  }

  fillDeadModules(es); // gets the dead module from config file or DB.
  fillLorentzAngle(es); // gets the Lorentz angle from the config file or DB.
  fillMapandGeom(es);  //gets the map and geometry from the DB (to kill ROCs)  
}

void SiPixelDigitizerAlgorithm::make_digis ( )

private

Definition at line 1065 of file SiPixelDigitizerAlgorithm.cc.

References _signal, ecalMGPA::adc(), PixelDigi::channelToPixel(), doMissCalibrate, i, internal_coll, link_coll, LogDebug, min, missCalibrate(), simi, theAdcFullScale, theElectronPerADC, thePixelThresholdInE, theThresholdInE_BPix, and theThresholdInE_FPix.

Referenced by digitize().

                                           {
  internal_coll.reserve(50); internal_coll.clear();

#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") << " make digis "<<" "
                   << " pixel threshold FPix" << theThresholdInE_FPix << " " 
                               << " pixel threshold BPix" << theThresholdInE_BPix << " "
                   << " List pixels passing threshold ";
#endif  
  
  // Loop over hit pixels

  for ( signal_map_iterator i = _signal.begin(); i != _signal.end(); i++) {
    
        float signalInElectrons = (*i).second ;   // signal in electrons

    // Do the miss calibration for calibration studies only.
    //if(doMissCalibrate) signalInElectrons = missCalibrate(signalInElectrons)

    // Do only for pixels above threshold

       if ( signalInElectrons >= thePixelThresholdInE) { // check threshold
             
      int chan =  (*i).first;  // channel number
      pair<int,int> ip = PixelDigi::channelToPixel(chan);
      int adc=0;  // ADC count as integer

      // Do the miss calibration for calibration studies only.
      if(doMissCalibrate) {
    int row = ip.first;  // X in row
    int col = ip.second; // Y is in col
    adc = int(missCalibrate(col,row,signalInElectrons)); //full misscalib.
      } else { // Just do a simple electron->adc conversion
    adc = int( signalInElectrons / theElectronPerADC ); // calibrate gain
      }
      adc = min(adc, theAdcFullScale); // Check maximum value

#ifdef TP_DEBUG
      LogDebug ("Pixel Digitizer") 
    << (*i).first << " " << (*i).second << " " << signalInElectrons 
    << " " << adc << ip.first << " " << ip.second ;
#endif
           
      // Load digis
      internal_coll.push_back( PixelDigi( ip.first, ip.second, adc));
  
      //digilink     
      if((*i).second.hits().size()>0){
    simi.clear();
    unsigned int il=0;
    for( vector<const PSimHit*>::const_iterator ihit = (*i).second.hits().begin();
         ihit != (*i).second.hits().end(); ihit++) {
      simi[(**ihit).trackId()].push_back((*i).second.individualampl()[il]);
      il++;
    }
    
    //sum the contribution of the same trackid 
    for( simlink_map::iterator simiiter=simi.begin();
         simiiter!=simi.end();
         simiiter++){
      
      float sum_samechannel=0;
      for (unsigned int iii=0;iii<(*simiiter).second.size();iii++){
        sum_samechannel+=(*simiiter).second[iii];
      }
      float fraction=sum_samechannel/(*i).second;
      if (fraction>1.) fraction=1.;
      link_coll.push_back(PixelDigiSimLink((*i).first,(*simiiter).first,((*i).second.hits().front())->eventId(),fraction));
    }
    
      }
    }
    
  }
}

std::vector<PixelDigiSimLink> SiPixelDigitizerAlgorithm::make_link ( )

inline

Definition at line 70 of file SiPixelDigitizerAlgorithm.h.

References link_coll.

Referenced by cms::SiPixelDigitizer::produce().

                                          {
    return link_coll; }

float SiPixelDigitizerAlgorithm::missCalibrate ( int  col, int  row, float  amp  ) const

private

Definition at line 1363 of file SiPixelDigitizerAlgorithm.cc.

References BPix_p0, BPix_p1, BPix_p2, BPix_p3, detID, electronsPerVCAL, electronsPerVCAL_Offset, FPix_p0, FPix_p1, FPix_p2, FPix_p3, p1, p2, p3, PixelSubdetector::PixelBarrel, and DetId::subdetId().

Referenced by make_digis().

                                                      {

  // Central values
  //const float p0=0.00352, p1=0.868, p2=112., p3=113.; // pix(0,0,0)
  //  const float p0=0.00382, p1=0.886, p2=112.7, p3=113.0; // average roc=0
  //const float p0=0.00492, p1=1.998, p2=90.6, p3=134.1; // average roc=6
  // Smeared (rms)
  //const float s0=0.00020, s1=0.051, s2=5.4, s3=4.4; // average roc=0
  //const float s0=0.00015, s1=0.043, s2=3.2, s3=3.1; // col average roc=0

  // Make 2 sets of parameters for Fpix and BPIx:

  float p0=0.0;
  float p1=0.0;
  float p2=0.0;
  float p3=0.0;

  unsigned int Sub_detid=DetId(detID).subdetId();

    if (Sub_detid == PixelSubdetector::PixelBarrel){// barrel layers
      p0 = BPix_p0;
      p1 = BPix_p1;
      p2 = BPix_p2;
      p3 = BPix_p3;
    } else {// forward disks
      p0 = FPix_p0;
      p1 = FPix_p1;
      p2 = FPix_p2;
      p3 = FPix_p3;
    }

  //  const float electronsPerVCAL = 65.5; // our present VCAL calibration (feb 2009)
  //  const float electronsPerVCAL_Offset = -414.0; // our present VCAL calibration (feb 2009)
  float newAmp = 0.; //Modified signal

  // Convert electrons to VCAL units
  float signal = (signalInElectrons-electronsPerVCAL_Offset)/electronsPerVCAL;

  // Simulate the analog response with fixed parametrization
  newAmp = p3 + p2 * tanh(p0*signal - p1);


  // Use the pixel-by-pixel calibrations
  //transform to ROC index coordinates
  //int chipIndex=0, colROC=0, rowROC=0;
  //pIndexConverter->transformToROC(col,row,chipIndex,colROC,rowROC);

  // Use calibration from a file
  //int chanROC = PixelIndices::pixelToChannelROC(rowROC,colROC); // use ROC coordinates
  //float pp0=0, pp1=0,pp2=0,pp3=0;
  //map<int,CalParameters,less<int> >::const_iterator it=calmap.find(chanROC);
  //CalParameters y  = (*it).second;
  //pp0 = y.p0;
  //pp1 = y.p1;
  //pp2 = y.p2;
  //pp3 = y.p3;

  //
  // Use random smearing 
  // Randomize the pixel response
  //float pp0  = RandGaussQ::shoot(p0,s0);
  //float pp1  = RandGaussQ::shoot(p1,s1);
  //float pp2  = RandGaussQ::shoot(p2,s2);
  //float pp3  = RandGaussQ::shoot(p3,s3);

  //newAmp = pp3 + pp2 * tanh(pp0*signal - pp1); // Final signal

  //cout<<" misscalibrate "<<col<<" "<<row<<" "<<chipIndex<<" "<<colROC<<" "
  //  <<rowROC<<" "<<signalInElectrons<<" "<<signal<<" "<<newAmp<<" "
  //  <<(signalInElectrons/theElectronPerADC)<<endl;

  return newAmp;
}  

void SiPixelDigitizerAlgorithm::module_killing_conf ( void  )

private

Definition at line 1552 of file SiPixelDigitizerAlgorithm.cc.

References _signal, PixelDigi::channelToPixel(), DeadModules, cond::rpcobgas::detid, detID, i, Module, and use_module_killing_.

Referenced by digitize().

                                                       {
  if(!use_module_killing_)
    return;
  
  bool isbad=false;
  int detid = detID;
  
  Parameters::iterator itDeadModules=DeadModules.begin();
  
  for(; itDeadModules != DeadModules.end(); ++itDeadModules){
    int Dead_detID = itDeadModules->getParameter<int>("Dead_detID"); 
    if(detid==Dead_detID){
      isbad=true;
      break;
    }
  }
  
  if(!isbad)
    return;
  
  std::string Module = itDeadModules->getParameter<std::string>("Module");
  
  if(Module=="whole"){
    for(signal_map_iterator i = _signal.begin();i != _signal.end(); i++) {    
      i->second.set(0.); // reset amplitude
    }
  }
  
  for(signal_map_iterator i = _signal.begin();i != _signal.end(); i++) {    
    pair<int,int> ip = PixelDigi::channelToPixel(i->first);//get pixel pos
    
    if(Module=="tbmA" && ip.first>=80 && ip.first<=159){
      i->second.set(0.);
    }
    
    if( Module=="tbmB" && ip.first<=79){
      i->second.set(0.);

    }

  }
}

void SiPixelDigitizerAlgorithm::module_killing_DB ( void  )

private

Definition at line 1597 of file SiPixelDigitizerAlgorithm.cc.

References _signal, PixelDigi::channelToPixel(), cond::rpcobgas::detid, detID, SiPixelQuality::disabledModuleType::errorType, i, sipixelobjects::PixelROC::idInDetUnit(), j, map_, path(), edm::ESHandle< class >::product(), SiPixelBadModule_, sipixelobjects::PixelROC::toGlobal(), and use_module_killing_.

Referenced by digitize().

                                                     {
  if(!use_module_killing_)
    return;
  
  bool isbad=false;
  uint32_t detid = detID;
  
  std::vector<SiPixelQuality::disabledModuleType>disabledModules = SiPixelBadModule_->getBadComponentList();

  SiPixelQuality::disabledModuleType badmodule;

  for (size_t id=0;id<disabledModules.size();id++)
    {
      if(detid==disabledModules[id].DetID){
    isbad=true;
        badmodule = disabledModules[id];
    break;
      }
    }
  
  if(!isbad)
    return;

  if(badmodule.errorType == 0){ // this is a whole dead module.

    for(signal_map_iterator i = _signal.begin();i != _signal.end(); i++) {
           i->second.set(0.); // reset amplitude
    }
  }
  else { // all other module types: half-modules and single ROCs.
       // Get Bad ROC position:
       //follow the example of getBadRocPositions in CondFormats/SiPixelObjects/src/SiPixelQuality.cc 
  std::vector<GlobalPixel> badrocpositions (0);
  std::pair<uint8_t, uint8_t> coord(1,1);
   for(unsigned int j = 0; j < 16; j++){
     if (SiPixelBadModule_->IsRocBad(detid, j) == true){

    std::vector<CablingPathToDetUnit> path = map_.product()->pathToDetUnit(detid);
    typedef  std::vector<CablingPathToDetUnit>::const_iterator IT;
      for  (IT it = path.begin(); it != path.end(); ++it) {
          const PixelROC* myroc = map_.product()->findItem(*it);
          if( myroc->idInDetUnit() == j) {
                      LocalPixel::RocRowCol  local = { 39, 25};   //corresponding to center of ROC row, col
              GlobalPixel global = myroc->toGlobal( LocalPixel(local) );
          badrocpositions.push_back(global);
         break;
       }
      }
     }
   }// end of getBadRocPositions

  for(signal_map_iterator i = _signal.begin();i != _signal.end(); i++) {    
    pair<int,int> ip = PixelDigi::channelToPixel(i->first);//get pixel pos

  for(std::vector<GlobalPixel>::const_iterator it = badrocpositions.begin(); it != badrocpositions.end(); ++it){
    //  std::cout<<"getBadRocPositions, detidnumber: " << detid << "row is "<<it->row<<" col is "<<it->col<<std::endl;
  if(it->row==120){
  if((ip.first>=(it->row)-40  && ip.first<=(it->row)+39) && (ip.second>=(it->col)-25 && ip.second<=(it->col)+26)){
             i->second.set(0.);
  }
  }
  else {
  if((ip.first>=(it->row)-39  && ip.first<=(it->row)+40) && (ip.second>=(it->col)-25 && ip.second<=(it->col)+26)){
             i->second.set(0.);
       }
          }
  }
  }
  }
}

void SiPixelDigitizerAlgorithm::pixel_inefficiency ( )

private

Definition at line 1246 of file SiPixelDigitizerAlgorithm.cc.

References _signal, PixelDigi::channelToPixel(), PixelIndices::DColumn(), PixelIndices::DColumnInModule(), detID, flatDistribution_, i, PXBDetId::layer(), LogDebug, numColumns, numRows, pIndexConverter, PixelSubdetector::PixelBarrel, rand(), DetId::subdetId(), thePixelChipEfficiency, thePixelColEfficiency, thePixelEfficiency, and PixelIndices::transformToROC().

Referenced by digitize().

                                                   {


  // Predefined efficiencies
  float pixelEfficiency  = 1.0;
  float columnEfficiency = 1.0;
  float chipEfficiency   = 1.0;

  // setup the chip indices conversion
  unsigned int Subid=DetId(detID).subdetId();
  if    (Subid==  PixelSubdetector::PixelBarrel){// barrel layers
    int layerIndex=PXBDetId(detID).layer();
    pixelEfficiency  = thePixelEfficiency[layerIndex-1];
    columnEfficiency = thePixelColEfficiency[layerIndex-1];
    chipEfficiency   = thePixelChipEfficiency[layerIndex-1];
 
    // This should never happen
    if(numColumns>416)  LogWarning ("Pixel Geometry") <<" wrong columns in barrel "<<numColumns;
    if(numRows>160)  LogWarning ("Pixel Geometry") <<" wrong rows in barrel "<<numRows;
    
  } else {                // forward disks
   
    // For endcaps take same for each endcap
    pixelEfficiency  = thePixelEfficiency[3];
    columnEfficiency = thePixelColEfficiency[3];
    chipEfficiency   = thePixelChipEfficiency[3];
 
    // Sometimes the forward pixels have wrong size, 
    // this crashes the index conversion, so exit.
    if(numColumns>260 || numRows>160) {
      if(numColumns>260)  LogWarning ("Pixel Geometry") <<" wrong columns in endcaps "<<numColumns;
      if(numRows>160)  LogWarning ("Pixel Geometry") <<" wrong rows in endcaps "<<numRows;
      return;
    }
  } // if barrel/forward

#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") << " enter pixel_inefficiency " << pixelEfficiency << " " 
                   << columnEfficiency << " " << chipEfficiency;
#endif
  
  // Initilize the index converter
  //PixelIndices indexConverter(numColumns,numRows);
  int chipIndex = 0;
  int rowROC = 0;
  int colROC = 0;
  map<int, int, less<int> >chips, columns;
  map<int, int, less<int> >::iterator iter;
  
  // Find out the number of columns and rocs hits
  // Loop over hit pixels, amplitude in electrons, channel = coded row,col
  for (signal_map_iterator i = _signal.begin();i != _signal.end();i++) {
    
    int chan = i->first;
    pair<int,int> ip = PixelDigi::channelToPixel(chan);
    int row = ip.first;  // X in row
    int col = ip.second; // Y is in col
    //transform to ROC index coordinates   
    pIndexConverter->transformToROC(col,row,chipIndex,colROC,rowROC);
    int dColInChip = pIndexConverter->DColumn(colROC); // get ROC dcol from ROC col 
    //dcol in mod
    int dColInDet = pIndexConverter->DColumnInModule(dColInChip,chipIndex); 
      
    chips[chipIndex]++;
    columns[dColInDet]++;
  }
  
  // Delete some ROC hits.
  for ( iter = chips.begin(); iter != chips.end() ; iter++ ) {
    //float rand  = RandFlat::shoot();
    float rand  = flatDistribution_->fire();
    if( rand > chipEfficiency ) chips[iter->first]=0;
  }

  // Delete some Dcol hits.
  for ( iter = columns.begin(); iter != columns.end() ; iter++ ) {
    //float rand  = RandFlat::shoot();
    float rand  = flatDistribution_->fire();
    if( rand > columnEfficiency ) columns[iter->first]=0;
  }
  
  // Now loop again over pixels to kill some of them.
  // Loop over hit pixels, amplitude in electrons, channel = coded row,col
  for(signal_map_iterator i = _signal.begin();i != _signal.end(); i++) {    

    //    int chan = i->first;
    pair<int,int> ip = PixelDigi::channelToPixel(i->first);//get pixel pos
    int row = ip.first;  // X in row
    int col = ip.second; // Y is in col
    //transform to ROC index coordinates   
    pIndexConverter->transformToROC(col,row,chipIndex,colROC,rowROC);
    int dColInChip = pIndexConverter->DColumn(colROC); //get ROC dcol from ROC col 
    //dcol in mod
    int dColInDet = pIndexConverter->DColumnInModule(dColInChip,chipIndex); 


    //float rand  = RandFlat::shoot();
    float rand  = flatDistribution_->fire();
    if( chips[chipIndex]==0 || columns[dColInDet]==0 
    || rand>pixelEfficiency ) {
      // make pixel amplitude =0, pixel will be lost at clusterization    
      i->second.set(0.); // reset amplitude, 
    } // end if

  } // end pixel loop
  
} // end pixel_indefficiency

void SiPixelDigitizerAlgorithm::pixel_inefficiency_db ( void  )

private

Definition at line 1528 of file SiPixelDigitizerAlgorithm.cc.

References _signal, PixelDigi::channelToPixel(), cond::rpcobgas::detid, detID, i, SiPixelGainCalibrationOfflineSimService::isDead(), theSiPixelGainCalibrationService_, and use_ineff_from_db_.

Referenced by digitize().

                                                         {
  if(!use_ineff_from_db_)
    return;
  
  // Loop over hit pixels, amplitude in electrons, channel = coded row,col
  for(signal_map_iterator i = _signal.begin();i != _signal.end(); i++) {    
    
    //    int chan = i->first;
    pair<int,int> ip = PixelDigi::channelToPixel(i->first);//get pixel pos
    int row = ip.first;  // X in row
    int col = ip.second; // Y is in col
    uint32_t detid = detID;
    //transform to ROC index coordinates   
    if(theSiPixelGainCalibrationService_->isDead(detid, col, row)){
      //      std::cout << "now in isdead check, row " << detid << " " << col << "," << row << std::endl;
      // make pixel amplitude =0, pixel will be lost at clusterization    
      i->second.set(0.); // reset amplitude, 
    } // end if
  } // end pixel loop
} // end pixel_indefficiency

void SiPixelDigitizerAlgorithm::primary_ionization ( const PSimHit &  hit )

private

Definition at line 590 of file SiPixelDigitizerAlgorithm.cc.

References _ionization_points, relval_parameters_module::energy, PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), fluctuateCharge, fluctuateEloss(), GeVperElectron, i, LogDebug, PV3DBase< T, PVType, FrameType >::mag(), NumberOfSegments, PSimHit::pabs(), PSimHit::particleType(), evf::utils::pid, point, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by digitize().

                                                                     {

// Straight line approximation for trajectory inside active media

  const float SegmentLength = 0.0010; //10microns in cm
  float energy;
  
  // Get the 3D segment direction vector 
  LocalVector direction = hit.exitPoint() - hit.entryPoint(); 
  
  float eLoss = hit.energyLoss();  // Eloss in GeV
  float length = direction.mag();  // Track length in Silicon
  
  NumberOfSegments = int ( length / SegmentLength); // Number of segments
  if(NumberOfSegments < 1) NumberOfSegments = 1;
  
#ifdef TP_DEBUG
  LogDebug ("Pixel Digitizer") 
    << " enter primary_ionzation " << NumberOfSegments 
    << " shift = " 
    << (hit.exitPoint().x()-hit.entryPoint().x()) << " " 
    << (hit.exitPoint().y()-hit.entryPoint().y()) << " " 
    << (hit.exitPoint().z()-hit.entryPoint().z()) << " "
    << hit.particleType() <<" "<< hit.pabs() ; 
#endif  

  float* elossVector = new float[NumberOfSegments];  // Eloss vector

  if( fluctuateCharge ) {
    //MP DA RIMUOVERE ASSOLUTAMENTE
    int pid = hit.particleType();
    //int pid=211;  // assume it is a pion
  
    float momentum = hit.pabs();
    // Generate fluctuated charge points
    fluctuateEloss(pid, momentum, eLoss, length, NumberOfSegments, 
           elossVector);
  }
  
  _ionization_points.resize( NumberOfSegments); // set size

  // loop over segments
  for ( int i = 0; i != NumberOfSegments; i++) {
    // Divide the segment into equal length subsegments 
    Local3DPoint point = hit.entryPoint() + 
      float((i+0.5)/NumberOfSegments) * direction;

    if( fluctuateCharge ) 
      energy = elossVector[i]/GeVperElectron; // Convert charge to elec.
    else
      energy = hit.energyLoss()/GeVperElectron/float(NumberOfSegments);
    
    EnergyDepositUnit edu( energy, point); //define position,energy point
    _ionization_points[i] = edu; // save
    
#ifdef TP_DEBUG
    LogDebug ("Pixel Digitizer") 
      << i << " " << _ionization_points[i].x() << " " 
      << _ionization_points[i].y() << " " 
      << _ionization_points[i].z() << " " 
      << _ionization_points[i].energy();
#endif
    
  }  // end for loop 

  delete[] elossVector;

}

edm::DetSet< PixelDigi >::collection_type SiPixelDigitizerAlgorithm::run	(	const std::vector< PSimHit > &	input,
		PixelGeomDetUnit *	pixdet,
		GlobalVector	bfield
	)

Definition at line 449 of file SiPixelDigitizerAlgorithm.cc.

References _bfield, _detp, _PixelHits, _signal, detID, digitize(), GeomDet::geographicalId(), collect_tpl::input, link_coll, LogDebug, and DetId::rawId().

Referenced by cms::SiPixelDigitizer::produce().

                                        {
  
  _detp = pixdet; //cache the PixGeomDetUnit
  _PixelHits=input; //cache the SimHit
  _bfield=bfield; //cache the drift direction

  // Pixel Efficiency moved from the constructor to the method run because
  // the information of the det are not available in the constructor
  // Effciency parameters. 0 - no inefficiency, 1-low lumi, 10-high lumi
  
  detID= _detp->geographicalId().rawId();
  
  _signal.clear();
  
  // initalization  of pixeldigisimlinks
  link_coll.clear();
  
  //Digitization of the SimHits of a given pixdet
  vector<PixelDigi> collector =digitize(pixdet);

  // edm::DetSet<PixelDigi> collector;
  
#ifdef TP_DEBUG
  LogDebug ("PixelDigitizer") << "[SiPixelDigitizerAlgorithm] converted " << collector.size() << " PixelDigis in DetUnit" << detID; 
#endif
  
  return collector;
}

Member Data Documentation

GlobalVector SiPixelDigitizerAlgorithm::_bfield

private

Definition at line 316 of file SiPixelDigitizerAlgorithm.h.

Referenced by DriftDirection(), and run().

std::vector<SignalPoint> SiPixelDigitizerAlgorithm::_collection_points

private

Definition at line 351 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), drift(), and induce_signal().

const PixelGeomDetUnit* SiPixelDigitizerAlgorithm::_detp

private

Definition at line 306 of file SiPixelDigitizerAlgorithm.h.

Referenced by DriftDirection(), and run().

std::vector<EnergyDepositUnit> SiPixelDigitizerAlgorithm::_ionization_points

private

Definition at line 350 of file SiPixelDigitizerAlgorithm.h.

Referenced by drift(), and primary_ionization().

std::vector<PSimHit> SiPixelDigitizerAlgorithm::_PixelHits

private

Definition at line 310 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and run().

signal_map_type SiPixelDigitizerAlgorithm::_signal

private

Definition at line 354 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), digitize(), induce_signal(), make_digis(), module_killing_conf(), module_killing_DB(), pixel_inefficiency(), pixel_inefficiency_db(), and run().

bool SiPixelDigitizerAlgorithm::addChargeVCALSmearing

private

Definition at line 289 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), and SiPixelDigitizerAlgorithm().

bool SiPixelDigitizerAlgorithm::addNoise

private

Definition at line 288 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

bool SiPixelDigitizerAlgorithm::addNoisyPixels

private

Definition at line 290 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), digitize(), and SiPixelDigitizerAlgorithm().

bool SiPixelDigitizerAlgorithm::addPixelInefficiency

private

Definition at line 292 of file SiPixelDigitizerAlgorithm.h.

bool SiPixelDigitizerAlgorithm::addThresholdSmearing

private

Definition at line 297 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

bool SiPixelDigitizerAlgorithm::alpha2Order

private

Definition at line 244 of file SiPixelDigitizerAlgorithm.h.

Referenced by drift(), DriftDirection(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::BPix_p0

private

Definition at line 281 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::BPix_p1

private

Definition at line 282 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::BPix_p2

private

Definition at line 283 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::BPix_p3

private

Definition at line 284 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

std::map<int,CalParameters,std::less<int> > SiPixelDigitizerAlgorithm::calmap

private

Definition at line 357 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::ClusterWidth

private

Definition at line 249 of file SiPixelDigitizerAlgorithm.h.

Referenced by induce_signal(), and SiPixelDigitizerAlgorithm().

edm::ParameterSet SiPixelDigitizerAlgorithm::conf_

private

Definition at line 234 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

Parameters SiPixelDigitizerAlgorithm::DeadModules

private

Definition at line 91 of file SiPixelDigitizerAlgorithm.h.

Referenced by module_killing_conf().

uint32_t SiPixelDigitizerAlgorithm::detID

private

Definition at line 307 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), DriftDirection(), missCalibrate(), module_killing_conf(), module_killing_DB(), pixel_inefficiency(), pixel_inefficiency_db(), and run().

float SiPixelDigitizerAlgorithm::Dist300

private

Definition at line 243 of file SiPixelDigitizerAlgorithm.h.

Referenced by drift(), and SiPixelDigitizerAlgorithm().

bool SiPixelDigitizerAlgorithm::doMissCalibrate

private

Definition at line 329 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), make_digis(), and SiPixelDigitizerAlgorithm().

double SiPixelDigitizerAlgorithm::electronsPerVCAL

private

Definition at line 269 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

double SiPixelDigitizerAlgorithm::electronsPerVCAL_Offset

private

Definition at line 270 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

CLHEP::RandFlat* SiPixelDigitizerAlgorithm::flatDistribution_

private

Definition at line 390 of file SiPixelDigitizerAlgorithm.h.

Referenced by pixel_inefficiency(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

SiG4UniversalFluctuation* SiPixelDigitizerAlgorithm::fluctuate

private

Definition at line 342 of file SiPixelDigitizerAlgorithm.h.

Referenced by fluctuateEloss(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

bool SiPixelDigitizerAlgorithm::fluctuateCharge

private

Definition at line 291 of file SiPixelDigitizerAlgorithm.h.

Referenced by primary_ionization(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::FPix_p0

private

Definition at line 277 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::FPix_p1

private

Definition at line 278 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::FPix_p2

private

Definition at line 279 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::FPix_p3

private

Definition at line 280 of file SiPixelDigitizerAlgorithm.h.

Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().

CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::gaussDistribution_

private

Definition at line 391 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::gaussDistributionVCALNoise_

private

Definition at line 392 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

edm::ESHandle<TrackerGeometry> SiPixelDigitizerAlgorithm::geom_

private

Definition at line 88 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::GeVperElectron

private

Definition at line 239 of file SiPixelDigitizerAlgorithm.h.

Referenced by primary_ionization(), and SiPixelDigitizerAlgorithm().

std::vector<PixelDigi> SiPixelDigitizerAlgorithm::internal_coll

private

Definition at line 313 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and make_digis().

std::vector<PixelDigiSimLink> SiPixelDigitizerAlgorithm::link_coll

private

Definition at line 315 of file SiPixelDigitizerAlgorithm.h.

Referenced by make_digis(), make_link(), and run().

edm::ESHandle<SiPixelFedCablingMap> SiPixelDigitizerAlgorithm::map_

private

Definition at line 87 of file SiPixelDigitizerAlgorithm.h.

Referenced by module_killing_DB().

float SiPixelDigitizerAlgorithm::moduleThickness

private

Definition at line 304 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and drift().

int SiPixelDigitizerAlgorithm::NumberOfSegments

private

Definition at line 237 of file SiPixelDigitizerAlgorithm.h.

Referenced by primary_ionization(), and SiPixelDigitizerAlgorithm().

int SiPixelDigitizerAlgorithm::numColumns

private

Definition at line 302 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), digitize(), induce_signal(), and pixel_inefficiency().

int SiPixelDigitizerAlgorithm::numRows

private

Definition at line 303 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), digitize(), induce_signal(), and pixel_inefficiency().

PixelIndices* SiPixelDigitizerAlgorithm::pIndexConverter

private

Definition at line 348 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and pixel_inefficiency().

float SiPixelDigitizerAlgorithm::PixelChipEff

private

Definition at line 320 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::PixelChipEfficiency

private

Definition at line 323 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::PixelColEff

private

Definition at line 319 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::PixelColEfficiency

private

Definition at line 322 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::PixelEff

private

Definition at line 318 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::PixelEfficiency

private

Definition at line 321 of file SiPixelDigitizerAlgorithm.h.

bool SiPixelDigitizerAlgorithm::pixelInefficiency

private

Definition at line 294 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and SiPixelDigitizerAlgorithm().

CLHEP::HepRandomEngine& SiPixelDigitizerAlgorithm::rndEngine

private

Definition at line 388 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::Sigma0

private

Definition at line 242 of file SiPixelDigitizerAlgorithm.h.

Referenced by drift(), and SiPixelDigitizerAlgorithm().

simlink_map SiPixelDigitizerAlgorithm::simi

private

Definition at line 353 of file SiPixelDigitizerAlgorithm.h.

Referenced by make_digis().

edm::ESHandle<SiPixelQuality> SiPixelDigitizerAlgorithm::SiPixelBadModule_

private

Definition at line 84 of file SiPixelDigitizerAlgorithm.h.

Referenced by module_killing_DB().

edm::ESHandle<SiPixelLorentzAngle> SiPixelDigitizerAlgorithm::SiPixelLorentzAngle_

private

Definition at line 81 of file SiPixelDigitizerAlgorithm.h.

Referenced by DriftDirection().

CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::smearedChargeDistribution_

private

Definition at line 399 of file SiPixelDigitizerAlgorithm.h.

CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::smearedThreshold_BPix_

private

Definition at line 397 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::smearedThreshold_FPix_

private

Definition at line 396 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::tanLorentzAnglePerTesla_BPix

private

Definition at line 275 of file SiPixelDigitizerAlgorithm.h.

Referenced by DriftDirection().

float SiPixelDigitizerAlgorithm::tanLorentzAnglePerTesla_FPix

private

Definition at line 274 of file SiPixelDigitizerAlgorithm.h.

Referenced by DriftDirection().

int SiPixelDigitizerAlgorithm::theAdcFullScale

private

Definition at line 252 of file SiPixelDigitizerAlgorithm.h.

Referenced by make_digis(), and SiPixelDigitizerAlgorithm().

int SiPixelDigitizerAlgorithm::theColsInChip

private

Definition at line 299 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::theElectronPerADC

private

Definition at line 251 of file SiPixelDigitizerAlgorithm.h.

Referenced by make_digis(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::theGainSmearing

private

Definition at line 330 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::theGaussianDistribution

private

Definition at line 402 of file SiPixelDigitizerAlgorithm.h.

float SiPixelDigitizerAlgorithm::theNoiseInElectrons

private

Definition at line 253 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), digitize(), and SiPixelDigitizerAlgorithm().

GaussianTailNoiseGenerator* SiPixelDigitizerAlgorithm::theNoiser

private

Definition at line 343 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::theOffsetSmearing

private

Definition at line 331 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::thePixelChipEfficiency[6]

private

Definition at line 326 of file SiPixelDigitizerAlgorithm.h.

Referenced by pixel_inefficiency(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::thePixelColEfficiency[6]

private

Definition at line 325 of file SiPixelDigitizerAlgorithm.h.

Referenced by pixel_inefficiency(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::thePixelEfficiency[6]

private

Definition at line 324 of file SiPixelDigitizerAlgorithm.h.

Referenced by pixel_inefficiency(), and SiPixelDigitizerAlgorithm().

int SiPixelDigitizerAlgorithm::thePixelLuminosity

private

Definition at line 295 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::thePixelThreshold

private

Definition at line 259 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), and digitize().

float SiPixelDigitizerAlgorithm::thePixelThresholdInE

private

Definition at line 261 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and make_digis().

float SiPixelDigitizerAlgorithm::theReadoutNoise

private

Definition at line 254 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

int SiPixelDigitizerAlgorithm::theRowsInChip

private

Definition at line 300 of file SiPixelDigitizerAlgorithm.h.

SiPixelGainCalibrationOfflineSimService* SiPixelDigitizerAlgorithm::theSiPixelGainCalibrationService_

private

Definition at line 380 of file SiPixelDigitizerAlgorithm.h.

Referenced by pixel_inefficiency_db(), and ~SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::theSmearedChargeRMS

private

Definition at line 257 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise().

float SiPixelDigitizerAlgorithm::theThresholdInE_BPix

private

Definition at line 264 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), digitize(), make_digis(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::theThresholdInE_FPix

private

Definition at line 263 of file SiPixelDigitizerAlgorithm.h.

Referenced by add_noise(), digitize(), make_digis(), and SiPixelDigitizerAlgorithm().

double SiPixelDigitizerAlgorithm::theThresholdSmearing_BPix

private

Definition at line 267 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

double SiPixelDigitizerAlgorithm::theThresholdSmearing_FPix

private

Definition at line 266 of file SiPixelDigitizerAlgorithm.h.

Referenced by SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::theTofLowerCut

private

Definition at line 272 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and SiPixelDigitizerAlgorithm().

float SiPixelDigitizerAlgorithm::theTofUpperCut

private

Definition at line 273 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and SiPixelDigitizerAlgorithm().

double SiPixelDigitizerAlgorithm::tMax

private

Definition at line 339 of file SiPixelDigitizerAlgorithm.h.

Referenced by fluctuateEloss(), and SiPixelDigitizerAlgorithm().

const PixelTopology* SiPixelDigitizerAlgorithm::topol

private

Definition at line 311 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and induce_signal().

bool SiPixelDigitizerAlgorithm::use_deadmodule_DB_

private

Definition at line 375 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize().

bool SiPixelDigitizerAlgorithm::use_ineff_from_db_

private

Definition at line 372 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), and pixel_inefficiency_db().

bool SiPixelDigitizerAlgorithm::use_LorentzAngle_DB_

private

Definition at line 376 of file SiPixelDigitizerAlgorithm.h.

Referenced by DriftDirection().

bool SiPixelDigitizerAlgorithm::use_module_killing_

private

Definition at line 374 of file SiPixelDigitizerAlgorithm.h.

Referenced by digitize(), module_killing_conf(), and module_killing_DB().