MatacqProducer::PosEstimator Class Reference

Public Member Functions
int	eventLength () const
int	firstOrbit () const
void	init (MatacqProducer *mp)
bool	invalid () const
int64_t	pos (int orb) const
	PosEstimator ()
void	verbosity (int verb)

Private Attributes
int	eventLength_
int	firstOrbit_
bool	invalid_
int	orbitStepMean_
int	verbosity_

Detailed Description

Estimates matacq event position in a file from its orbit id. This estimator requires that every event in the file has the same length. A linear extrapolation of pos=f(orbit) function from first and last event is performed. It gives only a rough estimate, relevant only to initiliaze the event search.

Definition at line 99 of file MatacqProducer.h.

Constructor & Destructor Documentation

MatacqProducer::PosEstimator::PosEstimator ( )

inline

Definition at line 104 of file MatacqProducer.h.

References init.

                  :eventLength_(0), orbitStepMean_(0), firstOrbit_(0),
           invalid_(true), verbosity_(0) { }

Member Function Documentation

int MatacqProducer::PosEstimator::eventLength ( ) const

inline

Definition at line 109 of file MatacqProducer.h.

Referenced by MatacqProducer::getMatacqEvent().

{ return eventLength_; }

int MatacqProducer::PosEstimator::firstOrbit ( ) const

inline

Definition at line 110 of file MatacqProducer.h.

{ return firstOrbit_; }

void MatacqProducer::PosEstimator::init

(

MatacqProducer *

mp

)

Definition at line 738 of file MatacqProducer.cc.

References gather_cfg::cout, data, MatacqRawEvent::getDccLen(), MatacqRawEvent::getOrbitId(), sistrip::invalid_, plotBeamSpotDB::last, MatacqProducer::mread(), MatacqProducer::mrewind(), MatacqProducer::mseek(), MatacqProducer::msize(), nEvents, now(), alignCSCRings::s, and MatacqProducer::verbosity_.

Referenced by MatacqProducer::getMatacqFile().

                                                       {
  mp->mrewind();

  const size_t headerSize = 8*8;
  unsigned char data[headerSize];
  if(!mp->mread((char*)data, headerSize)){
    if(verbosity_) cout << "[Matacq " << now() << "] reached end of file!\n"; 
    firstOrbit_ = eventLength_ = orbitStepMean_ = 0;
    return;
  } else{
    firstOrbit_ = MatacqRawEvent::getOrbitId(data, headerSize);
    eventLength_ = MatacqRawEvent::getDccLen(data, headerSize);
    if(verbosity_>1) cout << "[Matacq " << now() << "] First event orbit: " << firstOrbit_
                          << " event length: " << eventLength_
                          << "*8 byte\n";
  }

  mp->mrewind();
  
  if(eventLength_==0){    
    if(verbosity_) cout << "[Matacq " << now() << "] event length is null!" << endl; 
    return;
  }

  filepos_t s;
  mp->msize(s);

  //number of complete events:
  const unsigned nEvents = s/eventLength_/8;

  if(nEvents==0){
    if(verbosity_) cout << "[Matacq " << now() << "] File is empty!" << endl;
    orbitStepMean_ = 0;
    return;
  }

  if(verbosity_>1) cout << "[Matacq " << now() << "] File size: " << s
                        << " Number of events: " << nEvents << endl;
  
  //position of last complete events:
  off_t last = (nEvents-1)*(off_t)eventLength_*8;
  mp->mseek(last, SEEK_SET, "Moving to beginning of last complete "
        "matacq event");
  if(!mp->mread((char*) data, headerSize, "Reading matacq header", true)){
    LogWarning("Matacq") << "Fast matacq event retrieval failure. "
      "Falling back to safe retrieval mode.";
    orbitStepMean_ = 0;
  }
  
  int32_t lastOrb = MatacqRawEvent::getOrbitId(data, headerSize);
  int32_t lastLen = MatacqRawEvent::getDccLen(data, headerSize);

  if(verbosity_>1) cout << "[Matacq " << now() << "] Last event orbit: " << lastOrb
                        << " last event length: " << lastLen << endl;
  
  //some consistency check
  if(lastLen!=eventLength_){
    LogWarning("Matacq")
      //throw cms::Exception("Matacq")
      << "Fast matacq event retrieval failure: it looks like "
      "the matacq file contains events of different sizes.";
      //      " Falling back to safe retrieval mode.";
    invalid_ = false; //true;
    orbitStepMean_ = 112; //0;
    return;
  }

  orbitStepMean_ = (lastOrb - firstOrbit_)/nEvents;
  
  if(verbosity_>1) cout << "[Matacq " << now() << "] Orbit step mean: " << orbitStepMean_
                        << "\n";

  invalid_ = false;
}

bool MatacqProducer::PosEstimator::invalid ( ) const

inline

Definition at line 107 of file MatacqProducer.h.

References sistrip::invalid_.

Referenced by MatacqProducer::getMatacqEvent().

{ return invalid_; }

int64_t MatacqProducer::PosEstimator::pos

(

int

orb

)

const

Definition at line 813 of file MatacqProducer.cc.

References gather_cfg::cout, now(), alignCSCRings::r, and MatacqProducer::verbosity_.

Referenced by Vispa.Gui.PortWidget.PortWidget::connectionPoint(), Vispa.Gui.VispaWidget.VispaWidget::dragWidget(), MatacqProducer::getMatacqEvent(), and Vispa.Gui.VispaWidget.VispaWidget::paintEvent().

                                                    {
  if(orb<firstOrbit_) return -1;
  uint64_t r = orbitStepMean_!=0?
    (((uint64_t)(orb-firstOrbit_))/orbitStepMean_)*eventLength_*8
    :0;
  if(verbosity_>2) cout << "[Matacq " << now() << "] Estimated Position for orbit  " << orb
                        << ": " << r << endl;
  return r;
}

void MatacqProducer::PosEstimator::verbosity ( int  verb )

inline

Definition at line 111 of file MatacqProducer.h.

Referenced by MatacqProducer::MatacqProducer().

{ verbosity_ = verb; }

Member Data Documentation

int MatacqProducer::PosEstimator::eventLength_

private

Definition at line 113 of file MatacqProducer.h.

int MatacqProducer::PosEstimator::firstOrbit_

private

Definition at line 115 of file MatacqProducer.h.

bool MatacqProducer::PosEstimator::invalid_

private

Definition at line 116 of file MatacqProducer.h.

int MatacqProducer::PosEstimator::orbitStepMean_

private

Definition at line 114 of file MatacqProducer.h.

int MatacqProducer::PosEstimator::verbosity_

private

Definition at line 117 of file MatacqProducer.h.