Changeset 5879 for trunk/MagicSoft/Mars/mtemp/mifae/programs
- Timestamp:
- 01/18/05 10:53:10 (20 years ago)
- Location:
- trunk/MagicSoft/Mars/mtemp/mifae/programs
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mtemp/mifae/programs/Makefile
r5775 r5879 18 18 include ../../../Makefile.conf.general 19 19 20 PROGRAMS = makeHillas 20 PROGRAMS = makeHillas controlPlot 21 21 #psffit falseSource srcPos optimizeCuts controlPlot 22 22 SOLIB = $(MARSSYS)/libmars.so -
trunk/MagicSoft/Mars/mtemp/mifae/programs/controlPlot.cc
r4395 r5879 64 64 //----------------------------------------------------------------------------- 65 65 66 const TString defaultcard=" srcpos.datacard";66 const TString defaultcard="controlplot.datacard"; 67 67 const Float_t conver = 189./0.6; // conversion factor degrees to mm 68 68 … … 112 112 // include containers in parameter list 113 113 plist.AddToList(&tlist); 114 plist.AddToList(&islands);115 114 116 115 MF cut(filter); … … 127 126 MControlPlots controlplots(outputFile); 128 127 controlplots.SetProduceFile(kFALSE); 128 controlplots.SetMode(MControlPlots::kOn); 129 129 130 130 MContinue applycutS(&cutS); … … 157 157 if(widthcut.Length()) 158 158 tlist.AddToList(&applycutW); 159 if(filter.Length()) 160 tlist.AddToList(&applycut); 159 161 if(centercut.Length()) 160 162 tlist.AddToList(&applycutC); 161 if(filter.Length())162 tlist.AddToList(&applycut);163 163 tlist.AddToList(&controlplots); 164 164 … … 175 175 return; 176 176 177 178 // containers 179 MParList plist2; 180 MTaskList tlist2; 181 MIslands islands2; 182 183 // include containers in parameter list 184 plist2.AddToList(&tlist2); 185 186 // tasks 177 187 MReadTree read2("Parameters", offFile); 178 188 read2.DisableAutoScheme(); 179 tlist.AddToListBefore(&read2, &read, "All"); 180 tlist.RemoveFromList(&read); 189 tlist2.AddToList(&read2); 181 190 182 191 controlplots.SetMode(MControlPlots::kOff); 183 192 controlplots.SetProduceFile(kTRUE); 184 185 if (!evtloop.Eventloop(nmaxevents)) 193 controlplots.Reset(); 194 195 // include tasks in task list 196 if(sizecut.Length()) 197 tlist2.AddToList(&applycutS); 198 if(distcut.Length()) 199 tlist2.AddToList(&applycutD); 200 if(lengthcut.Length()) 201 tlist2.AddToList(&applycutL); 202 if(widthcut.Length()) 203 tlist2.AddToList(&applycutW); 204 if(filter.Length()) 205 tlist2.AddToList(&applycut); 206 if(centercut.Length()) 207 tlist2.AddToList(&applycutC); 208 tlist2.AddToList(&controlplots); 209 210 // Eventloop 211 MEvtLoop evtloop2; 212 evtloop2.SetParList(&plist2); 213 if (!evtloop2.Eventloop(nmaxevents)) 186 214 return; 187 188 tlist.PrintStatistics(); 215 216 tlist2.PrintStatistics(); 217 218 // do off-data if input file was specified 219 if(!offFile.Length()) 220 return; 221 189 222 } 223 190 224 //----------------------------------------------------------------------- 191 225 … … 268 302 ifun >> ch; 269 303 distcut+=ch; 304 distcut+="*"; 305 distcut+=conver; 270 306 distcut+=") && (MHillasSrc.fDist<"; 271 307 ifun >> ch; 272 308 distcut+=ch; 309 distcut+="*"; 310 distcut+=conver; 273 311 distcut+=")"; 274 312 } … … 281 319 ifun >> ch; 282 320 widthcut+=ch; 321 widthcut+="*"; 322 widthcut+=conver; 283 323 widthcut+=") && (MHillas.fWidth<"; 284 324 ifun >> ch; 285 325 widthcut+=ch; 326 widthcut+="*"; 327 widthcut+=conver; 286 328 widthcut+=")"; 287 329 } … … 294 336 ifun >> ch; 295 337 lengthcut+=ch; 338 lengthcut+="*"; 339 lengthcut+=conver; 296 340 lengthcut+=") && (MHillas.fLength<"; 297 341 ifun >> ch; 298 342 lengthcut+=ch; 343 lengthcut+="*"; 344 lengthcut+=conver; 299 345 lengthcut+=")"; 300 346 } … … 304 350 cout << "readDataCards Warning: overriding existing center cut" << endl; 305 351 char ch[10]; 306 centercut="sqrt(MHillas.fMeanX*MHillas.fMeanX+MHillas.fMeanY*MHillas.fMeanY) < 0.00317*";352 centercut="sqrt(MHillas.fMeanX*MHillas.fMeanX+MHillas.fMeanY*MHillas.fMeanY)<"; 307 353 ifun >> ch; 308 354 centercut+=ch; 355 centercut+="*"; 356 centercut+=conver; 309 357 } 310 358 } -
trunk/MagicSoft/Mars/mtemp/mifae/programs/controlplot.datacard
r4395 r5879 1 1 2 2 // Maximun number of on and off events to be processed) 3 NEVENTS 9999 99993 NEVENTS 9999 4 4 5 5 // On-data file name pattern 6 ONFILES ./srcPosPrueba.root6 ONFILES /mnt/users/jrico/magic/hillasCrab_200411_new.root 7 7 8 8 // Off-data file name pattern 9 OFFFILES ./srcPosOffPrueba.root9 OFFFILES /mnt/users/jrico/magic/hillasCrab_200411_new.root 10 10 11 11 // output file name 12 HITSFILE ./ rotateprueba2.ps12 HITSFILE ./prueba.ps 13 13 14 14 ////////// … … 16 16 ////////// 17 17 18 // Size cut (lower and upper) in # of photons19 SIZECUT 2000 999999918 // Size cut (lower and upper) 19 SIZECUT 0 9999999 20 20 21 21 // Dist cut (lower and upper) in degrees 22 DISTCUT 0.2 1.122 DISTCUT 0.2 1.2 23 23 24 24 // Width cut (lower and upper) in degrees 25 WIDTHCUT 0 0.1 225 WIDTHCUT 0 0.15 26 26 27 27 // Length cut (lower and upper) in degrees 28 LENGTHCUT 0 0.2 628 LENGTHCUT 0 0.25 29 29 30 30 // maxX and maxY upper cut in degrees 31 CENTERCUT 1. 131 CENTERCUT 1.2 32 32 33 33 // Other Additional cut
Note:
See TracChangeset
for help on using the changeset viewer.