Changeset 5827
- Timestamp:
- 01/13/05 19:28:53 (20 years ago)
- Location:
- trunk/MagicSoft/Mars/mtemp/mifae/programs
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mtemp/mifae/programs/makeHillas.cc
r5826 r5827 61 61 #include "MContinue.h" 62 62 #include "MReportDrive.h" 63 #include "MTopology.h" 64 #include "MTopologyCalc.h" 63 65 64 66 #include "TApplication.h" … … 93 95 Float_t ltail = 1.5; 94 96 Int_t islflag = 0; 97 Int_t topflag = 0; 95 98 Float_t lnew = 40; 96 99 Int_t kmethod = 1; … … 207 210 pedloop2.SetNoStorage(); 208 211 pedloop2.SetExtractor(extractor); 209 pedloop2.SetExtractorResolution();210 212 pedloop2.SetPedestals(pedloop1.GetPedestalCam()); 211 213 pedloop2.SetInput(&pedcaliter); … … 374 376 MIslands isl2; 375 377 MIslands isl3; 378 MTopology topology; 376 379 377 380 isl.SetName("MIslands"); … … 401 404 plist5.AddToList(&runhead); 402 405 plist5.AddToList(&reportdrive); 406 plist5.AddToList(&topology); 403 407 404 408 // cuts … … 435 439 MHillasCalc hcalc; 436 440 MHillasSrcCalc csrc1; 441 442 MTopologyCalc topcalc; 437 443 438 444 MContinue applycut(&cut); … … 453 459 write.AddContainer("MIslands" , "Parameters"); 454 460 write.AddContainer("MReportDrive" , "Parameters"); 461 write.AddContainer("MTopology", "Parameters"); 455 462 456 463 if (islflag == 2) … … 494 501 if(filter.Length()) 495 502 tlist5.AddToList(&applycut); 503 if(topflag>0) 504 tlist5.AddToList(&topcalc); 496 505 tlist5.AddToList(&write); 497 506 if(display) … … 503 512 tlist5.AddToList(disphillas); 504 513 } 514 505 515 506 516 // Create and setup the eventloop … … 635 645 ifun >> sext >> hifirst >> hilast >> lofirst >> lolast >> wsize; 636 646 647 if(strcmp(word.Data(),"TOPFLAG")==0) 648 ifun >> topflag; 649 637 650 if(strcmp(word.Data(),"ISLFLAG")==0) 638 651 { … … 674 687 cout << "Maximum number of events: " << nmaxevents << endl; 675 688 if(filter.Length()) 676 cout << "Applying rejection cut: " << filter << endl;689 cout << "Applying selection cut: " << filter << endl; 677 690 cout << "Output file name: " << outname << endl; 678 691 if(display) … … 692 705 if (islflag == 1 || islflag == 2) 693 706 cout << "Island calcultation..." << "using algorithm #" << kalgorithm <<endl; 707 if (topflag > 0) 708 cout << "Topology parameters will be computed" << endl; 694 709 if (islflag == 2) 695 710 { -
trunk/MagicSoft/Mars/mtemp/mifae/programs/makehillas.datacard
r5775 r5827 11 11 12 12 // Maximun number of (data) events to be processed) 13 NEVENTS 500013 NEVENTS 999999 14 14 15 15 // data file directory … … 30 30 // Enclose all conditions between brakets, like: "(x<y) && (z<5)" 31 31 // (see MF description class for more details) 32 // FILTER (MHillas.fLength<100) && (MHillas.fLength>50)32 FILTER (MNewImagePar.fNumUsedPixels<400) 33 33 34 34 // Display flag … … 61 61 62 62 // Cleaning level (tail cut, boundary cut, number of rings, cleaning method - 1=standard, 2=democratic) 63 CLEANLEVEL 2.5 2.01 163 CLEANLEVEL 4.0 3.5 1 1 64 64 65 65 //Island calculations … … 74 74 ISLANDCLEAN 1 40 75 75 76 // topology flag 77 // 0: do not calculate topology parameters 78 // 1: calculate topology parameters 79 TOPFLAG 1 80
Note:
See TracChangeset
for help on using the changeset viewer.