X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FAliMUONCheck.cxx;h=af096200aa020cc2c1b2786a01181c2b07a0b96b;hb=e6f3d431a97b8c7454ebfe50c664deb9beffb020;hp=d83fbf986de552bc072c06c4e6cfa67fbc942c7c;hpb=a55f49a01992f9b256e43a101357e2575981a0b3;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/AliMUONCheck.cxx b/MUON/AliMUONCheck.cxx index d83fbf986de..af096200aa0 100644 --- a/MUON/AliMUONCheck.cxx +++ b/MUON/AliMUONCheck.cxx @@ -32,9 +32,6 @@ #include "AliMUONCheck.h" #include "AliMUONConstants.h" -#include "AliMUONTrack.h" -#include "AliMUONTrackParam.h" -#include "AliMUONTrackExtrap.h" #include "AliMUONMCDataInterface.h" #include "AliMUONDataInterface.h" #include "AliMpCDB.h" @@ -52,7 +49,7 @@ #include "AliESDEvent.h" #include "AliESDMuonTrack.h" #include "AliESDVertex.h" -#include "AliMagFMaps.h" +#include "AliMagF.h" #include "AliLog.h" #include @@ -66,21 +63,26 @@ ClassImp(AliMUONCheck) /// \endcond +//_____________________________________________________________________________ +const TString& AliMUONCheck::GetDefaultOutFileName() +{ + /// Default output file name + static const TString kDefaultOutFileName = "output.txt"; + return kDefaultOutFileName; +} + +//_____________________________________________________________________________ AliMUONCheck::AliMUONCheck(const char* galiceFile, const char* esdFile,Int_t firstEvent, Int_t lastEvent,const char* outDir) : TObject(), fFileName(galiceFile), fFileNameSim(), fesdFileName(esdFile), -foutDir(outDir), +fkOutDir(outDir), +fOutFileName(GetDefaultOutFileName()), fFirstEvent(firstEvent), fLastEvent(lastEvent) { - /// ctor - - // Load mapping - if ( ! AliMpCDB::LoadMpSegmentation() ) { - AliFatal("Could not access mapping from OCDB !"); - } + /// ctor } //_____________________________________________________________________________ @@ -91,16 +93,12 @@ AliMUONCheck::AliMUONCheck(const char* galiceFile, const char* galiceFileSim, fFileName(galiceFile), fFileNameSim(galiceFileSim), fesdFileName(esdFile), -foutDir(outDir), +fkOutDir(outDir), +fOutFileName(GetDefaultOutFileName()), fFirstEvent(firstEvent), fLastEvent(lastEvent) { /// ctor - - // Load mapping - if ( ! AliMpCDB::LoadMpSegmentation() ) { - AliFatal("Could not access mapping from OCDB !"); - } } //_____________________________________________________________________________ @@ -169,7 +167,7 @@ AliMUONCheck::CheckESD(Bool_t pdc06TriggerResponse) if (!tree) { Error("CheckESD", "no ESD tree found"); - AliError(Form("CheckESD", "no ESD tree found")); + AliError(Form("no ESD tree found")); return ; } fESD->ReadFromTree(tree); @@ -184,8 +182,7 @@ AliMUONCheck::CheckESD(Bool_t pdc06TriggerResponse) for (ievent = fFirstEvent; ievent < endOfLoop; ++ievent ) { nev++; - tree->GetEvent(ievent); - if (!fESD) + if (tree->GetEvent(ievent) <= 0) { Error("CheckESD", "no ESD object found for event %d", ievent); return ; @@ -277,6 +274,10 @@ AliMUONCheck::CheckESD(Bool_t pdc06TriggerResponse) for ( Int_t iTrack1 = 0; iTrack1GetMuonTrack(iTrack1); + + // skip fake tracks (ghosts) + if (!muonTrack->ContainTrackerData()) continue; + ftracktot++; thetaX = muonTrack->GetThetaX(); @@ -344,10 +345,9 @@ AliMUONCheck::CheckESD(Bool_t pdc06TriggerResponse) }//if(pdc06TriggerResponse) - gSystem->cd(foutDir); + gSystem->cd(fkOutDir); - FILE *outtxt=fopen("output.txt","a"); - freopen("output.txt","a",outtxt); + FILE *outtxt=fopen(fOutFileName.Data(),"a"); if(pdc06TriggerResponse){ fprintf(outtxt," \n"); @@ -460,7 +460,7 @@ AliMUONCheck::CheckKine() p->Print(""); Int_t pdg=p->GetPdgCode(); - if (abs(pdg) == 13) + if (TMath::Abs(pdg) == 13) { ++nmu2; } @@ -527,9 +527,8 @@ AliMUONCheck::CheckKine() printf("*** Leaving MuonKine() *** \n"); printf("**************************************************************** \n"); - gSystem->cd(foutDir); - FILE *outtxt=fopen("output.txt","a"); - freopen("output.txt","a",outtxt); + gSystem->cd(fkOutDir); + FILE *outtxt=fopen(fOutFileName.Data(),"a"); fprintf(outtxt," \n"); fprintf(outtxt,"=================================================================\n"); fprintf(outtxt,"================ MUONkine SUMMARY ================\n"); @@ -672,7 +671,7 @@ AliMUONCheck::CheckTrackRef() }//evt loop - gSystem->cd(foutDir); + gSystem->cd(fkOutDir); TCanvas *c6 = new TCanvas("c6","TOF",400,10,600,700); c6->Divide(1,2); c6->cd(1); @@ -769,7 +768,7 @@ AliMUONCheck::CheckOccupancy(Bool_t perDetEle) const ndigits = digitStore->GetSize(detele,1-cathode); dEchannelsNonBending[ichamber][idetele] = nchannels; - dEoccupancyBending[ichamber][idetele] = ndigits; + dEoccupancyNonBending[ichamber][idetele] = ndigits; cHchannelsNonBending[ichamber] += nchannels; cHoccupancyNonBending[ichamber] += ndigits; totalchannelsNonBending += nchannels; @@ -789,16 +788,18 @@ AliMUONCheck::CheckOccupancy(Bool_t perDetEle) const // Output values - for ( Int_t ichamber = 0; ichamber < 14; ++ichamber ) + for ( Int_t ichamber = 0; ichamber < nchambers; ++ichamber ) { printf(">>> Chamber %2d nChannels Bending %5d nChannels NonBending %5d \n", ichamber+1, cHoccupancyBending[ichamber], - cHoccupancyNonBending[ichamber]); - printf(">>> Chamber %2d Occupancy Bending %5.2f %% Occupancy NonBending %5.2f %% \n", - ichamber+1, - 100.*((Float_t) cHoccupancyBending[ichamber])/((Float_t) cHchannelsBending[ichamber]), - 100.*((Float_t) cHoccupancyNonBending[ichamber])/((Float_t) cHchannelsBending[ichamber]) ); + cHoccupancyNonBending[ichamber]); + if ( cHchannelsBending[ichamber] != 0 && cHchannelsBending[ichamber] ) { + printf(">>> Chamber %2d Occupancy Bending %5.2f %% Occupancy NonBending %5.2f %% \n", + ichamber+1, + 100.*((Float_t) cHoccupancyBending[ichamber])/((Float_t) cHchannelsBending[ichamber]), + 100.*((Float_t) cHoccupancyNonBending[ichamber])/((Float_t) cHchannelsBending[ichamber])); + } if ( perDetEle ) { @@ -811,65 +812,23 @@ AliMUONCheck::CheckOccupancy(Bool_t perDetEle) const idetele+100*(ichamber+1), dEoccupancyBending[ichamber][idetele], dEoccupancyNonBending[ichamber][idetele]); - printf(">>> DetEle %4d Occupancy Bending %5.2f %% Occupancy NonBending %5.2f %% \n", - idetele+100*(ichamber+1), - 100.*((Float_t) dEoccupancyBending[ichamber][idetele])/((Float_t) dEchannelsBending[ichamber][idetele]), - 100.*((Float_t) dEoccupancyNonBending[ichamber][idetele])/((Float_t) dEchannelsBending[ichamber][idetele])); + + if ( dEchannelsBending[ichamber][idetele] != 0 && dEchannelsBending[ichamber][idetele] !=0 ) { + printf(">>> DetEle %4d Occupancy Bending %5.2f %% Occupancy NonBending %5.2f %% \n", + idetele+100*(ichamber+1), + 100.*((Float_t) dEoccupancyBending[ichamber][idetele])/((Float_t) dEchannelsBending[ichamber][idetele]), + 100.*((Float_t) dEoccupancyNonBending[ichamber][idetele])/((Float_t) dEchannelsBending[ichamber][idetele])); + } } } } } - printf(">>> Muon Spectrometer Occupancy Bending %5.2f %% Occupancy NonBending %5.2f %% \n", - 100.*((Float_t) totaloccupancyBending)/((Float_t) totalchannelsBending), - 100.*((Float_t) totaloccupancyNonBending)/((Float_t) totalchannelsNonBending) ); - -} - -//_____________________________________________________________________________ -void -AliMUONCheck::CheckRecTracks () const -{ - /// Reads and dumps rec tracks objects - - AliWarning("Reimplement me ? or use AliMUONDumper simply ?"); - - // // waiting for mag field in CDB - // AliInfoStream() << "Loading field map...\n"; - // if (!AliTracker::GetFieldMap()) { - // AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG); - // AliTracker::SetFieldMap(field, kFALSE); - // } - // - // // Loading data - // fLoader->LoadTracks("READ"); - // - // Int_t endOfLoop = fLastEvent+1; - // if ( fLastEvent == -1 ) endOfLoop = fRunLoader->GetNumberOfEvents(); - // - // for (Int_t ievent=fFirstEvent; ieventGetEvent(ievent); - // - // fRecData->SetTreeAddress("RT"); - // fRecData->GetRecTracks(); - // TClonesArray* recTracks = fRecData->RecTracks(); - // - // Int_t nrectracks = (Int_t) recTracks->GetEntriesFast(); // - // printf(">>> Event %d, Number of Recconstructed tracks %d \n",ievent, nrectracks); - // - // // Set the magnetic field for track extrapolations - // AliMUONTrackExtrap::SetField(AliTracker::GetFieldMap()); - // - // // Loop over tracks - // for (Int_t iRecTracks = 0; iRecTracks < nrectracks; iRecTracks++) { - // AliMUONTrack* recTrack = (AliMUONTrack*) recTracks->At(iRecTracks); - // AliMUONTrackParam* trackParam = (AliMUONTrackParam*) (recTrack->GetTrackParamAtHit())->First(); - // AliMUONTrackExtrap::ExtrapToZ(trackParam,0.); - // recTrack->Print("full"); - // } - // fRecData->ResetRecTracks(); - // } - // fLoader->UnloadTracks(); + if ( totalchannelsBending != 0 && totalchannelsNonBending != 0 ) { + printf(">>> Muon Spectrometer Occupancy Bending %5.2f %% Occupancy NonBending %5.2f %% \n", + 100.*((Float_t) totaloccupancyBending)/((Float_t) totalchannelsBending), + 100.*((Float_t) totaloccupancyNonBending)/((Float_t) totalchannelsNonBending)); + } } //_____________________________________________________________________________