ATLAS Offline Software
Loading...
Searching...
No Matches
Functions
TRTCalibrator.cxx File Reference
#include "TRTCalibrator.h"
#include "Calibrator.h"
#include "TRT_CalibData/TrackInfo.h"
#include "TrkTrack/Track.h"
#include "TrkParameters/TrackParameters.h"
#include "TrkEventPrimitives/LocalParameters.h"
#include "EventPrimitives/EventPrimitives.h"
#include "TrkRIO_OnTrack/RIO_OnTrack.h"
#include "InDetRIO_OnTrack/TRT_DriftCircleOnTrack.h"
#include "InDetPrepRawData/TRT_DriftCircle.h"
#include "TrkPrepRawData/PrepRawData.h"
#include "AtlasDetDescr/AtlasDetectorID.h"
#include "InDetIdentifier/TRT_ID.h"
#include "TRT_ConditionsData/RtRelation.h"
#include "TRT_ConditionsData/BasicRtRelation.h"
#include "TRT_ConditionsData/DinesRtRelation.h"
#include "TRT_ConditionsData/BinnedRtRelation.h"
#include "CommissionEvent/ComTime.h"
#include "TRT_ReadoutGeometry/TRT_DetectorManager.h"
#include <TNtuple.h>
#include <TFile.h>
#include <fstream>

Go to the source code of this file.

Functions

bool TRTCalibrator::calibrate ATLAS_NOT_THREAD_SAFE ()
 Install fatal handler with default options.

Function Documentation

◆ ATLAS_NOT_THREAD_SAFE()

bool TRTCalibrator::calibrate ATLAS_NOT_THREAD_SAFE ( )
inline

Install fatal handler with default options.

This is meant to be easy to call from python via ctypes.

Definition at line 309 of file TRTCalibrator.cxx.

309 {
310
311 float sid;
312 Identifier ident;
313 const InDetDD::TRT_BaseElement* strawelement;
314 const InDetDD::TRT_BarrelElement* barrelelement;
315 const InDetDD::TRT_EndcapElement* endcapelement;
316 int nTRThist=0, ndethist=0, nlayhist=0, nmodhist=0, nbrdhist=0, nchphist=0, nstwhist=0;
317 int nTRThistAr=0, ndethistAr=0, nlayhistAr=0, nmodhistAr=0, nbrdhistAr=0, nchphistAr=0, nstwhistAr=0;
318 databundle hitdata;
319
320 std::string infile=m_hittuple;
321 std::string outfile="calibout";
322
323 std::string detname="Detector";
324 std::string layname="Layer";
325 if (!m_SplitBarrel) {
326 detname="WholeBarrel";
327 layname="WholeBarrelLayer";
328 }
329
330 //create Calibrator objects for each sub-level
331 Calibrator TRT(0,"TRT",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
332 Calibrator Detector(1,detname.data(),m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
333 Calibrator Layer(2,layname.data(),m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
334 Calibrator Module(3,"Module",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
335 Calibrator Board(4,"Board",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
336 Calibrator Chip(5,"Chip",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
337 Calibrator Straw(6,"Straw",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
338
339
340 // make a map of calibrators (so that they can be looped over)
341 std::map<std::string,Calibrator*> calibrators;
342 calibrators["TRT"]=&TRT;
343 calibrators["Detector"]=&Detector;
344 calibrators["Layer"]=&Layer;
345 calibrators["Module"]=&Module;
346 calibrators["Board"]=&Board;
347 calibrators["Chip"]=&Chip;
348 calibrators["Straw"]=&Straw;
349
350 //use default configurations unless "user" option is used
351 if(m_calsub!="user"){
352 m_doRt=m_config[m_calsub].CalibrateRt;
353 m_doT0=m_config[m_calsub].CalibrateT0;
354 m_doRes=m_config[m_calsub].FitResidual;
355 m_beQuiet=m_config[m_calsub].NoHistograms;
356 m_doOutPrint=m_config[m_calsub].PrintT0Out;
357 m_doRtPrint=m_config[m_calsub].PrintRtOut;
358 m_doLogPrint=m_config[m_calsub].PrintLog;
359 }
360
361 //configure the calibrators
362 for (const std::string& s : m_doRt) calibrators[s]->dort=true;
363 for (const std::string& s : m_doT0) calibrators[s]->dot0=true;
364 for (const std::string& s : m_doRes) calibrators[s]->dores=true;
365 for (const std::string& s : m_beQuiet) calibrators[s]->bequiet=true;
366 for (const std::string& s : m_doOutPrint) calibrators[s]->printt0=true;
367 for (const std::string& s : m_doRtPrint) calibrators[s]->printrt=true;
368 for (const std::string& s : m_doLogPrint) calibrators[s]->printlog=true;
369 for (const std::string& s : m_useBoardRef) calibrators[s]->usebref=true;
370 for (std::pair<const std::string, Calibrator*>& p : calibrators) {
371 Calibrator* calib = p.second;
372 calib->usep0=m_useP0;
373 calib->floatp3=m_floatP3;
374 calib->useshortstraws=m_DoShortStrawCorrection;
375
376 if (m_doRt.size()==0) calib->nort=true;
377 if (m_doT0.size()==0) calib->not0=true;
378
379 std::set<int> selset;
380 GetSubLevels(m_selstring,calib->level+1,&selset); //get the set of selections for level "level+1" from the m_selstring ...
381 for (std::set<int>::iterator imod=selset.begin(); imod!=selset.end(); ++imod){
382 calib->selection = selset;//... and configure the calibrator with them
383 }
384 }
385
386 // FOR ARGON:
387 Calibrator TRT_Ar (0,"TRT_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
388 Calibrator Detector_Ar(1,"Detector_Ar",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
389 Calibrator Layer_Ar (2,"Layer_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
390 Calibrator Module_Ar (3,"Module_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
391 Calibrator Board_Ar (4,"Board_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
392 Calibrator Chip_Ar (5,"Chip_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
393 Calibrator Straw_Ar (6,"Straw_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
394
395 // make a map of calibrators (so that they can be looped over)
396 std::map<std::string,Calibrator*> calibratorsAr;
397 calibratorsAr["TRT"] =&TRT_Ar;
398 calibratorsAr["Detector"] =&Detector_Ar;
399 calibratorsAr["Layer"] =&Layer_Ar;
400 calibratorsAr["Module"] =&Module_Ar;
401 calibratorsAr["Board"] =&Board_Ar;
402 calibratorsAr["Chip"] =&Chip_Ar;
403 calibratorsAr["Straw"] =&Straw_Ar;
404
405 //configure the calibrators
406 for (const std::string& s : m_doRt) calibratorsAr[s]->dort=true;
407 for (const std::string& s : m_doT0) calibratorsAr[s]->dot0=true;
408 for (const std::string& s : m_doRes) calibratorsAr[s]->dores=true;
409 for (const std::string& s : m_beQuiet) calibratorsAr[s]->bequiet=true;
410 for (const std::string& s : m_doOutPrint) calibratorsAr[s]->printt0=true;
411 for (const std::string& s : m_doRtPrint) calibratorsAr[s]->printrt=true;
412 for (const std::string& s : m_doLogPrint) calibratorsAr[s]->printlog=true;
413 for (const std::string& s : m_useBoardRef) calibratorsAr[s]->usebref=true;
414 for (std::pair<const std::string, Calibrator*>& p : calibratorsAr) {
415 Calibrator* calib = p.second;
416 calib->usep0=m_useP0;
417 calib->floatp3=m_floatP3;
418 calib->useshortstraws=m_DoShortStrawCorrection;
419
420 if (m_doRt.size()==0) calib->nort=true;
421 if (m_doT0.size()==0) calib->not0=true;
422
423 std::set<int> selset;
424 GetSubLevels(m_selstring,calib->level+1,&selset); //get the set of selections for level "level+1" from the m_selstring ...
425 for (std::set<int>::iterator imod=selset.begin(); imod!=selset.end(); ++imod){
426 calib->selection = selset;//... and configure the calibrator with them
427 }
428 }
429
430
431
432 m_options = '_' + TRT.GetOptString() + '_' + Detector.GetOptString() + '_' + Layer.GetOptString() + '_' + Module.GetOptString() + '_' + Board.GetOptString() + '_' + Chip.GetOptString() + '_' + Straw.GetOptString();
433 if(m_calsub!="user") m_selstring = '_' + TRT.GetSelString() + '_' + Detector.GetSelString() + '_' + Layer.GetSelString() + '_' + Module.GetSelString() + '_' + Board.GetSelString() + '_' + Chip.GetSelString() + '_' + Straw.GetSelString();
434
435 //print some info
436 ATH_MSG_INFO(" ");
437 ATH_MSG_INFO( "INPUT FILE : " << infile );
438 ATH_MSG_INFO( "OUTPUT FILE : " << outfile );
439 ATH_MSG_INFO( "SELECTION STRING : " << m_selstring );
440 ATH_MSG_INFO( "OPTION STRING : " << m_options );
441 ATH_MSG_INFO( "RT RELATION : " << m_rtrel );
442 ATH_MSG_INFO( "MIN STATISTICS : RT=" << m_minrt << ", T0=" << m_mint0 );
443 ATH_MSG_INFO( "T0 OFFSET : " << m_t0offset );
444 ATH_MSG_INFO(" ");
445 ATH_MSG_INFO( TRT.PrintInfo() );
446 ATH_MSG_INFO( Detector.PrintInfo() );
447 ATH_MSG_INFO( Layer.PrintInfo() );
448 ATH_MSG_INFO( Module.PrintInfo() );
449 ATH_MSG_INFO( Board.PrintInfo() );
450 ATH_MSG_INFO( Chip.PrintInfo() );
451 ATH_MSG_INFO( Straw.PrintInfo() );
452 ATH_MSG_INFO(" ");
453
454 if (m_DoArXenonSep){
455 //print some info
456 ATH_MSG_INFO(" ");
457 ATH_MSG_INFO( "Setup for the ARGON Straws: " );
458 ATH_MSG_INFO( "INPUT FILE : " << infile );
459 ATH_MSG_INFO( "OUTPUT FILE : " << outfile );
460 ATH_MSG_INFO( "SELECTION STRING : " << m_selstring );
461 ATH_MSG_INFO( "OPTION STRING : " << m_options );
462 ATH_MSG_INFO( "RT RELATION : " << m_rtrel );
463 ATH_MSG_INFO( "MIN STATISTICS : RT=" << m_minrt << ", T0=" << m_mint0 );
464 ATH_MSG_INFO( "T0 OFFSET : " << m_t0offset );
465 ATH_MSG_INFO( " ");
466 ATH_MSG_INFO( TRT_Ar.PrintInfo() );
467 ATH_MSG_INFO( Detector_Ar.PrintInfo() );
468 ATH_MSG_INFO( Layer_Ar.PrintInfo() );
469 ATH_MSG_INFO( Module_Ar.PrintInfo() );
470 ATH_MSG_INFO( Board_Ar.PrintInfo() );
471 ATH_MSG_INFO( Chip_Ar.PrintInfo() );
472 ATH_MSG_INFO( Straw_Ar.PrintInfo() );
473 ATH_MSG_INFO(" ");
474 }
475 //read the chip reference t0 values from finedelays.txt
476 int rbrd, rchp, rdet, dum;
477 float rt0;
478 std::map<std::string,float> reft0map;
479 std::ifstream t0ref("finedelays.txt",std::ios::in);
480 if(t0ref.is_open()) ATH_MSG_INFO(" Opened finedelays.txt ");
481 for(int iref=0;iref<208;iref++){
482 t0ref >> rbrd >> rchp >> rdet >> dum >> rt0;
483 reft0map[std::string(Form("_%i_%i_%i",rdet,rbrd,rchp))]=rt0;
484 }
485
486 t0ref.close();
487
488
489 //set type of rt-relation
490 bool isdines = (m_rtrel.value()).find("dines")!=std::string::npos;
491 bool isbinned = (m_rtrel.value()).find("binned")!=std::string::npos;
492 int rtint;
493 if (isdines) rtint=2;
494 else if (isbinned) rtint=1;
495 else rtint=0;
496 ATH_MSG_INFO(" Rt relation is: " << rtint);
497
498 int npop,isid;
499 std::map<std::string,epdata> ephasemap;
500
501 // open the output histogram file
502 std::unique_ptr<TFile> histfile(TFile::Open("calibout.root","RECREATE"));
503
504 if(!histfile || histfile->IsZombie()) {
505 ATH_MSG_ERROR ("Failed to open calibout.root ");
506 return false;
507 } else {
508 ATH_MSG_INFO ("Opened calibout.root ");
509 }
510
511 std::ifstream myFile (infile.data(), std::ios::in | std::ios::binary);
512 ATH_MSG_INFO( " Opened " << infile << " as binary histogram file " );
513 int ihist=0;
514 int ihistAr=0;
515
516 while(true){
517
518 //read a binary histogram
519 //dangerous cast reproduces pre-existing c style cast
520 myFile.read (reinterpret_cast<char*>(&npop),sizeof(int)); //number of populated bins
521 if (myFile.eof()) break;
522 std::vector<int> chist(2 * npop + 2);
523 //dangerous cast reproduces pre-existing c style cast
524 if (npop>0) myFile.read (reinterpret_cast<char*>(chist.data()+2), sizeof(int)*2*npop);
525 //dangerous cast reproduces pre-existing c style cast
526 myFile.read (reinterpret_cast<char*>(&isid),sizeof(int)); //the straw id
527 sid = (float)isid;
528 if(sid<0) continue;
529
530 chist[0]=npop;
531 chist[1]=isid;
532 //get the straw address (barrel/ec, layer, module, ...) based on the straw identifier
533 ident=(Identifier)(isid);
534 int chip=0,board=-1;
535 m_neighbourSvc->getChip(ident,chip);
536 if(abs(m_TRTID->barrel_ec(ident))<2){
537 board=m_neighbourSvc->chipToBoardBarrel(chip,m_TRTID->layer_or_wheel(ident));
538 }
539 else if (chip<12) {
540 board=0;
541 }
542 else {
543 chip=chip-20;
544 board=1;
545 }
546 if (m_SplitBarrel) hitdata.det=(int)m_TRTID->barrel_ec(ident);
547 else hitdata.det=abs((int)m_TRTID->barrel_ec(ident));
548 if(hitdata.det!=1 && hitdata.det!=-1 && hitdata.det!=2 && hitdata.det!=-2) {
549 ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " detector decoded to be : " << hitdata.det);
550 continue;
551 }
552 hitdata.lay=(int)m_TRTID->layer_or_wheel(ident);
553 if( (hitdata.det==1 || hitdata.det==-1) && (hitdata.lay<0 || hitdata.lay>2) ) {
554 ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " barrel layer decoded to be : " << hitdata.lay);
555 continue;
556 } else if( (hitdata.det==2 || hitdata.det==-2) && (hitdata.lay<0 || hitdata.lay>13) ) {
557 ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " endcap layer decoded to be : " << hitdata.lay);
558 continue;
559 }
560 hitdata.mod=(int)m_TRTID->phi_module(ident);
561 if( hitdata.mod < 0 || hitdata.mod > 31 ) {
562 ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " phi-module decoded to be : " << hitdata.lay);
563 continue;
564 }
565
566 hitdata.brd=board;
567 hitdata.chp=chip;
568 hitdata.stl=(int)m_TRTID->straw_layer(ident);
569 if( hitdata.stl < 0 || hitdata.stl > m_TRTID->straw_layer_max(ident) ) {
570 ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " straw-layer decoded to be : " << hitdata.stl);
571 continue;
572 }
573
574 hitdata.stw=(int)m_TRTID->straw(ident);
575 if( hitdata.stw < 0 || hitdata.stw > m_TRTID->straw_max(ident) ) {
576 ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " straw decoded to be : " << hitdata.stw);
577 continue;
578 }
579
580 hitdata.sid=isid;
581
582 //get the old rt parameters based on straw identifier
583 const float* pcal ;
584 std::vector<float> rvalues ;
585 const float defaultpcal[] = {0,0,0,0} ;
586 if (isdines){
587 const TRTCond::DinesRtRelation* rtr = dynamic_cast<const TRTCond::DinesRtRelation*>(m_trtcaldbTool->getRtRelation(ident)) ;
588 pcal = rtr ? rtr->cal() : defaultpcal ;
589 }
590
591 else {
592 const TRTCond::BasicRtRelation* rtr = dynamic_cast<const TRTCond::BasicRtRelation*>(m_trtcaldbTool->getRtRelation(ident)) ;
593 pcal = rtr ? rtr->cal() : defaultpcal ;
594 }
595 hitdata.rtpar[0]=pcal[0];
596 hitdata.rtpar[1]=pcal[1];
597 hitdata.rtpar[2]=pcal[2];
598 hitdata.rtpar[3]=pcal[3];
599
600
601 //build map keys
602 MakeBDKeys(hitdata.det, hitdata.lay, hitdata.mod, hitdata.brd, hitdata.chp, hitdata.sid);
603
604 //make the level hierachy dictionaries
605 m_trt.t[m_Tkey].d[m_Dkey].l[m_Lkey].m[m_Mkey].b[m_Bkey].c[m_Ckey].s[m_Skey].z=0;
606 m_trt_acc.t[m_Tkey].d[m_Dkey_acc].l[m_Lkey_acc].m[m_Mkey].b[m_Bkey].c[m_Ckey].s[m_Skey].z=0;
607
608 //populate the hit data structure to be added
609
610
611 hitdata.ievt=0;
612 hitdata.tres=0;
613 hitdata.weight=1.0;
614 hitdata.res=0;
615 hitdata.t=0;
616 hitdata.r=0;
617 hitdata.t0=m_trtcaldbTool->getT0(ident);
618
619 strawelement = m_trtmanager->getElement(ident);
620 if(hitdata.det==1 || hitdata.det==-1) {
621 barrelelement=static_cast<const InDetDD::TRT_BarrelElement*>(strawelement);
622 hitdata.x=(barrelelement->center(ident)).x();
623 hitdata.y=(barrelelement->center(ident)).y();
624 hitdata.z=(barrelelement->center(ident)).z();
625 }else{
626 endcapelement=static_cast<const InDetDD::TRT_EndcapElement*>(strawelement);
627 hitdata.x=(endcapelement->center(ident)).x();
628 hitdata.y=(endcapelement->center(ident)).y();
629 hitdata.z=(endcapelement->center(ident)).z();
630 }
631
632 //in the short straws corrections, autodetect if it was applied on the previous step
633 if ( m_DoShortStrawCorrection && hitdata.lay==0 && hitdata.stl<9){
634 //If correction was done in ctes in db (readed), undo the correction: This fixes problem on t0 averaging....
635 double t0test1=m_trtcaldbTool->getT0((Identifier)301998432);
636 double t0test2=m_trtcaldbTool->getT0((Identifier)302001504);
637 if (t0test1 != t0test2) hitdata.t0+=0.75; //short straw compensation
638 }
639
640 //the reference t0 from finedelays
641 hitdata.rt0=reft0map[std::string(Form("_%i_%i_%i",hitdata.det,hitdata.brd,hitdata.chp))];
642
643 // Prepare for Xe-Ar mixed conditions:
644 int isArgonStraw = 0;
645 if (m_TRTStrawSummaryTool->getStatusHT(ident, Gaudi::Hive::currentContext()) != TRTCond::StrawStatus::Good) {
646 isArgonStraw = 1;
647 }
648
649
650 //add histogram to the Calibrators (A and C side separated)
651 if(!m_DoArXenonSep){
652
653 nTRThist += TRT.AddHit(m_Tkey,hitdata,chist.data(),true);
654 if (Detector.CheckSelection(hitdata.det)) { //only add the histogram if it is in the selection
655 if (m_SplitBarrel) ndethist += Detector.AddHit(m_Dkey,hitdata,chist.data(),true);
656 else ndethist += Detector.AddHit(m_Dkey_acc,hitdata,chist.data(),true);
657 if (Layer.CheckSelection(hitdata.lay)) {
658 if (m_SplitBarrel) nlayhist += Layer.AddHit(m_Lkey,hitdata,chist.data(),true);
659 else { nlayhist += Layer.AddHit(m_Lkey_acc,hitdata,chist.data(),true); continue;}
660 if (Module.CheckSelection(hitdata.mod)) {
661 nmodhist += Module.AddHit(m_Mkey,hitdata,chist.data(),true);
662 if (Board.CheckSelection(hitdata.brd)) {
663 nbrdhist += Board.AddHit(m_Bkey,hitdata,chist.data(),true);
664 if (Chip.CheckSelection(hitdata.chp)) {
665 nchphist += Chip.AddHit(m_Ckey,hitdata,chist.data(),true);
666 if (Straw.CheckSelection(hitdata.stw))
667 nstwhist += Straw.AddHit(m_Skey,hitdata,chist.data(),true);
668 }
669 }
670 }
671 }
672 } // Here it closes
673 ihist++;
674 } else
675 if(isArgonStraw==0){ // Separate Ar and Xe in endcaps. Here Xe
676 nTRThist += TRT.AddHit(m_Tkey,hitdata,chist.data(),true);
677 if (Detector.CheckSelection(hitdata.det)) { //only add the histogram if it is in the selection
678 if (m_SplitBarrel) ndethist += Detector.AddHit(m_Dkey,hitdata,chist.data(),true);
679 else ndethist += Detector.AddHit(m_Dkey_acc,hitdata,chist.data(),true);
680 if (Layer.CheckSelection(hitdata.lay)) {
681 if (m_SplitBarrel) nlayhist += Layer.AddHit(m_Lkey,hitdata,chist.data(),true);
682 else { nlayhist += Layer.AddHit(m_Lkey_acc,hitdata,chist.data(),true); continue;}
683 if (Module.CheckSelection(hitdata.mod)) {
684 nmodhist += Module.AddHit(m_Mkey,hitdata,chist.data(),true);
685 if (Board.CheckSelection(hitdata.brd)) {
686 nbrdhist += Board.AddHit(m_Bkey,hitdata,chist.data(),true);
687 if (Chip.CheckSelection(hitdata.chp)) {
688 nchphist += Chip.AddHit(m_Ckey,hitdata,chist.data(),true);
689 if (Straw.CheckSelection(hitdata.stw)) {
690 nstwhist += Straw.AddHit(m_Skey,hitdata,chist.data(),true);
691 }
692 }
693 }
694 }
695 }
696 } // Here it closes
697 ihist++;
698 } else { // ARGON HITS
699 nTRThistAr += TRT_Ar.AddHit(m_Tkey,hitdata,chist.data(),true);
700 if (Detector_Ar.CheckSelection(hitdata.det)) { //only add the histogram if it is in the selection
701 if (m_SplitBarrel) ndethistAr += Detector_Ar.AddHit(m_Dkey,hitdata,chist.data(),true);
702 else ndethistAr += Detector_Ar.AddHit(m_Dkey_acc,hitdata,chist.data(),true);
703 if (Layer_Ar.CheckSelection(hitdata.lay)) {
704 if (m_SplitBarrel) nlayhistAr += Layer_Ar.AddHit(m_Lkey,hitdata,chist.data(),true);
705 else { nlayhistAr += Layer_Ar.AddHit(m_Lkey_acc,hitdata,chist.data(),true); continue;}
706 if (Module_Ar.CheckSelection(hitdata.mod)) {
707 nmodhistAr += Module_Ar.AddHit(m_Mkey,hitdata,chist.data(),true);
708 if (Board_Ar.CheckSelection(hitdata.brd)) {
709 nbrdhistAr += Board_Ar.AddHit(m_Bkey,hitdata,chist.data(),true);
710 if (Chip_Ar.CheckSelection(hitdata.chp)) {
711 nchphistAr += Chip_Ar.AddHit(m_Ckey,hitdata,chist.data(),true);
712 if (Straw_Ar.CheckSelection(hitdata.stw)) {
713 nstwhistAr += Straw_Ar.AddHit(m_Skey,hitdata,chist.data(),true);
714 }
715 }
716 }
717 }
718 }
719 } // Here it closes
720 ihistAr++;
721 } // END ARGON HITS
722
723
724 if ((ihist%10000==9999) | (ihist==m_nevents-1)){
725 ATH_MSG_INFO( Form("%7i HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihist+1,nTRThist, ndethist, nlayhist, nmodhist, nbrdhist, nchphist, nstwhist) );
726 if(m_DoArXenonSep) ATH_MSG_INFO( Form("%7i Ar HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihistAr+1,nTRThistAr, ndethistAr, nlayhistAr, nmodhistAr, nbrdhistAr, nchphistAr, nstwhistAr) );
727 }
728
729 ihist++;
730 } // Finish Straw Loop
731
732 ATH_MSG_INFO( Form("%7i HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihist+1,nTRThist, ndethist, nlayhist, nmodhist, nbrdhist, nchphist, nstwhist) );
733 if(m_DoArXenonSep) ATH_MSG_INFO( Form("%7i Ar HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihistAr+1,nTRThistAr, ndethistAr, nlayhistAr, nmodhistAr, nbrdhistAr, nchphistAr, nstwhistAr) );
734
735 TRT.UpdateOldConstants();
736 Detector.UpdateOldConstants();
737 Layer.UpdateOldConstants();
738 Module.UpdateOldConstants();
739 Board.UpdateOldConstants();
740 Chip.UpdateOldConstants();
741
742 if (ihistAr>0 ){
743 TRT_Ar.UpdateOldConstants();
744 Detector_Ar.UpdateOldConstants();
745 Layer_Ar.UpdateOldConstants();
746 Module_Ar.UpdateOldConstants();
747 Board_Ar.UpdateOldConstants();
748 Chip_Ar.UpdateOldConstants();
749 }
750
751 myFile.close();
752
753
754 ATH_MSG_INFO(" ");
755 ATH_MSG_INFO( TRT.PrintStat() );
756 ATH_MSG_INFO( Detector.PrintStat() );
757 ATH_MSG_INFO( Layer.PrintStat() );
758 ATH_MSG_INFO( Module.PrintStat() );
759 ATH_MSG_INFO( Board.PrintStat() );
760 ATH_MSG_INFO( Chip.PrintStat() );
761 ATH_MSG_INFO( Straw.PrintStat() );
762 ATH_MSG_INFO(" ");
763
764
765 if(m_DoArXenonSep){
766 ATH_MSG_INFO(" ");
767 ATH_MSG_INFO( TRT_Ar.PrintStat() );
768 ATH_MSG_INFO( Detector_Ar.PrintStat() );
769 ATH_MSG_INFO( Layer_Ar.PrintStat() );
770 ATH_MSG_INFO( Module_Ar.PrintStat() );
771 ATH_MSG_INFO( Board_Ar.PrintStat() );
772 ATH_MSG_INFO( Chip_Ar.PrintStat() );
773 ATH_MSG_INFO( Straw_Ar.PrintStat() );
774 ATH_MSG_INFO(" ");
775 }
776
777
778
779 //caldata startdata(true,tbins,rbins);
780 caldata startdata(true,55,100);
781 //caldata startdata(true,64,64);
782 startdata.t0=5.0;
783 std::map<std::string,TDirectory*> dirmap;
784 TDirectory* trtdir=gDirectory;
785 TDirectory* detdir=gDirectory;
786 TDirectory* laydir=gDirectory;
787 TDirectory* moddir=gDirectory;
788 TDirectory* brddir=gDirectory;
789 TDirectory* chpdir=gDirectory;
790
791 std::ofstream rtcalfile(Form("%s_rt.txt",outfile.data()),std::ios::out);
792 std::ofstream binrtcalfile(Form("%s_binrt.txt",outfile.data()),std::ios::out);
793 std::ofstream t0calfile(Form("%s_t0.txt",outfile.data()),std::ios::out);
794
795
796 //if (m_SplitBarrel){
797
798 for (std::pair<const std::string, BDdetector>& pt : m_trt.t) {
799
800 if (TRT.Skip()) break;
801 if (TRT.HasKey(pt.first)) {
802 trtdir = TRT.Calibrate(histfile.get(),pt.first,SubLev(m_options,1),&startdata);
803 if (TRT.printt0) t0calfile << Form("-3 -1 -1 -1 -1 : %e %e",TRT.data[pt.first].t0,TRT.data[pt.first].t0err) << '\n';
804 if (TRT.printrt) rtcalfile << Form("-3 -1 -1 -1 -1 : %i %e %e %e %e",rtint,TRT.data[pt.first].rtpar[0],TRT.data[pt.first].rtpar[1],TRT.data[pt.first].rtpar[2],TRT.data[pt.first].rtpar[3]) << '\n';
805 }
806 for (std::pair<const std::string, BDlayer>& pd : pt.second.d) {
807
808 if(Detector.Skip()) break;
809 if(Detector.HasKey(pd.first)){
810 detdir = Detector.Calibrate(trtdir,pd.first,SubLev(m_options,2),&TRT.data[pt.first]);
811 }
812 for (std::pair<const std::string, BDmodule>& pl : pd.second.l) {
813
814 if(Layer.Skip()) break;
815 if(Layer.HasKey(pl.first)){
816 laydir = Layer.Calibrate(detdir,pl.first,SubLev(m_options,3),&Detector.data[pd.first]);
817 if (Layer.printt0) t0calfile << Form("%i %i -1 -1 -1 : %e %e",Layer.data[pl.first].det,Layer.data[pl.first].lay,Layer.data[pl.first].t0,Layer.data[pl.first].t0err) << '\n';
818 if (Layer.printrt) rtcalfile << Form("%i %i -1 -1 -1 : %i %e %e %e %e",Layer.data[pl.first].det,Layer.data[pl.first].lay,rtint,Layer.data[pl.first].rtpar[0],Layer.data[pl.first].rtpar[1],Layer.data[pl.first].rtpar[2],Layer.data[pl.first].rtpar[3]) << '\n';
819 if (!m_SplitBarrel) {
820 if (Layer.printt0) t0calfile << Form("%i %i -1 -1 -1 : %e %e",-Layer.data[pl.first].det,Layer.data[pl.first].lay,Layer.data[pl.first].t0,Layer.data[pl.first].t0err) << '\n';
821 if (Layer.printrt) rtcalfile << Form("%i %i -1 -1 -1 : %i %e %e %e %e",-Layer.data[pl.first].det,Layer.data[pl.first].lay,rtint,Layer.data[pl.first].rtpar[0],Layer.data[pl.first].rtpar[1],Layer.data[pl.first].rtpar[2],Layer.data[pl.first].rtpar[3]) << '\n';
822 }
823 }
824 for (std::pair<const std::string, BDboard>& pm : pl.second.m) {
825
826 if(Module.Skip()) break;
827 if(Module.HasKey(pm.first)){
828 moddir = Module.Calibrate(laydir,pm.first,SubLev(m_options,4),&Layer.data[pl.first]);
829 if (Module.printt0) t0calfile << Form("%i %i %i -1 -1 : %e %e",Module.data[pm.first].det,Module.data[pm.first].lay,Module.data[pm.first].mod,Module.data[pm.first].t0,Module.data[pm.first].t0err) << '\n';
830 if (Module.printrt) rtcalfile << Form("%i %i %i -1 -1 : %i %e %e %e %e",Module.data[pm.first].det,Module.data[pm.first].lay,Module.data[pm.first].mod,rtint,Module.data[pm.first].rtpar[0],Module.data[pm.first].rtpar[1],Module.data[pm.first].rtpar[2],Module.data[pm.first].rtpar[3]) << '\n';
831 }
832 for (std::pair<const std::string, BDchip>& pb : pm.second.b) {
833
834 if(Board.Skip()) break;
835 if(Board.HasKey(pb.first)){
836 brddir = Board.Calibrate(moddir,pb.first,SubLev(m_options,5),&Module.data[pm.first]);
837 }
838 for (std::pair<const std::string, BDstraw>& pc : pb.second.c) {
839
840 if(Chip.Skip()) break;
841 if(Chip.HasKey(pc.first)){
842 chpdir = Chip.Calibrate(brddir,pc.first,SubLev(m_options,6),&Board.data[pb.first]);
843 }
844 for (std::pair<const std::string, BDzero>& ps : pc.second.s) {
845
846 if(Straw.Skip()) break;
847 if(Straw.HasKey(ps.first)){
848 Straw.Calibrate(chpdir,ps.first,SubLev(m_options,7),&Chip.data[pc.first]);
849 if (Straw.printt0) t0calfile << Form("%i %i %i %i %i : %e %e",Straw.data[ps.first].det,Straw.data[ps.first].lay,Straw.data[ps.first].mod,Straw.data[ps.first].stl,Straw.data[ps.first].stw,Straw.data[ps.first].t0,Straw.data[ps.first].t0err) << '\n';
850 if (Straw.printrt) rtcalfile << Form("%i %i %i %i %i : %i %e %e %e %e",Straw.data[ps.first].det,Straw.data[ps.first].lay,Straw.data[ps.first].mod,Straw.data[ps.first].stl,Straw.data[ps.first].stw,rtint,Straw.data[ps.first].rtpar[0],Straw.data[ps.first].rtpar[1],Straw.data[ps.first].rtpar[2],Straw.data[ps.first].rtpar[3]) << '\n';
851
852
853 }
854 }
855 }
856 }
857 }
858 }
859 }
860 }
861
862 std::map<std::string,TDirectory*> dirArmap;
863 TDirectory* trtdirAr=gDirectory;
864 TDirectory* detdirAr=gDirectory;
865 TDirectory* laydirAr=gDirectory;
866 TDirectory* moddirAr=gDirectory;
867 TDirectory* brddirAr=gDirectory;
868 TDirectory* chpdirAr=gDirectory;
869
870 if (m_DoArXenonSep){
871
872 for (std::pair<const std::string, BDdetector>& pt : m_trt.t) {
873
874 if (TRT_Ar.Skip()) break;
875 if (TRT_Ar.HasKey(pt.first)) {
876 trtdirAr = TRT_Ar.Calibrate(histfile.get(),pt.first,SubLev(m_options,1),&startdata);
877 }
878 for (std::pair<const std::string, BDlayer>& pd : pt.second.d) {
879
880 if(Detector_Ar.Skip()) break;
881 if(Detector_Ar.HasKey(pd.first)){
882 detdirAr = Detector_Ar.Calibrate(trtdirAr,pd.first,SubLev(m_options,2),&TRT_Ar.data[pt.first]);
883 }
884 for (std::pair<const std::string, BDmodule>& pl : pd.second.l) {
885
886 if(Layer_Ar.Skip()) break;
887 if(Layer_Ar.HasKey(pl.first)){
888 laydirAr = Layer_Ar.Calibrate(detdirAr,pl.first,SubLev(m_options,3),&Detector_Ar.data[pd.first]);
889
890 }
891 for (std::pair<const std::string, BDboard>& pm : pl.second.m) {
892
893 if(Module_Ar.Skip()) break;
894 if(Module_Ar.HasKey(pm.first)){
895 moddirAr = Module_Ar.Calibrate(laydirAr,pm.first,SubLev(m_options,4),&Layer_Ar.data[pl.first]);
896 if (Module_Ar.printt0) t0calfile << Form("%i %i %i -1 -1 : %e %e",Module_Ar.data[pm.first].det,Module_Ar.data[pm.first].lay,Module_Ar.data[pm.first].mod,Module_Ar.data[pm.first].t0,Module_Ar.data[pm.first].t0err) << '\n';
897 if (Layer_Ar.printrt) rtcalfile << Form("%i %i %i -1 -1 : %i %e %e %e %e",Module_Ar.data[pm.first].det,Module_Ar.data[pm.first].lay,Module_Ar.data[pm.first].mod,rtint,Module_Ar.data[pm.first].rtpar[0],Module_Ar.data[pm.first].rtpar[1],Module_Ar.data[pm.first].rtpar[2],Module_Ar.data[pm.first].rtpar[3]) << '\n';
898 }
899 for (std::pair<const std::string, BDchip>& pb : pm.second.b) {
900 if(Board_Ar.Skip()) break;
901 if(Board_Ar.HasKey(pb.first)){
902 brddirAr = Board_Ar.Calibrate(moddirAr,pb.first,SubLev(m_options,5),&Module_Ar.data[pm.first]);
903 }
904 for (std::pair<const std::string, BDstraw>& pc : pb.second.c) {
905
906 if(Chip_Ar.Skip()) break;
907 if(Chip_Ar.HasKey(pc.first)){
908 chpdirAr = Chip_Ar.Calibrate(brddirAr,pc.first,SubLev(m_options,6),&Board_Ar.data[pb.first]);
909 }
910 for (std::pair<const std::string, BDzero>& ps : pc.second.s) {
911
912 if(Straw_Ar.Skip()) break;
913 if(Straw_Ar.HasKey(ps.first)){
914 Straw_Ar.Calibrate(chpdirAr,ps.first,SubLev(m_options,7),&Chip_Ar.data[pc.first]);
915 if (Straw_Ar.printt0) t0calfile << Form("%i %i %i %i %i : %e %e",Straw_Ar.data[ps.first].det,Straw_Ar.data[ps.first].lay,Straw_Ar.data[ps.first].mod,Straw_Ar.data[ps.first].stl,Straw_Ar.data[ps.first].stw,Straw_Ar.data[ps.first].t0,Straw_Ar.data[ps.first].t0err) << '\n';
916 if (Straw_Ar.printrt) rtcalfile << Form("%i %i %i %i %i : %i %e %e %e %e",Straw_Ar.data[ps.first].det,Straw_Ar.data[ps.first].lay,Straw_Ar.data[ps.first].mod,Straw_Ar.data[ps.first].stl,Straw_Ar.data[ps.first].stw,rtint,Straw_Ar.data[ps.first].rtpar[0],Straw_Ar.data[ps.first].rtpar[1],Straw_Ar.data[ps.first].rtpar[2],Straw_Ar.data[ps.first].rtpar[3]) << '\n';
917 }
918 }
919 }
920 }
921 }
922 }
923 }
924 }
925
926}// Close for ARGON
927
928
929
930
931 //make ntuples for each calibrator
932 if (!TRT.bequiet) TRT.WriteStat(histfile.get());
933 if (!Detector.bequiet) Detector.WriteStat(histfile.get());
934 if (!Layer.bequiet) Layer.WriteStat(histfile.get());
935 if (!Module.bequiet) Module.WriteStat(histfile.get());
936 if (!Board.bequiet) Board.WriteStat(histfile.get());
937 Chip.WriteStat(histfile.get());
938 Straw.WriteStat(histfile.get());
939
940 TRT.DumpConstants();
941 Detector.DumpConstants();
942 Layer.DumpConstants();
943 Module.DumpConstants();
944 Board.DumpConstants();
945 Chip.DumpConstants();
946 Straw.DumpConstants();
947
948 rtcalfile.close();
949 binrtcalfile.close();
950 t0calfile.close();
951
952 //AR
953 //make ntuples for each calibrator
954 if(m_DoArXenonSep){
955 if (!TRT_Ar.bequiet) TRT_Ar.WriteStat(histfile.get());
956 if (!Detector_Ar.bequiet) Detector_Ar.WriteStat(histfile.get());
957 if (!Layer_Ar.bequiet) Layer_Ar.WriteStat(histfile.get());
958 if (!Module_Ar.bequiet) Module_Ar.WriteStat(histfile.get());
959 if (!Board_Ar.bequiet) Board_Ar.WriteStat(histfile.get());
960 Chip_Ar.WriteStat(histfile.get());
961 Straw_Ar.WriteStat(histfile.get());
962
963 TRT_Ar.DumpConstants();
964 Detector_Ar.DumpConstants();
965 Layer_Ar.DumpConstants();
966 Module_Ar.DumpConstants();
967 Board_Ar.DumpConstants();
968 Chip_Ar.DumpConstants();
969 Straw_Ar.DumpConstants();
970 }
971
972 histfile->ls();
973 ATH_MSG_INFO( "writing out calibout.root");
974 histfile->Write();
975
976 return true;
977
978}
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition AthMsgStreamMacros.h:33
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition AthMsgStreamMacros.h:31
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition AthMsgStreamMacros.h:32
Calibrator
A structure to contain data associated with the calibration of a certain sub-module.
Definition Calibrator.h:178
InDetDD::TRT_BarrelElement
Extended TRT_BaseElement to describe a TRT readout element, this is a planar layer with n ( order of ...
Definition TRT_BarrelElement.h:43
InDetDD::TRT_BaseElement
Virtual base class of TRT readout elements.
Definition TRT_BaseElement.h:56
InDetDD::TRT_BaseElement::center
virtual const Amg::Vector3D & center() const override final
Element Surface: center of a straw layer.
InDetDD::TRT_EndcapElement
Extended class of a TRT_BaseElement to describe a readout elment in the endcap.
Definition TRT_EndcapElement.h:43
TRTCond::BasicRtRelation
Default rt-relation class in the TRT: a 3rd degree polynomial.
Definition BasicRtRelation.h:26
TRTCond::BasicRtRelation::cal
const float * cal() const
return to calibration constants
Definition BasicRtRelation.h:76
TRTCond::DinesRtRelation
Default rt-relation class in the TRT: a 3rd degree polynomial.
Definition DinesRtRelation.h:26
TRTCond::DinesRtRelation::cal
const float * cal() const
return to calibration constants
Definition DinesRtRelation.h:76
caldata
A structure to contain data associated with the calibration of a certain sub-module.
Definition Calibrator.h:104
databundle
A structure to contain hit data.
Definition Calibrator.h:77
databundle::res
float res
space residual
Definition Calibrator.h:90
databundle::det
int det
detector (barrel -1 or 1, or end-cap -2 or 2)
Definition Calibrator.h:79
databundle::rt0
float rt0
reference t0 (offset from board mean)
Definition Calibrator.h:94
databundle::ievt
int ievt
event number
Definition Calibrator.h:87
databundle::chp
int chp
chip
Definition Calibrator.h:83
databundle::t
float t
raw time
Definition Calibrator.h:91
databundle::r
float r
drift radius from r(t) relation
Definition Calibrator.h:92
databundle::brd
int brd
board
Definition Calibrator.h:82
databundle::tres
float tres
time residual
Definition Calibrator.h:88
databundle::rtpar
std::array< float, 4 > rtpar
rt-parameters used in the reconstruction
Definition Calibrator.h:95
databundle::t0
float t0
t0 value used in reconstruction
Definition Calibrator.h:93
databundle::stw
int stw
straw number (within the strawlayer)
Definition Calibrator.h:85
databundle::stl
int stl
straw-layer
Definition Calibrator.h:84
databundle::z
float z
straw z position
Definition Calibrator.h:98
databundle::x
float x
straw x position
Definition Calibrator.h:96
databundle::lay
int lay
layer
Definition Calibrator.h:80
databundle::y
float y
straw y position
Definition Calibrator.h:97
databundle::weight
float weight
histogram weight (currently set to 1)
Definition Calibrator.h:89
databundle::mod
int mod
phi module
Definition Calibrator.h:81
databundle::sid
int sid
straw ID
Definition Calibrator.h:86
TRT::Hit::ident
@ ident
Definition HitInfo.h:77
TRT
Definition HitInfo.h:33
keylayer_zslicemap.pb
dict pb
Definition keylayer_zslicemap.py:188
makeDTCalibBlob_pickPhase.pd
pd
Definition makeDTCalibBlob_pickPhase.py:342
module_driven_slicing.histfile
histfile
Definition module_driven_slicing.py:570
test_pyathena.pt
pt
Definition test_pyathena.py:11
Generated on for ATLAS Offline Software by doxygen 1.17.0