DeviationsFromFileSensor2D Class Reference

#include <DeviationsFromFileSensor2D.h>

Public Member Functions
	DeviationsFromFileSensor2D ()
std::pair< ALIdouble, ALIdouble >	getDevis (ALIdouble intersX, ALIdouble intersY)
void	readFile (ALIFileIn &ifdevi)
void	setOffset (ALIdouble offX, ALIdouble offY)
	~DeviationsFromFileSensor2D ()

Static Public Member Functions
static const ALIbool	apply ()
static void	setApply (ALIbool val)

Private Attributes
bool	firstScanDir
vvd	theDeviations
ALIuint	theNPoints
ALIdouble	theOffsetX
ALIdouble	theOffsetY
int	theScanSenseX
int	theScanSenseY
ALIint	verbose

Static Private Attributes
static ALIbool	theApply = true

Detailed Description

Definition at line 23 of file DeviationsFromFileSensor2D.h.

Constructor & Destructor Documentation

DeviationsFromFileSensor2D()

DeviationsFromFileSensor2D::DeviationsFromFileSensor2D ( )

inline

Definition at line 26 of file DeviationsFromFileSensor2D.h.

                               {
    theOffsetX = 0.;
    theOffsetY = 0.;
  };

References theOffsetX, and theOffsetY.

~DeviationsFromFileSensor2D()

DeviationsFromFileSensor2D::~DeviationsFromFileSensor2D ( )

inline

Definition at line 30 of file DeviationsFromFileSensor2D.h.

{};

Member Function Documentation

apply()

static const ALIbool DeviationsFromFileSensor2D::apply ( )

inlinestatic

Definition at line 45 of file DeviationsFromFileSensor2D.h.

{ return theApply; }

References theApply.

Referenced by heavyIonTools.ConfigureHeavyIons::__call__(), coreTools.RunOnData::__call__(), trackTools.MakeAODTrackCandidates::__call__(), runJetUncertainties.RunJetUncertainties::__call__(), metTools.AddMETCollection::__call__(), heavyIonTools.ProductionDefaults::__call__(), cmsswVersionTools.PickRelValInputFiles::__call__(), coreTools.RemoveMCMatching::__call__(), trackTools.MakePATTrackCandidates::__call__(), trigTools.SwitchOnTrigger::__call__(), runMETCorrectionsAndUncertainties.RunMETCorrectionsAndUncertainties::__call__(), heavyIonTools.SelectionDefaults::__call__(), heavyIonTools.DisbaleMonteCarloDeps::__call__(), trigTools.SwitchOnTriggerStandAlone::__call__(), tauTools.AddTauCollection::__call__(), trackTools.MakeTrackCandidates::__call__(), trigTools.SwitchOnTriggerMatching::__call__(), trigTools.SwitchOnTriggerMatchingStandAlone::__call__(), trigTools.SwitchOnTriggerMatchEmbedding::__call__(), jetTools.AddJetCollection::__call__(), jetTools.SwitchJetCollection::__call__(), jetTools.UpdateJetCollection::__call__(), jetTools.AddJetID::__call__(), jetTools.SetTagInfos::__call__(), OptOSensor2D::detailedTraversesLightRay(), and OptOSensor2D::fastTraversesLightRay().

getDevis()

std::pair< ALIdouble, ALIdouble > DeviationsFromFileSensor2D::getDevis	(	ALIdouble	intersX,
		ALIdouble	intersY
	)

Definition at line 146 of file DeviationsFromFileSensor2D.cc.

                                                                                                       {
  //intersX += 10000;
  //intersY += 10000;
  if (verbose >= 4)
    std::cout << " entering getdevis " << intersX << " " << intersY << " " << this << std::endl;
  vvd::iterator vvdite;
  vd::iterator vdite;
 
  // assume scan is in Y first, else revers intersection coordinates
  /*  if( !firstScanDirY ) {
    ALIdouble tt = intersY;
    intersY = intersX;
    intersX = tt;
  }
*/
 
  //---------- look which point in the deviation matrices correspond to intersX,Y
  //----- Look for each column, between which rows intersY is
  //assume first dir is Y
  auto yrows = std::make_unique<unsigned int[]>(theNPoints);
 
  unsigned int ii = 0;
  ALIbool insideMatrix = false;
  for (vvdite = theDeviations.begin(); vvdite != (theDeviations.end() - 1); ++vvdite) {
    for (vdite = (*vvdite).begin(); vdite != ((*vvdite).end() - 1); ++vdite) {
      if (verbose >= 5)
        std::cout << " check posy " << (*(vdite))->posY() << " " << (*(vdite + 1))->posY() << " " << (*(vdite))
                  << std::endl;
      // if posy is between this point and previous one
 
      //-     std::cout << "intersy" << intersY << " " <<  (*(vdite))->posY() << std::endl;
 
      if ((intersY - (*(vdite))->posY()) * theScanSenseY > 0 &&
          (intersY - (*(vdite + 1))->posY()) * theScanSenseY < 0) {
        //-std::cout << " ii " << ii << std::endl;
        yrows[ii] = vdite - (*vvdite).begin();
        if (verbose >= 3)
          std::cout << intersY << " DeviationsFromFileSensor2D yrows " << ii << " " << yrows[ii] << " : "
                    << (*(vdite))->posY() << std::endl;
        insideMatrix = true;
        break;
      }
    }
    ii++;
  }
  if (insideMatrix == 0) {
    std::cerr << "!!EXITING intersection in Y outside matrix of deviations from file " << intersY << std::endl;
    exit(1);
  }
  insideMatrix = false;
 
  vd thePoints;
  thePoints.clear();
  //----- For each row in 'yrows' look between which columns intersX is
  unsigned int rn;
  DeviationSensor2D *dev1, *dev2;
  for (ii = 0; ii < theNPoints - 1; ii++) {
    rn = yrows[ii];
    //-    std::cout << ii << " rn " << rn << std::endl;
    dev1 = (*(theDeviations.begin() + ii))[rn];  // column ii, row yrows[ii]
    rn = yrows[ii + 1];
    dev2 = (*(theDeviations.begin() + ii + 1))[rn];  // column ii+1, row yrows[ii+1]
    if ((intersX - dev1->posX()) * theScanSenseX > 0 && (intersX - dev2->posX()) * theScanSenseX < 0) {
      thePoints.push_back(dev1);
      thePoints.push_back(dev2);
      insideMatrix = true;
      if (verbose >= 3)
        std::cout << " column up " << ii << " " << dev1->posX() << " " << dev2->posX() << " : " << intersX << std::endl;
    }
 
    rn = yrows[ii] + 1;
    if (rn == theNPoints)
      rn = theNPoints - 1;
    dev1 = (*(theDeviations.begin() + ii))[rn];  // column ii, row yrows[ii]+1
    rn = yrows[ii + 1] + 1;
    if (rn == theNPoints)
      rn = theNPoints - 1;
    dev2 = (*(theDeviations.begin() + ii + 1))[rn];  // column ii+1, row yrows[ii+1]+1
    if ((intersX - dev1->posX()) * theScanSenseX > 0 && (intersX - dev2->posX()) * theScanSenseX < 0) {
      thePoints.push_back(dev1);
      thePoints.push_back(dev2);
      if (verbose >= 3)
        std::cout << " column down " << ii << " " << dev1->posX() << " " << dev2->posX() << " : " << intersX
                  << std::endl;
    }
 
    if (thePoints.size() == 4)
      break;
  }
 
  if (insideMatrix == 0) {
    std::cerr << "!!EXITING intersection in X outside matrix of deviations from file " << intersX << std::endl;
    exit(1);
  }
 
  //----------- If loop finished and not 4 points, point is outside scan bounds
 
  //----- calculate deviation in x and y interpolating between four points
  ALIdouble dist, disttot = 0, deviX = 0, deviY = 0;
 
  if (verbose >= 4)
    std::cout << " thepoints size " << thePoints.size() << std::endl;
 
  for (ii = 0; ii < 4; ii++) {
    dist = sqrt(pow(thePoints[ii]->posX() - intersX, 2) + pow(thePoints[ii]->posY() - intersY, 2));
    disttot += 1. / dist;
    deviX += thePoints[ii]->devX() / dist;
    deviY += thePoints[ii]->devY() / dist;
    if (verbose >= 4) {
      //t      std::cout << ii << " point " << *thePoints[ii] << std::endl;
      std::cout << ii << " distances: " << dist << " " << deviX << " " << deviY << " devX " << thePoints[ii]->devX()
                << " devY " << thePoints[ii]->devY() << std::endl;
    }
  }
  deviX /= disttot;
  deviY /= disttot;
 
  // add offset
  deviX += theOffsetX;
  deviY += theOffsetY;
  if (verbose >= 4) {
    std::cout << " devisX/Y: " << deviX << " " << deviY << " intersX/Y " << intersX << " " << intersY << std::endl;
  }
 
  //change sign!!!!!?!?!?!?!?!
  deviX *= -1;
  deviY *= -1;
  return std::pair<ALIdouble, ALIdouble>(deviX * 1.e-6, deviY * 1e-6);  // matrix is in microrad
}

References beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, MillePedeFileConverter_cfg::e, beamvalidation::exit(), cuy::ii, DeviationSensor2D::posX(), RecoTauValidation_cfi::posX, RecoTauValidation_cfi::posY, funct::pow(), mathSSE::sqrt(), theDeviations, theNPoints, theOffsetX, theOffsetY, theScanSenseX, theScanSenseY, and verbose.

readFile()

void DeviationsFromFileSensor2D::readFile

(

ALIFileIn &

ifdevi

)

Definition at line 20 of file DeviationsFromFileSensor2D.cc.

                                                           {
  verbose = ALIUtils::debug;
  //  verbose = 4;
  if (verbose >= 3)
    std::cout << "DeviationsFromFileSensor2D::readFile " << this << ifdevi.name() << std::endl;
 
  theScanSenseX = 0;
  theScanSenseY = 0;
 
  ALIuint nl = 1;
 
  ALIdouble oldposX = 0, oldposY = 0;
  vd vcol;
  std::vector<ALIstring> wl;
  /*  //------ first line with dimension factors //now with global option
  ifdevi.getWordsInLine( wl );
  if( wl[0] != "DIMFACTOR:" || wl.size() != 3) {
    std::cerr << "Error reading sensor2D deviation file " << ifdevi.name() << std::endl
     << " first line has to be DIMFACTOR: 'posDimFactor' 'angDimFactor' " << std::endl;
    ALIUtils::dumpVS( wl, " ");
    exit(1);
  }
  ALIdouble posDimFactor = atof(wl[1].c_str());
  ALIdouble angDimFactor = atof(wl[2].c_str());
  */
 
  for (;;) {
    ifdevi.getWordsInLine(wl);
    if (ifdevi.eof())
      break;
 
    DeviationSensor2D* dev = new DeviationSensor2D();
    dev->fillData(wl);
 
    if (verbose >= 5) {
      ALIUtils::dumpVS(wl, "deviation sensor2D", std::cout);
    }
 
    //--- line 2 of data
    if (nl == 2) {
      //--------- get if scan is first in Y or X
      firstScanDir = ydir;
      if (verbose >= 3)
        std::cout << "firstScanDir " << firstScanDir << " " << dev->posX() << " " << oldposX << " " << dev->posY()
                  << " " << oldposY << std::endl;
      if (std::abs(dev->posX() - oldposX) > std::abs(dev->posY() - oldposY)) {
        std::cerr << "!!!EXITING first scan direction has to be Y for the moment " << ifdevi.name() << std::endl;
        firstScanDir = xdir;
        exit(1);
      }
      //-------- get sense of first scan direction
      if (firstScanDir == ydir) {
        if (dev->posY() > oldposY) {
          theScanSenseY = +1;
        } else {
          theScanSenseY = -1;
        }
        if (verbose >= 3)
          std::cout << " theScanSenseY " << theScanSenseY << std::endl;
      } else {
        if (dev->posX() > oldposX) {
          theScanSenseX = +1;
        } else {
          theScanSenseX = -1;
        }
        if (verbose >= 3)
          std::cout << " theScanSenseX " << theScanSenseX << std::endl;
      }
    }
 
    //-      std::cout << "firstScanDir " << firstScanDir << " " <<  dev->posX() <<  " " << oldposX << " " <<  dev->posY() <<  " " << oldposY << std::endl;
 
    //------- check if change of row: clear current std::vector of column values
    if ((firstScanDir == ydir && (dev->posY() - oldposY) * theScanSenseY < 0) ||
        (firstScanDir == xdir && (dev->posX() - oldposX) * theScanSenseX < 0)) {
      theDeviations.push_back(vcol);
      vcol.clear();
 
      //-      std::cout << " theDevi size " << theDeviations.size() << " " << ifdevi.name()  << std::endl;
      //-------- get sense of second scan direction
      if (theScanSenseY == 0) {
        if (dev->posY() > oldposY) {
          theScanSenseY = +1;
        } else {
          theScanSenseY = -1;
        }
      }
      if (theScanSenseX == 0) {
        if (dev->posX() > oldposX) {
          theScanSenseX = +1;
        } else {
          theScanSenseX = -1;
        }
      }
      if (verbose >= 3)
        std::cout << " theScanSenseX " << theScanSenseX << " theScanSenseY " << theScanSenseY << std::endl;
    }
 
    oldposX = dev->posX();
    oldposY = dev->posY();
 
    //--- fill deviation std::vectors
    vcol.push_back(dev);
    nl++;
  }
  theDeviations.push_back(vcol);
 
  //----- Calculate std::vector size
  ALIdouble dvsiz = (sqrt(ALIdouble(nl - 1)));
  //----- Check that it is a square of points (<=> vsiz is integer)
  ALIuint vsiz = ALIuint(dvsiz);
  if (vsiz != dvsiz) {
    if (ALIUtils::debug >= 0)
      std::cerr << "!!WARNING: error reading deviation file: number of X points <> number of Y points : Number of "
                   "points in X "
                << dvsiz << " nl " << nl - 1 << " file " << ifdevi.name() << std::endl;
    //    exit(1);
  }
  theNPoints = vsiz;
 
  if (verbose >= 4) {
    std::cout << " Filling deviation from file: " << ifdevi.name() << " theNPoints " << theNPoints << std::endl;
  }
}

References funct::abs(), beam_dqm_sourceclient-live_cfg::cerr, gather_cfg::cout, ALIUtils::debug, ALIUtils::dumpVS(), ALIFileIn::eof(), beamvalidation::exit(), DeviationSensor2D::fillData(), firstScanDir, ALIFileIn::getWordsInLine(), ALIFileIn::name(), DeviationSensor2D::posX(), DeviationSensor2D::posY(), mathSSE::sqrt(), theDeviations, theNPoints, theScanSenseX, theScanSenseY, verbose, LaserClient_cfi::wl, xdir, and ydir.

setApply()

static void DeviationsFromFileSensor2D::setApply ( ALIbool  val )

inlinestatic

Definition at line 48 of file DeviationsFromFileSensor2D.h.

{ theApply = val; }

References theApply, and heppy_batch::val.

Referenced by Fit::calculateSimulatedMeasurementsWithOriginalValues(), and OptOSensor2D::detailedTraversesLightRay().

setOffset()

void DeviationsFromFileSensor2D::setOffset ( ALIdouble  offX, ALIdouble  offY  )

inline

Definition at line 39 of file DeviationsFromFileSensor2D.h.

                                                 {
    theOffsetX = offX;
    theOffsetY = offY;
  }

References theOffsetX, and theOffsetY.

Member Data Documentation

firstScanDir

bool DeviationsFromFileSensor2D::firstScanDir

private

Definition at line 51 of file DeviationsFromFileSensor2D.h.

Referenced by readFile().

theApply

ALIbool DeviationsFromFileSensor2D::theApply = true

staticprivate

Definition at line 56 of file DeviationsFromFileSensor2D.h.

Referenced by apply(), and setApply().

theDeviations

vvd DeviationsFromFileSensor2D::theDeviations

private

Definition at line 55 of file DeviationsFromFileSensor2D.h.

Referenced by getDevis(), and readFile().

theNPoints

ALIuint DeviationsFromFileSensor2D::theNPoints

private

Definition at line 54 of file DeviationsFromFileSensor2D.h.

Referenced by getDevis(), and readFile().

theOffsetX

ALIdouble DeviationsFromFileSensor2D::theOffsetX

private

Definition at line 60 of file DeviationsFromFileSensor2D.h.

Referenced by DeviationsFromFileSensor2D(), getDevis(), and setOffset().

theOffsetY

ALIdouble DeviationsFromFileSensor2D::theOffsetY

private

Definition at line 60 of file DeviationsFromFileSensor2D.h.

Referenced by DeviationsFromFileSensor2D(), getDevis(), and setOffset().

theScanSenseX

int DeviationsFromFileSensor2D::theScanSenseX

private

Definition at line 52 of file DeviationsFromFileSensor2D.h.

Referenced by getDevis(), and readFile().

theScanSenseY

int DeviationsFromFileSensor2D::theScanSenseY

private

Definition at line 52 of file DeviationsFromFileSensor2D.h.

Referenced by getDevis(), and readFile().

verbose

ALIint DeviationsFromFileSensor2D::verbose

private

Definition at line 58 of file DeviationsFromFileSensor2D.h.

Referenced by getDevis(), python.diff_provenance.difference::list_diff(), python.diffProv.difference::list_diff(), python.diff_provenance.difference::module_diff(), python.diffProv.difference::module_diff(), python.diff_provenance.difference::onefilemodules(), python.diffProv.difference::onefilemodules(), core.TriggerMatchAnalyzer.TriggerMatchAnalyzer::process(), core.SkimAnalyzerCount.SkimAnalyzerCount::process(), objects.VertexAnalyzer.VertexAnalyzer::process(), confdbOfflineConverter.OfflineConverter::query(), and readFile().