]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MUON/MUONCheck.C
From Laurent:
[u/mrichter/AliRoot.git] / MUON / MUONCheck.C
index 81dceaa346564783760c92bb26680aa4b6fd2653..f22401b12ba3e724a71ce155735e144b4a86a980 100644 (file)
@@ -168,9 +168,8 @@ void MUONdigits(char * filename="galice.root", Int_t event2Check=0)
     Int_t icathode, ncathodes;
     ncathodes=2;
     //Loop on cathodes 
-    for(icathode=0; icathode<ncathodes; icathode++) {
-      printf(">>> Cathode %d\n",icathode);
-      muondata.GetCathode(icathode);
+    //    for(icathode=0; icathode<ncathodes; icathode++) {
+      muondata.GetDigits();
       // Loop on chambers
       for( ichamber=0; ichamber<nchambers; ichamber++) {
        printf(">>> Chamber %d\n",ichamber+1);
@@ -193,14 +192,15 @@ void MUONdigits(char * filename="galice.root", Int_t event2Check=0)
          Int_t TCharges1 = mDigit->TrackCharge(1);
          Int_t TCharges2 = mDigit->TrackCharge(2);
          Int_t idDE = mDigit->DetElemId();
-                 
+         //      printf(">>> Cathode %d\n",Cathode);
+         
          printf(">>>IdDE %d Digit %4d cathode %1d hit %4d PadX %3d PadY %3d Signal %4d Physics %4d Track0 %4d TrackCharge0 %4d Track1 %4d TrackCharge1 %4d Track2 %4d TrackCharge2 %4d \n",
                 idDE, idigit, Cathode,Hit, PadX, PadY, Signal, Physics, Track0, 
                 TCharges0, Track1, TCharges1, Track2, TCharges2);
        } // end digit loop
       } // end chamber loop
       muondata.ResetDigits();
-    } // end cathode loop
+      //    } // end cathode loop
     if (event2Check!=0) ievent=nevents;
   }  // end event loop
   MUONLoader->UnloadDigits();
@@ -239,10 +239,9 @@ void MUONrecpoints(char * filename="galice.root", Int_t event2Check=0) {
       printf(">>> Chamber %d\n",ichamber);
       sprintf(branchname,"MUONRawClusters%d",ichamber+1);
       //printf(">>>  branchname %s\n",branchname);
-  
       Int_t irecpoint, nrecpoints;
       nrecpoints = (Int_t) muondata.RawClusters(ichamber)->GetEntriesFast();
-      
+      printf("number of recpoints = %6d \n",nrecpoints);
       for(irecpoint=0; irecpoint<nrecpoints; irecpoint++) {
        mRecPoint = static_cast<AliMUONRawCluster*>(muondata.RawClusters(ichamber)->At(irecpoint));
 //     Int_t       fTracks[3];        //labels of overlapped tracks
@@ -278,9 +277,9 @@ void MUONrecpoints(char * filename="galice.root", Int_t event2Check=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);
+       Int_t de = mRecPoint->GetDetElemId();
+       printf(">>> RecPoint %4d  DetElem %2d 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,de,x0,x1,y0,y1,z0,z1,Q0,Q1,Track0, Track1, Track2, chi2_0, chi2_1);
       } // end recpoint loop
     } // end chamber loop
     muondata.ResetRawClusters();