|
|
|
#include <SiPixelDigitizerAlgorithm.h>
Definition at line 59 of file SiPixelDigitizerAlgorithm.h.
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.
| SiPixelDigitizerAlgorithm::SiPixelDigitizerAlgorithm | ( | const edm::ParameterSet & | conf, |
| CLHEP::HepRandomEngine & | eng | ||
| ) |
Definition at line 141 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, AlCaHLTBitMon_ParallelJobs::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 if (thePixelLuminosity==-20){ thePixelColEfficiency[0] = conf_.getParameter<double>("thePixelColEfficiency_BPix1"); thePixelColEfficiency[1] = conf_.getParameter<double>("thePixelColEfficiency_BPix2"); thePixelColEfficiency[2] = conf_.getParameter<double>("thePixelColEfficiency_BPix3"); thePixelColEfficiency[3] = conf_.getParameter<double>("thePixelColEfficiency_FPix1"); thePixelColEfficiency[4] = conf_.getParameter<double>("thePixelColEfficiency_FPix2"); // Not used, but leave it in in case we want use it to later cout<<"\nReading in custom Pixel efficiencies "<<thePixelColEfficiency[0]<<" , "<<thePixelColEfficiency[1]<<" , " <<thePixelColEfficiency[2]<<" , "<<thePixelColEfficiency[3]<<" , "<<thePixelColEfficiency[4]<<"\n"; if (thePixelColEfficiency[0]<=0.5) {cout <<"\n\nDid you mean to set the Pixel efficiency at "<<thePixelColEfficiency[0] <<", or did you mean for this to be the inefficiency?\n\n\n";} } // 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 437 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;
}
| void SiPixelDigitizerAlgorithm::add_noise | ( | ) | [private] |
Definition at line 1159 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 494 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 729 of file SiPixelDigitizerAlgorithm.cc.
References _collection_points, _ionization_points, alpha2Order, PSimHit::detUnitId(), Dist300, DriftDirection(), i, LogDebug, lorentzAngle, moduleThickness, FitTarget::Sigma, Sigma0, SiPixelLorentzAngle_, 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
//float LorentzAng = SiPixelLorentzAngle_->getLorentzAngle(hit.detUnitId());
lorentzAngle = SiPixelLorentzAngle_->getLorentzAngle(hit.detUnitId());
//LocalVector driftDir=DriftDirection(LorentzAng); // get the charge drift direction
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 1461 of file SiPixelDigitizerAlgorithm.cc.
References _bfield, _detp, alpha2Order, detID, LogDebug, lorentzAngle, PixelSubdetector::PixelBarrel, GloballyPositioned< T >::position(), GloballyPositioned< T >::rotation(), 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_){
alpha2 = lorentzAngle * lorentzAngle;
//std::cout << "detID is: "<< it->first <<"The LA per tesla is: "<< it->second << std::endl;
dir_x = -( lorentzAngle * Bfield.y() + alpha2 * Bfield.z()* Bfield.x() );
dir_y = +( lorentzAngle * 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 110 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 120 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 134 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 674 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 842 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, PV3DBase< T, PVType, FrameType >::x(), PV2DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), PV2DBase< T, PVType, FrameType >::y(), and detailsBasic3DVector::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 97 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 1080 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 1378 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 1554 of file SiPixelDigitizerAlgorithm.cc.
References _signal, PixelDigi::channelToPixel(), DeadModules, detID, cond::rpcobgas::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 1596 of file SiPixelDigitizerAlgorithm.cc.
References _signal, PixelDigi::channelToPixel(), detID, cond::rpcobgas::detid, SiPixelQuality::disabledModuleType::errorType, i, sipixelobjects::PixelROC::idInDetUnit(), j, map_, path(), edm::ESHandle< T >::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;
//cout<<"Hit in: "<< detid <<" errorType "<< badmodule.errorType<<" BadRocs="<<hex<<SiPixelBadModule_->getBadRocs(detid)<<dec<<" "<<endl;
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);
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){
if(it->row >= 80 && ip.first >= 80 ){
if((fabs(ip.second - it->col) < 26) ) {i->second.set(0.);}
else if(it->row==120 && ip.second-it->col==26){i->second.set(0.);}
else if(it->row==119 && it->col-ip.second==26){i->second.set(0.);}
}
else if(it->row < 80 && ip.first < 80 ){
if((fabs(ip.second - it->col) < 26) ){i->second.set(0.);}
else if(it->row==40 && ip.second-it->col==26){i->second.set(0.);}
else if(it->row==39 && it->col-ip.second==26){i->second.set(0.);}
}
}
}
}
}
| void SiPixelDigitizerAlgorithm::pixel_inefficiency | ( | ) | [private] |
Definition at line 1261 of file SiPixelDigitizerAlgorithm.cc.
References _signal, PixelDigi::channelToPixel(), PixelIndices::DColumn(), PixelIndices::DColumnInModule(), detID, flatDistribution_, i, PXBDetId::layer(), LogDebug, numColumns, numRows, pIndexConverter, PixelSubdetector::PixelBarrel, 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 1530 of file SiPixelDigitizerAlgorithm.cc.
References _signal, PixelDigi::channelToPixel(), detID, cond::rpcobgas::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 602 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 461 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;
}
Definition at line 317 of file SiPixelDigitizerAlgorithm.h.
Referenced by DriftDirection(), and run().
std::vector<SignalPoint> SiPixelDigitizerAlgorithm::_collection_points [private] |
Definition at line 352 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), drift(), and induce_signal().
const PixelGeomDetUnit* SiPixelDigitizerAlgorithm::_detp [private] |
Definition at line 307 of file SiPixelDigitizerAlgorithm.h.
Referenced by DriftDirection(), and run().
std::vector<EnergyDepositUnit> SiPixelDigitizerAlgorithm::_ionization_points [private] |
Definition at line 351 of file SiPixelDigitizerAlgorithm.h.
Referenced by drift(), and primary_ionization().
std::vector<PSimHit> SiPixelDigitizerAlgorithm::_PixelHits [private] |
Definition at line 311 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), and run().
Definition at line 355 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 290 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), and SiPixelDigitizerAlgorithm().
bool SiPixelDigitizerAlgorithm::addNoise [private] |
Definition at line 289 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
bool SiPixelDigitizerAlgorithm::addNoisyPixels [private] |
Definition at line 291 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), digitize(), and SiPixelDigitizerAlgorithm().
bool SiPixelDigitizerAlgorithm::addPixelInefficiency [private] |
Definition at line 293 of file SiPixelDigitizerAlgorithm.h.
bool SiPixelDigitizerAlgorithm::addThresholdSmearing [private] |
Definition at line 298 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 282 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::BPix_p1 [private] |
Definition at line 283 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::BPix_p2 [private] |
Definition at line 284 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::BPix_p3 [private] |
Definition at line 285 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
std::map<int,CalParameters,std::less<int> > SiPixelDigitizerAlgorithm::calmap [private] |
Definition at line 358 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().
Definition at line 234 of file SiPixelDigitizerAlgorithm.h.
Referenced by SiPixelDigitizerAlgorithm().
Definition at line 91 of file SiPixelDigitizerAlgorithm.h.
Referenced by module_killing_conf().
uint32_t SiPixelDigitizerAlgorithm::detID [private] |
Definition at line 308 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 330 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 391 of file SiPixelDigitizerAlgorithm.h.
Referenced by pixel_inefficiency(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
Definition at line 343 of file SiPixelDigitizerAlgorithm.h.
Referenced by fluctuateEloss(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
bool SiPixelDigitizerAlgorithm::fluctuateCharge [private] |
Definition at line 292 of file SiPixelDigitizerAlgorithm.h.
Referenced by primary_ionization(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::FPix_p0 [private] |
Definition at line 278 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::FPix_p1 [private] |
Definition at line 279 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::FPix_p2 [private] |
Definition at line 280 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::FPix_p3 [private] |
Definition at line 281 of file SiPixelDigitizerAlgorithm.h.
Referenced by missCalibrate(), and SiPixelDigitizerAlgorithm().
CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::gaussDistribution_ [private] |
Definition at line 392 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::gaussDistributionVCALNoise_ [private] |
Definition at line 393 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
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 314 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), and make_digis().
std::vector<PixelDigiSimLink> SiPixelDigitizerAlgorithm::link_coll [private] |
Definition at line 316 of file SiPixelDigitizerAlgorithm.h.
Referenced by make_digis(), make_link(), and run().
float SiPixelDigitizerAlgorithm::lorentzAngle [private] |
Definition at line 276 of file SiPixelDigitizerAlgorithm.h.
Referenced by drift(), and DriftDirection().
Definition at line 87 of file SiPixelDigitizerAlgorithm.h.
Referenced by module_killing_DB().
float SiPixelDigitizerAlgorithm::moduleThickness [private] |
Definition at line 305 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 303 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), digitize(), induce_signal(), and pixel_inefficiency().
int SiPixelDigitizerAlgorithm::numRows [private] |
Definition at line 304 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), digitize(), induce_signal(), and pixel_inefficiency().
Definition at line 349 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), and pixel_inefficiency().
float SiPixelDigitizerAlgorithm::PixelChipEff [private] |
Definition at line 321 of file SiPixelDigitizerAlgorithm.h.
float SiPixelDigitizerAlgorithm::PixelChipEfficiency [private] |
Definition at line 324 of file SiPixelDigitizerAlgorithm.h.
float SiPixelDigitizerAlgorithm::PixelColEff [private] |
Definition at line 320 of file SiPixelDigitizerAlgorithm.h.
float SiPixelDigitizerAlgorithm::PixelColEfficiency [private] |
Definition at line 323 of file SiPixelDigitizerAlgorithm.h.
float SiPixelDigitizerAlgorithm::PixelEff [private] |
Definition at line 319 of file SiPixelDigitizerAlgorithm.h.
float SiPixelDigitizerAlgorithm::PixelEfficiency [private] |
Definition at line 322 of file SiPixelDigitizerAlgorithm.h.
bool SiPixelDigitizerAlgorithm::pixelInefficiency [private] |
Definition at line 295 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), and SiPixelDigitizerAlgorithm().
CLHEP::HepRandomEngine& SiPixelDigitizerAlgorithm::rndEngine [private] |
Definition at line 389 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 354 of file SiPixelDigitizerAlgorithm.h.
Referenced by make_digis().
Definition at line 84 of file SiPixelDigitizerAlgorithm.h.
Referenced by module_killing_DB().
Definition at line 81 of file SiPixelDigitizerAlgorithm.h.
Referenced by drift().
CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::smearedChargeDistribution_ [private] |
Definition at line 400 of file SiPixelDigitizerAlgorithm.h.
CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::smearedThreshold_BPix_ [private] |
Definition at line 398 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::smearedThreshold_FPix_ [private] |
Definition at line 397 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 300 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 331 of file SiPixelDigitizerAlgorithm.h.
Referenced by SiPixelDigitizerAlgorithm().
CLHEP::RandGaussQ* SiPixelDigitizerAlgorithm::theGaussianDistribution [private] |
Definition at line 403 of file SiPixelDigitizerAlgorithm.h.
float SiPixelDigitizerAlgorithm::theNoiseInElectrons [private] |
Definition at line 253 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), digitize(), and SiPixelDigitizerAlgorithm().
Definition at line 344 of file SiPixelDigitizerAlgorithm.h.
Referenced by add_noise(), SiPixelDigitizerAlgorithm(), and ~SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::theOffsetSmearing [private] |
Definition at line 332 of file SiPixelDigitizerAlgorithm.h.
Referenced by SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::thePixelChipEfficiency[6] [private] |
Definition at line 327 of file SiPixelDigitizerAlgorithm.h.
Referenced by pixel_inefficiency(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::thePixelColEfficiency[6] [private] |
Definition at line 326 of file SiPixelDigitizerAlgorithm.h.
Referenced by pixel_inefficiency(), and SiPixelDigitizerAlgorithm().
float SiPixelDigitizerAlgorithm::thePixelEfficiency[6] [private] |
Definition at line 325 of file SiPixelDigitizerAlgorithm.h.
Referenced by pixel_inefficiency(), and SiPixelDigitizerAlgorithm().
int SiPixelDigitizerAlgorithm::thePixelLuminosity [private] |
Definition at line 296 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 301 of file SiPixelDigitizerAlgorithm.h.
SiPixelGainCalibrationOfflineSimService* SiPixelDigitizerAlgorithm::theSiPixelGainCalibrationService_ [private] |
Definition at line 381 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 340 of file SiPixelDigitizerAlgorithm.h.
Referenced by fluctuateEloss(), and SiPixelDigitizerAlgorithm().
const PixelTopology* SiPixelDigitizerAlgorithm::topol [private] |
Definition at line 312 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), and induce_signal().
bool SiPixelDigitizerAlgorithm::use_deadmodule_DB_ [private] |
Definition at line 376 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize().
bool SiPixelDigitizerAlgorithm::use_ineff_from_db_ [private] |
Definition at line 373 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), and pixel_inefficiency_db().
bool SiPixelDigitizerAlgorithm::use_LorentzAngle_DB_ [private] |
Definition at line 377 of file SiPixelDigitizerAlgorithm.h.
Referenced by DriftDirection().
bool SiPixelDigitizerAlgorithm::use_module_killing_ [private] |
Definition at line 375 of file SiPixelDigitizerAlgorithm.h.
Referenced by digitize(), module_killing_conf(), and module_killing_DB().
1.7.3