X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FMUONCheck.C;h=2c6476a6077e182894b2107d3fbc68c5fbbfa7fa;hb=191aa23aea5b4c8309bb08a03cc985f36da6ceaa;hp=76db4f6a597f81eb4a3571ea48551805ee8bb293;hpb=d8d3b5b80718dff35317a286ea44232f566f42a8;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/MUONCheck.C b/MUON/MUONCheck.C index 76db4f6a597..2c6476a6077 100644 --- a/MUON/MUONCheck.C +++ b/MUON/MUONCheck.C @@ -19,9 +19,8 @@ // Macro for checking aliroot output and associated files contents // Gines Martinez, Subatech June 2003 // - -// ROOT includes #if !defined(__CINT__) || defined(__MAKECINT__) +// ROOT includes #include "TBranch.h" #include "TClonesArray.h" #include "TFile.h" @@ -46,11 +45,12 @@ #include "AliMUONGlobalTrigger.h" #include "AliMUONLocalTrigger.h" #include "AliMUONTrack.h" -#ENDIF +#endif -void MUONkine(char * filename="galice.root") + +void MUONkine(char * filename="galice.root",Int_t event2Check=0) { - TClonesArray * ListOfParticles = new TClonesArray("TParticle",1000); + // TClonesArray * ListOfParticles = new TClonesArray("TParticle",1000); TParticle * particle = new TParticle(); // Creating Run Loader and openning file containing Hits AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONFolder","READ"); @@ -63,22 +63,24 @@ void MUONkine(char * filename="galice.root") nevents = RunLoader->GetNumberOfEvents(); for(ievent=0; ieventGetEvent(ievent); RunLoader->TreeK()->GetBranch("Particles")->SetAddress(&particle); - nparticles = RunLoader->TreeK()->GetEntries(); + nparticles = (Int_t)RunLoader->TreeK()->GetEntries(); printf(">>> Event %d, Number of particles is %d \n",ievent, nparticles); for(iparticle=0; iparticleTreeK()->GetEvent(iparticle); - particle->Print(""); + particle->Print(""); } + if (event2Check!=0) ievent=nevents; } RunLoader->UnloadKinematics(); } -void MUONhits(char * filename="galice.root") +void MUONhits(char * filename="galice.root", Int_t event2Check=0) { // Creating Run Loader and openning file containing Hits AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONFolder","READ"); @@ -94,6 +96,7 @@ void MUONhits(char * filename="galice.root") nevents = RunLoader->GetNumberOfEvents(); for(ievent=0; ievent>> Event %d \n",ievent); // Getting event ievent RunLoader->GetEvent(ievent); @@ -113,6 +116,7 @@ void MUONhits(char * filename="galice.root") for(ihit=0; ihit(muondata.Hits()->At(ihit)); Int_t Nch = mHit->Chamber(); // chamber number + Int_t detele = mHit-> DetElemId(); // Detection element if defined Int_t hittrack = mHit->Track(); Float_t x = mHit->X(); Float_t y = mHit->Y(); @@ -121,17 +125,18 @@ void MUONhits(char * filename="galice.root") Float_t theta = mHit->Theta(); Float_t phi = mHit->Phi(); Float_t momentum = mHit->Momentum(); - printf(">>> Hit %2d Chamber %2d Track %4d x %6.3f y %6.3f z %7.3f elos %g theta %6.3f phi %5.3f momentum %5.3f\n", - ihit, Nch,hittrack,x,y,z,elos,theta,phi, momentum); + printf(">>> Hit %2d Chamber %2d DetEle %4d Track %4d x %6.3f y %6.3f z %7.3f elos %g momentum %5.3f\n", + ihit, Nch, detele, hittrack,x,y,z,elos,momentum); } muondata.ResetHits(); } // end track loop + if (event2Check!=0) ievent=nevents; } // end event loop MUONLoader->UnloadHits(); } -void MUONdigits(char * filename="galice.root") +void MUONdigits(char * filename="galice.root", Int_t event2Check=0) { // Creating Run Loader and openning file containing Hits AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONFolder","READ"); @@ -150,6 +155,7 @@ void MUONdigits(char * filename="galice.root") AliMUONDigit * mDigit; for(ievent=0; ievent>> Event %d \n",ievent); RunLoader->GetEvent(ievent); @@ -157,7 +163,7 @@ void MUONdigits(char * filename="galice.root") Int_t ichamber, nchambers; nchambers = AliMUONConstants::NCh(); ; muondata.SetTreeAddress("D"); - char branchname[30]; + // char branchname[30]; Int_t icathode, ncathodes; ncathodes=2; @@ -187,18 +193,18 @@ void MUONdigits(char * filename="galice.root") Int_t TCharges1 = mDigit->TrackCharge(1); Int_t TCharges2 = mDigit->TrackCharge(2); - printf(">>> Digit %4d cathode %1d hit %4d PadX %3d PadY %3d Signal %4d Physics %4d Track0 %4d TrackCharge0 %4d Track1 %'d - TrackCharge1 %4d Track2 %4d TrackCharge2 %4d \n",idigit, Cathode,Hit, PadX, PadY, Signal, Physics, + printf(">>> Digit %4d cathode %1d hit %4d PadX %3d PadY %3d Signal %4d Physics %4d Track0 %4d TrackCharge0 %4d Track1 %'d TrackCharge1 %4d Track2 %4d TrackCharge2 %4d \n",idigit, Cathode,Hit, PadX, PadY, Signal, Physics, Track0, TCharges0, Track1, TCharges1, Track2, TCharges2); } // end digit loop } // end chamber loop muondata.ResetDigits(); } // end cathode loop + if (event2Check!=0) ievent=nevents; } // end event loop MUONLoader->UnloadDigits(); } -void MUONrecpoints(char * filename="galice.root") { +void MUONrecpoints(char * filename="galice.root", Int_t event2Check=0) { // Creating Run Loader and openning file containing Hits AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONFolder","READ"); @@ -214,9 +220,10 @@ void MUONrecpoints(char * filename="galice.root") { Int_t ievent, nevents; nevents = RunLoader->GetNumberOfEvents(); - AliMUONRawCluster * mRecPoint; + AliMUONRawCluster * mRecPoint = 0; for(ievent=0; ievent>> Event %d \n",ievent); RunLoader->GetEvent(ievent); // Addressing @@ -256,30 +263,31 @@ void MUONrecpoints(char * filename="galice.root") { // // charge chi2 compatibility // // 2 none give satisfactory chi2 - Int_t Track0 = mRecPoint->fTracks[0]; - Int_t Track1 = mRecPoint->fTracks[1]; - Int_t Track2 = mRecPoint->fTracks[2]; - Int_t Q0 = mRecPoint->fQ[0]; - Int_t Q1 = mRecPoint->fQ[1]; - Float_t x0 = mRecPoint->fX[0]; - Float_t x1 = mRecPoint->fX[1]; - Float_t y0 = mRecPoint->fY[0]; - Float_t y1 = mRecPoint->fY[1]; - Float_t z0 = mRecPoint->fZ[0]; - Float_t z1 = mRecPoint->fZ[1]; - Float_t chi2_0 = mRecPoint->fChi2[0]; - Float_t chi2_1 = mRecPoint->fChi2[1]; + Int_t Track0 = mRecPoint->GetTrack(0); + Int_t Track1 = mRecPoint->GetTrack(1); + Int_t Track2 = mRecPoint->GetTrack(2); + Int_t Q0 = mRecPoint->GetCharge(0); + Int_t Q1 = mRecPoint->GetCharge(1); + Float_t x0 = mRecPoint->GetX(0); + Float_t x1 = mRecPoint->GetX(1); + Float_t y0 = mRecPoint->GetY(0); + Float_t y1 = mRecPoint->GetY(1); + Float_t z0 = mRecPoint->GetZ(0); + Float_t z1 = mRecPoint->GetZ(1); + Float_t chi2_0 = mRecPoint->GetChi2(0); + Float_t chi2_1 = mRecPoint->GetChi2(1); printf(">>> RecPoint %4d x %6.3f %6.3f y %6.3f %6.3f z %6.3f %6.3f Q0 %4d Q1 %4d Hit %4d Track1 %4d Track2 %4d Chi2 %6.3f %6.3f \n", irecpoint, x0, x1, y0, y1, z0, z1, Q0, Q1, Track0, Track1, Track2, chi2_0, chi2_1); } // end recpoint loop } // end chamber loop muondata.ResetRawClusters(); + if (event2Check!=0) ievent=nevents; } // end event loop MUONLoader->UnloadRecPoints(); } -void MUONTestTrigger (char * filename="galice.root"){ +void MUONTestTrigger (char * filename="galice.root", Int_t event2Check=0){ // reads and dumps trigger objects from MUON.RecPoints.root TClonesArray * globalTrigger; TClonesArray * localTrigger; @@ -304,6 +312,7 @@ void MUONTestTrigger (char * filename="galice.root"){ AliMUONLocalTrigger *locTrg; for (ievent=0; ieventGetEvent(ievent); muondata.SetTreeAddress("GLT"); @@ -365,13 +374,14 @@ void MUONTestTrigger (char * filename="galice.root"){ } // end of loop on Local Trigger muondata.ResetTrigger(); + if (event2Check!=0) ievent=nevents; } // end loop on event MUONLoader->UnloadRecPoints(); } -void MUONRecTracks (char * filename="galice.root"){ +void MUONRecTracks (char * filename="galice.root", Int_t event2Check=0 ){ // reads and dumps trigger objects from MUON.RecPoints.root TClonesArray * RecTracks; @@ -390,9 +400,10 @@ void MUONRecTracks (char * filename="galice.root"){ Int_t ievent, nevents; nevents = RunLoader->GetNumberOfEvents(); - AliMUONTrack * rectrack; + // AliMUONTrack * rectrack; for (ievent=0; ieventGetEvent(ievent); muondata.SetTreeAddress("RT"); @@ -405,6 +416,7 @@ void MUONRecTracks (char * filename="galice.root"){ printf(">>> Event %d Number of Recconstructed tracks %d \n",ievent, nrectracks); muondata.ResetRecTracks(); + if (event2Check!=0) ievent=nevents; } // end loop on event MUONLoader->UnloadTracks(); }