Changeset 1852 for trunk/MagicSoft
- Timestamp:
- 03/21/03 11:10:36 (22 years ago)
- Location:
- trunk/MagicSoft/Mars
- Files:
-
- 1 added
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/Changelog
r1851 r1852 1 1 -*-*- END OF LINE -*-*- 2 3 2003/03/21: Thomas Bretz 4 * manalysis/MEnergyEstParam.[h,cc] 5 - Added StopMapping and Print functions. 6 7 2003/03/21: Abelardo Moralejo 8 9 * mhist/MHMatrix.[h,cc]: 10 - Added third argument (a filter) to the second instantiation 11 of the Fill procedure. 12 13 * macros/CT1EnergyEst.C: 14 - Example of the parameter calculation and use of the energy 15 estimation method for CT1. 16 2 17 3 18 2003/03/21: Thomas Bretz -
trunk/MagicSoft/Mars/manalysis/MEnergyEstParam.cc
r1844 r1852 177 177 // -------------------------------------------------------------------------- 178 178 // 179 // Set the f ourcoefficients for the estimation of the impact parameter.179 // Set the five coefficients for the estimation of the impact parameter. 180 180 // Argument must ba a TArrayD of size 5. 181 181 // … … 193 193 // -------------------------------------------------------------------------- 194 194 // 195 // Set the fourcoefficients for the estimation of the energy.195 // Set the seven coefficients for the estimation of the energy. 196 196 // Argument must ba a TArrayD of size 7. 197 197 // … … 209 209 // -------------------------------------------------------------------------- 210 210 // 211 // Set the fourcoefficients for the estimation of impact and energy.211 // Set the twelve coefficients for the estimation of impact and energy. 212 212 // Argument must ba a TArrayD of size 12. 213 213 // … … 261 261 } 262 262 263 void MEnergyEstParam::StopMapping() 264 { 265 fMatrix = NULL; 266 fPairs->Clear(); 267 fHillasSrc->Clear(); 268 fEnergy->Clear(); 269 } 270 263 271 // -------------------------------------------------------------------------- 264 272 // … … 272 280 273 281 AddToBranchList(hillas+".fDist"); 282 } 283 284 void MEnergyEstParam::Print(Option_t *opt) 285 { 286 for (int i=0; i<fA.GetSize(); i++) 287 *fLog << all << "fA[" << i << "]=" << fA[i] << endl; 288 for (int i=0; i<fB.GetSize(); i++) 289 *fLog << all << "fB[" << i << "]=" << fB[i] << endl; 274 290 } 275 291 … … 295 311 296 312 /* MY PARAM */ 297 const Double_t e1 = fB[0] + fB[1]*size + fB[3]*length + fB[4]*size*width;298 const Double_t e2 = fB[2] + fB[5]*size*width;313 //const Double_t e1 = fB[0] + fB[1]*size + fB[3]*length + fB[4]*size*width; 314 //const Double_t e2 = fB[2] + fB[5]*size*width; 299 315 300 316 /* MARCOS */ 301 //const Double_t e1 = fB[0] + fB[1]*size + fB[3]*length + fB[4]*size*length;302 //const Double_t e2 = fB[2] + fB[5]*length;317 const Double_t e1 = fB[0] + fB[1]*size + fB[3]*length + fB[4]*size*length; 318 const Double_t e2 = fB[2] + fB[5]*length; 303 319 304 320 TIter NextH(fHillasSrc); … … 317 333 318 334 /* MARCOS */ 319 //const Double_t ir = i0 * (i1 + fA[1]*dist); // [cm]320 /*const*/ // Double_t er = e0 * (e1 + e2*ir); // [GeV]335 const Double_t ir = i0 * (i1 + fA[1]*dist); // [cm] 336 /*const*/ Double_t er = e0 * (e1 + e2*0/*ir*/); // [GeV] 321 337 322 338 /* MY PARAM */ 323 339 // if (width==0) return kCONTINUE; 324 const Double_t ir = i0 * (i1 + dist*(fA[1]/width + fA[4]/log10(size))); // [cm]325 Double_t er = e0 * (e1 + e2*ir); // [GeV]340 //const Double_t ir = i0 * (i1 + dist*(fA[1]/width + fA[4]/log10(size))); // [cm] 341 //Double_t er = e0 * (e1 + e2*ir); // [GeV] 326 342 327 343 /* MKA */ … … 339 355 } 340 356 341 357 est->SetEnergy(er); 342 358 est->SetImpact(ir); 343 359 est->SetReadyToSave(); -
trunk/MagicSoft/Mars/manalysis/MEnergyEstParam.h
r1847 r1852 50 50 51 51 void InitMapping(MHMatrix *mat); 52 void StopMapping(); 52 53 53 54 Int_t GetNumCoeff() const { return fA.GetSize()+fB.GetSize(); } … … 57 58 void SetCoeffB(const TArrayD &arr); 58 59 60 void Print(Option_t *o=NULL); 61 59 62 ClassDef(MEnergyEstParam, 0) // Task to estimate the energy 60 63 }; -
trunk/MagicSoft/Mars/mhist/MHMatrix.cc
r1831 r1852 68 68 #include "MData.h" 69 69 #include "MDataArray.h" 70 #include "MF.h" 71 70 72 71 73 ClassImp(MHMatrix); … … 600 602 // -------------------------------------------------------------------------- 601 603 // 602 Bool_t MHMatrix::Fill(MParList *plist, MTask *read )604 Bool_t MHMatrix::Fill(MParList *plist, MTask *read, MF *filter) 603 605 { 604 606 // … … 614 616 615 617 tlist.AddToList(read); 618 619 620 if (filter) 621 { 622 tlist.AddToList(filter); 623 fillh.SetFilter(filter); 624 } 625 616 626 tlist.AddToList(&fillh); 627 617 628 618 629 MEvtLoop evtloop; -
trunk/MagicSoft/Mars/mhist/MHMatrix.h
r1829 r1852 22 22 class MParList; 23 23 class MDataArray; 24 class MF; 24 25 25 26 class MHMatrix : public MH … … 92 93 Double_t operator[](Int_t col) { return fM(fRow, col); } 93 94 94 Bool_t Fill(MParList *plist, MTask *read );95 Bool_t Fill(MParList *plist, MTask *read, MF *filter=0); 95 96 96 97 TString GetDataMember() const; -
trunk/MagicSoft/Mars/mhist/Makefile
r1821 r1852 23 23 # 24 24 INCLUDES = -I. -I../mbase -I../mraw -I../manalysis -I../mmc \ 25 -I../mgui -I../mgeom -I../mdata 25 -I../mgui -I../mgeom -I../mdata -I../mfilter 26 26 27 27 #------------------------------------------------------------------------------
Note:
See TracChangeset
for help on using the changeset viewer.