#include <string.h>
#include <iostream>
#include <TEnv.h>
#include "EdbLog.h"
#include "EdbScanProc.h"
#include "EdbScanTracking.h"
#include "EdbMosaicIO.h"

Include dependency graph for tagtra.cpp:

Functions
int	GetPattern (EdbPattern &pat, EdbID id, int side)

int	main (int argc, char *argv[])

void	print_help_message ()

void	set_default (TEnv &cenv)

void	TrackSetTag (EdbID idset, TEnv &env, EdbScanProc &sproc)

Function Documentation

◆ GetPattern()

int GetPattern	(	EdbPattern &	pat,
		EdbID	id,
		int	side
	)

 {
   int n=0;
   EdbMosaicIO mio;
   TString file;
   file.Form("p%3.3d/%d.%d.%d.%d.tag.root",
         id.ePlate, id.eBrick, id.ePlate, id.eMajor, id.eMinor);
   TFile f( file.Data() );
   EdbPattern *p = (EdbPattern*)(f.Get( Form("gpat%d_%d_0",id.ePlate, side) ) );
   if(p) {
     n = p->N();
     //p->SetSegmentsW(1.);                        //TODO set correct weights peak weights in mostag
     pat.AddPattern(*p);
     delete p;
   }
   //  p->SetScanID(id); p->SetSide(side);
   printf( "read %d segments from %s\n",n, file.Data() );
   return n;
 }

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

 {
   if (argc < 2)   { print_help_message();  return 0; }
   
   TEnv cenv("trackenv");
   set_default(cenv);
   gEDBDEBUGLEVEL        = cenv.GetValue("tagtra.EdbDebugLevel" , 1);
   const char *env       = cenv.GetValue("tagtra.env"            , "tagtra.rootrc");
   const char *outdir    = cenv.GetValue("tagtra.outdir"         , "..");
   
   bool      do_set              = false;
   bool      do_pred             = false;
   Int_t     pred_plate  = 0, to_plate=0;
   Int_t     brick=0, plate=0, major=0, minor=0;
   
   for(int i=1; i<argc; i++ ) {
     char *key  = argv[i];
     
     if(!strncmp(key,"-set=",5))
     {
       if(strlen(key)>5) sscanf(key+5,"%d.%d.%d.%d",&brick,&plate,&major,&minor);
       do_set=true;
     }
     else if(!strncmp(key,"-o=",3)) 
     {
       if(strlen(key)>3) outdir=key+3;
     }
     else if(!strncmp(key,"-v=",3))
     {
       if(strlen(key)>3) gEDBDEBUGLEVEL = atoi(key+3);
     }
   }
   
   if(!do_set)   { print_help_message(); return 0; }
   
   cenv.SetValue("tagtra.env"                  , env);
   cenv.ReadFile( cenv.GetValue("tagtra.env"   , "tagtra.rootrc") ,kEnvLocal);
   cenv.SetValue("tagtra.outdir"               , outdir);
   cenv.WriteFile("tagtra.save.rootrc");
   
   if(do_set) {
     EdbScanProc sproc;
     sproc.eProcDirClient=outdir;
     printf("\n----------------------------------------------------------------------------\n");
     printf("tracking set %d.%d.%d.%d\n", brick,plate, major,minor);
     printf("----------------------------------------------------------------------------\n\n");
     
     EdbID id(brick,plate,major,minor);
     EdbScanSet *ss = sproc.ReadScanSet(id);
     ss->Brick().SetID(brick);
     TrackSetTag(id,cenv,sproc);
   }
   cenv.WriteFile("tagtra.save.rootrc");
   
   return 1;
 }

◆ print_help_message()

void print_help_message ( )

 {
   cout<< "\nUsage: \n\t  tagtra -set=ID [ -o=DATA_DIRECTORY -v=DEBUG] \n\n";
   cout<< "\t\t  ID              - id of the dataset formed as BRICK.PLATE.MAJOR.MINOR \n";
   cout<< "\t\t  DEBUG           - verbosity level: 0-print nothing, 1-errors only, 2-normal, 3-print all messages\n";
   cout<< "\t\t  -o              - the data directory\n";
   cout<< "\nExample: \n";
   cout<< "\t  tagtra -id=4554.10.1.0 -o=/scratch/BRICKS \n";
   cout<< "\n If the data location directory if not explicitly defined\n";
   cout<< " the current directory will be assumed to be the brick directory \n";
   cout<<endl;
 }

◆ set_default()

void set_default ( TEnv & cenv )

 {
   // default parameters for tracking
   cenv.SetValue("fedra.track.TrZmap", "4000 0 200000   4000 0 200000   30" );
   cenv.SetValue("fedra.track.npass"     ,   1 );
   cenv.SetValue("fedra.track.minPlate"  ,-999 );
   cenv.SetValue("fedra.track.maxPlate"  , 999 );
   cenv.SetValue("fedra.track.refPlate"  , 999 );
   cenv.SetValue("fedra.track.nsegmin"   ,   2 );
   cenv.SetValue("fedra.track.ngapmax"   ,   4 );
   cenv.SetValue("fedra.track.DZGapMax"  , 5000. );
   cenv.SetValue("fedra.track.DRmax"     , 100. );
   cenv.SetValue("fedra.track.DTmax"     , 0.07 );
   
   cenv.SetValue("fedra.track.Sigma0" , "50 50 0.005 0.005");
   cenv.SetValue("fedra.track.PulsRamp0"      , "0.1 0.1");
   cenv.SetValue("fedra.track.PulsRamp04"     , "0.1 0.1");
   cenv.SetValue("fedra.track.Degrad"         ,  4   );
   
   cenv.SetValue("fedra.track.probmin"   , 0.001 );
   cenv.SetValue("fedra.track.momentum"  , 2 );
   cenv.SetValue("fedra.track.mass"      , 0.14 );
   cenv.SetValue("fedra.track.do_use_mcs", 0 );
   cenv.SetValue("fedra.track.RadX0"     , 5810.);
   
   cenv.SetValue("fedra.track.do_shtag"       , false );
   cenv.SetValue("fedra.track.do_misalign"    , false );
   cenv.SetValue("fedra.track.misalign_offset", 5000.);
   cenv.SetValue("fedra.track.do_comb"        ,    0 );
   cenv.SetValue("fedra.track.NsegMin"        ,    2 );
   cenv.SetValue("tagtra.outdir"          , "..");
   cenv.SetValue("tagtra.env"             , "tagtra.rootrc");
   cenv.SetValue("tagtra.EdbDebugLevel"   , 1);
 }

◆ TrackSetTag()

void TrackSetTag	(	EdbID	idset,
		TEnv &	env,
		EdbScanProc &	sproc
	)