EdbEDAIO Class Reference

#include <EdbEDA.h>

Inheritance diagram for EdbEDAIO:

Public Types
enum	OutputFileMode { kBern , kOtherLabs }

Public Member Functions
	EdbEDAIO ()
EdbPattern *	GetPatternIPL (int pid)
int	IdMuon (const char *filename="../cs_info.txt")
void	OpenTextEditor (const char *filename)
EdbPattern *	ReadCouples (int ipl, EdbPattern *pat=NULL)
EdbPattern *	ReadCouplesPID (int pid, EdbPattern *pat=NULL)
void	ReadFeedbackFile (const char *filename=NULL)
EdbVertex *	ReadFeedbackFile2008 (const char *filename=NULL)
EdbPVRec *	ReadFilteredText (const char *filename=NULL)
void	WriteFeedbackFile (const char *filename=NULL)
void	WriteFeedbackFile2008 (EdbVertex *v=NULL, const char *filename=NULL)
void	WriteFilteredText (const char *filename=NULL)
virtual	~EdbEDAIO ()

Public Attributes
TCut	eCutCP
int	eOutputFileMode

Member Enumeration Documentation

OutputFileMode

enum EdbEDAIO::OutputFileMode

Enumerator
kBern
kOtherLabs

                       {
        kBern,
        kOtherLabs
    };

Constructor & Destructor Documentation

EdbEDAIO()

EdbEDAIO::EdbEDAIO ( )

inline

:eOutputFileMode(kBern), eCutCP("eN1==1&&eN2==1&&eCHI2P<s.eW*0.12-1"){}

~EdbEDAIO()

virtual EdbEDAIO::~EdbEDAIO ( )

inlinevirtual

{}

Member Function Documentation

GetPatternIPL()

EdbPattern * EdbEDAIO::GetPatternIPL

(

int

pid

)

                                           {
    // return Pattern with Couples.
    // if Couples are not loaded yet, load it in PVRec for TS.
    
    EdbEDATrackSet *set = gEDA->GetTrackSet("TS");
    EdbPVRec *pvr = set->GetPVRec();
 
    EdbPattern *pat=NULL;
    
    for(int i=0;i<pvr->Npatterns();i++){
        EdbPattern *p = pvr->GetPattern(i);
        if(p==NULL) continue;
        if(p->N()>0){
            if(p->GetSegment(0)->Plate()==ipl) {
                pat=p;
                break;
            }
        }
    }
    
    if(pat==NULL) {
        printf("No pattern in TS, pl=%d.\n", ipl);
        return NULL;
    }
    
    // Check if couples data were already loaded.
    if( IsIncludeCouples(pat)==0 ) {
        
        if     (set->GetDataSet()!=NULL) gEDA->ReadCouples(ipl,pat);
        else if(set->GetScanSet()!=NULL) {
            printf("ReadCouples SS ipl=%d\n", ipl);
            EdbID id = *(EdbID *)set->GetScanSet()->eIDS.At(0);
            id.ePlate=ipl;
            id.Print();
            float z = pat->Z();
            gEDA->ScanProc()->ReadPatCP(*pat,id,eCutCP);
            pat->SetZ(z);
            pat->SetSegmentsZ();
            printf("z = %lf\n", pat->Z());
            pat->SetSegmentsPlate(ipl);
        }
    }
    return pat;
}

IdMuon()

int EdbEDAIO::IdMuon

(

const char *

filename = "../cs_info.txt"

)

                                        {
    // code by GL
  char            buf[256];
  char            key[256];
  char m1[12], m2[10], m3[10],m4[10],m5[10],m6[10];
  static int val=-2;
  char mu[3];
  
  if(val>=0) return val;
  if(val==-1) return val;
 
  FILE *fp=fopen(filename,"r");
  if (fp==NULL)   {
    printf("ERROR open file: %s \n", filename);
    val = -1;
    return val;
  }
  
  // printf( "Read affine transformation from file: %s\n", fname );
  while( fgets(buf,256,fp)!=NULL ) {
 
    for( int i=0; i<(int)strlen(buf); i++ ) 
      if( buf[i]=='#' )  {
        buf[i]='0';                       // cut out comments starting from #
        break;
      }
    
   if( sscanf(buf,"%s",key)!=1 ) continue;
   if      ( !strcmp(key,"Prediction")   )
     {
       sscanf(buf+strlen(key),"%s %s %s %s %s %s ",m1,m2,m3,m4,m5,m6);
       mu[0] = m2[3];
       mu[1] = m2[4];
       mu[2] = m2[5];
       val = atol(mu);
       //       cout << m2[3] << " " << m2[4] << " " << m2[5]<< " " << val << endl;
 
     }
   if      ( !strcmp(key,"No")   )
     {
       val=-1;
     }
 
  }
  fclose(fp);   
 
  return val;
}

OpenTextEditor()

void EdbEDAIO::OpenTextEditor

(

const char *

filename

)

                                                 {
    TEveBrowser* browser = gEve->GetBrowser();
    browser->StartEmbedding(1);
    new TGTextEditor(filename, gClient->GetRoot());
    browser->StopEmbedding();
    
    unsigned int i,k=0;
    for(i=0;i<strlen(filename);i++){
        if(filename[i]=='\\'||filename[i]=='/') k=i+1;
    }
    browser->SetTabTitle(filename+k, 1);
}

ReadCouples()

EdbPattern * EdbEDAIO::ReadCouples	(	int	ipl,
		EdbPattern *	pat = NULL
	)

                                                           {
    printf("ReadCouples DS ipl=%d\n", ipl);
    return ReadCouplesPID(gEDA->GetPID(ipl), pat0);
}

ReadCouplesPID()

EdbPattern * EdbEDAIO::ReadCouplesPID	(	int	pid,
		EdbPattern *	pat = NULL
	)

                                                              {
            
    EdbDataSet *dset = gEDA->GetDataSet();
    EdbDataPiece *piece = dset->GetPiece(pid);
    EdbLayer *l = piece->GetLayer(0);
    EdbPattern *pat = new EdbPattern( 0.,0., l->Z());
    pat->SetPID(pid);
    piece->InitCouplesTree("READ");
    piece->GetCPData_new(pat);
    piece->CloseCPData();
    pat->SetSegmentsZ();
    pat->SetID(pid);
    pat->SetSegmentsPID();
    pat->Transform(    l->GetAffineXY()   );
    pat->TransformA(   l->GetAffineTXTY() );
    pat->TransformShr( l->Shr()  );
    
    if(pat0==0) return pat;
    
    EdbScanCond cond;
    
    // copy segments.
    for(int i=0;i<pat->N();i++){
        EdbSegP *s = pat->GetSegment(i);
        if(s->Plate()==0) s->SetPlate(gEDA->GetIPL(pid));
    
        // fill COV for vertexing
        s->SetErrors();
        cond.FillErrorsCov(s->TX(), s->TY(), s->COV());
    
        pat0->AddSegment(*s);
    }
    delete pat;
    return pat0;
}

ReadFeedbackFile()

void EdbEDAIO::ReadFeedbackFile

(

const char *

filename = NULL

)

                                                   {
    // Read feedback file format (ver 2009 Oct).
    
    gEDA->Reset();
    
    
    char FlagPart[][10]={"", "mu", "charm", "e", "e-pair"};
    char FlagCS[][10]   ={"", "cs"};
    
    // if filename is not given, open file browser.
    if(filename==NULL){
        TGFileInfo *fi=new TGFileInfo;
        fi->fIniDir    = StrDup(".");
        const char *filetypes[] = { "Feedback File", "*.feedback", "All files","*",0,0};
        fi->fFileTypes=filetypes;
        new TGFileDialog(gClient->GetRoot(), gEve?gEve->GetMainWindow():0, kFDOpen, fi);
        filename = fi->fFilename;
    }
    
    if(NULL==filename) return;
    
    FILE *fp = fopen(filename,"rt");
    if(fp==NULL) {
        ErrorMessage(Form("File open error : %s . stop.\n", filename));
        return;
    }
    
    EdbPVRec *pvr = new EdbPVRec;
    EdbScanCond cond;
    EdbTrackP *t=0;
    EdbEDAText *tx = 0;
    
    char buf[256];
    
    char asetname[10];
    EdbEDAAreaSet *aset=NULL;
    
    while(fgets(buf,sizeof(buf),fp)){
        TString line=buf;
        printf("%s", line.Data());
 
        // track comment = ""
        if(line.Contains("track comment")){
            int istart = line.Index("\"");
            int iend   = line.Index("\"",1,istart+1, TString::kExact);
            if(istart!=kNPOS && iend!=kNPOS && istart+1!=iend){
                char buf[256];
                strncpy(buf, line.Data()+istart+1, iend-istart-1);
                buf[iend-istart-1]='\0';
                tx->AddText(buf);
            }
        }
        
        if(line.Contains("Area information :")){
            aset = new EdbEDAAreaSet;
            aset->SetDraw(kTRUE);
            sscanf( line.Data()+line.Index(":")+1, "%s", asetname);
            printf("Area infor %s\n", asetname);
        }
        if(aset){
            if(line.Contains(Form("area : %s", asetname))){
                int ipl;
                double xmin, xmax, ymin, ymax, z;
                sscanf(line.Data() + line.Index(Form("area : %s", asetname)) + strlen(Form("area : %s", asetname)),
                    "%d %lf %lf %lf %lf %lf", &ipl, &xmin, &xmax, &ymin, &ymax, &z);
                aset->AddArea(ipl, xmin, xmax, ymin, ymax, z);
            }
        }
        
        // Remove comments
        line.ReplaceAll("\t"," ");
        int iposcmt = line.Index("//");
        if(iposcmt!=kNPOS) line.Remove(iposcmt, line.Length()-iposcmt);
        
        // Check number of columns.
        TObjArray *tokens = line.Tokenize(" ");
        int ntokens = tokens->GetEntries();
        delete tokens;
        
        
        if(ntokens==0) continue;
        
        else if( ntokens == 10 ){
            // Vertices
            float x,y,z; int id,isprimary,istau, ischarm;
            sscanf(line.Data(),"%d %f %f %f %d %d %d", &id, &x, &y, &z, &isprimary, &ischarm, &istau);
            EdbVertex *v = new EdbVertex();
            v->SetXYZ(x,y,z); v->SetID(id);
            v->SetFlag(isprimary);
            pvr->AddVertex(v);
            printf("Vertex %d %f %f %f\n",v->ID(), v->X(), v->Y(), v->Z());
        }
 
        else if( ntokens == 27 ){
            // Tracks
            float x,y,z,ax,ay, ip_upvtx, ip_downvtx,  p,pmin,pmax;
            int id_track, id_upvtx, id_downvtx,  manual, particle_id, scanback, darkness;
            int OutOfBrick, LastPlate, nseg, iplRmax1, iplRmax2, result;
            float RmaxT, RmaxL, rmst, rmsl;
            
            int n,nn;
            sscanf(line.Data(),         "%d %d %d %f %f %f %f %f %f%n", &id_track, &id_upvtx, &id_downvtx, &x, &y, &z, &ax, &ay, &ip_upvtx, &nn);
            sscanf(line.Data()+(n=nn),  "%f %f %f %f %d %d %d %d %d%n", &ip_downvtx, &p,&pmin,&pmax, &manual, &particle_id, &scanback, &darkness, &OutOfBrick, &nn);
            sscanf(line.Data()+(n+=nn), "%d %d %f %f %f %f %d %d %d",   &LastPlate, &nseg, &RmaxT, &RmaxL, &rmst, &rmsl, &iplRmax1, &iplRmax2, &result);
            
            // Create Track. "t" is defined out of main loop.
            t = new EdbTrackP;
            t->Set(id_track, x, y, ax, ay, 0, 0);
            t->SetZ(z);
            t->SetTrack(id_track);
            pvr->AddTrack(t);
            
            tx = new EdbEDAText(t->X(), t->Y(), t->Z(), "[");
            if(0<=particle_id && particle_id<=4) { if(tx->N()>1)tx->AddText(", "); tx->AddText(FlagPart[particle_id]);}
            if(scanback==1)                      { if(tx->N()>1)tx->AddText(", "); tx->AddText(FlagCS[scanback]);}
            tx->AddText("]");
            if(tx->N()==2) tx->SetText("");
            tx->SetReference(tx);
            gEDA->AddDrawObject(tx);
            
            // fill COV for vertexing
            t->SetErrors();
            cond.FillErrorsCov(t->TX(), t->TY(), t->COV());
            
            // associate to vertex.
            for(int i=0; i<pvr->Nvtx(); i++){
                EdbVertex *v = pvr->GetVertex(i);
                if(id_upvtx==v->ID()||id_downvtx==v->ID()){
                    EdbVTA *vta = new EdbVTA(t, v);
                    vta->SetZpos( t->Z()>v->Z() ? 1 : 0);
                    vta->SetFlag(2);
                    vta->SetImp( id_upvtx==v->ID()? ip_upvtx : ip_downvtx);
                    v->AddVTA(vta);
                }
            }
        }
        
        else if( ntokens ==  9 ){
            // Segments
            int ipl, type, irec, ph;
            float x,y,z,ax,ay;
            sscanf(line.Data(),"%d %f %f %f %f %f %d %d %d", &ipl, &x, &y, &z, &ax, &ay, &type, &irec, &ph);
            
            EdbSegP *s = new EdbSegP(t->ID(),x,y,ax,ay,0,0);
            s->SetZ(z);
            s->SetPID(ipl);
            s->SetPlate(ipl);
            s->SetW(ph);
            s->SetTrack(t->ID());
            s->SetSide(type);
 
            // fill COV for vertexing
            s->SetErrors();
            cond.FillErrorsCov(s->TX(), s->TY(), s->COV());
            
            // Add segment in PVRec and Track, keeping consistency of pointer in them.
            EdbSegP *s_in_pattern = pvr->AddSegment(*s);
            t->AddSegment( s_in_pattern);
            
            if(NULL!=tx){
                EdbSegP seg(*s);
                seg.PropagateTo(seg.Z()+1000);
                tx->SetXYZ(seg.X(),seg.Y(),seg.Z());
            }
        }
        
        
    }
    
    for(int i=0;i<pvr->Npatterns(); i++) pvr->GetPattern(i)->SetPID(i);
    for(int i=0;i<pvr->Npatterns(); i++) pvr->GetPattern(i)->SetSegmentsPID();
    for(int i=0;i<pvr->Ntracks(); i++) pvr->GetTrack(i)->SetCounters();
    
    printf("\nEdbPVRec summary. %d vertices, %d tracks.\n", pvr->Nvtx(), pvr->Ntracks());
    pvr->Print();
    
    fclose(fp);
 
    if(pvr==NULL) return;
    for(int i=0;i<pvr->Npatterns(); i++){
        EdbPattern *pat = pvr->GetPattern(i);
        pat->SetPID(gEDA->GetPID(pat->ID()));
    }
 
    gEDA->SetPVR(pvr);
    gEDA->GetTrackSet("TS")->AddTracksPVR(pvr);
    gEDA->GetVertexSet()->AddVertices(pvr->eVTX);
    
    if(aset){
        EdbEDATrackSet *set = gEDA->GetTrackSet(asetname);
        if(set) set->SetAreaSet(aset);
    }
    gEDA->Redraw();
    gEDA->UpdateScene();
}

ReadFeedbackFile2008()

EdbVertex * EdbEDAIO::ReadFeedbackFile2008

(

const char *

filename = NULL

)

                                                              {
    using namespace std;
 
    // if filename is not given, open file browser.
    if(filename==NULL){
        TGFileInfo *fi=new TGFileInfo;
        fi->fIniDir    = StrDup(".");
        const char *filetypes[] = { "Feedback File", "*.feedback", "All files","*",0,0};
        fi->fFileTypes=filetypes;
        new TGFileDialog(gClient->GetRoot(), gEve->GetMainWindow(), kFDOpen, fi);
        filename = fi->fFilename;
    }
    
    if(NULL==filename) return NULL;
        
    // check the file.
    ifstream fin(filename);
    if(!fin) {
        printf("can not open file : %s\n", filename);
        return NULL;
    } printf("read feedback %s\n", filename);
    
 
    gEDA->ClearVertices();
    
    EdbEDATrackSet *set = gEDA->GetTrackSet("TS");
    set->ClearTracks();
    
    // read file line by line.
    char buf[512];
    while( fin && fin.getline(buf,512) ){
        string line(buf);
        
        // ignore comment
        if(line.find("//")!=line.npos) line[line.find("//")]='\0';
        
        // check the number of columnline.
        istringstream ss(line);
        string s2;
        int ncolumn;
        for(ncolumn=0;ss>>s2;ncolumn++){
            //cout <<s2 << endl;
        }
        
        // ncolumn == 7 --> vertex info.  ncolumn == 17 --> track info
        if(ncolumn==7){
            double x,y,z;
            int id,isprimary,istau, ischarm;
            sscanf(line.c_str(),"%d %lf %lf %lf %d %d %d", &id, &x, &y, &z, &isprimary, &ischarm, &istau);
            EdbVertex *v = new EdbVertex();
            v->SetXYZ(x,y,z);
            v->SetID(id);
            gEDA->AddVertex(v);
        }
        else if (ncolumn==17){
            double x,y,z,ax,ay, ip_upvtx, ip_downvtx;
            double p,pmin,pmax;
            int id_track, id_upvtx, id_downvtx;
            int manual, particle_id, scanback, darkness;
            
            istringstream ss(line);
            
            ss  >> id_track >> id_upvtx >> id_downvtx >> x >> y >> z >> ax >> ay >> ip_upvtx >> ip_downvtx 
                >> p >> pmin >> pmax >> manual >> particle_id >> scanback >> darkness;
            
            EdbSegP s(id_track, x, y, ax, ay);
            s.SetZ(z);
            
            EdbTrackP *t = set->FindTrack(&s);
            
            if(t) {
                printf("%s track %d\n", line.c_str(), t->ID());
                set->SetTrack(t);
                for(int i=0; i<gEDA->GetVertexSet()->N(); i++){
                    EdbVertex *v = gEDA->GetVertexSet()->GetVertex(i);
                    if(id_upvtx!=v->ID()) continue;
                    EdbVTA *vta = new EdbVTA(t, v);
                    vta->SetZpos(1);
                    vta->SetFlag(2);
                    vta->SetImp(EdbEDAUtil::CalcIP(&s,v));
                    v->AddVTA(vta);
                }
            }
            else  printf("track not found. id = %d\n", id_track);
        }
    }
    gEDA->Redraw();
    if( gEDA->GetVertexSet()->N()) return gEDA->GetVertexSet()->GetVertex(0);
    else return NULL;
}

ReadFilteredText()

EdbPVRec * EdbEDAIO::ReadFilteredText

(

const char *

filename = NULL

)

                                                         {
    int i,j;
    
    if(filename==NULL||strlen(filename)<1){ // File dialog
        TGFileInfo *fi=new TGFileInfo;
        fi->fIniDir    = StrDup(".");
        const char *filetypes[] = { "Filtered text", "*.txt", "All files","*",0,0};
        fi->fFileTypes = filetypes;
        new TGFileDialog(gClient->GetRoot(), gEve->GetMainWindow(), kFDOpen, fi);
        if(fi->fFilename!=NULL) {
            filename = fi->fFilename;
        }
    }
    if(filename==NULL) return NULL;
    
    // check the file.
    ifstream fin(filename);
    if(!fin) {
        printf("can not open file : %s\n", filename);
        return NULL;
    } 
    else printf("Read Filtered Text File : %s\n", filename);
    
 
    TObjArray *vtxarr = new TObjArray;
    TObjArray *trkarr = new TObjArray;
    TObjArray *segarr = new TObjArray;  
    
    EdbPVRec *PVR = new EdbPVRec;
    
    EdbDataSet *dset = gEDA->GetDataSet();
 
    for(i=0;i<dset->N();i++){
        EdbDataPiece *piece = dset->GetPiece(i);
        EdbLayer* l = piece->GetLayer(0);
        EdbPattern *pat = new EdbPattern( 0.0, 0.0, l->Z());
        pat->SetPID(i);
        PVR->AddPattern(pat);
    }
    
    char buf[512];
    while( fin && fin.getline(buf,512) ){
        string line(buf);
 
        std::istringstream ssline(line);
        std::string column1;
        ssline >> column1;
        
        std::cout <<line <<std::endl;
        if(column1.find("//")!=column1.npos) continue;
 
        if(column1=="vertex"){
            int id;
            float x, y, z;
            ssline >> id >> x >> y >> z;
            EdbVertex *v = new EdbVertex;
            v->SetXYZ(x,y,z);
            v->SetID(id);
            vtxarr->Add(v);
            PVR->AddVertex(v);
        }
        
        if(column1=="track"){
            int id, flag, ipl1, ipl2, nseg, ivtx1, ivtx2;
            double tx,ty, ip_1ry, ip1, ip2;
            ssline >> id >> flag >> tx >> ty >> ipl1 >> ipl2 >> nseg >> ip_1ry >> ivtx1 >> ip1 >> ivtx2 >> ip2;
            EdbTrackP *t = new EdbTrackP;
            t->SetID(id);
            t->SetFlag(flag);
            for(i=0;i<vtxarr->GetEntriesFast();i++){
                EdbVertex *v = (EdbVertex *) vtxarr->At(i);
                if(v->ID()==ivtx1||v->ID()==ivtx2){
                    EdbVTA *vta = new EdbVTA(t, v);
                    vta->SetFlag(2);
                    v->AddVTA(vta);
                }
            }
            trkarr->Add(t);
            PVR->AddTrack(t);
        }
        
        if(column1=="segment"){
            int itrk, ipl, iseg, sflag, tflag;
            float x, y, z, ax, ay, chi2, w;
            ssline >> itrk >> tflag >> ipl >> iseg >> x >> y >> z >> ax >> ay >> w >> chi2 >> sflag;
            
            EdbSegP *s = new EdbSegP (iseg, x, y, ax, ay, w, sflag);
            s->SetPID(gEDA->GetPID(ipl));
            s->SetZ(z);
            s->SetDZ(300);
            s->SetTrack(itrk);
            s->SetChi2(chi2);
            
            for(i=0;i<trkarr->GetEntriesFast();i++){
                EdbTrackP *t = (EdbTrackP *) trkarr->At(i);
                if(t->ID()==itrk&&t->Flag()==tflag) {
                    t->AddSegment(s);
                }
            }
            
            EdbPattern *pat = PVR->GetPattern( gEDA->GetPID(ipl));
            pat->AddSegment(*s);
            segarr->Add(s);
            
        }
        
        if(column1=="area"){
            int ipl;
            double xmin,xmax,ymin,ymax,z;
            ssline >> ipl >> xmin >> xmax >> ymin >> ymax >> z;
            EdbEDAArea *a = new EdbEDAArea(ipl, xmin, xmax, ymin, ymax, z);
            gEDA->GetTrackSet("TS")->GetAreaSet()->AddArea(a);
        }
        
        
    }
    fin.close();
 
    printf("Vertices %4d\n", vtxarr->GetEntriesFast());
    printf("Tracks   %4d\n", trkarr->GetEntriesFast());
    printf("Segments %4d\n", segarr->GetEntriesFast());
    
    
    // Clear all track set
    for(i=0; i<gEDA->NTrackSets();i++){
        gEDA->GetTrackSet(i)->Clear();
    }
    
    for(i=0;i<trkarr->GetEntriesFast();i++){
        EdbTrackP *t = (EdbTrackP *) trkarr->At(i);
        t->SetCounters();
        EdbSegP *s = t->GetSegmentFirst();
        t->Set(t->ID(), s->X(), s->Y(), s->TX(), s->TY(), 0, t->Flag());
        t->SetZ(s->Z());
        
        if( 0 <= t->Flag()) gEDA->GetTrackSet("TS")->AddTrack(t);
        if( -100 < t->Flag() && t->Flag() <     0) gEDA->GetTrackSet("MN")->AddTrack(t);
        if( -200 < t->Flag() && t->Flag() <=  -99) gEDA->GetTrackSet("SB")->AddTrack(t);
        if(-1000 < t->Flag() && t->Flag() <= -200) gEDA->GetTrackSet("SF")->AddTrack(t);
    }
    
    EdbEDAVertexSet *vset = gEDA->GetVertexSet();
    gEDA->ClearSelectedVertex();
    vset->Clear();
    for(i=0;i<vtxarr->GetEntriesFast();i++){
        EdbVertex *v = (EdbVertex *) vtxarr->At(i);
        
        for(j=0;j<v->N();j++){
            EdbVTA    *vta = v->GetVTa(j);
            EdbTrackP *t   = vta->GetTrack();
            EdbSegP   *s   = t->GetSegmentFirst();
            vta->SetZpos( s->Z()>v->Z() ? 1 : 0 );
            vta->SetImp(CalcIP(s,v));
        }
        vset->AddVertex(v);
    }
    
    return PVR;
}

WriteFeedbackFile()

void EdbEDAIO::WriteFeedbackFile

(

const char *

filename = NULL

)

                                                    {
    // write feedback file
    // if eOutputFileMode == 0 (kBern), output comments with "//" escape sequence (default).
    //                                  also track id is not 1,2,3,4, but id of track.
    // if eOutputFileMode == 1 (kOtherLabs). doesn't output comments.
    //
    // to convert Bern format, use
    // gawk -f ConvertFeedback.awk tmp.feedback
    //
    // -- ConvertFeedback.awk --
    // BEGIN{ itrk=1;}
    // !/\/\//&&NF==27{print itrk++,$2,$3,$4,$5,$6,$7,$8,$9,$10,$11,$12,$13,$14,$15,$16,$18,$19,$20,$21,$22,$23,$24,$25,$26,$27;}
    // !/\/\//&&NF==9{print $1,$2,$3,$4,$5,$6,$7,$8,$9;}
    // !/\/\//&&NF==10{print $1,$2,$3,$4,$5,$6,$7,$8,$9,$10;}
    
    
    if( gEDA->GetVertexSet()->N() == 0 ){
        ErrorMessage("No vertex is found. stop.");
        return;
    }
    
    if(filename==NULL){ // File dialog
        TGFileInfo *fi=new TGFileInfo;
        fi->fIniDir    = StrDup(".");
        const char *filetypes[] = { "Feedback file", "*.feedback", "All files","*",0,0};
        fi->fFileTypes = filetypes;
        new TGFileDialog(gClient->GetRoot(), gEve->GetMainWindow(), kFDSave, fi);
        if(fi->fFilename!=NULL) filename = fi->fFilename;
    }
    
    FILE *fp;
    if(NULL==filename||strlen(filename)==0) {
        return;
    }
    else {
        fp=fopen(filename,"wt");
        if(NULL==fp){
            printf("Couldn't open file : %s\n", filename);
            return;
        }
    }
 
    printf("--------- Feedback format---------------\n");
    printf("VID        X        Y        Z 1ry charm tau Ndw Nup OutOfBrick\n");
    printf("trkid v1 v2        x        y        z      tx      ty   ip1   ip2 ");
    printf("    p  pmin  pmax mn pa sb dk ou ps ");
    printf(" n RmaxT RmaxL   rmsT   rmsL pl1 pl2 res\n");
    printf("        ipl        x        y        z      tx      ty type irec ngrains\n");
    
    if(eOutputFileMode==kBern){
        fprintf(fp,"// --------- Feedback format---------------\n");
        fprintf(fp,"// VID X Y Z 1ry charm tau Ndw Nup OutOfBrick\n");
        fprintf(fp,"//   OutOfBrick   : 1= vertex in dead material, 0= others\n");
        fprintf(fp,"// \n");
        fprintf(fp,"// Trkid v1 v2 x y z tx ty ip1 ip2 p pmin pmax mn pa sb dk ou ps n RmaxT RmaxL rmsT rmsL pl1 pl2 res\n");
        fprintf(fp,"//   Manual       : 0= auto, 1= manually added\n");
        fprintf(fp,"//   Particle     : 1= muon, 2= charm, 3= electron, 4= e+e- pair\n");
        fprintf(fp,"//   Scanback     : 1= scanback, 0=not scanback\n");
        fprintf(fp,"//   Darkness     : 0= MIP, 1= BLACK, 2= GRAY\n");
        fprintf(fp,"//   OutOfBrick   : 1= pass through up, 2= edge out, 0 others\n");
        fprintf(fp,"//   Last plate   : 1= pass through down, n= edge out last plate number, 0 others\n");
        fprintf(fp,"//   Out_flag(res): 1= 1ry vertex track, 2= e+e- pair, 3= Low momentum (DTheta_RMS >20 mrad),\n");
        fprintf(fp,"//                  4= Scan Forth to be done, 5= Scan Forth  done\n");
        fprintf(fp,"// \n");
        fprintf(fp,"// Ipl x y z tx ty type irec ngrains\n");
        fprintf(fp,"//   Type         : 0 BT, 1 microtrack top, 2 microtrack bottom\n");
        fprintf(fp,"//   Irec         : 0 Automatic, 1 SB, 2 Manual\n");
        fprintf(fp,"// \n");
        fprintf(fp,"///////////////////////////////////////////// \n");
        fprintf(fp,"// ---------- Vertex Data ---------------- // \n");
        fprintf(fp,"///////////////////////////////////////////// \n");
    }
    // Vertices
    
    TObjArray *vertices = (gEDA->GetVertexSet())->GetVertices();
    
    printf("%d vertices\n", vertices->GetEntriesFast());
    EdbVertex *v_1ry=NULL;
    double z=1e9;
    for(int i=0;i<vertices->GetEntriesFast();i++){
        EdbVertex *v = (EdbVertex *) vertices->At(i);
        if(v->Z()<z) {
            z = v->Z();
            v_1ry=v;
        }
    }
    int vtxpl = gEDA->GetIPLZ(v_1ry->Z());
    printf("Vertex plate = %d. z=%7.1lf\n", vtxpl, gEDA->GetZ(vtxpl));
    if(eOutputFileMode==kBern) fprintf(fp, "// primary vertex plate = %d. dz = %f\n", vtxpl, gEDA->GetZ(vtxpl)-v_1ry->Z());
    
    if(eOutputFileMode==kBern) fprintf(fp, "// VID     X        Y        Z 1ry charm tau Ndw Nup OutOfBrick\n");
    
    for(int i=0; i<vertices->GetEntriesFast(); i++){
        EdbVertex *v = (EdbVertex *) vertices->At(i);
        v->SetID(i+1);
        int ISPRIMARY = v==v_1ry ? 1 : 0;
        int ISCHARM = 0;
        int ISTAU = 0;
        int Nup=0, Ndown=0;
        for(int j=0;j<v->N();j++){
            EdbTrackP *t = v->GetTrack(j);
            if(t->GetSegmentFirst()->Z()>v->Z()) Ndown++;
            if(t->GetSegmentLast()->Z()<v->Z()) Nup++;
        }
        int OutOfBrick = v->Z()<gEDA->GetZ(1)? 1:0; // 1: if vertex is outside of the brick. 0: others.
        
        if(eOutputFileMode==kBern){
            fprintf(fp,"%3d %8.1lf %8.1lf %8.1lf %3d %5d %3d %3d %3d %3d\n", 
                    v->ID(), v->X(), v->Y(), v->Z(),ISPRIMARY,ISCHARM,ISTAU, Ndown, Nup, OutOfBrick);
        } else {
            fprintf(fp,"%d %8.1lf %8.1lf %8.1lf %3d %5d %3d %3d %3d %3d\n", 
                    v->ID(), v->X(), v->Y(), v->Z(),ISPRIMARY,ISCHARM,ISTAU, Ndown, Nup, OutOfBrick);
        }
    }
    
    // Tracks + Segments
    int itrk=1;
    for(int i=0;i<gEDA->NTrackSets();i++){
        EdbEDATrackSet *set = gEDA->GetTrackSet(i);
        if(set->N()==0) continue;
        if(set->GetDraw()==0) continue;
        printf("TrackSet %s : %3d tracks\n", set->GetName(), set->N()); 
        if(eOutputFileMode==kBern){
            EdbID& id=set->GetID();
            fprintf(fp,"////////////////////////////////////////////////////////// \n");
            fprintf(fp,"// ---------- TrackSet : %s %6d 0 %d %4d ----------- //\n", set->GetName(), id.eBrick, id.eMajor, id.eMinor);
            fprintf(fp,"////////////////////////////////////////////////////////// \n");
        }
        for(int j=0;j<set->N();j++){
            EdbTrackP *t = set->GetTrack(j);
            EdbSegP *sf = t->GetSegmentFirst();
            EdbSegP *sl = t->GetSegmentLast();
            EdbVertex * v1 = NULL;
            EdbVertex * v2 = NULL;
            for(int k=0;k<vertices->GetEntriesFast();k++){
                EdbVertex *v = (EdbVertex *) vertices->At(k);
                for(int l=0;l<v->N();l++){
                    EdbTrackP *tv = v->GetTrack(l);
                    if(t==tv){
                        if(NULL==v1) v1 = v;
                        else {
                            if(v1->Z()>v->Z()) { // v1 should 1ry.
                                v2 = v1;
                                v1 = v;
                            }
                            else v2 = v;
            }   }   }   }
            if(v1==NULL) v1=v_1ry; // all tracks are attached to 1ry vertex at this moment.
            
            // t  : pointer to the track
            // v1 : upstream vertex
            // v2 : downstream vertex
            // sf : first segment to be used for vertex calculation
            // fp : file pointer to print data.
            
            int ivtx1 = v1 ? v1->ID() : -1;
            int ivtx2 = v2 ? v2->ID() : -1;
            
            float ip1= v1? CalcIP(sf,v1) : -1.0; // v1->DistSeg(sf): -1.0;
            float ip2= v2? CalcIP(sl,v2) : -1.0; // v2->DistSeg(sl): -1.0;
            
            double p,pmin,pmax;
            CalcP(t, p, pmin, pmax); // pmin = pmin, pmax = pmax
            
            // Flags.... same as Feedback file.
            int MAN = 0;
            int SCANBACK=0, idscanback=-1;
            int MUON    =0;
            int DARKNESS=0; // no function to define.
            
            if(gEDA->GetTrackSet("MN")->FindTrack(t)) MAN=1;
            
            EdbTrackP * sbt = set==gEDA->GetTrackSet("SB")? t : gEDA->CheckScanback(t); // check if the track is a scanback track
            if(sbt!=NULL) idscanback = sbt->ID();
            if (idscanback!=-1) SCANBACK = 1;
            if (idscanback==gEDA->IdMuon()&&gEDA->IdMuon()!=-1) MUON =1;
            
            int OutOfBrick=0; // OutOfBrick: 1= pass through, 2=edge out, 0 others
            int LastPlate=0;  // Last plate: 1= pass though, n=edge out track last plate, 0 others 
            
            int nseg = t->N();
            
            // Small Kink search
 
            // remove Micro-tracks.
            double RmaxT=0.0, RmaxL=0.0, Rmax=0.0, dtmaxt=0.0, dtmaxl=0.0, rmstmax=-1.0, rmslmax=-1.0;
            int iplRmax1=0, iplRmax2=0;
            
            EdbTrackP *tt = new EdbTrackP;
            tt->Copy(*t);
            tt->SetCounters();
            
            for(int k=0;k<tt->N()&&tt->N()>1;k++){
                EdbSegP *s = tt->GetSegment(k);
                if(s->Side()!=0) tt->RemoveSegment(s); // if it's microtrack, remove.
            }
            
            int n = tt->N();
            
            if(n<=1) {
                // cannot calculate rms.
                RmaxT=RmaxL=-1.;
                iplRmax1=iplRmax2=-1;
            }
            else if( n<4 ){
                // if nseg <=3, cannot calculate kink, just calculate rms but not R.
                RmaxT=RmaxL=-1.;
                iplRmax1=iplRmax2=-1;
            }
            else {
                // Calculate rms and R.
                for(int k=0; k<n-1;k++){
                    EdbSegP *ss1 = tt->GetSegment(k);
                    EdbSegP *ss2 = tt->GetSegment(k+1);
                    
                    int ipl1 = ss1->Plate();
                    int ipl2 = ss2->Plate();
 
                    // if ss2 is on pl 57, skip.
                    if(ipl2>=gEDA->GetLastPlate()) continue;
                    // if ss1 is more downstream than 3 plates from vertex, skip. search only kink between 1-x,2-x,3-x.
                    if(ipl1>vtxpl+2) continue;
                    
                    double dtt,dtl;
                    CalcDTTransLongi(ss1, ss2, &dtt, &dtl);
                    dtt=fabs(dtt);
                    dtl=fabs(dtl);
                    
                    double r,rt,rl;
                    int ndata;
                    DTRMSTLGiven1Kink(tt, k, &r, &rt, &rl, &ndata);
 
                    double RT =dtt/rt/sqrt((double)ipl2-ipl1);
                    double RL =dtl/rl/sqrt((double)ipl2-ipl1);
                    double R=sqrt(RT*RT+RL*RL);
                    printf("ipl %2d-%2d dtt=%7.4lf dtl=%7.4lf rmst=%7.4lf rmsl=%7.4lf RT=%5.2lf RL=%5.2lf n=%d\n", 
                            ipl1, ipl2, dtt, dtl, rt, rl, RT,RL,ndata);
                    if(Rmax<R){
                        Rmax  = R;
                        RmaxT = RT;
                        RmaxL = RL;
                        iplRmax1=ipl1;
                        iplRmax2=ipl2;
                        dtmaxt = dtt;
                        dtmaxl = dtl;
                        rmstmax = rt;
                        rmslmax = rl;
                    }
                }
            }
            // nseg in 5 plate.
            int nin5pl = 0;
            for(int k=0;k<n;k++){
                EdbSegP *s = tt->GetSegment(k);
                if(s->Plate()>=vtxpl && s->Plate()<vtxpl+5) nin5pl++;
            }
            if( nin5pl<=2 ){
                // nseg <=2 in first 5 plate is out of criteria. just calculate RMS but not R.
                RmaxT=RmaxL=-1.;
                iplRmax1=iplRmax2=-1;
            }
            
            double rms, rmst, rmsl;
            if(n<=1) rms=rmst=rmsl=-1.;
            else DTRMSTL(tt, &rms, &rmst, &rmsl);
            
            delete tt;
            
            int result=1;
            
            const char *comment = set->GetComment(t);
            if(eOutputFileMode==kBern){
                // Bern format
                if(eOutputFileMode==kBern){
                    fprintf(fp,"// tid v1 v2        x        y        z      tx      ty    ip1    ip2 ");
                    fprintf(fp,"    p  pmin  pmax mn pa sb dk of ps ");
                    fprintf(fp," n RmaxT RmaxL    rmsT    rmsL pl1 pl2 res\n");
                }
                fprintf(fp, "%6d %2d %2d %8.1f %8.1f %8.1f %7.4f %7.4f %6.1f %6.1f ",
                    t->ID(), ivtx1, ivtx2, sf->X(), sf->Y(), sf->Z(), sf->TX(), sf->TY(), ip1, ip2);
                
            } else {
                // Official
                fprintf(fp, "%d %2d %2d %8.1f %8.1f %8.1f %7.4f %7.4f %6.1f %6.1f ",
                    itrk++, ivtx1, ivtx2, sf->X(), sf->Y(), sf->Z(), sf->TX(), sf->TY(), ip1, ip2);
            }
            fprintf(fp, "%5.1lf %5.1lf %5.1lf %2d %2d %2d %2d %2d %2d ", 
                p,pmin,pmax, MAN, MUON, SCANBACK, DARKNESS, OutOfBrick, LastPlate);
            fprintf(fp, "%2d %5.2lf %5.2lf %7.4lf %7.4lf %3d %3d %3d\n",
                nseg, RmaxT, RmaxL, rmst, rmsl, iplRmax1, iplRmax2, result);
 
            // print on the shell.
            fprintf(stdout, "%5d %2d %2d %8.1f %8.1f %8.1f %7.4f %7.4f %6.1f %6.1f ",
                t->ID(), ivtx1, ivtx2, sf->X(), sf->Y(), sf->Z(), sf->TX(), sf->TY(), ip1, ip2);
            fprintf(stdout, "%5.1lf %5.1lf %5.1lf %2d %2d %2d %2d %2d %2d ", 
                p,pmin,pmax, MAN, MUON, SCANBACK, DARKNESS, OutOfBrick, LastPlate);
            fprintf(stdout, "%2d %5.2lf %5.2lf %7.4lf %7.4lf %3d %3d %3d\n",
                nseg, RmaxT, RmaxL, rmst, rmsl, iplRmax1, iplRmax2, result);
            
            // Check if the rms is smaller than angular resolution or not.
            // if it's smaller, then print a note in the feedback file.
            
            // minimum value of rms (angular resolution).
            double angle = sqrt(t->TX()*t->TX()+t->TY()*t->TY());
            double resT  = 0.0015*sqrt(2.0);
            double resL  = 0.0015*(1+angle*5)*sqrt(2.0);
            double RresT = RmaxT>0 ? dtmaxt/resT : -1.;
            double RresL = RmaxL>0 ? dtmaxl/resL : -1.;
 
            if( eOutputFileMode==kBern ){
                if( ( rmst<resT || rmsl<resL ) && rmst>0)
                     fprintf(fp, "//     rms for R = %7.4lf %7.4lf , note : rms is smaller than angular resolution.  "
                                "RresT RresL resT resL = %5.2lf %5.2lf %7.4lf %7.4lf\n", 
                                rmstmax, rmslmax, RresT, RresL, resT, resL);
                else fprintf(fp, "//     rms for R = %7.4lf %7.4lf                                                   "
                                "RresT RresL resT resL = %5.2lf %5.2lf %7.4lf %7.4lf\n", 
                                rmstmax, rmslmax, RresT, RresL, resT, resL);
                
                if( nin5pl<3 ) fprintf(fp, "//     nBT = %d. nBT in first five plate = %d.\n", n, nin5pl);
                fprintf(fp,      "//     track comment = \"%s\"\n", comment);
            }
            
            // Fill microtrack information
            
            int irec = 0;    // 0:TS, 1:PredScan, 2: Manual
            TString name = set->GetName();
            if( name.BeginsWith("TS") ) irec=0;
            if( name.BeginsWith("SB") ) irec=1;
            if( name.BeginsWith("SF") ) irec=1;
            if( name.BeginsWith("MN") ) irec=2;
            
            for(int k=0;k<t->N();k++){
                EdbSegP *s = t->GetSegment(k);
 
                int type = s->Side(); // 0:BT, 1:MTtop, 2:MTbottom
                EdbID idseg=s->ScanID();
                EdbEDATrackSet *setseg = gEDA->GetTrackSet( idseg);
                if(setseg){
                    TString name = setseg->GetName();
                    if( name.BeginsWith("TS") ) irec=0;
                    if( name.BeginsWith("SB") ) irec=1;
                    if( name.BeginsWith("SF") ) irec=1;
                    if( name.BeginsWith("MN") ) irec=2;
                }
                
                if(eOutputFileMode==kBern){
                    // Bern format
                    if(k==0) fprintf(fp,"//       ipl        x        y        z      tx      ty type irec ngrains\n");
                    fprintf(fp,"%12d %8.1f %8.1f %8.1f %7.4f %7.4f %4d %4d %4d\n",
                        s->Plate(), s->X(), s->Y(), s->Z(), s->TX(), s->TY(), type, irec, (int) s->W());
                } else {
                    // Official
                    fprintf(fp,"%d %8.1f %8.1f %8.1f %7.4f %7.4f %4d %4d %4d\n",
                        s->Plate(), s->X(), s->Y(), s->Z(), s->TX(), s->TY(), type, irec, (int) s->W());
                }
            }
            
        }
    }
    
    // Scanning area
    
    EdbEDATrackSet *set = gEDA->GetTrackSet("TS");
    EdbEDAAreaSet *aset = set->GetAreaSet();
    if(aset->N()!=0){
        fprintf(fp,"//////////////////////////////////////////\n");
        fprintf(fp,"// ------ Area information : %s ------- //\n", set->GetName());
        fprintf(fp,"//////////////////////////////////////////\n");
        for(int i=0;i<aset->N();i++){
            EdbEDAArea *a = aset->GetArea(i);
            fprintf(fp,"// area : TS %2d %lf %lf %lf %lf %lf\n", 
                a->Plate(), a->Xmin(), a->Xmax(),
                a->Ymin(), a->Ymax(), a->Z());
        }
    }
    
    fclose(fp);
 
    printf("file %s was written.\n", filename);
 
}

WriteFeedbackFile2008()

void EdbEDAIO::WriteFeedbackFile2008	(	EdbVertex *	v = NULL,
		const char *	filename = NULL
	)

                                                                      {
 
// write feedback file. original code from Frederic.
// the flags are correct if you work in the correct directory.
 
    // extract vertex position and multiplicity
    // extract vertex tracks infos
    if(v==NULL) v = gEDA->GetSelectedVertex();
    if(filename==NULL){ // File dialog
        TGFileInfo *fi=new TGFileInfo;
        fi->fIniDir    = StrDup(".");
        const char *filetypes[] = { "Feedback file", "*.feedback", "All files","*",0,0};
        fi->fFileTypes = filetypes;
        new TGFileDialog(gClient->GetRoot(), gEve->GetMainWindow(), kFDSave, fi);
        if(fi->fFilename!=NULL) filename = fi->fFilename;
    }
    
    if(NULL==filename||strlen(filename)==0) { printf("Cancel to write feedback\n"); return;}
    FILE *fp=fopen(filename,"wt");
    if(NULL==fp){
        printf("Couldn't open file : %s\n", filename);
        return;
    }
 
 
    // put first the primary vertex
 
    fprintf(fp,"%i %.1f %.1f %.1f %i %i %i \n",1,v->X(),v->Y(),v->Z(),1,0,0);
 
    // tracks from primary vertex
    for(int j=0; j<v->N(); j++){
        int MANUAL=0, MUON=0, SCANBACK=0, DARKNESS=0;
        EdbTrackP *t = v->GetTrack(j);
        EdbSegP *s = t->GetSegmentFirst();
 
        // check if the track is a scanback track
        int idscanback = -1;
        EdbTrackP *tsb = gEDA->CheckScanback(t);
        if(tsb) idscanback = tsb->ID();
 
        int idmuon     = IdMuon();
        if (idscanback!=-1) SCANBACK = 1;
        if (idscanback==idmuon&&idmuon!=-1) {
            MUON =1;
            printf("      Track: %d is muon\n", t->ID());
        }
 
        double p, pmin, pmax;
        CalcP(t, p, pmin, pmax); // pmin = pmin, pmax = pmax
 
        fprintf(fp,"%4i %i %i %8.1f %8.1f %8.1f %7.4f %7.4f %6.2f %3.1f %6.2lf %5.2lf %8.1lf %i %i %i %i \n",
        t->ID(),1,-1,s->X(),s->Y(),s->Z(),s->TX(), s->TY(), v->Impact(j),-1.,p,pmin,pmax,MANUAL,MUON,SCANBACK,DARKNESS);
 
    }
    printf("-- FeedbackFile : %s--\n", filename);
    fclose(fp);
 
#ifdef _WINDOWS
    char buf[256];
    sprintf(buf, "type %s",filename);
    for(int i=0;i<strlen(buf);i++) if(buf[i]=='/') buf[i]='\\';
    gSystem->Exec(buf);
#else
    gSystem->Exec(Form("cat %s",filename));
#endif
}

WriteFilteredText()

void EdbEDAIO::WriteFilteredText

(

const char *

filename = NULL

)

                                                    {
    int i,j,k,l;
    FILE *fp;
    
    if(filename==NULL){ // File dialog
        TGFileInfo *fi=new TGFileInfo;
        fi->fIniDir    = StrDup(".");
        const char *filetypes[] = { "filtered tracks", "*.txt", "All files","*",0,0};
        fi->fFileTypes = filetypes;
        new TGFileDialog(gClient->GetRoot(), gEve->GetMainWindow(), kFDSave, fi);
        if(fi->fFilename!=NULL) {
            filename = fi->fFilename;
        }
    }
    if(filename==NULL) return;
 
    fp=fopen(filename,"wt");
    
    if(fp==NULL) {
        printf("File open error : %s\n", filename);
        return;
    }
    
    printf("Write Filtered Text File : %s\n", filename);
 
    // Vertices
    
    TObjArray *vertices = (gEDA->GetVertexSet())->GetVertices();
 
    EdbVertex *v_1ry=NULL;
    double z=1e9;
    for(i=0;i<vertices->GetEntriesFast();i++){
        EdbVertex *v = (EdbVertex *) vertices->At(i);
        if(v->Z()<z) {
            z = v->Z();
            v_1ry=v;
        }
    }
 
    fprintf(fp,"// vtx  id        X        Y        Z 1ry charm tau\n");
    for( i=0; i<vertices->GetEntriesFast(); i++){
        EdbVertex *v = (EdbVertex *) vertices->At(i);
        v->SetID(i);
        int ISPRIMARY = v==v_1ry ? 1 : 0;
        int ISCHARM = 0;
        int ISTAU = 0;
        fprintf(fp,"vertex %3d %8.1lf %8.1lf %8.1lf %3d %5d %3d\n", 
                    v->ID(), v->X(), v->Y(), v->Z(),ISPRIMARY,ISCHARM,ISTAU);
    }
    
    printf("vertices %d\n", vertices->GetEntriesFast());
    
//  ID_Track ID_UpVtx ID_DownVtx X Y Z SX SY IP_UpVtx IP_DownVtx P Pmin Pmax Manual Mu/e Scanback Darknes
 
    fprintf(fp,"// trk    id flag      tx      ty pl1 pl2 nseg   ip1ry vtx1   ip1 vtx2   ip2      p   pmin       pmax man part sb dark\n");
 
    // Tracks + Segments
    
    for(i=0;i<gEDA->NTrackSets();i++){
        EdbEDATrackSet *set = gEDA->GetTrackSet(i);
        if(set->GetDraw()==0) continue;
        for(j=0;j<set->N();j++){
            EdbTrackP *t = set->GetTrack(j);
            EdbSegP *sf = t->GetSegmentFirst();
            EdbSegP *sl = t->GetSegmentLast();
            EdbVertex * vtx1st = NULL;
            EdbVertex * vtx2nd = NULL;
            for(k=0;k<vertices->GetEntriesFast();k++){
                EdbVertex *v = (EdbVertex *) vertices->At(k);
                for(l=0;l<v->N();l++){
                    EdbTrackP *tv = v->GetTrack(l);
                    if(t==tv){
                        if(NULL==vtx1st) vtx1st = v;
                        else {
                            if(vtx1st->Z()>v->Z()) { // vtx1st should 1ry.
                                vtx2nd = vtx1st;
                                vtx1st = v;
                            }
                            else vtx2nd = v;
            }   }   }   }
 
            // Flags.... same as Feedback file.
            int MAN = 0;
            int SCANBACK=0, idscanback=-1;
            int MUON    =0;
            int DARKNESS=0; // no function to define.
 
            if(-99<=t->Flag()&&t->Flag()<0) MAN=1;
            
            EdbTrackP * sbt = gEDA->CheckScanback(t); // check if the track is a scanback track
            if(sbt!=NULL) idscanback = sbt->ID();
            if (idscanback!=-1) SCANBACK = 1;
            if (idscanback==IdMuon()&&IdMuon()!=-1) MUON =1;
 
            double p,pmin,pmax;
            CalcP(t, p, pmin, pmax); // pmin = pmin, pmax = pmax
            int ivtx1 = vtx1st ? vtx1st->ID() : -1;
            int ivtx2 = vtx2nd ? vtx2nd->ID() : -1;
            float ip1=-1.0, ip2=-1.0;
            if(vtx1st){
                for(k=0;k<vtx1st->N();k++){
                    EdbVTA *vta = vtx1st->GetVTa(k);
                    if(t==vta->GetTrack()) ip1 = vta->Imp();
                }
            }
            if(vtx2nd){
                for(k=0;k<vtx2nd->N();k++){
                    EdbVTA *vta = vtx2nd->GetVTa(k);
                    if(t==vta->GetTrack()) ip2 = vta->Imp();
                }
            }
            double ip_1ry = CalcIP(sf, v_1ry);
            
            fprintf(fp,"track %6d %4d %7.4f %7.4f %3d %3d %4d %7.1lf %4d %5.1f %4d %5.1f ",
                t->ID(), t->Flag(), t->TX(), t->TY(), sf->Plate(), sl->Plate(), t->N(), 
                ip_1ry, ivtx1, ip1, ivtx2, ip2);
            fprintf(fp,"%6.2lf %6.2lf %10.2lf %3d %4d %2d %4d\n", p,pmin,pmax, MAN, MUON, SCANBACK, DARKNESS);
        }
        printf("TrackSet : %s %d tracks\n", set->GetName(), set->N());
    }
 
    fprintf(fp,"// seg   trkid tflag ipl  segid        x        y        z      tx      ty  w   chi2 flag\n");
    
    for(i=0;i<gEDA->NTrackSets();i++){
        EdbEDATrackSet *set = gEDA->GetTrackSet(i);
        if(set->GetDraw()==0) continue;
        for(j=0;j<set->N();j++){
            EdbTrackP *t = set->GetTrack(j);
            for(k=0;k<t->N();k++){
                EdbSegP *s = t->GetSegment(k);
                fprintf(fp,"segment %6d %5d %3d %6d %8.1f %8.1f %8.1f %7.4f %7.4f %2d %6.4f %4d\n",
                    t->ID(), t->Flag(), s->Plate(), s->ID(), s->X(),s->Y(),s->Z(),s->TX(),s->TY(),
                    (int)s->W(),s->Chi2(),s->Flag());
            }
        }
    }
    
    
    // Scanning area
    fprintf(fp,"// area ipl xmin xmax ymin ymax z\n");
    EdbEDAAreaSet *aset = gEDA->GetTrackSet("TS")->GetAreaSet();
    for(i=0; i<aset->N(); i++){
        EdbEDAArea *a = aset->GetArea(i);
        fprintf(fp,"area %2d %8.1lf %8.1lf %8.1lf %8.1lf %8.1lf\n",
            a->Plate(), a->Xmin(), a->Xmax(), a->Ymin(), a->Ymax(), a->Z());
    }
    
    
    
    fclose(fp);
}

Member Data Documentation

eCutCP

TCut EdbEDAIO::eCutCP

eOutputFileMode

int EdbEDAIO::eOutputFileMode

---

The documentation for this class was generated from the following files:

• src/libEDA/EdbEDA.h
• src/libEDA/EdbEDA.C