X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=PHOS%2FAliPHOSvImpacts.cxx;h=56ffb673c4e0c3749b6fcfc2e61ae74ccd5c0b99;hp=3cae449e928328487768036804534d90e82ce10c;hb=6350d0c03d9aed9b63ea3a0ef57f3aa3252efced;hpb=88cb7938ca21d4a80991d4e7aa564008c29340f7 diff --git a/PHOS/AliPHOSvImpacts.cxx b/PHOS/AliPHOSvImpacts.cxx index 3cae449e928..56ffb673c4e 100644 --- a/PHOS/AliPHOSvImpacts.cxx +++ b/PHOS/AliPHOSvImpacts.cxx @@ -16,6 +16,23 @@ /* $Id$ */ +/* History of cvs commits: + * + * $Log$ + * Revision 1.24 2006/11/14 17:11:15 hristov + * Removing inheritances from TAttLine, TAttMarker and AliRndm in AliModule. The copy constructor and assignment operators are moved to the private part of the class and not implemented. The corresponding changes are propagated to the detectors + * + * Revision 1.23 2006/09/13 07:31:01 kharlov + * Effective C++ corrections (T.Pocheptsov) + * + * Revision 1.22 2005/06/17 07:39:07 hristov + * Removing GetDebug and SetDebug from AliRun and AliModule. Using AliLog for the messages + * + * Revision 1.21 2005/05/28 14:19:05 schutz + * Compilation warnings fixed by T.P. + * + */ + //_________________________________________________________________________ // Implementation version vImpacts of PHOS Manager class. // This class inherits from v1 and adds impacts storing. @@ -29,8 +46,10 @@ // --- ROOT system --- -#include +//#include +#include #include +#include // --- Standard library --- @@ -40,21 +59,33 @@ #include "AliPHOSImpact.h" #include "AliPHOSvImpacts.h" #include "AliRun.h" +#include "AliLoader.h" +#include "AliMC.h" +#include "AliLog.h" ClassImp(AliPHOSvImpacts) //____________________________________________________________________________ -AliPHOSvImpacts::AliPHOSvImpacts():AliPHOSv1() +AliPHOSvImpacts::AliPHOSvImpacts(): + fEMCImpacts(new TList), + fCPVImpacts(new TList), + fPPSDImpacts(new TList), + fNEMCImpacts(), + fNCPVImpacts(), + fNPPSDImpacts() { // ctor - fEMCImpacts = 0 ; - fPPSDImpacts = 0 ; - fCPVImpacts = 0 ; } //____________________________________________________________________________ AliPHOSvImpacts::AliPHOSvImpacts(const char *name, const char *title): -AliPHOSv1(name,title) + AliPHOSv1(name,title), + fEMCImpacts(new TList), + fCPVImpacts(new TList), + fPPSDImpacts(0), + fNEMCImpacts(), + fNCPVImpacts(), + fNPPSDImpacts() { // ctor : title is used to identify the layout // @@ -68,9 +99,6 @@ AliPHOSv1(name,title) // - fEMCImpacts, fCPVImpacts which are // TList of EMC and CPV modules respectively, each // modules contains TClonesArray of AliPHOSImpacts - - fEMCImpacts = new TList(); - fCPVImpacts = new TList(); Int_t nPHOSModules = GetGeometry()->GetNModules(); Int_t nCPVModules = GetGeometry()->GetNModules(); @@ -116,7 +144,7 @@ AliPHOSvImpacts::~AliPHOSvImpacts() } //____________________________________________________________________________ -void AliPHOSvImpacts::AddImpact( char* det, Int_t shunt, Int_t primary, Int_t track, Int_t module, +void AliPHOSvImpacts::AddImpact(const char* det, Int_t shunt, Int_t primary, Int_t track, Int_t module, Int_t pid, TLorentzVector p, Float_t *xyz) { // Add an impact to the impact list. @@ -134,13 +162,14 @@ void AliPHOSvImpacts::AddImpact( char* det, Int_t shunt, Int_t primary, Int_t tr nImpacts= fNCPVImpacts[module]; fNCPVImpacts[module]++ ; } + else + AliFatal(Form("Wrong PHOS configuration: det=%s",det)); new((*impacts)[nImpacts]) AliPHOSImpact(shunt,primary,track,pid,p,xyz) ; - if (fDebug==1) { - printf("Module %d %s: ",module,det); - (dynamic_cast((impacts->At(nImpacts))))->Print(); - } + AliDebugClass(1,Form("Module %d %s: ",module,det)); + if (AliLog::GetGlobalDebugLevel()>0) + (static_cast((impacts->At(nImpacts))))->Print(); } //____________________________________________________________________________ @@ -153,8 +182,8 @@ void AliPHOSvImpacts::MakeBranch(Option_t *opt) Int_t bufferSize = 32000 ; Int_t splitlevel = 0 ; - TreeH()->Branch("PHOSEmcImpacts" , "TList", &fEMCImpacts , bufferSize, splitlevel); - TreeH()->Branch("PHOSCpvImpacts" , "TList", &fCPVImpacts , bufferSize, splitlevel); + fLoader->TreeH()->Branch("PHOSEmcImpacts" , "TList", &fEMCImpacts , bufferSize, splitlevel); + fLoader->TreeH()->Branch("PHOSCpvImpacts" , "TList", &fCPVImpacts , bufferSize, splitlevel); } @@ -167,12 +196,12 @@ void AliPHOSvImpacts::ResetHits() Int_t i; for (i=0; iGetNModules(); i++) { - (dynamic_cast(fEMCImpacts->At(i))) -> Clear(); + (static_cast(fEMCImpacts->At(i))) -> Clear(); fNEMCImpacts[i] = 0 ; } for (i=0; iGetNModules(); i++) { - (dynamic_cast(fCPVImpacts->At(i))) -> Clear(); + (static_cast(fCPVImpacts->At(i))) -> Clear(); fNCPVImpacts[i] = 0 ; } @@ -191,13 +220,14 @@ void AliPHOSvImpacts::StepManager(void) TLorentzVector pos ; // Lorentz vector of the track current position Int_t copy ; - Int_t tracknumber = gAlice->CurrentTrack() ; - Int_t primary = gAlice->GetPrimary( gAlice->CurrentTrack() ); + Int_t tracknumber = gAlice->GetMCApp()->GetCurrentTrackNumber() ; + Int_t primary = gAlice->GetMCApp()->GetPrimary( gAlice->GetMCApp()->GetCurrentTrackNumber() ); TString name = GetGeometry()->GetName() ; // Add impact to EMC - if( gMC->CurrentVolID(copy) == gMC->VolId("PXTL") && + static Int_t idPXTL = gMC->VolId("PXTL"); + if( gMC->CurrentVolID(copy) == idPXTL && gMC->IsTrackEntering() ) { gMC->TrackMomentum(pmom); gMC->TrackPosition(pos) ; @@ -213,8 +243,8 @@ void AliPHOSvImpacts::StepManager(void) gMC -> Gmtod (pm, pd, 2); // transform 3-momentum from master to daughter system // Select tracks coming to the crystal from up or down sides - if (pd[1]<0 && xyzd[1] > GetGeometry()->GetCrystalSize(1)/2-0.1 || - pd[1]>0 && xyzd[1] < -GetGeometry()->GetCrystalSize(1)/2+0.1) { + if ((pd[1]<0 && xyzd[1] > GetGeometry()->GetCrystalSize(1)/2-0.1) || + (pd[1]>0 && xyzd[1] < -GetGeometry()->GetCrystalSize(1)/2+0.1)) { // Select tracks coming to the crystal from up or down sides Int_t pid = gMC->TrackPid(); Int_t module; @@ -226,7 +256,8 @@ void AliPHOSvImpacts::StepManager(void) // Add impact to CPV - if( gMC->CurrentVolID(copy) == gMC->VolId("PCPQ") && + static Int_t idPCPQ = gMC->VolId("PCPQ"); + if( gMC->CurrentVolID(copy) == idPCPQ && gMC->IsTrackEntering() ) { gMC->TrackMomentum(pmom); gMC->TrackPosition(pos) ;