FEDRA emulsion software from the OPERA Collaboration
Namespaces | Enumerations | Functions | Variables
tlg2couples.C File Reference
#include "Riostream.h"
#include "TFile.h"
#include "EdbPVRec.h"
#include <iostream>
Include dependency graph for tlg2couples.C:

Namespaces

 File
 
 File::Format
 
 File::Info
 
 File::Section
 
 Header
 
 Header::Center
 
 Header::Extents
 
 Header::Identifier
 
 Header::Transform
 
 MIPBaseTrack
 
 MIPEmulsionTrack
 
 TrackIndexEntry
 
 View
 

Enumerations

enum  { File::Info::Track = 1 , File::Info::BaseTrack = 2 , File::Info::Field = 3 }
 
enum  { File::Section::Data = 0x20 , File::Section::Header = 0x40 }
 
enum  {
  File::Format::Old = 0x08 , File::Format::Old2 = 0x02 , File::Format::NoExtents = 0x01 , File::Format::Normal = 0x03 ,
  File::Format::NormalWithIndex = 0x04 , File::Format::NormalDouble_tWithIndex = 0x05 , File::Format::Detailed = 0x06 , File::Format::MultiSection = 0x07
}
 
enum  en_side { TOP , BOT , TOP , BOT }
 

Functions

void Add_MIPBaseTrack (int i, float project_dz=-105.)
 
void Add_MIPEmulsionTrack (en_side side, int i)
 
int main (int argc, char *argv[])
 
void print_Header ()
 
void print_MIPBaseTrack (int i)
 
void print_MIPEmulsionTrack (int s, int i)
 
void print_View (int side, int i)
 
void printrackIndexEntry (int s, int i)
 
TClonesArray segbot ("EdbSegP")
 
TClonesArray segtop ("EdbSegP")
 
void tlg2couples (const char *input_file, char *output_file=0, float project_dz=-105.)
 

Variables

TTree * couples
 
EdbSegCouple * cp = 0
 
Int_t Pid1 =0
 
Int_t Pid2 =0
 
EdbSegP * s = 0
 
EdbSegP * s1 = 0
 
EdbSegP * s2 = 0
 
Float_t Xv =0
 
Float_t Yv =0
 

Enumeration Type Documentation

◆ en_side

enum en_side
Enumerator
TOP 
BOT 
TOP 
BOT 
36 { TOP , BOT } ;
TOP
@ TOP
Definition: tlg2couples.C:36
BOT
@ BOT
Definition: tlg2couples.C:36

Function Documentation

◆ Add_MIPBaseTrack()

void Add_MIPBaseTrack ( int  i,
float  project_dz = -105. 
)
235 {
236  using namespace MIPBaseTrack ;
237 
238  EdbSegP* t = (EdbSegP*) segtop[ TopId ] ;
239 
240  EdbSegP* b = (EdbSegP*) segbot[ BottomId ] ;
241 
242  float x_ = X + SX*( project_dz ) ;
243  float y_ = Y + SY*( project_dz ) ;
244 
245  float z_ = Z + project_dz ;
246 
247  s = new EdbSegP( i , x_, y_, SX, SY, Count ) ;
248  s->SetChi2( Sigma ) ;
249  s->SetVolume( AreaSum ) ;
250  s->SetZ( z_ ) ;
251 
252  s->SetAid ( -90 , -90 );
253  s->SetFlag( -90 );
254  s->SetDZ( b->Z() - t->Z() );
255  s->SetDZem( 0 ) ;
256 
257  Pid1=1;
258  Pid2=2;
259  Xv=0;
260  Yv=0;
261  cp= new EdbSegCouple( TopId, BottomId );
262  cp->SetN1(1);
263  cp->SetN2(1);
264  cp->SetN1tot(1);
265  cp->SetN2tot(1);
266  cp->SetCHI2P( Sigma ) ;
267 
268  s1 = new EdbSegP( *t );
269 
270  s2 = new EdbSegP( *b );
271 
272  couples->Fill();
273 }
t
TTree * t
Definition: check_shower.C:4
EdbSegCouple
Definition: EdbSegCouple.h:14
EdbSegCouple::SetN1
void SetN1(int n1)
Definition: EdbSegCouple.h:40
EdbSegCouple::SetN2
void SetN2(int n2)
Definition: EdbSegCouple.h:41
EdbSegCouple::SetN1tot
void SetN1tot(int n)
Definition: EdbSegCouple.h:42
EdbSegCouple::SetN2tot
void SetN2tot(int n)
Definition: EdbSegCouple.h:43
EdbSegCouple::SetCHI2P
void SetCHI2P(float chi2)
Definition: EdbSegCouple.h:45
EdbSegP
Definition: EdbSegP.h:18
EdbSegP::SetVolume
void SetVolume(float w)
Definition: EdbSegP.h:133
EdbSegP::SetDZem
void SetDZem(float dz)
Definition: EdbSegP.h:124
EdbSegP::Z
Float_t Z() const
Definition: EdbSegP.h:150
EdbSegP::SetZ
void SetZ(float z)
Definition: EdbSegP.h:122
EdbSegP::SetChi2
void SetChi2(float chi2)
Definition: EdbSegP.h:132
EdbSegP::SetDZ
void SetDZ(float dz)
Definition: EdbSegP.h:123
EdbSegP::SetFlag
void SetFlag(int flag)
Definition: EdbSegP.h:127
EdbSegP::SetAid
void SetAid(int a, int v, int side=0)
Definition: EdbSegP.h:135
MIPBaseTrack
Definition: tlg2pattern.C:117
MIPBaseTrack::TopId
Int_t TopId
Definition: tlg2pattern.C:128
MIPBaseTrack::Count
UInt_t Count
Definition: tlg2pattern.C:120
MIPBaseTrack::Y
Double_t Y
Definition: tlg2pattern.C:122
MIPBaseTrack::Z
Double_t Z
Definition: tlg2pattern.C:123
MIPBaseTrack::SX
Double_t SX
Definition: tlg2pattern.C:124
MIPBaseTrack::Sigma
Double_t Sigma
Definition: tlg2pattern.C:127
MIPBaseTrack::AreaSum
UInt_t AreaSum
Definition: tlg2pattern.C:119
MIPBaseTrack::X
Double_t X
Definition: tlg2pattern.C:121
MIPBaseTrack::SY
Double_t SY
Definition: tlg2pattern.C:125
MIPBaseTrack::BottomId
Int_t BottomId
Definition: tlg2pattern.C:129
Xv
Float_t Xv
Definition: tlg2couples.C:27
segtop
TClonesArray segtop("EdbSegP")
s1
EdbSegP * s1
Definition: tlg2couples.C:29
Yv
Float_t Yv
Definition: tlg2couples.C:27
couples
TTree * couples
Definition: tlg2couples.C:25
Pid2
Int_t Pid2
Definition: tlg2couples.C:26
s
EdbSegP * s
Definition: tlg2couples.C:31
segbot
TClonesArray segbot("EdbSegP")
s2
EdbSegP * s2
Definition: tlg2couples.C:30
Pid1
Int_t Pid1
Definition: tlg2couples.C:26
cp
EdbSegCouple * cp
Definition: tlg2couples.C:28

◆ Add_MIPEmulsionTrack()

void Add_MIPEmulsionTrack ( en_side  side,
int  i 
)
203 {
204  using namespace MIPEmulsionTrack ;
205 
206  if ( side == TOP) {
207  new( segtop[i] ) EdbSegP( i , X, Y, SX, SY, Count ) ;
208  ((EdbSegP*) segtop[i])->SetZ( Z );
209  //((EdbSegP*) segtop[i])->SetZ( Header::Top_BottomZ );
210  ((EdbSegP*) segtop[i])->SetChi2( Sigma );
211  ((EdbSegP*) segtop[i])->SetVolume( AreaSum ) ;
212 
213  ((EdbSegP*) segtop[i])->SetAid ( -91 , viewid );
214  ((EdbSegP*) segtop[i])->SetFlag( -91 );
215  ((EdbSegP*) segtop[i])->SetDZ( TopZ-BottomZ );
216  ((EdbSegP*) segtop[i])->SetDZem( 0 );
217  }
218 
219  if ( side == BOT) {
220  new( segbot[i] ) EdbSegP( i , X, Y, SX, SY, Count ) ;
221  ((EdbSegP*) segbot[i])->SetZ( Z );
222  //(EdbSegP*) segtop[i])->SetZ( Header::Bottom_TopZ );
223  ((EdbSegP*) segbot[i])->SetChi2( Sigma );
224  ((EdbSegP*) segbot[i])->SetVolume( AreaSum ) ;
225 
226  ((EdbSegP*) segbot[i])->SetAid ( -92 , viewid );
227  ((EdbSegP*) segbot[i])->SetFlag( -92 );
228  ((EdbSegP*) segbot[i])->SetDZ( TopZ-BottomZ );
229  ((EdbSegP*) segbot[i])->SetDZem( 0 );
230  }
231 }
MIPEmulsionTrack
Definition: tlg2pattern.C:99
MIPEmulsionTrack::Z
Double_t Z
Definition: tlg2pattern.C:105
MIPEmulsionTrack::BottomZ
Double_t BottomZ
Definition: tlg2pattern.C:111
MIPEmulsionTrack::Sigma
Double_t Sigma
Definition: tlg2pattern.C:109
MIPEmulsionTrack::X
Double_t X
Definition: tlg2pattern.C:103
MIPEmulsionTrack::AreaSum
UInt_t AreaSum
Definition: tlg2pattern.C:101
MIPEmulsionTrack::Count
UInt_t Count
Definition: tlg2pattern.C:102
MIPEmulsionTrack::SY
Double_t SY
Definition: tlg2pattern.C:107
MIPEmulsionTrack::SX
Double_t SX
Definition: tlg2pattern.C:106
MIPEmulsionTrack::viewid
Int_t viewid
Definition: tlg2pattern.C:113
MIPEmulsionTrack::TopZ
Double_t TopZ
Definition: tlg2pattern.C:110
MIPEmulsionTrack::Y
Double_t Y
Definition: tlg2pattern.C:104

◆ main()

int main ( int  argc,
char *  argv[] 
)
794 {
795  if(argc==2)
796  tlg2couples( argv[1], 0, -105. );
797  else if(argc==3)
798  tlg2couples( argv[1], argv[2], -105. );
799  else if(argc==4) {
800  float project_dz ;
801  sscanf( argv[3], "%f", &project_dz );
802  tlg2couples( argv[1], argv[2], project_dz );
803  } else {
804  cout << "tlg2couples.exe <input_file> [output_file] [project_dz (default= -105.)]\n";
805  return 0 ;
806  }
807 
808  return 1 ;
809 }
tlg2couples
void tlg2couples(const char *input_file, char *output_file=0, float project_dz=-105.)
Definition: tlg2couples.C:276

◆ print_Header()

void print_Header ( )
142 {
143 
144  using namespace Header ;
145  cout << "Header"<< endl;
146  cout << (Int_t) infotype << endl;
147  cout << headerformat << endl;
148 
149  using namespace Header::Identifier;
150  cout << Part0 <<" "<< Part1 <<" "<< Part2 <<" "<< Part3 << endl;
151  using namespace Header::Center;
152  cout << X <<" "<< Y << endl;
153  using namespace Header::Extents ;
154  cout << MinX <<" "<< MaxX <<" "<< MinY <<" "<< MaxY << endl;
155  using namespace Header::Transform ;
156  cout << MXX <<" "<< MXY <<" "<< MYX <<" "<< MYY << endl;
157  cout << TX <<" "<< TY <<" "<< RX <<" "<< RY << endl;
158 
159  using namespace Header;
160  cout << TopViewsLength <<" "<< BottomViewsLength << endl ;
161  cout << TopTracksLength <<" "<< BottomTracksLength <<" "<< TracksLength << endl;
162  cout << Top_TopZ <<" "<< Top_BottomZ <<" "<< Bottom_TopZ <<" "<< Bottom_BottomZ << endl;
163 
164 }
Header::Center
Definition: tlg2pattern.C:77
Header::Center::X
Double_t X
Definition: tlg2pattern.C:77
Header::Center::Y
Double_t Y
Definition: tlg2pattern.C:77
Header::Extents
Definition: tlg2pattern.C:78
Header::Extents::MaxX
Double_t MaxX
Definition: tlg2pattern.C:78
Header::Extents::MinX
Double_t MinX
Definition: tlg2pattern.C:78
Header::Extents::MinY
Double_t MinY
Definition: tlg2pattern.C:78
Header::Extents::MaxY
Double_t MaxY
Definition: tlg2pattern.C:78
Header::Identifier
Definition: tlg2pattern.C:76
Header::Identifier::Part1
Int_t Part1
Definition: tlg2pattern.C:76
Header::Identifier::Part3
Int_t Part3
Definition: tlg2pattern.C:76
Header::Identifier::Part0
Int_t Part0
Definition: tlg2pattern.C:76
Header::Identifier::Part2
Int_t Part2
Definition: tlg2pattern.C:76
Header::Transform
Definition: tlg2pattern.C:79
Header::Transform::RX
Double_t RX
Definition: tlg2pattern.C:79
Header::Transform::MYX
Double_t MYX
Definition: tlg2pattern.C:79
Header::Transform::TY
Double_t TY
Definition: tlg2pattern.C:79
Header::Transform::MXX
Double_t MXX
Definition: tlg2pattern.C:79
Header::Transform::TX
Double_t TX
Definition: tlg2pattern.C:79
Header::Transform::RY
Double_t RY
Definition: tlg2pattern.C:79
Header::Transform::MYY
Double_t MYY
Definition: tlg2pattern.C:79
Header::Transform::MXY
Double_t MXY
Definition: tlg2pattern.C:79
Header
Definition: tlg2pattern.C:66
Header::infotype
Char_t infotype
Definition: tlg2pattern.C:69
Header::Bottom_TopZ
Double_t Bottom_TopZ
Definition: tlg2pattern.C:74
Header::BottomViewsLength
Int_t BottomViewsLength
Definition: tlg2pattern.C:82
Header::headerformat
UShort_t headerformat
Definition: tlg2pattern.C:70
Header::TopTracksLength
Int_t TopTracksLength
Definition: tlg2pattern.C:84
Header::TracksLength
Int_t TracksLength
Definition: tlg2pattern.C:86
Header::Top_TopZ
Double_t Top_TopZ
Definition: tlg2pattern.C:72
Header::Bottom_BottomZ
Double_t Bottom_BottomZ
Definition: tlg2pattern.C:75
Header::Top_BottomZ
Double_t Top_BottomZ
Definition: tlg2pattern.C:73
Header::BottomTracksLength
Int_t BottomTracksLength
Definition: tlg2pattern.C:85
Header::TopViewsLength
Int_t TopViewsLength
Definition: tlg2pattern.C:81

◆ print_MIPBaseTrack()

void print_MIPBaseTrack ( int  i)
184 {
185  using namespace MIPBaseTrack ;
186  cout << "MIPBaseTrack " << i <<" "<< id <<" " ;
187  cout << AreaSum <<" "<< Count <<" "
188  << X <<" "<< Y <<" "<< Z <<" "
189  << SX <<" "<< SY <<" "<< SZ <<" "
190  << Sigma <<" "<< TopId <<" "<< BottomId << endl;
191 } ;
MIPBaseTrack::SZ
Double_t SZ
Definition: tlg2pattern.C:126

◆ print_MIPEmulsionTrack()

void print_MIPEmulsionTrack ( int  s,
int  i 
)
174 {
175  using namespace MIPEmulsionTrack ;
176  cout << "MIPEmulsionTrack " << s <<" "<< i <<" " ;
177  cout << Field <<" "<< AreaSum <<" "<< Count <<" "
178  << X <<" "<< Y <<" "<< Z <<" "
179  << SX <<" "<< SY <<" "<< SZ <<" "
180  << Sigma <<" "<< TopZ <<" "<< BottomZ <<" "<< viewid << endl;
181 } ;
MIPEmulsionTrack::SZ
Double_t SZ
Definition: tlg2pattern.C:108
MIPEmulsionTrack::Field
UInt_t Field
Definition: tlg2pattern.C:100

◆ print_View()

void print_View ( int  side,
int  i 
)
167 {
168  using namespace View ;
169  //if (i<2 || i>848)
170  cout << "View " << side <<" "<< id <<" "<< topz <<" "<< bottomz <<" "<< px <<" "<< py << endl;
171 } ;
View::py
Double_t py
Definition: tlg2pattern.C:95
View::px
Double_t px
Definition: tlg2pattern.C:94
View::topz
Double_t topz
Definition: tlg2pattern.C:92
View::bottomz
Double_t bottomz
Definition: tlg2pattern.C:93
View
Definition: tlg2pattern.C:90

◆ printrackIndexEntry()

void printrackIndexEntry ( int  s,
int  i 
)
194 {
195  using namespace TrackIndexEntry ;
196  cout << "TrackIndexEntry " << s <<" "<< i <<" "<< TrackIndexEntry::Fragment <<" "
197  << TrackIndexEntry::View <<" "<< TrackIndexEntry::Track << endl;
198 } ;
TrackIndexEntry
Definition: tlg2pattern.C:133
TrackIndexEntry::View
Int_t View
Definition: tlg2pattern.C:136
TrackIndexEntry::Fragment
Int_t Fragment
Definition: tlg2pattern.C:135
TrackIndexEntry::Track
Int_t Track
Definition: tlg2pattern.C:137

◆ segbot()

TClonesArray segbot ( "EdbSegP"  )

◆ segtop()

TClonesArray segtop ( "EdbSegP"  )

◆ tlg2couples()

void tlg2couples ( const char *  input_file,
char *  output_file = 0,
float  project_dz = -105. 
)
279 {
280  if (!output_file) {
281  output_file = new char[256];
282  sprintf(output_file, "%s", input_file );
283  char* ext = strrchr(output_file,'.') ;
284  sprintf(ext,".cp.root\0");
285  }
286 
287  // CREATING OUTPUT FILE
288  TFile *o = new TFile( output_file, "RECREATE");
289 
290  couples = new TTree("couples","couples");
291 
292  couples->Branch( "pid1" , &Pid1 , "pid1/I" ) ;
293  couples->Branch( "pid2" , &Pid2 , "pid2/I" ) ;
294  couples->Branch( "xv" , &Xv , "xv/F" ) ;
295  couples->Branch( "yv" , &Yv , "yv/F" ) ;
296 
297  couples->Branch( "cp" , "EdbSegCouple" , &cp ,32000 ,99 ) ;
298  couples->Branch( "s1." , "EdbSegP" , &s1 ,32000 ,99 ) ;
299  couples->Branch( "s2." , "EdbSegP" , &s2 ,32000 ,99 ) ;
300  couples->Branch( "s." , "EdbSegP" , &s ,32000 ,99 ) ;
301 
302 
303  // CREATE TEMPORARY ARRAY to store top and Bottom
304 
305  //base = new TClonesArray("EdbSegP", 100) ;
306  if ( segtop.GetEntriesFast() ) segtop.Delete() ;
307  if ( segbot.GetEntriesFast() ) segbot.Delete() ;
308 
309  // create dummy variables to read unused field
310  Char_t dummych ;
311  UShort_t dummyus ;
312  UInt_t dummyui ;
313  Float_t dummyf ;
314  Double_t dummyd ;
315 
316  // Initialize the transformation
317  Header::Transform::MXX = 1 ;
318  Header::Transform::MXY = 0 ;
319  Header::Transform::MYX = 0 ;
320  Header::Transform::MYY = 1 ;
321  Header::Transform::TX = 0 ;
322  Header::Transform::TY = 0 ;
323  Header::Transform::RX = 0 ;
324  Header::Transform::RY = 0 ;
325 
326  // Read File
327  //fstream infile(input_file, ios::in | ios::binary);
328  fstream infile(input_file, ios::in | ios::binary);
329 
330  using namespace Header;
331  infile.read( (Char_t*) &infotype , sizeof(Char_t) );
332  infile.read( (Char_t*) &headerformat , sizeof(UShort_t) );
333 
334  if (infotype == ((Char_t) File::Info::Track | (Char_t)File::Section::Header))
335  switch (headerformat)
336  {
337  case (UShort_t) File::Format::MultiSection: break;
338  case (UShort_t) File::Format::Detailed : break;
339  case (UShort_t) File::Format::NormalDouble_tWithIndex : break;
340  case (UShort_t) File::Format::NormalWithIndex : break;
341  case (UShort_t) File::Format::Normal : break;
342  case (UShort_t) File::Format::NoExtents : break;
343  case (UShort_t) File::Format::Old : break;
344  case (UShort_t) File::Format::Old2 : break;
345 
346  default: { cout << "Unknown format\n"; return ; }
347  };
348 
349  if (headerformat == (UShort_t)File::Format::MultiSection)
350  {
351  infile.read( (Char_t*) &dummych , sizeof(Char_t) );
352  if (dummych != SectionTag)
353  {
354  cout << "The first section in a TLG file must contain tracks!" << endl;
355  return;
356  }
357  infile.read( (Char_t*) &dummyd , sizeof(Double_t) );
358  }
359 
360  int i;
361 
362  Header::Transform::MXX = Header::Transform::MYY = 1.0;
363  Header::Transform::MXY = Header::Transform::MYX = 0.0;
364  Header::Transform::TX = Header::Transform::TY = 0.0;
365  Header::Transform::RX = Header::Transform::RY= 0.0;
366 
367  if (headerformat == (UShort_t) File::Format::Old)
368  {
369  infile.read( (Char_t*) &Identifier::Part3 , sizeof(Int_t) );
370  Identifier::Part0 = Identifier::Part1 = Identifier::Part2 = 0;
371  }
372  else
373  {
374  infile.read( (Char_t*) &Identifier::Part0 , sizeof(Int_t) );
375  infile.read( (Char_t*) &Identifier::Part1 , sizeof(Int_t) );
376  infile.read( (Char_t*) &Identifier::Part2 , sizeof(Int_t) );
377  infile.read( (Char_t*) &Identifier::Part3 , sizeof(Int_t) );
378  }
379 
380  if (headerformat == (UShort_t) File::Format::Detailed || headerformat == (UShort_t)File::Format::MultiSection)
381  {
382  infile.read( (Char_t*) &Center::X , sizeof(Double_t) );
383  infile.read( (Char_t*) &Center::Y , sizeof(Double_t) );
384  infile.read( (Char_t*) &Extents::MinX , sizeof(Double_t) );
385  infile.read( (Char_t*) &Extents::MaxX , sizeof(Double_t) );
386  infile.read( (Char_t*) &Extents::MinX , sizeof(Double_t) );
387  infile.read( (Char_t*) &Extents::MaxY , sizeof(Double_t) );
388  infile.read( (Char_t*) &Transform::MXX , sizeof(Double_t) );
389  infile.read( (Char_t*) &Transform::MXY , sizeof(Double_t) );
390  infile.read( (Char_t*) &Transform::MYX , sizeof(Double_t) );
391  infile.read( (Char_t*) &Transform::MYY , sizeof(Double_t) );
392  infile.read( (Char_t*) &Transform::TX , sizeof(Double_t) );
393  infile.read( (Char_t*) &Transform::TY , sizeof(Double_t) );
394  infile.read( (Char_t*) &Transform::RX , sizeof(Double_t) );
395  infile.read( (Char_t*) &Transform::RY , sizeof(Double_t) );
396 
397 
398  infile.read( (Char_t*) &TopViewsLength , sizeof(Int_t) ); // n. of View
399  infile.read( (Char_t*) &BottomViewsLength , sizeof(Int_t) ); // n. of View
400  infile.read( (Char_t*) &Top_TopZ , sizeof(Double_t) ); // topext
401  infile.read( (Char_t*) &Top_BottomZ , sizeof(Double_t) ); // topInt_t
402  infile.read( (Char_t*) &Bottom_TopZ , sizeof(Double_t) ); // bottomInt_t
403  infile.read( (Char_t*) &Bottom_BottomZ , sizeof(Double_t) ); // bottomext
404 
405  using namespace View ;
406  for (int s=0; s<2; s++ ) {
407  int nviews = ( s==0 ? TopViewsLength : BottomViewsLength );
408  for (int i=0; i<nviews; i++ ) {
409  infile.read( (Char_t*) &id , sizeof(Int_t) );
410  infile.read( (Char_t*) &px , sizeof(Double_t) );
411  infile.read( (Char_t*) &py , sizeof(Double_t) );
412  infile.read( (Char_t*) &topz , sizeof(Double_t) );
413  infile.read( (Char_t*) &bottomz , sizeof(Double_t) );
414  // print_View( s, i ) ;
415  }
416  }
417 
418  infile.read( (Char_t*) &TopTracksLength , sizeof(Int_t) );
419  infile.read( (Char_t*) &BottomTracksLength , sizeof(Int_t) );
420  infile.read( (Char_t*) &TracksLength , sizeof(Int_t) );
421 
422  print_Header();
423 
424  using namespace MIPEmulsionTrack ;
425  for (int s=0; s<2; s++ ) {
426  en_side side = ( s==0? TOP : BOT ) ;
427  int ntracks = ( s==0 ? TopTracksLength : BottomTracksLength );
428  for (int i=0; i<ntracks; i++ ) {
429  // 84 bytes
430  infile.read( (Char_t*) &Field , sizeof(UInt_t) );
431  infile.read( (Char_t*) &AreaSum , sizeof(UInt_t) );
432  infile.read( (Char_t*) &Count , sizeof(UInt_t) );
433  infile.read( (Char_t*) &X , sizeof(Double_t) );
434  infile.read( (Char_t*) &Y , sizeof(Double_t) );
435  infile.read( (Char_t*) &Z , sizeof(Double_t) );
436  infile.read( (Char_t*) &SX , sizeof(Double_t) );
437  infile.read( (Char_t*) &SY , sizeof(Double_t) );
438  infile.read( (Char_t*) &SZ , sizeof(Double_t) );
439  infile.read( (Char_t*) &Sigma , sizeof(Double_t) );
440  infile.read( (Char_t*) &TopZ , sizeof(Double_t) );
441  infile.read( (Char_t*) &BottomZ , sizeof(Double_t) );
442  infile.read( (Char_t*) &viewid , sizeof(Int_t) );
443 
444  // print_MIPEmulsionTrack( s, i ) ;
445  Add_MIPEmulsionTrack ( side, i ) ;
446  }
447  }
448 
449  using namespace MIPBaseTrack ;
450  for (int i=0; i<TracksLength; i++ ) {
451  // .. bytes
452  infile.read( (Char_t*) &MIPBaseTrack::AreaSum , sizeof(UInt_t) );
453  infile.read( (Char_t*) &MIPBaseTrack::Count , sizeof(UInt_t) );
454  infile.read( (Char_t*) &MIPBaseTrack::X , sizeof(Double_t) );
455  infile.read( (Char_t*) &MIPBaseTrack::Y , sizeof(Double_t) );
456  infile.read( (Char_t*) &MIPBaseTrack::Z , sizeof(Double_t) );
457  infile.read( (Char_t*) &MIPBaseTrack::SX , sizeof(Double_t) );
458  infile.read( (Char_t*) &MIPBaseTrack::SY , sizeof(Double_t) );
459  infile.read( (Char_t*) &MIPBaseTrack::SZ , sizeof(Double_t) );
460  infile.read( (Char_t*) &MIPBaseTrack::Sigma , sizeof(Double_t) );
461  infile.read( (Char_t*) &MIPBaseTrack::TopId , sizeof(Int_t) );
462  infile.read( (Char_t*) &MIPBaseTrack::BottomId , sizeof(Int_t) );
463 
464  // print_MIPBaseTrack( i ) ;
465  Add_MIPBaseTrack( i , project_dz) ;
466  }
467 
468  using namespace TrackIndexEntry;
469  for (int s=0; s<2; s++ ) {
470  int ntracks = ( s==0 ? TopTracksLength : BottomTracksLength );
471  for (int i=0; i<ntracks; i++ )
472  {
473  infile.read( (Char_t*) &TrackIndexEntry::Fragment , sizeof(Int_t) );
474  infile.read( (Char_t*) &TrackIndexEntry::View , sizeof(Int_t) );
475  infile.read( (Char_t*) &TrackIndexEntry::Track , sizeof(Int_t) );
476 
477  // printrackIndexEntry( s , i ) ;
478  }
479  }
480  }
481  else if (headerformat == (UShort_t) File::Format::NormalDouble_tWithIndex )
482  {
483  using namespace Header;
484  infile.read( (Char_t*) &Center::X , sizeof(Double_t) );
485  infile.read( (Char_t*) &Center::Y , sizeof(Double_t) );
486  infile.read( (Char_t*) &Extents::MinX , sizeof(Double_t) );
487  infile.read( (Char_t*) &Extents::MaxX , sizeof(Double_t) );
488  infile.read( (Char_t*) &Extents::MinX , sizeof(Double_t) );
489  infile.read( (Char_t*) &Extents::MaxY , sizeof(Double_t) );
490  infile.read( (Char_t*) &dummyf , sizeof(Float_t) ); // for format compliance
491 
492  // 16 bytes
493  infile.read( (Char_t*) &TopTracksLength , sizeof(Int_t) ); // n. of MIPIndexedEmulsionTrack
494  infile.read( (Char_t*) &BottomTracksLength , sizeof(Int_t) ); // n. of MIPIndexedEmulsionTrack
495  infile.read( (Char_t*) &TracksLength , sizeof(Int_t) ); // n. of MIPBaseTrack
496  infile.read( (Char_t*) &dummyui , sizeof(UInt_t) ); // count of fields
497 
498  // 34 bytes
499  infile.read( (Char_t*) &Top_TopZ , sizeof(Double_t) ); // topext
500  infile.read( (Char_t*) &Top_BottomZ , sizeof(Double_t) ); // topInt_t
501  infile.read( (Char_t*) &Bottom_TopZ , sizeof(Double_t) ); // bottomInt_t
502  infile.read( (Char_t*) &Bottom_BottomZ , sizeof(Double_t) ); // bottomext
503  infile.read( (Char_t*) &dummych , sizeof(Char_t) ); // skip fields
504  infile.read( (Char_t*) &dummych , sizeof(Char_t) ); // skip fields
505 
506  print_Header();
507 
508  using namespace MIPEmulsionTrack ;
509  for (int s=0; s<2; s++ ) {
510  en_side side = ( s==0? TOP : BOT ) ;
511  int ntracks = ( s==0 ? TopTracksLength : BottomTracksLength );
512  for (int i=0; i<ntracks; i++ ) {
513  // 84 bytes
514  infile.read( (Char_t*) &Field , sizeof(UInt_t) );
515  infile.read( (Char_t*) &AreaSum , sizeof(UInt_t) );
516  infile.read( (Char_t*) &Count , sizeof(UInt_t) );
517  infile.read( (Char_t*) &X , sizeof(Double_t) );
518  infile.read( (Char_t*) &Y , sizeof(Double_t) );
519  infile.read( (Char_t*) &Z , sizeof(Double_t) );
520  infile.read( (Char_t*) &SX , sizeof(Double_t) );
521  infile.read( (Char_t*) &SY , sizeof(Double_t) );
522  infile.read( (Char_t*) &SZ , sizeof(Double_t) );
523  infile.read( (Char_t*) &Sigma , sizeof(Double_t) );
524  infile.read( (Char_t*) &TopZ , sizeof(Double_t) );
525  infile.read( (Char_t*) &BottomZ , sizeof(Double_t) );
526 
527  //print_MIPEmulsionTrack( s, i ) ;
528  Add_MIPEmulsionTrack ( side , i );
529  }
530  }
531 
532  using namespace MIPBaseTrack ;
533  for (int i=0; i<TracksLength; i++ ) {
534  // .. bytes
535  infile.read( (Char_t*) &MIPBaseTrack::AreaSum , sizeof(UInt_t) );
536  infile.read( (Char_t*) &MIPBaseTrack::Count , sizeof(UInt_t) );
537  infile.read( (Char_t*) &MIPBaseTrack::X , sizeof(Double_t) );
538  infile.read( (Char_t*) &MIPBaseTrack::Y , sizeof(Double_t) );
539  infile.read( (Char_t*) &MIPBaseTrack::Z , sizeof(Double_t) );
540  infile.read( (Char_t*) &MIPBaseTrack::SX , sizeof(Double_t) );
541  infile.read( (Char_t*) &MIPBaseTrack::SY , sizeof(Double_t) );
542  infile.read( (Char_t*) &MIPBaseTrack::SZ , sizeof(Double_t) );
543  infile.read( (Char_t*) &MIPBaseTrack::Sigma , sizeof(Double_t) );
544  if (headerformat == (UShort_t)File::Format::NoExtents) {
545  // Int_tercept errors and slope errors
546  for (int j=0; j<6; j++)
547  infile.read( (Char_t*) &dummyd , sizeof(Double_t) );
548  }
549  infile.read( (Char_t*) &MIPBaseTrack::TopId , sizeof(Int_t) );
550  infile.read( (Char_t*) &MIPBaseTrack::BottomId , sizeof(Int_t) );
551 
552  //print_MIPBaseTrack( i ) ;
553  Add_MIPBaseTrack( i , project_dz ) ;
554  }
555 
556  using namespace TrackIndexEntry ;
557  for (int s=0; s<2; s++ ) {
558  int ntracks = ( s==0 ? TopTracksLength : BottomTracksLength );
559  for (int i=0; i<ntracks; i++ )
560  {
561  infile.read( (Char_t*) &TrackIndexEntry::Fragment , sizeof(Int_t) );
562  infile.read( (Char_t*) &TrackIndexEntry::View , sizeof(Int_t) );
563  infile.read( (Char_t*) &TrackIndexEntry::Track , sizeof(Int_t) );
564 
565  //printrackIndexEntry( s , i ) ;
566  }
567  }
568  }
569  else
570  {
571 
572  /*
573  if (headerformat == (ushort)File.Format.Normal || headerformat == (ushort)File.Format.NormalWithIndex)
574  {
575  m_Extents.MinX = r.ReadSingle();
576  m_Extents.MaxX = r.ReadSingle();
577  m_Extents.MinY = r.ReadSingle();
578  m_Extents.MaxY = r.ReadSingle();
579  r.ReadSingle(); //skip a meaningless Float_t...
580  }
581  else
582  {
583  r.ReadSingle(); r.ReadSingle();
584  r.ReadSingle(); r.ReadSingle(); r.ReadSingle();
585  m_Extents.MinX = m_Extents.MaxX = m_Center.X;
586  m_Extents.MinY = m_Extents.MaxY = m_Center.Y;
587  }
588 
589  View topview = new View();
590  View bottomview = new View();
591  MIPIndexedEmulsionTrack [] toptracks = new MIPIndexedEmulsionTrack[r.ReadUInt_t32()];
592  MIPIndexedEmulsionTrack [] bottomtracks = new MIPIndexedEmulsionTrack[r.ReadUInt_t32()];
593  m_Tracks = new MIPBaseTrack[r.ReadUInt_t32()];
594 
595  i = r.ReadInt_t32(); // count of fields
596 
597  Double_t topext, topInt_t, bottomInt_t, bottomext;
598  topext = r.ReadSingle();
599  topInt_t = r.ReadSingle();
600  bottomInt_t = r.ReadSingle();
601  bottomext = r.ReadSingle();
602 
603  r.ReadBytes(2 * i); // skip fields
604 
605  MIPIndexedEmulsionTrack [] tkarr = toptracks;
606  View vw = topview;
607  do
608  {
609  if (headerformat == (ushort)File.Format.Normal || headerformat == (ushort)File.Format.NormalWithIndex)
610  {
611  for (i = 0; i < tkarr.Length; i++)
612  {
613  MIPEmulsionTrackInfo info = new MIPEmulsionTrackInfo();
614  info.Field = r.ReadUInt_t32();
615  info.AreaSum = r.ReadUInt_t32();
616  info.Count = (ushort)r.ReadUInt_t32();
617  info.Int_tercept.X = r.ReadSingle();
618  info.Int_tercept.Y = r.ReadSingle();
619  info.Int_tercept.Z = r.ReadSingle();
620  info.Slope.X = r.ReadSingle();
621  info.Slope.Y = r.ReadSingle();
622  info.Slope.Z = r.ReadSingle();
623  info.Sigma = r.ReadSingle();
624  info.TopZ = r.ReadSingle();
625  info.BottomZ = r.ReadSingle();
626  tkarr[i] = new MIPIndexedEmulsionTrack(info, i, vw);
627  }
628  }
629  else
630  {
631  if (headerformat == (ushort)File.Format.NoExtents)
632  for (i = 0; i < tkarr.Length; i++)
633  {
634  MIPEmulsionTrackInfo info = new MIPEmulsionTrackInfo();
635  info.AreaSum = 0;
636  info.Field = r.ReadUInt_t32();
637  info.Count = (ushort)r.ReadUInt_t32();
638  info.Int_tercept.X = r.ReadSingle();
639  info.Int_tercept.Y = r.ReadSingle();
640  info.Int_tercept.Z = r.ReadSingle();
641  info.Slope.X = r.ReadSingle();
642  info.Slope.Y = r.ReadSingle();
643  info.Slope.Z = r.ReadSingle();
644  info.Sigma = r.ReadSingle();
645  info.TopZ = r.ReadSingle();
646  info.BottomZ = r.ReadSingle();
647  r.ReadSingle(); r.ReadSingle(); r.ReadSingle(); // Int_tercept errors
648  r.ReadSingle(); r.ReadSingle(); r.ReadSingle(); // slope errors
649  tkarr[i] = new MIPIndexedEmulsionTrack(info, i, vw);
650  }
651  else
652  for (i = 0; i < tkarr.Length; i++)
653  {
654  MIPEmulsionTrackInfo info = new MIPEmulsionTrackInfo();
655  info.AreaSum = 0;
656  info.Field = r.ReadUInt_t32();
657  info.Count = (ushort)r.ReadUInt_t32();
658  r.ReadSingle(); r.ReadSingle(); r.ReadSingle();
659  info.Int_tercept.X = r.ReadSingle();
660  info.Int_tercept.Y = r.ReadSingle();
661  info.Int_tercept.Z = r.ReadSingle();
662  info.Slope.X = r.ReadSingle();
663  info.Slope.Y = r.ReadSingle();
664  info.Slope.Z = r.ReadSingle();
665  info.Sigma = r.ReadSingle();
666  info.TopZ = r.ReadSingle();
667  info.BottomZ = r.ReadSingle();
668  tkarr[i] = new MIPIndexedEmulsionTrack(info, i, vw);
669  }
670  }
671  tkarr = (tkarr == toptracks) ? bottomtracks : null;
672  vw = (vw == topview) ? bottomview : null;
673  }
674  while (tkarr != null);
675 
676  m_Top = new Side(topext, topInt_t);
677  m_Bottom = new Side(bottomInt_t, bottomext);
678  ((Side)m_Top).SetViews(new View[1] { topview });
679  ((Side)m_Bottom).SetViews(new View[1] { bottomview });
680  ((Side)m_Top).SetTracks(toptracks);
681  ((Side)m_Bottom).SetTracks(bottomtracks);
682  topview.m_Side = (Side)m_Top;
683  bottomview.m_Side = (Side)m_Bottom;
684  topview.m_Tracks = toptracks;
685  bottomview.m_Tracks = bottomtracks;
686  topview.m_Id = 0;
687  bottomview.m_Id = 0;
688  topview.m_Position.X = bottomview.m_Position.X = m_Center.X;
689  topview.m_Position.Y = bottomview.m_Position.Y = m_Center.Y;
690  topview.m_TopZ = m_Top.TopZ;
691  topview.m_BottomZ = m_Top.BottomZ;
692  bottomview.m_TopZ = m_Bottom.TopZ;
693  bottomview.m_BottomZ = m_Bottom.BottomZ;
694 
695  if (headerformat == (ushort)File.Format.Normal || headerformat == (ushort)File.Format.NormalWithIndex)
696  {
697  for (i = 0; i < m_Tracks.Length; i++)
698  {
699  MIPEmulsionTrackInfo info = new MIPEmulsionTrackInfo();
700  info.AreaSum = (ushort)r.ReadUInt_t32();
701  info.Count = (ushort)r.ReadUInt_t32();
702  info.Int_tercept.X = r.ReadSingle();
703  info.Int_tercept.Y = r.ReadSingle();
704  info.Int_tercept.Z = r.ReadSingle();
705  info.Slope.X = r.ReadSingle();
706  info.Slope.Y = r.ReadSingle();
707  info.Slope.Z = r.ReadSingle();
708  info.Sigma = r.ReadSingle();
709  if (headerformat == (ushort)File.Format.NoExtents)
710  {
711  r.ReadSingle(); r.ReadSingle(); r.ReadSingle(); // Int_tercept errors
712  r.ReadSingle(); r.ReadSingle(); r.ReadSingle(); // slope errors
713  }
714  MIPIndexedEmulsionTrack toptk = (MIPIndexedEmulsionTrack)m_Top.m_Tracks[r.ReadUInt_t32()];
715  MIPIndexedEmulsionTrack bottomtk = (MIPIndexedEmulsionTrack)m_Bottom.m_Tracks[r.ReadUInt_t32()];
716  info.TopZ = MIPIndexedEmulsionTrack.GetInfo(toptk).TopZ;
717  info.BottomZ = MIPIndexedEmulsionTrack.GetInfo(bottomtk).BottomZ;
718  m_Tracks[i] = new MIPBaseTrack(info, toptk, bottomtk, i);
719  }
720  }
721  else
722  {
723  if (headerformat == (ushort)File.Format.NoExtents)
724  {
725  for (i = 0; i < m_Tracks.Length; i++)
726  {
727  MIPEmulsionTrackInfo info = new MIPEmulsionTrackInfo();
728  info.AreaSum = 0;
729  info.Count = (ushort)r.ReadUInt_t32();
730  info.Int_tercept.X = r.ReadSingle();
731  info.Int_tercept.Y = r.ReadSingle();
732  info.Int_tercept.Z = r.ReadSingle();
733  info.Slope.X = r.ReadSingle();
734  info.Slope.Y = r.ReadSingle();
735  info.Slope.Z = r.ReadSingle();
736  info.Sigma = r.ReadSingle();
737  r.ReadSingle(); r.ReadSingle(); r.ReadSingle(); // Int_tercept errors
738  r.ReadSingle(); r.ReadSingle(); r.ReadSingle(); // slope errors
739  MIPIndexedEmulsionTrack toptk = (MIPIndexedEmulsionTrack)m_Top.m_Tracks[r.ReadUInt_t32()];
740  MIPIndexedEmulsionTrack bottomtk = (MIPIndexedEmulsionTrack)m_Bottom.m_Tracks[r.ReadUInt_t32()];
741  info.TopZ = MIPIndexedEmulsionTrack.GetInfo(toptk).TopZ;
742  info.BottomZ = MIPIndexedEmulsionTrack.GetInfo(bottomtk).BottomZ;
743  m_Tracks[i] = new MIPBaseTrack(info, toptk, bottomtk, i);
744  }
745  }
746  else
747  {
748  for (i = 0; i < m_Tracks.Length; i++)
749  {
750  MIPEmulsionTrackInfo info = new MIPEmulsionTrackInfo();
751  info.Count = (ushort)r.ReadUInt_t32();
752  info.Int_tercept.X = r.ReadSingle();
753  info.Int_tercept.Y = r.ReadSingle();
754  info.Int_tercept.Z = r.ReadSingle();
755  info.Slope.X = r.ReadSingle();
756  info.Slope.Y = r.ReadSingle();
757  info.Slope.Z = r.ReadSingle();
758  info.Sigma = r.ReadSingle();
759  MIPIndexedEmulsionTrack toptk = (MIPIndexedEmulsionTrack)m_Top.m_Tracks[r.ReadUInt_t32()];
760  MIPIndexedEmulsionTrack bottomtk = (MIPIndexedEmulsionTrack)m_Bottom.m_Tracks[r.ReadUInt_t32()];
761  info.AreaSum = MIPIndexedEmulsionTrack.GetInfo(toptk).AreaSum + MIPIndexedEmulsionTrack.GetInfo(bottomtk).AreaSum;
762  info.TopZ = MIPIndexedEmulsionTrack.GetInfo(toptk).TopZ;
763  info.BottomZ = MIPIndexedEmulsionTrack.GetInfo(bottomtk).BottomZ;
764  m_Tracks[i] = new MIPBaseTrack(info, toptk, bottomtk, i);
765  }
766  }
767  }
768  if (headerformat == (ushort)File.Format.NormalWithIndex)
769  {
770  tkarr = (MIPIndexedEmulsionTrack [])m_Top.m_Tracks;
771  while (true)
772  {
773  for (i = 0; i < tkarr.Length; i++)
774  {
775  tkarr[i].m_OriginalRawData.Fragment = r.ReadInt_t32();
776  tkarr[i].m_OriginalRawData.View = r.ReadInt_t32();
777  tkarr[i].m_OriginalRawData.Track = r.ReadInt_t32();
778  }
779  if (tkarr == (MIPIndexedEmulsionTrack [])m_Top.m_Tracks) tkarr = (MIPIndexedEmulsionTrack [])m_Bottom.m_Tracks; else break;
780  }
781  }
782  */
783  }
784 
785  infile.close() ;
786 
787  o->Write();
788  o->Close();
789 
790 }
Z
struct @8 Z
return
return
Definition: energy.C:49
o
AcqOdyssey * o
Definition: hwinit.C:2
dummyf
Float_t dummyf
Definition: merge_Energy_SytematicSources_Electron.C:9
File::Format::NormalWithIndex
@ NormalWithIndex
Definition: tlg2pattern.C:59
File::Format::NormalDouble_tWithIndex
@ NormalDouble_tWithIndex
Definition: tlg2pattern.C:60
File::Format::NoExtents
@ NoExtents
Definition: tlg2pattern.C:57
File::Format::Old
@ Old
Definition: tlg2pattern.C:55
File::Format::Detailed
@ Detailed
Definition: tlg2pattern.C:61
File::Format::Old2
@ Old2
Definition: tlg2pattern.C:56
File::Format::Normal
@ Normal
Definition: tlg2pattern.C:58
File::Format::MultiSection
@ MultiSection
Definition: tlg2pattern.C:62
File::Info::Track
@ Track
Definition: tlg2pattern.C:53
File::Info::Field
@ Field
Definition: tlg2pattern.C:53
File::Section::Header
@ Header
Definition: tlg2pattern.C:54
Header::SectionTag
const Char_t SectionTag
Definition: tlg2pattern.C:67
MIPBaseTrack::id
UInt_t id
Definition: tlg2pattern.C:118
Add_MIPBaseTrack
void Add_MIPBaseTrack(int i, float project_dz=-105.)
Definition: tlg2couples.C:234
print_Header
void print_Header()
Definition: tlg2couples.C:141
Add_MIPEmulsionTrack
void Add_MIPEmulsionTrack(en_side side, int i)
Definition: tlg2couples.C:202
en_side
en_side
Definition: tlg2pattern.C:37

Variable Documentation

◆ couples

TTree* couples

◆ cp

EdbSegCouple* cp = 0

◆ Pid1

Int_t Pid1 =0

◆ Pid2

Int_t Pid2 =0

◆ s

EdbSegP* s = 0

◆ s1

EdbSegP* s1 = 0

◆ s2

EdbSegP* s2 = 0

◆ Xv

Float_t Xv =0

◆ Yv

Float_t Yv =0