|
|
|
Classes | |
| struct | RawHits |
Functions | |
| double | genericAverage (const reco::DeDxHitCollection &, float expo=1.) |
| bool | shapeSelection (const std::vector< uint8_t > &ls) |
| void | trajectoryRawHits (const edm::Ref< std::vector< Trajectory > > &trajectory, vector< RawHits > &hits, bool usePixel, bool useStrip) |
| void | trajectoryRawHits (const edm::Ref< std::vector< Trajectory > > &trajectory, std::vector< RawHits > &hits, bool usePixel, bool useStrip) |
| double DeDxTools::genericAverage | ( | const reco::DeDxHitCollection & | hits, |
| float | expo = 1. |
||
| ) |
Definition at line 129 of file DeDxTools.cc.
References DeDxDiscriminatorTools::charge(), i, n, funct::pow(), and query::result.
Referenced by GenericAverageDeDxEstimator::dedx().
| bool DeDxTools::shapeSelection | ( | const std::vector< uint8_t > & | ampls | ) |
Definition at line 141 of file DeDxTools.cc.
References i.
Referenced by DeDxEstimatorProducer::produce(), and DeDxDiscriminatorProducer::produce().
{
// ---------------- COMPTAGE DU NOMBRE DE MAXIMA --------------------------
//----------------------------------------------------------------------------
// printf("ShapeTest \n");
Int_t NofMax=0; Int_t recur255=1; Int_t recur254=1;
bool MaxOnStart=false;bool MaxInMiddle=false, MaxOnEnd =false;
Int_t MaxPos=0;
// Début avec max
if(ampls.size()!=1 && ((ampls[0]>ampls[1])
|| (ampls.size()>2 && ampls[0]==ampls[1] && ampls[1]>ampls[2] && ampls[0]!=254 && ampls[0]!=255)
|| (ampls.size()==2 && ampls[0]==ampls[1] && ampls[0]!=254 && ampls[0]!=255)) ){
NofMax=NofMax+1; MaxOnStart=true; }
// Maximum entouré
if(ampls.size()>2){
for (unsigned int i =1; i < ampls.size()-1; i++) {
if( (ampls[i]>ampls[i-1] && ampls[i]>ampls[i+1])
|| (ampls.size()>3 && i>0 && i<ampls.size()-2 && ampls[i]==ampls[i+1] && ampls[i]>ampls[i-1] && ampls[i]>ampls[i+2] && ampls[i]!=254 && ampls[i]!=255) ){
NofMax=NofMax+1; MaxInMiddle=true; MaxPos=i;
}
if(ampls[i]==255 && ampls[i]==ampls[i-1]) {
recur255=recur255+1;
MaxPos=i-(recur255/2);
if(ampls[i]>ampls[i+1]){NofMax=NofMax+1;MaxInMiddle=true;}
}
if(ampls[i]==254 && ampls[i]==ampls[i-1]) {
recur254=recur254+1;
MaxPos=i-(recur254/2);
if(ampls[i]>ampls[i+1]){NofMax=NofMax+1;MaxInMiddle=true;}
}
}
}
// Fin avec un max
if(ampls.size()>1){
if(ampls[ampls.size()-1]>ampls[ampls.size()-2]
|| (ampls.size()>2 && ampls[ampls.size()-1]==ampls[ampls.size()-2] && ampls[ampls.size()-2]>ampls[ampls.size()-3] )
|| ampls[ampls.size()-1]==255){
NofMax=NofMax+1; MaxOnEnd=true; }
}
// Si une seule strip touchée
if(ampls.size()==1){ NofMax=1;}
// --- SELECTION EN FONCTION DE LA FORME POUR LES MAXIMA UNIQUES ---------
//------------------------------------------------------------------------
/*
____
| |____
____| | |
| | | |____
____| | | | |
| | | | | |____
__|____|____|____|____|____|____|__
C_Mnn C_Mn C_M C_D C_Dn C_Dnn
*/
// bool shapetest=true;
bool shapecdtn=false;
// Float_t C_M; Float_t C_D; Float_t C_Mn; Float_t C_Dn; Float_t C_Mnn; Float_t C_Dnn;
Float_t C_M=0.0; Float_t C_D=0.0; Float_t C_Mn=10000; Float_t C_Dn=10000; Float_t C_Mnn=10000; Float_t C_Dnn=10000;
Int_t CDPos;
Float_t coeff1=1.7; Float_t coeff2=2.0;
Float_t coeffn=0.10; Float_t coeffnn=0.02; Float_t noise=4.0;
if(NofMax==1){
if(MaxOnStart==true){
C_M=(Float_t)ampls[0]; C_D=(Float_t)ampls[1];
if(ampls.size()<3) shapecdtn=true ;
else if(ampls.size()==3){C_Dn=(Float_t)ampls[2] ; if(C_Dn<=coeff1*coeffn*C_D+coeff2*coeffnn*C_M+2*noise || C_D==255) shapecdtn=true;}
else if(ampls.size()>3){ C_Dn=(Float_t)ampls[2]; C_Dnn=(Float_t)ampls[3] ;
if((C_Dn<=coeff1*coeffn*C_D+coeff2*coeffnn*C_M+2*noise || C_D==255)
&& C_Dnn<=coeff1*coeffn*C_Dn+coeff2*coeffnn*C_D+2*noise){
shapecdtn=true;}
}
}
if(MaxOnEnd==true){
C_M=(Float_t)ampls[ampls.size()-1]; C_D=(Float_t)ampls[ampls.size()-2];
if(ampls.size()<3) shapecdtn=true ;
else if(ampls.size()==3){C_Dn=(Float_t)ampls[0] ; if(C_Dn<=coeff1*coeffn*C_D+coeff2*coeffnn*C_M+2*noise || C_D==255) shapecdtn=true;}
else if(ampls.size()>3){C_Dn=(Float_t)ampls[ampls.size()-3] ; C_Dnn=(Float_t)ampls[ampls.size()-4] ;
if((C_Dn<=coeff1*coeffn*C_D+coeff2*coeffnn*C_M+2*noise || C_D==255)
&& C_Dnn<=coeff1*coeffn*C_Dn+coeff2*coeffnn*C_D+2*noise){
shapecdtn=true;}
}
}
if(MaxInMiddle==true){
C_M=(Float_t)ampls[MaxPos];
int LeftOfMaxPos=MaxPos-1;if(LeftOfMaxPos<=0)LeftOfMaxPos=0;
int RightOfMaxPos=MaxPos+1;if(RightOfMaxPos>=(int)ampls.size())RightOfMaxPos=ampls.size()-1;
//int after = RightOfMaxPos; int before = LeftOfMaxPos; if (after>=(int)ampls.size() || before<0) std::cout<<"invalid read MaxPos:"<<MaxPos <<"size:"<<ampls.size() <<std::endl;
if(ampls[LeftOfMaxPos]<ampls[RightOfMaxPos]){ C_D=(Float_t)ampls[RightOfMaxPos]; C_Mn=(Float_t)ampls[LeftOfMaxPos];CDPos=RightOfMaxPos;} else{ C_D=(Float_t)ampls[LeftOfMaxPos]; C_Mn=(Float_t)ampls[RightOfMaxPos];CDPos=LeftOfMaxPos;}
if(C_Mn<coeff1*coeffn*C_M+coeff2*coeffnn*C_D+2*noise || C_M==255){
if(ampls.size()==3) shapecdtn=true ;
else if(ampls.size()>3){
if(CDPos>MaxPos){
if(ampls.size()-CDPos-1==0){
C_Dn=0; C_Dnn=0;
}
if(ampls.size()-CDPos-1==1){
C_Dn=(Float_t)ampls[CDPos+1];
C_Dnn=0;
}
if(ampls.size()-CDPos-1>1){
C_Dn=(Float_t)ampls[CDPos+1];
C_Dnn=(Float_t)ampls[CDPos+2];
}
if(MaxPos>=2){
C_Mnn=(Float_t)ampls[MaxPos-2];
}
else if(MaxPos<2) C_Mnn=0;
}
if(CDPos<MaxPos){
if(CDPos==0){
C_Dn=0; C_Dnn=0;
}
if(CDPos==1){
C_Dn=(Float_t)ampls[0];
C_Dnn=0;
}
if(CDPos>1){
C_Dn=(Float_t)ampls[CDPos-1];
C_Dnn=(Float_t)ampls[CDPos-2];
}
if(ampls.size()-LeftOfMaxPos>1 && MaxPos+2<(int)(ampls.size())-1){
C_Mnn=(Float_t)ampls[MaxPos+2];
}else C_Mnn=0;
}
if((C_Dn<=coeff1*coeffn*C_D+coeff2*coeffnn*C_M+2*noise || C_D==255)
&& C_Mnn<=coeff1*coeffn*C_Mn+coeff2*coeffnn*C_M+2*noise
&& C_Dnn<=coeff1*coeffn*C_Dn+coeff2*coeffnn*C_D+2*noise) {
shapecdtn=true;
}
}
}
}
}
if(ampls.size()==1){shapecdtn=true;}
return shapecdtn;
}
| void DeDxTools::trajectoryRawHits | ( | const edm::Ref< std::vector< Trajectory > > & | trajectory, |
| vector< RawHits > & | hits, | ||
| bool | usePixel, | ||
| bool | useStrip | ||
| ) |
Definition at line 17 of file DeDxTools.cc.
References DeDxTools::RawHits::angleCosine, DeDxTools::RawHits::charge, SiStripRecHit2D::cluster(), DeDxTools::RawHits::detId, TrackingRecHit::geographicalId(), i, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), DeDxTools::RawHits::NSaturating, and DeDxTools::RawHits::trajectoryMeasurement.
{
// vector<RawHits> hits;
const vector<TrajectoryMeasurement> & measurements = trajectory->measurements();
for(vector<TrajectoryMeasurement>::const_iterator it = measurements.begin(); it!=measurements.end(); it++){
//FIXME: check that "updated" State is the best one (wrt state in the middle)
TrajectoryStateOnSurface trajState=it->updatedState();
if( !trajState.isValid()) continue;
const TrackingRecHit * recHit=(*it->recHit()).hit();
LocalVector trackDirection = trajState.localDirection();
double cosine = trackDirection.z()/trackDirection.mag();
if(const SiStripMatchedRecHit2D* matchedHit=dynamic_cast<const SiStripMatchedRecHit2D*>(recHit)){
if(!useStrip) continue;
RawHits mono,stereo;
mono.trajectoryMeasurement = &(*it);
stereo.trajectoryMeasurement = &(*it);
mono.angleCosine = cosine;
stereo.angleCosine = cosine;
const std::vector<uint8_t> & amplitudes = matchedHit->monoHit()->cluster()->amplitudes();
mono.charge = accumulate(amplitudes.begin(), amplitudes.end(), 0);
mono.NSaturating =0;
for(unsigned int i=0;i<amplitudes.size();i++){if(amplitudes[i]>=254)mono.NSaturating++;}
const std::vector<uint8_t> & amplitudesSt = matchedHit->stereoHit()->cluster()->amplitudes();
stereo.charge = accumulate(amplitudesSt.begin(), amplitudesSt.end(), 0);
stereo.NSaturating =0;
for(unsigned int i=0;i<amplitudes.size();i++){if(amplitudes[i]>=254)stereo.NSaturating++;}
mono.detId= matchedHit->monoHit()->geographicalId();
stereo.detId= matchedHit->stereoHit()->geographicalId();
hits.push_back(mono);
hits.push_back(stereo);
}else if(const ProjectedSiStripRecHit2D* projectedHit=dynamic_cast<const ProjectedSiStripRecHit2D*>(recHit)) {
if(!useStrip) continue;
const SiStripRecHit2D* singleHit=&(projectedHit->originalHit());
RawHits mono;
mono.trajectoryMeasurement = &(*it);
mono.angleCosine = cosine;
const std::vector<uint8_t> & amplitudes = singleHit->cluster()->amplitudes();
mono.charge = accumulate(amplitudes.begin(), amplitudes.end(), 0);
mono.NSaturating =0;
for(unsigned int i=0;i<amplitudes.size();i++){if(amplitudes[i]>=254)mono.NSaturating++;}
mono.detId= singleHit->geographicalId();
hits.push_back(mono);
}else if(const SiStripRecHit2D* singleHit=dynamic_cast<const SiStripRecHit2D*>(recHit)){
if(!useStrip) continue;
RawHits mono;
mono.trajectoryMeasurement = &(*it);
mono.angleCosine = cosine;
const std::vector<uint8_t> & amplitudes = singleHit->cluster()->amplitudes();
mono.charge = accumulate(amplitudes.begin(), amplitudes.end(), 0);
mono.NSaturating =0;
for(unsigned int i=0;i<amplitudes.size();i++){if(amplitudes[i]>=254)mono.NSaturating++;}
mono.detId= singleHit->geographicalId();
hits.push_back(mono);
}else if(const SiStripRecHit1D* single1DHit=dynamic_cast<const SiStripRecHit1D*>(recHit)){
if(!useStrip) continue;
RawHits mono;
mono.trajectoryMeasurement = &(*it);
mono.angleCosine = cosine;
const std::vector<uint8_t> & amplitudes = single1DHit->cluster()->amplitudes();
mono.charge = accumulate(amplitudes.begin(), amplitudes.end(), 0);
mono.NSaturating =0;
for(unsigned int i=0;i<amplitudes.size();i++){if(amplitudes[i]>=254)mono.NSaturating++;}
mono.detId= single1DHit->geographicalId();
hits.push_back(mono);
}else if(const SiPixelRecHit* pixelHit=dynamic_cast<const SiPixelRecHit*>(recHit)){
if(!usePixel) continue;
RawHits pixel;
pixel.trajectoryMeasurement = &(*it);
pixel.angleCosine = cosine;
pixel.charge = pixelHit->cluster()->charge();
pixel.NSaturating=-1;
pixel.detId= pixelHit->geographicalId();
hits.push_back(pixel);
}
}
// return hits;
}
| void DeDxTools::trajectoryRawHits | ( | const edm::Ref< std::vector< Trajectory > > & | trajectory, |
| std::vector< RawHits > & | hits, | ||
| bool | usePixel, | ||
| bool | useStrip | ||
| ) |
1.7.3