|
|
|
#include <ShallowGainCalibration.h>
Public Member Functions | |
| ShallowGainCalibration (const edm::ParameterSet &) | |
Private Member Functions | |
| bool | IsFarFromBorder (TrajectoryStateOnSurface *trajState, const uint32_t detid, const edm::EventSetup *iSetup) |
| void | produce (edm::Event &, const edm::EventSetup &) |
| double | thickness (DetId id) |
Private Attributes | |
| std::map< DetId, double > | m_thicknessMap |
| const TrackerGeometry * | m_tracker |
| std::string | Prefix |
| std::string | Suffix |
| edm::InputTag | theTracksLabel |
Definition at line 63 of file ShallowGainCalibration.h.
| ShallowGainCalibration::ShallowGainCalibration | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 12 of file ShallowGainCalibration.cc.
References Prefix, and Suffix.
: theTracksLabel( iConfig.getParameter<edm::InputTag>("Tracks") ), Suffix ( iConfig.getParameter<std::string>("Suffix") ), Prefix ( iConfig.getParameter<std::string>("Prefix") ) { produces <std::vector<int> > ( Prefix + "trackindex" + Suffix ); produces <std::vector<unsigned int> > ( Prefix + "rawid" + Suffix ); produces <std::vector<float> > ( Prefix + "localdirx" + Suffix ); produces <std::vector<float> > ( Prefix + "localdiry" + Suffix ); produces <std::vector<float> > ( Prefix + "localdirz" + Suffix ); produces <std::vector<unsigned short> > ( Prefix + "firststrip" + Suffix ); produces <std::vector<unsigned short> > ( Prefix + "nstrips" + Suffix ); produces <std::vector<bool> > ( Prefix + "saturation" + Suffix ); produces <std::vector<bool> > ( Prefix + "overlapping" + Suffix ); produces <std::vector<bool> > ( Prefix + "farfromedge" + Suffix ); produces <std::vector<unsigned int> > ( Prefix + "charge" + Suffix ); produces <std::vector<float> > ( Prefix + "path" + Suffix ); produces <std::vector<float> > ( Prefix + "chargeoverpath" + Suffix ); produces <std::vector<unsigned char> > ( Prefix + "amplitude" + Suffix ); produces <std::vector<double> > ( Prefix + "gainused" + Suffix ); }
| bool ShallowGainCalibration::IsFarFromBorder | ( | TrajectoryStateOnSurface * | trajState, |
| const uint32_t | detid, | ||
| const edm::EventSetup * | iSetup | ||
| ) | [private] |
Definition at line 227 of file ShallowGainCalibration.cc.
References BoundSurface::bounds(), gather_cfg::cout, edm::EventSetup::get(), Bounds::length(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), Parameters::parameters, LocalTrajectoryError::positionError(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().
Referenced by produce().
{
edm::ESHandle<TrackerGeometry> tkGeom; iSetup->get<TrackerDigiGeometryRecord>().get( tkGeom );
LocalPoint HitLocalPos = trajState->localPosition();
LocalError HitLocalError = trajState->localError().positionError() ;
const GeomDetUnit* it = tkGeom->idToDetUnit(DetId(detid));
if (dynamic_cast<const StripGeomDetUnit*>(it)==0 && dynamic_cast<const PixelGeomDetUnit*>(it)==0) {
std::cout << "this detID doesn't seem to belong to the Tracker" << std::endl;
return false;
}
const BoundPlane plane = it->surface();
const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));
double DistFromBorder = 1.0;
double HalfLength = it->surface().bounds().length() /2.0;
if(trapezoidalBounds)
{
std::vector<float> const & parameters = (*trapezoidalBounds).parameters();
HalfLength = parameters[3];
}else if(rectangularBounds){
HalfLength = it->surface().bounds().length() /2.0;
}else{return false;}
if (fabs(HitLocalPos.y())+HitLocalError.yy() >= (HalfLength - DistFromBorder) ) return false;
return true;
}
| void ShallowGainCalibration::produce | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDProducer.
Definition at line 35 of file ShallowGainCalibration.cc.
References a, SiStripCluster::amplitudes(), DeDxDiscriminatorTools::charge(), shallow::findTrackIndex(), SiStripCluster::firstStrip(), SiStripCluster::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), h, IsFarFromBorder(), TrajectoryStateOnSurface::isValid(), edm::ESHandleBase::isValid(), TrajectoryStateOnSurface::localDirection(), m_tracker, PV3DBase< T, PVType, FrameType >::mag(), Trajectory::measurements(), NULL, scaleCards::path, Prefix, edm::Event::put(), Suffix, theTracksLabel, thickness(), testEve_cfg::tracks, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
std::auto_ptr<std::vector<int> > trackindex ( new std::vector<int> );
std::auto_ptr<std::vector<unsigned int> > rawid ( new std::vector<unsigned int> );
std::auto_ptr<std::vector<float> > localdirx ( new std::vector<float> );
std::auto_ptr<std::vector<float> > localdiry ( new std::vector<float> );
std::auto_ptr<std::vector<float> > localdirz ( new std::vector<float> );
std::auto_ptr<std::vector<unsigned short> > firststrip ( new std::vector<unsigned short> );
std::auto_ptr<std::vector<unsigned short> > nstrips ( new std::vector<unsigned short> );
std::auto_ptr<std::vector<bool> > saturation ( new std::vector<bool> );
std::auto_ptr<std::vector<bool> > overlapping ( new std::vector<bool> );
std::auto_ptr<std::vector<bool> > farfromedge ( new std::vector<bool> );
std::auto_ptr<std::vector<unsigned int> > charge ( new std::vector<unsigned int> );
std::auto_ptr<std::vector<float> > path ( new std::vector<float> );
std::auto_ptr<std::vector<float> > chargeoverpath( new std::vector<float> );
std::auto_ptr<std::vector<unsigned char> > amplitude ( new std::vector<unsigned char> );
std::auto_ptr<std::vector<double> > gainused ( new std::vector<double> );
edm::ESHandle<TrackerGeometry> theTrackerGeometry; iSetup.get<TrackerDigiGeometryRecord>().get( theTrackerGeometry );
m_tracker=&(* theTrackerGeometry );
edm::ESHandle<SiStripGain> gainHandle; iSetup.get<SiStripGainRcd>().get(gainHandle);
edm::Handle<edm::View<reco::Track> > tracks; iEvent.getByLabel(theTracksLabel, tracks);
edm::Handle<TrajTrackAssociationCollection> associations; iEvent.getByLabel(theTracksLabel, associations);
for( TrajTrackAssociationCollection::const_iterator association = associations->begin(); association != associations->end(); association++) {
const Trajectory* traj = association->key.get();
const reco::Track* track = association->val.get();
vector<TrajectoryMeasurement> measurements = traj->measurements();
for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
TrajectoryStateOnSurface trajState = measurement_it->updatedState();
if( !trajState.isValid() ) continue;
const TrackingRecHit* hit = (*measurement_it->recHit()).hit();
const SiStripRecHit1D* sistripsimple1dhit = dynamic_cast<const SiStripRecHit1D*>(hit);
const SiStripRecHit2D* sistripsimplehit = dynamic_cast<const SiStripRecHit2D*>(hit);
const SiStripMatchedRecHit2D* sistripmatchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
const SiStripCluster* Cluster = NULL;
uint32_t DetId = 0;
for(unsigned int h=0;h<2;h++){
if(!sistripmatchedhit && h==1){
continue;
}else if(sistripmatchedhit && h==0){
Cluster = &sistripmatchedhit->monoCluster();
DetId = sistripmatchedhit->monoId();
}else if(sistripmatchedhit && h==1){
Cluster = &sistripmatchedhit->stereoCluster();;
DetId = sistripmatchedhit->stereoId();
}else if(sistripsimplehit){
Cluster = (sistripsimplehit->cluster()).get();
DetId = sistripsimplehit->geographicalId().rawId();
}else if(sistripsimple1dhit){
Cluster = (sistripsimple1dhit->cluster()).get();
DetId = sistripsimple1dhit->geographicalId().rawId();
}else{
continue;
}
LocalVector trackDirection = trajState.localDirection();
double cosine = trackDirection.z()/trackDirection.mag();
const vector<uint8_t>& Ampls = Cluster->amplitudes();
int FirstStrip = Cluster->firstStrip();
int APVId = FirstStrip/128;
bool Saturation = false;
bool Overlapping = false;
unsigned int Charge = 0;
double Path = (10.0*thickness(DetId))/fabs(cosine);
double PrevGain = -1;
if(gainHandle.isValid()){
SiStripApvGain::Range detGainRange = gainHandle->getRange(DetId);
PrevGain = *(detGainRange.first + APVId);
}
for(unsigned int a=0;a<Ampls.size();a++){
Charge+=Ampls[a];
if(Ampls[a] >=254)Saturation =true;
amplitude->push_back( Ampls[a] );
}
double ChargeOverPath = (double)Charge / Path ;
if(FirstStrip==0 )Overlapping=true;
if(FirstStrip==128 )Overlapping=true;
if(FirstStrip==256 )Overlapping=true;
if(FirstStrip==384 )Overlapping=true;
if(FirstStrip==512 )Overlapping=true;
if(FirstStrip==640 )Overlapping=true;
if(FirstStrip<=127 && FirstStrip+Ampls.size()>127)Overlapping=true;
if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlapping=true;
if(FirstStrip<=383 && FirstStrip+Ampls.size()>383)Overlapping=true;
if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlapping=true;
if(FirstStrip<=639 && FirstStrip+Ampls.size()>639)Overlapping=true;
if(FirstStrip+Ampls.size()==127 )Overlapping=true;
if(FirstStrip+Ampls.size()==255 )Overlapping=true;
if(FirstStrip+Ampls.size()==383 )Overlapping=true;
if(FirstStrip+Ampls.size()==511 )Overlapping=true;
if(FirstStrip+Ampls.size()==639 )Overlapping=true;
if(FirstStrip+Ampls.size()==767 )Overlapping=true;
trackindex ->push_back( shallow::findTrackIndex(tracks, track) );
rawid ->push_back( DetId );
localdirx ->push_back( trackDirection.x() );
localdiry ->push_back( trackDirection.y() );
localdirz ->push_back( trackDirection.z() );
firststrip ->push_back( FirstStrip );
nstrips ->push_back( Ampls.size() );
saturation ->push_back( Saturation );
overlapping ->push_back( Overlapping );
farfromedge ->push_back( IsFarFromBorder(&trajState,DetId, &iSetup) );
charge ->push_back( Charge );
path ->push_back( Path );
chargeoverpath->push_back( ChargeOverPath );
gainused ->push_back( PrevGain );
}
}
}
iEvent.put(trackindex, Prefix + "trackindex" + Suffix );
iEvent.put(rawid , Prefix + "rawid" + Suffix );
iEvent.put(localdirx , Prefix + "localdirx" + Suffix );
iEvent.put(localdiry , Prefix + "localdiry" + Suffix );
iEvent.put(localdirz , Prefix + "localdirz" + Suffix );
iEvent.put(firststrip, Prefix + "firststrip" + Suffix );
iEvent.put(nstrips, Prefix + "nstrips" + Suffix );
iEvent.put(saturation, Prefix + "saturation" + Suffix );
iEvent.put(overlapping, Prefix + "overlapping" + Suffix );
iEvent.put(farfromedge, Prefix + "farfromedge" + Suffix );
iEvent.put(charge, Prefix + "charge" + Suffix );
iEvent.put(path, Prefix + "path" + Suffix );
iEvent.put(chargeoverpath,Prefix + "chargeoverpath"+ Suffix );
iEvent.put(amplitude, Prefix + "amplitude" + Suffix );
iEvent.put(gainused, Prefix + "gainused" + Suffix );
}
| double ShallowGainCalibration::thickness | ( | DetId | id | ) | [private] |
Definition at line 261 of file ShallowGainCalibration.cc.
References errorMatrix2Lands_multiChannel::id, TrackerGeometry::idToDetUnit(), m_thicknessMap, and m_tracker.
Referenced by produce().
{
map<DetId,double>::iterator th=m_thicknessMap.find(id);
if(th!=m_thicknessMap.end())
return (*th).second;
else {
double detThickness=1.;
//compute thickness normalization
const GeomDetUnit* it = m_tracker->idToDetUnit(DetId(id));
bool isPixel = dynamic_cast<const PixelGeomDetUnit*>(it)!=0;
bool isStrip = dynamic_cast<const StripGeomDetUnit*>(it)!=0;
if (!isPixel && ! isStrip) {
//FIXME throw exception
edm::LogWarning("DeDxHitsProducer") << "\t\t this detID doesn't seem to belong to the Tracker";
detThickness = 1.;
}else{
detThickness = it->surface().bounds().thickness();
}
m_thicknessMap[id]=detThickness;//computed value
return detThickness;
}
}
std::map<DetId,double> ShallowGainCalibration::m_thicknessMap [private] |
Definition at line 78 of file ShallowGainCalibration.h.
Referenced by thickness().
const TrackerGeometry* ShallowGainCalibration::m_tracker [private] |
Definition at line 77 of file ShallowGainCalibration.h.
Referenced by produce(), and thickness().
std::string ShallowGainCalibration::Prefix [private] |
Definition at line 69 of file ShallowGainCalibration.h.
Referenced by produce(), and ShallowGainCalibration().
std::string ShallowGainCalibration::Suffix [private] |
Definition at line 68 of file ShallowGainCalibration.h.
Referenced by produce(), and ShallowGainCalibration().
Definition at line 67 of file ShallowGainCalibration.h.
Referenced by produce().
1.7.3