|
|
|
#include <OptOSensor2D.h>
Public Member Functions | |
| void | constructSolidShape () |
| ALIdouble * | convertPointToLocalCoordinates (const CLHEP::Hep3Vector &point) |
| virtual void | defaultBehaviour (LightRay &lightray, Measurement &meas) |
| virtual void | detailedTraversesLightRay (LightRay &lightray) |
| virtual void | fastTraversesLightRay (LightRay &lightray) |
| virtual void | fillExtraEntry (std::vector< ALIstring > &wordlist) |
| virtual void | makeMeasurement (LightRay &lightray, Measurement &meas) |
| OptOSensor2D () | |
| OptOSensor2D (OpticalObject *parent, const ALIstring &type, const ALIstring &name, const ALIbool copy_data) | |
| ~OptOSensor2D () | |
Private Attributes | |
| DeviationsFromFileSensor2D * | deviFromFile |
| ALIbool | fdevi_from_file |
Definition at line 19 of file OptOSensor2D.h.
| OptOSensor2D::OptOSensor2D | ( | ) | [inline] |
Definition at line 24 of file OptOSensor2D.h.
{ };
| OptOSensor2D::OptOSensor2D | ( | OpticalObject * | parent, |
| const ALIstring & | type, | ||
| const ALIstring & | name, | ||
| const ALIbool | copy_data | ||
| ) | [inline] |
Definition at line 25 of file OptOSensor2D.h.
: OpticalObject( parent, type, name, copy_data), fdevi_from_file(0){ };
| OptOSensor2D::~OptOSensor2D | ( | ) | [inline] |
Definition at line 27 of file OptOSensor2D.h.
{ };
| void OptOSensor2D::constructSolidShape | ( | ) | [virtual] |
Reimplemented from OpticalObject.
Definition at line 454 of file OptOSensor2D.cc.
References GlobalOptionMgr::getGlobalOptionValue(), GlobalOptionMgr::getInstance(), and m.
{
ALIdouble go;
GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
gomgr->getGlobalOptionValue("VisScale", go );
theSolidShape = new CocoaSolidShapeBox( "Box", go*4.*cm/m, go*4.*cm/m, go*1.*cm/m ); //COCOA internal units are meters
}
| ALIdouble * OptOSensor2D::convertPointToLocalCoordinates | ( | const CLHEP::Hep3Vector & | point | ) |
Definition at line 399 of file OptOSensor2D.cc.
References gather_cfg::cout, ALIUtils::debug, ALIUtils::dump3v(), and ALIUtils::dumprm().
{
ALIdouble* interslc = new ALIdouble[2];
//----- X value
CLHEP::HepRotation rmt = rmGlob();
CLHEP::Hep3Vector XAxism(1.,0.,0.);
XAxism*=rmt;
if( ALIUtils::debug >= 5) ALIUtils::dump3v( (this)->centreGlob(), "centre glob sensor2D" );
if( ALIUtils::debug >= 5) ALIUtils::dumprm( rmt, "rotation matrix sensor2D" );
//t ALIUtils::dump3v(point - (this)->centreGlob() , "inters - (this)->centreGlob()");
if( ALIUtils::debug >= 5) ALIUtils::dump3v(XAxism , "XAxism");
interslc[0] = (point - (this)->centreGlob() ) * XAxism;
//----- Y value
CLHEP::Hep3Vector YAxism(0.,1.,0.);
YAxism*=rmt;
if( ALIUtils::debug >= 5)
ALIUtils::dump3v(YAxism , "YAxism");
interslc[1] = (point - (this)->centreGlob() ) * YAxism;
if( ALIUtils::debug >= 5 ) {
std::cout << " intersection in local coordinates: X= " << interslc[0] << " Y= " << interslc[1] << std::endl;
ALIUtils::dump3v( point - (this)->centreGlob() , " inters - centre " );
}
return interslc;
}
| void OptOSensor2D::defaultBehaviour | ( | LightRay & | lightray, |
| Measurement & | meas | ||
| ) | [virtual] |
Reimplemented from OpticalObject.
Definition at line 33 of file OptOSensor2D.cc.
{
makeMeasurement( lightray, meas);
}
| void OptOSensor2D::detailedTraversesLightRay | ( | LightRay & | lightray | ) | [virtual] |
Reimplemented from OpticalObject.
Definition at line 247 of file OptOSensor2D.cc.
References DeviationsFromFileSensor2D::apply(), benchmark_cfg::cerr, gather_cfg::cout, ALIUtils::debug, LightRay::dumpData(), LightRay::intersect(), mergeVDriftHistosByStation::name, LightRay::point(), LightRay::refract(), DeviationsFromFileSensor2D::setApply(), LightRay::setPoint(), and tablePrinter::width.
{
if (ALIUtils::debug >= 4) std::cout << "%%% LR: DETAILED TRAVERSES SENSOR2D " << name() << std::endl;
if( DeviationsFromFileSensor2D::apply() && fdevi_from_file) {
DeviationsFromFileSensor2D::setApply( 0 );
//- std::cout << "fdeviFromFile" << fdevi_from_file << std::endl;
if(ALIUtils::debug >= 0 )std::cerr << "!!WARNING: sensor " << name() << " has read deviation from file and it will not be taken into account. Please use FAST TRAVERSES" << deviFromFile << std::endl;
}
//---------- If width is 0, just keep the intersection point
ALIdouble width = findExtraEntryValue("width");
if( width == 0 ) {
//---------- Get intersection
CLHEP::Hep3Vector ZAxis(0.,0,1.);
CLHEP::HepRotation rmt = rmGlob();
ZAxis = rmt * ZAxis;
lightray.intersect( ALIPlane(centreGlob(), ZAxis) );
CLHEP::Hep3Vector inters = lightray.point();
lightray.setPoint( inters );
if (ALIUtils::debug >= 2) {
lightray.dumpData("LightRay Sensor2D traversed: ");
}
return;
}
if (ALIUtils::debug >= 4) std::cout << std::endl << "$$$ LR: REFRACTION IN FORWARD PLATE " << std::endl;
//---------- Get forward plate
ALIPlane plate = getPlate(1, 1);
//---------- Refract while entering object
ALIdouble refra_ind1 = 1.;
ALIdouble refra_ind2 = findExtraEntryValueMustExist("refra_ind");
lightray.refract( plate, refra_ind1, refra_ind2 );
if (ALIUtils::debug >= 4) std::cout << std::endl << "$$$ LR: REFRACTION IN BACKWARD PLATE " << std::endl;
//---------- Get backward plate
plate = getPlate(0, 1);
//---------- Refract while exiting splitter
lightray.refract( plate, refra_ind2, refra_ind1 );
CLHEP::Hep3Vector inters = lightray.point();
lightray.setPoint( inters );
if (ALIUtils::debug >= 4) {
lightray.dumpData("LightRay Sensor2D traversed: ");
}
}
| void OptOSensor2D::fastTraversesLightRay | ( | LightRay & | lightray | ) | [virtual] |
Reimplemented from OpticalObject.
Definition at line 119 of file OptOSensor2D.cc.
References DeviationsFromFileSensor2D::apply(), benchmark_cfg::cerr, gather_cfg::cout, debug, ALIUtils::debug, LightRay::dumpData(), LightRay::intersect(), ALIUtils::LengthValueDimensionFactor(), Model::MeasurementList(), Entry::name(), mergeVDriftHistosByStation::name, LightRay::point(), LightRay::setPoint(), Entry::setValue(), edm::shift, LightRay::shiftAndDeviateWhileTraversing(), Measurement::value(), and Entry::value().
{
verbose = ALIUtils::debug;
if (ALIUtils::debug >= 2) std::cout << "LR: FAST TRAVERSES SENSOR2D " << name() << std::endl;
//---------- Shift and Deviate
ALIdouble shiftX = findExtraEntryValue("shiftX");
ALIdouble shiftY = findExtraEntryValue("shiftY");
ALIdouble shift;
ALIbool bb = findExtraEntryValueIfExists("shift", shift);
if( bb ) {
shiftX = shift;
shiftY = shift;
}
//---------- Get intersection
CLHEP::Hep3Vector ZAxis(0.,0,1.);
CLHEP::HepRotation rmt = rmGlob();
ZAxis = rmt * ZAxis;
lightray.intersect( ALIPlane(centreGlob(), ZAxis) );
CLHEP::Hep3Vector inters = lightray.point();
ALIdouble deviX, deviY, devi;
// if deviationsFromFile are applied and this sensors has them get the deviations that corresponds to the intersection point.
// Set this deviation as original one, as it will not be changed for derivatives
if( DeviationsFromFileSensor2D::apply() && fdevi_from_file) {
//- std::cout << " DeviationsFromFileSensor2D::apply() " << DeviationsFromFileSensor2D::apply() << std::endl;
if(ALIUtils::debug >= 4) {
std::cout << "fdeviFromFile" << fdevi_from_file << std::endl;
//- std::cout << "deviFromFile" << deviFromFile << std::endl;
}
//--------- get measurement value of the current sensor
std::vector< Measurement* >& measv = Model::MeasurementList();
unsigned int ii;
Measurement *omeas = 0;
for( ii = 0; ii < measv.size(); ii++ ) {
//- std::cout << " sensor2d finding meas " << measv[ii]->sensorName() << " " << name() << std::endl;
if( measv[ii]->sensorName() == name() ) {
omeas = measv[ii];
break;
}
}
if( omeas == 0 ) {
std::cerr << "!!!EXITING OptOSensor2D::fastTraversesLightRay: meas " << name() << " not found " << std::endl;
}
ALIdouble interslcx = omeas->value( 0 );
ALIdouble interslcy = omeas->value( 1 );
if(ALIUtils::debug >= 5) std::cout << " interslcx " << interslcx << " interslcy " << interslcy << std::endl;
//----- transform in milimeters and positive
//mum interslcx = interslcx*1.E6 + 10000.;
//mum interslcy = interslcy*1.E6 + 10000.;
ALIdouble df = ALIUtils::LengthValueDimensionFactor();
interslcx = interslcx/df + 0.010/df;
interslcy = interslcy/df + 0.010/df;
if(ALIUtils::debug >= 5) std::cout << " interslcx " << interslcx << " interslcy " << interslcy << std::endl;
//---------- Get deviations from file (they are in microrads)
std::pair< ALIdouble, ALIdouble> devis = deviFromFile->getDevis( interslcx, interslcy );
deviX = devis.second;
deviY = -devis.first;
// deviX = devis.first;
// deviY = devis.second;
//o deviX = *((deviFromFile->deviationsInX()).begin() + pointY*deviFromFile->nPoints() + pointX) * 1.E-6;;
//o deviY = *((deviFromFile->deviationsInY().begin()) + pointY*deviFromFile->nPoints() + pointX) * 1.E-6;
//---------- Set this deviation value as original one, as it will not be changed for derivatives (change the entry and also the ExtraEntryValueOriginalList())
ALIuint entryNo = extraEntryNo( "deviX" );
if( verbose >= 3 ) std::cout << "entrynox" << entryNo << name() << verbose << std::endl;
Entry* entryDeviX = *(ExtraEntryList().begin()+entryNo);
entryDeviX->setValue( deviX );
//- std::vector< ALIdouble >::const_iterator eevolite = static_cast<std::vector< ALIdouble >::iterator>( ExtraEntryValueOriginalList().begin() );
std::vector< ALIdouble > eevil = ExtraEntryValueOriginalList();
//- std::vector< ALIdouble >::const_iterator eevolite = ( ExtraEntryValueOriginalList().begin() );
std::vector< ALIdouble >::iterator eevolite = eevil.begin();
*(eevolite+entryNo) = deviX;
if( verbose >= 3 ) std::cout<< " entryDeviX name " << entryDeviX->name() << entryDeviX->value() << std::endl;
entryNo = extraEntryNo( "deviY" );
Entry* entryDeviY = *(ExtraEntryList().begin()+entryNo);
//- std::cout << "entrynoy" << entryNo << name() << std::endl;
entryDeviY->setValue( deviY );
*(eevolite+entryNo) = deviY;
//- std::cout<< entryDeviY << " entryDeviY name " << entryDeviY->name() << entryDeviY->value() << std::endl;
} else {
deviX = findExtraEntryValue("deviX");
deviY = findExtraEntryValue("deviY");
//??? why previous does not work??
if( fdevi_from_file ) {
if( ALIUtils::debug >= 5) std::cout << "fdeviFromFile" << fdevi_from_file << std::endl;
ALIuint entryNo = extraEntryNo( "deviX" );
Entry* entryDeviX = *(ExtraEntryList().begin()+entryNo);
if( verbose >= 3 ) std::cout<< entryDeviX << " entryDeviX name " << entryDeviX->name() << entryDeviX->value() << std::endl;
deviX = entryDeviX->value();
entryNo = extraEntryNo( "deviY" );
Entry* entryDeviY = *(ExtraEntryList().begin()+entryNo);
if( verbose >= 3 ) std::cout<< entryDeviY << " entryDeviY name " << entryDeviY->name() << entryDeviY->value() << std::endl;
deviY = entryDeviY->value();
} else {
bb = findExtraEntryValueIfExists("devi", devi);
if( bb ) {
deviX = devi;
deviY = devi;
}
}
}
if(ALIUtils::debug >= 4) {
std::cout << "devi " << devi << " devi x " << deviX << " devi y " << deviY << std::endl;
}
lightray.setPoint( inters );
// lightray.shiftAndDeviateWhileTraversing( this, shiftX, shiftY, 0., deviX, deviY, 0.);
lightray.shiftAndDeviateWhileTraversing( this, 'T' );
if (ALIUtils::debug >= 2) {
lightray.dumpData("Shifted and Deviated");
}
}
| void OptOSensor2D::fillExtraEntry | ( | std::vector< ALIstring > & | wordlist | ) | [virtual] |
Reimplemented from OpticalObject.
Definition at line 300 of file OptOSensor2D.cc.
References gather_cfg::cout, ALIUtils::debug, OpticalObject::fillExtraEntry(), ALIUtils::getFloat(), ALIFileIn::getInstance(), and ALIFileIn::getWordsInLine().
{
if(ALIUtils::debug >= 5) std::cout << "OptOSensor2D fillExtraEntry wordlist[1] " << wordlist[1] << std::endl;
//---------- check if it is deviation read from file
fdevi_from_file = 0;
//- std::cout << "WL " << wordlist[1]<< "WL " << wordlist[2]<< "WL " << wordlist[3] << std::endl;
if( wordlist[1] == ALIstring("devi") && wordlist[2] == ALIstring("from_file") ) {
//---------- Open file
ALIstring fnam;
if( wordlist.size() >= 4) {
fnam = wordlist[3];
} else {
//----- Build up file name if it does not exists
fnam = "dat/devi-";
fnam += shortName();
fnam += ".dat";
}
ALIFileIn& ifdevi = ALIFileIn::getInstance( fnam );
//----- Check that file exists
if(ALIUtils::debug >= 4) std::cout << "Opening deviation file: " << fnam << std::endl;
/*- if( !ifdevi ) {
std::cerr << " !!! Sensor2D Deviation file not found: " << fnam << " of object " << name() << std::endl;
exit(1);
}*/
deviFromFile = new DeviationsFromFileSensor2D();
fdevi_from_file = 1;
if(ALIUtils::debug >= 5 ) std::cout << "deviFromFile " << deviFromFile << std::endl;
//----- Read header
ALIstring sensor1_name, sensor2_name;
ALIint sensor1_side, sensor2_side;
ALIdouble sensor_dist;
ALIdouble prec_deviX,prec_deviY;
std::vector<ALIstring> wl;
ifdevi.getWordsInLine( wl );
sensor1_name = wl[0];
sensor2_name = wl[1];
sensor_dist = atof( wl[2].c_str() );
// 'c' means that light is traversing the glass
sensor1_side = 1;
if(sensor1_name[sensor1_name.size()-2] == 'c') {
sensor1_side = 2;
sensor1_name = sensor1_name.substr(0,sensor1_name.size()-1);
}
sensor2_side = 1;
if(sensor2_name[sensor2_name.size()-2] == 'c') {
sensor2_side = 2;
sensor2_name = sensor2_name.substr(0,sensor2_name.size()-1);
}
if(ALIUtils::debug >= 5) std::cout << "sensor1_name " << sensor1_name << " sensor2_name " << sensor2_name << " sensor_dist " << sensor_dist << " unknown " << wl[3] << std::endl;
ifdevi.getWordsInLine( wl );
prec_deviX = atof( wl[0].c_str() );
prec_deviY = atof( wl[1].c_str() );
if(ALIUtils::debug >= 5) std::cout << "prec_deviX " << prec_deviX << " prec_deviY " << prec_deviY << std::endl;
deviFromFile = new DeviationsFromFileSensor2D();
ALIdouble offsetX, offsetY;
if( wl.size() == 5 ) {
offsetX = ALIUtils::getFloat( wl[3] );
offsetY = ALIUtils::getFloat( wl[4] );
deviFromFile->setOffset( offsetX, offsetY );
}
deviFromFile->readFile( ifdevi );
fdevi_from_file = 1;
if(ALIUtils::debug >= 5 ) std::cout << "deviFromFile " << deviFromFile << std::endl;
//--- Fill extra entries 'deviX' & 'deviY' to compute derivatives
std::vector<ALIstring> wlo;
char chartmp[20];
wlo.push_back( wordlist[0] );
wlo.push_back("deviX");
wlo.push_back("0"); // value is set to 0 as it is on the file and the point of intersection is not computed yet
gcvt( prec_deviX, 10, chartmp );
wlo.push_back( ALIstring(chartmp) );
wlo.push_back("cal");
std::vector<ALIstring> wl2(wlo);
OpticalObject::fillExtraEntry( wlo );
wl2[1] = "deviY";
gcvt( prec_deviY, 10, chartmp );
wl2[3] = ALIstring( chartmp );
OpticalObject::fillExtraEntry( wl2 );
} else {
OpticalObject::fillExtraEntry( wordlist );
}
}
| void OptOSensor2D::makeMeasurement | ( | LightRay & | lightray, |
| Measurement & | meas | ||
| ) | [virtual] |
Reimplemented from OpticalObject.
Definition at line 42 of file OptOSensor2D.cc.
References gather_cfg::cout, ALIUtils::debug, LightRay::direction(), ALIUtils::dump3v(), LightRay::intersect(), mergeVDriftHistosByStation::name, dbtoconf::parent, LightRay::point(), Measurement::setLightRayDirection(), Measurement::setLightRayPosition(), Measurement::setValueSimulated(), Measurement::value(), Measurement::valueSimulated(), and Measurement::valueType().
{
//---------- check is last and not the only one
/* if( vocite - OptOList().begin() == 0 ||
OptOList().end() - vocite != 1) {
std::cerr << " last and not only one Optical Object should be 'sensor2D' (unless you specify (:T)raverse) " <<
OptOList().end() - vocite << " size " <<OptOList().size() << std::endl;
DumpBadOrderOptOs();
}*/
//---------- Get simulated value
//---------- Get intersection
CLHEP::Hep3Vector ZAxis(0.,0,1.);
CLHEP::HepRotation rmt = rmGlob();
ZAxis = rmt * ZAxis;
if( ALIUtils::debug >= 4 ) {
ALIUtils::dump3v( centreGlob(), " sensor2D centre Glob ");
ALIUtils::dump3v( ZAxis, " snsor2D normal ");
}
//- ALIUtils::dumprm( rmt, "rot global " );
lightray.intersect( ALIPlane(centreGlob(), ZAxis) );
CLHEP::Hep3Vector inters = lightray.point();
ALIdouble* interslc;
interslc = convertPointToLocalCoordinates( inters );
ALIdouble interslcx = interslc[0];
ALIdouble interslcy = interslc[1];
meas.setValueSimulated( 0, interslcx );
meas.setValueSimulated( 1, interslcy );
//----- Dump info
if (ALIUtils::debug >= 2) {
//--- Special for range studies
ALIstring chrg = "";
/*t if(Model::Ncmslinkrange >= 1 && Model::Ncmslinkrange <= 8 ) {
chrg = "RG";
} else {
chrg = "";
} t*/
CLHEP::Hep3Vector measvv( meas.value()[0], meas.value()[1], 0.);
measvv = rmt*measvv;
ALIUtils::dump3v( measvv, " $$$$$$MEAS IN LOCAL FRAME");
ALIUtils::dump3v( measvv+centreGlob(), " $$$$$$MEAS IN GLOBAL FRAME");
ALIdouble detH = 1000*meas.valueSimulated(0); if(fabs(detH) <= 1.e-9 ) detH = 0.;
ALIdouble detV = 1000*meas.valueSimulated(1); if(fabs(detV) <= 1.e-9 ) detV = 0.;
std::cout << "REAL value: " << chrg << meas.valueType(0) << ": " << 1000*meas.value()[0] << chrg << " " << meas.valueType(1) << ": " << 1000*meas.value()[1] << " (mm) " << (this)->name()
<< " DIFF= " << detH-1000*meas.value()[0] << " " << detV-1000*meas.value()[1] << std::endl;
std::cout << "SIMU value: " << chrg << " " << meas.valueType(0) << ": "
// << setprecision(3) << setw(4)
<< detH
<< chrg << " " << meas.valueType(1) << ": " << detV
<< " (mm) " << (this)->name() << std::endl;
/*- std::cout << "SIMU value: " << chrg << " " << meas.valueType(0) << ": "
// << setprecision(3) << setw(4)
<< detH / 0.3125
<< chrg << " " << meas.valueType(1) << ": " << detV / 0.3125
<< " STRIPS " << (this)->name() << std::endl; */
// << detH
// << chrg << " V: " << detV
// << " (mm) " << (this)->name() << std::endl;
ALIUtils::dump3v( 1000.*(inters - parent()->centreGlob()) , " $$$$$$SIMU inters - parent centre");
ALIUtils::dump3v( 1000.*(inters - centreGlob()) , " $$$$$$SIMU inters - centre");
}
//t delete &lightray;
// store the lightray position and direction
meas.setLightRayPosition( lightray.point() );
meas.setLightRayDirection( lightray.direction() );
delete[] interslc;
}
Definition at line 52 of file OptOSensor2D.h.
ALIbool OptOSensor2D::fdevi_from_file [private] |
Definition at line 53 of file OptOSensor2D.h.
1.7.3