FEDRA emulsion software from the OPERA Collaboration
Functions
run.C File Reference
#include <exception>
Include dependency graph for run.C:
This graph shows which files directly or indirectly include this file:

Functions

void Run (bool AbsCoord=false, float X0, float Y0, float DDX, float DDY, const char *File, bool Images=false)
 
int scanpreds (int PLATE, int MINOR)
 

Function Documentation

◆ Run()

void Run ( bool  AbsCoord = false,
float  X0,
float  Y0,
float  DDX,
float  DDY,
const char *  File,
bool  Images = false 
)
209{
210 // Scan emulsion zone with the use of faster async function TS->ScanAndFillView
211
212 // s->L->ON();
213 int res=0;
214 float X,Y;
215 int Vn=0;
216 int vid=0;
217 long start,end,cur;
218 EdbMark *xyz=new EdbMark();
219 xyz->SetX(X0);
220 xyz->SetY(Y0);
221 if(AbsCoord && FM->eFoundMarksSet) xyz->Transform(FM->eFoundMarksSet->Abs2Stage());
222 X0=xyz->X();
223 Y0=xyz->Y();
224 printf("Going to %f, %f\n",X0,Y0);
225 s->X->PosMove(X0,1e6,1e6);
226 s->Y->PosMove(Y0,1e6,1e6);
227 res=s->Y->WaitForMoveComplete(1e4);
228 res*=s->X->WaitForMoveComplete(1e4);
229 if(!res) return;
230 printf("Reached %f, %f\n",s->X->GetPos(),s->Y->GetPos());
231 start=atoi((gSystem->Now()).AsString());
232
233// if(!focus(Cluster_Threshold)) return;
234// finefocus(Cluster_Threshold);
235 end=atoi((gSystem->Now()).AsString());
236// printf("%d ms spent for focusing\n",end-start);
237
238 run=new EdbRun(File,"RECREATE");
239 run->SetTitle("Test of libACQ");
240 run->SetComment(" Test run: online data taking with libACQ. \n Parameters are selected to fit OPERA plates \n");
241 EdbRunHeader *header = run->GetHeader();
242 header->SetCCD( o->SizeX , // number of pixels
243 o->SizeY , //
244 -999 , // ccd pixel size in microns
245 -999 , // (is not important for analysis)
246 "MC1310",
247 "Mikrotron CMOS High Speed camera" );
248
249 header->SetObjective( 50. , // magnification nominal value
250 s->X->Invert*o->SizeX*o->PixelToMicronX, // pixel size X |
251 0. , // |
252 0. , // | CCD2Stage affine
253 s->Y->Invert*o->SizeY*o->PixelToMicronY, // pixel size Y | transformation:
254 0. , // | pixels -> microns
255 0. , // |
256 "NIKON",
257 "Dry Objective " );
258
259
260
261 header->SetArea( TS->Nx*TS->Ny , // number of views/area to scan
262 TS->ViewStepX , // X-step in microns
263 TS->ViewStepY , // Y-step in microns
264 TS->Layers , // NframesTop
265 TS->Layers , // NframesBot
266 0 ); //path
267
268 run->eMarks=FM->eFoundMarksSet;
269
270
271
272 header->SetPlate( 6 , // plate ID
273 s->EmTopLayer , // emulsion up
274 s->EmBase , // base
275 s->EmBotLayer , // emulsion down
276 1.0 , // shrinkage up
277 1.0 , // shrinkage down
278 "W2P6R",
279 "Welded brick2, plate 6" );
280
281// calculate fragment division
282int indFragX=0;
283int indFragY=0;
284int AreaID=0;
285float FragStepX = TS->ViewStepX * TS->Nx;
286float FragStepY = TS->ViewStepY * TS->Ny;
287int NumFragsX=int(DDX/FragStepX)+1;
288int NumFragsY=int(DDY/FragStepY)+1;
289
290for( indFragY=0; indFragY<NumFragsY; indFragY++)
291for( indFragX=0; indFragX<NumFragsX; indFragX++)
292{
293 X=X0+indFragX * FragStepX;
294 Y=Y0+indFragY * FragStepY;
295//finefocus(Cluster_Threshold,3);
296 if(!focus(Cluster_Threshold))
297 {
298 if(!finefocus(Cluster_Threshold,.4)) return;
299 if(!finefocus(Cluster_Threshold,1)) return;
300
301 }
303 o->Threshold=ThrTop;
304 TS->FastScanFragment(true, AreaID, AbsCoord, X,Y, run, Images);
306 o->Threshold=ThrBot;
307 TS->FastScanFragment(false, AreaID, AbsCoord, X,Y, run, Images);
308 AreaID++;
309
310}
311 start=atoi((gSystem->Now()).AsString());
312 printf("%d ms spent for scanning.\n",start-end);
313 run->Save();
314 run->Close();
315 delete xyz;
316
317}
X0
brick X0
Definition: RecDispMC.C:112
AcqMarksScanner::eFoundMarksSet
EdbMarksSet * eFoundMarksSet
Definition: AcqMarksScanner.h:41
AcqOdyssey::SizeX
int SizeX
Definition: AcqOdyssey.h:168
AcqOdyssey::PixelToMicronY
float PixelToMicronY
Definition: AcqOdyssey.h:179
AcqOdyssey::Threshold
int Threshold
Definition: AcqOdyssey.h:174
AcqOdyssey::SizeY
int SizeY
Definition: AcqOdyssey.h:168
AcqOdyssey::PixelToMicronX
float PixelToMicronX
Definition: AcqOdyssey.h:178
AcqTrackScanner::Layers
int Layers
Definition: AcqTrackScanner.h:30
AcqTrackScanner::Ny
int Ny
Definition: AcqTrackScanner.h:35
AcqTrackScanner::FastScanFragment
int FastScanFragment(bool IsTop, int AreaId, bool AbsCoord, float X0, float Y0, EdbRun *run, bool RawImages=false)
Definition: AcqTrackScanner.cxx:181
AcqTrackScanner::ViewStepY
float ViewStepY
Definition: AcqTrackScanner.h:34
AcqTrackScanner::Nx
int Nx
Definition: AcqTrackScanner.h:35
AcqTrackScanner::ViewStepX
float ViewStepX
Definition: AcqTrackScanner.h:34
EdbMark
Definition: EdbFiducial.h:21
EdbMark::Y
Float_t Y() const
Definition: EdbFiducial.h:43
EdbMark::SetY
void SetY(float y)
Definition: EdbFiducial.h:45
EdbMark::SetX
void SetX(float x)
Definition: EdbFiducial.h:44
EdbMark::X
Float_t X() const
Definition: EdbFiducial.h:42
EdbMarksSet::Abs2Stage
EdbAffine2D * Abs2Stage() const
Definition: EdbFiducial.cxx:271
EdbPoint2D::Transform
virtual void Transform(const EdbAffine2D *a)
Definition: EdbVirtual.cxx:107
EdbRunHeader
Definition: EdbRunHeader.h:95
EdbRunHeader::SetArea
void SetArea(int N, float stepx, float stepy, int ft, int fb, int path)
Definition: EdbRunHeader.h:209
EdbRunHeader::SetCCD
void SetCCD(int nx=0, int ny=0, float px=0, float py=0, char *name="", char *title="")
stage setting functions
Definition: EdbRunHeader.h:179
EdbRunHeader::SetPlate
void SetPlate(int id, float up, float base, float down, float shru, float shrd, char *name="", char *title="")
Definition: EdbRunHeader.h:200
EdbRunHeader::SetObjective
void SetObjective(float mag, float a, float b, float c, float d, float p, float q, char *name="", char *title="")
Definition: EdbRunHeader.h:191
EdbRun
Definition: EdbRun.h:75
EdbRun::eMarks
EdbMarksSet * eMarks
fiducial marks
Definition: EdbRun.h:101
EdbRun::Close
void Close()
Definition: EdbRun.cxx:439
EdbRun::Save
void Save()
Definition: EdbRun.cxx:424
EdbRun::SetTitle
void SetTitle(const char *tit)
Definition: EdbRun.h:133
EdbRun::SetComment
void SetComment(const char *com)
Definition: EdbRun.h:132
EdbRun::GetHeader
EdbRunHeader * GetHeader() const
Definition: EdbRun.h:138
ThrBot
int ThrBot
Definition: hwinit.C:12
FM
AcqMarksScanner * FM
Definition: hwinit.C:3
o
AcqOdyssey * o
Definition: hwinit.C:2
Cluster_Threshold
int Cluster_Threshold
Definition: hwinit.C:13
ThrTop
int ThrTop
Definition: hwinit.C:11
TS
AcqTrackScanner * TS
Definition: hwinit.C:4
focus
int focus(int clu_thre=300)
Definition: init.C:76
finefocus
int finefocus(int clu_thre=300, float fineness=4.)
Definition: init.C:109
run
EdbRun * run
Definition: check_raw.C:38
s
s
Definition: check_shower.C:55
File
Definition: tlg2couples.C:51
Header::Center::X
Double_t X
Definition: tlg2couples.C:76
Header::Center::Y
Double_t Y
Definition: tlg2couples.C:76

◆ scanpreds()

int scanpreds ( int  PLATE,
int  MINOR 
)
5{
6int BRICK=92;
7
8bool loadresult;
9try {loadresult=PC->LoadPlate(BRICK,PLATE,0,0); printf("Load plate returned %d.\n",loadresult);}
10catch (...) {printf("Exception caught in LoadPlate!\n"); return 0;}
11if(!loadresult) {printf("can't load plate! Exiting.\n"); return 0;}
12try { printf("Scanmarks is about to execute..\n"); loadresult=scanmarks();}
13catch (...) {cout << "Exception caught in scanmarks!" << endl; return 0;}
14if(!loadresult) return 0;
15 char cmd[128];
16// Initializing run
17 try
18 {
19 EdbMark *xyz=new EdbMark();
20 if(!FM->eFoundMarksSet) {printf("Marks are not scanned! Exiting.\n"); return 0;}
21 start=atoi((gSystem->Now()).AsString());
22 end=atoi((gSystem->Now()).AsString());
23 sprintf(cmd,"del .\\data\\r.root");
24 gSystem->Exec(cmd);
25 run=new EdbRun(".\\data\\r.root","RECREATE");
26 run->SetTitle("Test of prediction scan");
27 run->SetComment(" Test run: online data taking with libACQ. \n Parameters are selected to fit OPERA plates \n");
28 EdbRunHeader *header = run->GetHeader();
29 header->SetCCD( o->SizeX , // number of pixels
30 o->SizeY , //
31 -999 , // ccd pixel size in microns
32 -999 , // (is not important for analysis)
33 "MC1310",
34 "Mikrotron CMOS High Speed camera" );
35
36 header->SetObjective( 50. , // magnification nominal value
37 s->X->Invert*o->SizeX*o->PixelToMicronX, // pixel size X |
38 0. , // |
39 0. , // | CCD2Stage affine
40 s->Y->Invert*o->SizeY*o->PixelToMicronY, // pixel size Y | transformation:
41 0. , // | pixels -> microns
42 0. , // |
43 "NIKON",
44 "Dry Objective " );
45
46
47
48 header->SetArea( TS->Nx*TS->Ny , // number of views/area to scan
49 TS->ViewStepX , // X-step in microns
50 TS->ViewStepY , // Y-step in microns
51 TS->Layers , // NframesTop
52 TS->Layers , // NframesBot
53 0 ); //path
54
55 if(run->eMarks) delete (run->eMarks);
56 (run->eMarks)=FM->eFoundMarksSet;
57 printf("Markset written to file:\n");
58 (run->eMarks)->Print();
59
60
61
62 header->SetPlate( 6 , // plate ID
63 s->EmTopLayer , // emulsion up
64 s->EmBase , // base
65 s->EmBotLayer , // emulsion down
66 1.0 , // shrinkage up
67 1.0 , // shrinkage down
68 "W2P6R",
69 "Welded brick2, plate 6" );
70 }
71 catch (...) {cout << "Exception caught in Run initialization!" << endl; return 0;}
72
73// run initialized
74
75int scan_zone_id=98;
76int proc_zone_id=94;
77bool data_ready=false;
78int proctracks=0;
79int dummy;
80
81char fname[128];
82sprintf(fname,".\\data\\pred%d.%d.0.%d.txt",BRICK,PLATE,MINOR);
83FILE* fp=fopen(fname,"r");
84if(fp==NULL) {printf("Prediction file %s can't be opened.\n",fname); return 0;}
85char line[256];
86int id;
87float XX0,YY0,X,Y,TX,TY;
88int AreaID=0;
89
90try while(fgets(line,256,fp)!=NULL) //main loop
91 {
92// if(kbhit()){getch();printf("Paused. Continue? (y/n)"); if(getchar()=='n') break;}
93 sscanf(line,"%d %f %f %f %f",&id,&X,&Y,&TX,&TY);
94 printf("Predicted segment N %d: %f %f %f %f\n",id,X,Y,TX,TY);
95 xyz->SetX(X);
96 xyz->SetY(Y);
97 xyz->Transform(FM->eFoundMarksSet->Abs2Stage());
98 X=xyz->X();
99 Y=xyz->Y();
100 printf("Stage coordinates: %f %f \n",X,Y);
101 if(proctracks==0)
102 if(!focus(Cluster_Threshold))
103 {
104 if(!finefocus(Cluster_Threshold,.4)) return 0;
105 if(!finefocus(Cluster_Threshold,1)) return 0;
106
107 }
108
109 if(proctracks>0)
110 {
111 if(!focus(Cluster_Threshold))
112 {
113 if(!finefocus(Cluster_Threshold,.4)) {
114 s->X->PosMove(X,1e9,1e9);
115 s->Y->PosMove(Y,1e9,1e9);
116 s->X->WaitForMoveComplete(1e4);
117 s->Y->WaitForMoveComplete(1e4);
118 XX0=X; YY0=Y;
119 proctracks++;
120 continue;
121 }
122 if(!finefocus(Cluster_Threshold,1)) {
123 s->X->PosMove(X,1e9,1e9);
124 s->Y->PosMove(Y,1e9,1e9);
125 s->X->WaitForMoveComplete(1e4);
126 s->Y->WaitForMoveComplete(1e4);
127 XX0=X; YY0=Y;
128 proctracks++;
129 continue;
130 }
131
132 }
133
134 TS->Layers=Layers;
136 o->Threshold=ThrTop;
137 // TS->FocusScanPrediction(true, AreaID, true, XX0, YY0, run, false, 0, 0);
138 TS->FastScanPrediction(true, AreaID, true, XX0, YY0, run, false, 0, 0);
139// o->ShowImage(4);
140// AreaID++;
142 o->Threshold=ThrBot;
143// TS->FocusScanPrediction(false, AreaID, true, XX0, YY0, run, false, X, Y);
144 TS->FastScanPrediction(false, AreaID, true, XX0, YY0, run, false, X, Y);
145// o->ShowImage(4);
146 AreaID++;
147 }
148 XX0=X; YY0=Y; //to remember previous
149 proctracks++;
150 printf("Prepared to scan %d predictions\n",proctracks);
151 }
152
153 catch (...) {cout << "Exception caught in main loop!" << endl; return 0;}
154
155 // eND OF MAIN LOOP
156 try {
157 if(!focus(Cluster_Threshold))
158 {
159 if(!finefocus(Cluster_Threshold,.4)) return 0;
160 if(!finefocus(Cluster_Threshold,1)) return 0;
161
162 }
163
164 TS->Layers=Layers;
166 o->Threshold=ThrTop;
167// TS->FocusScanPrediction(true, AreaID, true, X, Y, run, false, 0, 0);
168 TS->FastScanPrediction(true, AreaID, true, X, Y, run, false, 0, 0);
169// o->ShowImage(4);
170// AreaID++;
172 o->Threshold=ThrBot;
173 // TS->FocusScanPrediction(false, AreaID, true, X, Y, run, false, 0, 0);
174 TS->FastScanPrediction(false, AreaID, true, X, Y, run, false, 0, 0);
175// o->ShowImage(4);
176 AreaID++;
177
178run->Save();
179run->Close();
180delete run;
181 }
182 catch (...) {cout << "Exception caught in Postprocessing!" << endl; return 0;}
183
184char rootfname[128];
185 sprintf(rootfname,"data\\tst%d.%d.0.%d.root",BRICK,PLATE,MINOR);
186//else sprintf(rootfname,"%s",ofile);
187printf("move data\\r.root %s\n",rootfname);
188sprintf(cmd,"move data\\r.root %s",rootfname);
189gSystem->Exec(cmd);
190/*Link1(BRICK,PLATE,98,false);
191 sprintf(cmd,"mv data/%02d_098.cp.root data/%02d_%03d.cp.root ",PLATE,PLATE,MINOR);
192 printf("mv data/%02d_098.cp.root data/%02d_%03d.cp.root\n ",PLATE,PLATE,MINOR);
193 gSystem->Exec(cmd);
194 sprintf(cmd,"mv data/tst%d.%d.0.%d.root data/tst%d.%d.0.%d.root",BRICK,PLATE,98,BRICK,PLATE,MINOR);
195 printf("mv data/tst%d.%d.0.%d.root data/tst%d.%d.0.%d.root\n",BRICK,PLATE,98,BRICK,PLATE,MINOR);
196 gSystem->Exec(cmd);
197*/
198
199fclose(fp);
200start=atoi((gSystem->Now()).AsString());
201printf("%d ms spent for scanning.\n",start-end);
202delete xyz;
203return 1;
204}
AcqSerialPlateChanger::LoadPlate
int LoadPlate(IdentifierS id)
Definition: AcqSerialPlateChanger.cxx:783
Layers
int Layers
Definition: hwinit.C:65
PC
AcqSerialPlateChanger * PC
Definition: hwinit.C:5
scanmarks
int scanmarks()
Definition: init.C:368
Print
Alg Print()
fname
const char * fname
Definition: mc2raw.cxx:41
dummy
Int_t dummy
Definition: merge_Energy_SytematicSources_Electron.C:7
fclose
fclose(pFile)
Header::Transform::TY
Double_t TY
Definition: tlg2couples.C:78
Header::Transform::TX
Double_t TX
Definition: tlg2couples.C:78
MIPBaseTrack::id
UInt_t id
Definition: tlg2couples.C:117
SHOWER_PAR::PLATE
int PLATE
Definition: shower_btr.C:18
NULL
#define NULL
Definition: nidaqmx.h:84
BRICK
Definition: RecDispMC.C:95