/************************************************************************** * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * * * Author: The ALICE Off-line Project. * * Contributors are mentioned in the code where appropriate. * * * * Permission to use, copy, modify and distribute this software and its * * documentation strictly for non-commercial purposes is hereby granted * * without fee, provided that the above copyright notice appears in all * * copies and that both the copyright notice and this permission notice * * appear in the supporting documentation. The authors make no claims * * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ /* $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. // Impacts stands for exact values of track coming to the detectors // EMC, CPV or PPSD. // Impacts are written to the same tree as hits are // but in separate branches. //--- //*-- Author: Yuri Kharlov (IHEP, Protvino/SUBATECH, Nantes) // --- ROOT system --- //#include #include #include #include // --- Standard library --- // --- AliRoot header files --- #include "AliPHOSGeometry.h" #include "AliPHOSImpact.h" #include "AliPHOSvImpacts.h" #include "AliRun.h" #include "AliLoader.h" #include "AliMC.h" #include "AliLog.h" ClassImp(AliPHOSvImpacts) //____________________________________________________________________________ AliPHOSvImpacts::AliPHOSvImpacts(): fEMCImpacts(new TList), fCPVImpacts(new TList), fPPSDImpacts(new TList), fNEMCImpacts(), fNCPVImpacts(), fNPPSDImpacts() { // ctor } //____________________________________________________________________________ AliPHOSvImpacts::AliPHOSvImpacts(const char *name, const char *title): AliPHOSv1(name,title), fEMCImpacts(new TList), fCPVImpacts(new TList), fPPSDImpacts(0), fNEMCImpacts(), fNCPVImpacts(), fNPPSDImpacts() { // ctor : title is used to identify the layout // // We store hits : // - fHits (the "normal" one), which retains the hits associated with // the current primary particle being tracked // (this array is reset after each primary has been tracked). // This part inherits from AliPHOSv1 // // We store impacts : // - fEMCImpacts, fCPVImpacts which are // TList of EMC and CPV modules respectively, each // modules contains TClonesArray of AliPHOSImpacts Int_t nPHOSModules = GetGeometry()->GetNModules(); Int_t nCPVModules = GetGeometry()->GetNModules(); Int_t iPHOSModule; TClonesArray * impacts; for (iPHOSModule=0; iPHOSModuleAdd(new TClonesArray("AliPHOSImpact",200)) ; fNEMCImpacts[iPHOSModule] = 0; impacts = dynamic_cast(fEMCImpacts->At(iPHOSModule)); } for (iPHOSModule=0; iPHOSModuleAdd(new TClonesArray("AliPHOSImpact",200)) ; fNCPVImpacts[iPHOSModule] = 0; impacts = dynamic_cast(fCPVImpacts->At(iPHOSModule)); } } //____________________________________________________________________________ AliPHOSvImpacts::~AliPHOSvImpacts() { // dtor // Delete hits if ( fHits ) { fHits->Delete() ; delete fHits ; fHits = 0 ; } // Delete impacts in EMC, CPV if ( fEMCImpacts ) { fEMCImpacts->Delete() ; delete fEMCImpacts ; fEMCImpacts = 0 ; } if ( fCPVImpacts ) { fCPVImpacts->Delete() ; delete fCPVImpacts ; fCPVImpacts = 0 ; } } //____________________________________________________________________________ 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. TClonesArray * impacts = 0; Int_t nImpacts = 0; if (strcmp(det,"EMC ")==0) { impacts = dynamic_cast(fEMCImpacts->At(module)); nImpacts= fNEMCImpacts[module]; fNEMCImpacts[module]++ ; } else if (strcmp(det,"CPV ")==0) { impacts = dynamic_cast(fCPVImpacts->At(module)); nImpacts= fNCPVImpacts[module]; fNCPVImpacts[module]++ ; } else AliFatal(Form("Wrong PHOS configuration: det=%s",det)); new((*impacts)[nImpacts]) AliPHOSImpact(shunt,primary,track,pid,p,xyz) ; AliDebugClass(1,Form("Module %d %s: ",module,det)); if (AliLog::GetGlobalDebugLevel()>0) (static_cast((impacts->At(nImpacts))))->Print(); } //____________________________________________________________________________ void AliPHOSvImpacts::MakeBranch(Option_t *opt) { // Create new branch in the current Hits Root Tree containing // a list of PHOS impacts (exact values of track coming to detector) AliDetector::MakeBranch(opt); Int_t bufferSize = 32000 ; Int_t splitlevel = 0 ; fLoader->TreeH()->Branch("PHOSEmcImpacts" , "TList", &fEMCImpacts , bufferSize, splitlevel); fLoader->TreeH()->Branch("PHOSCpvImpacts" , "TList", &fCPVImpacts , bufferSize, splitlevel); } //____________________________________________________________________________ void AliPHOSvImpacts::ResetHits() { // Reset impact branches for EMC, CPV and PPSD AliDetector::ResetHits(); Int_t i; for (i=0; iGetNModules(); i++) { (static_cast(fEMCImpacts->At(i))) -> Clear(); fNEMCImpacts[i] = 0 ; } for (i=0; iGetNModules(); i++) { (static_cast(fCPVImpacts->At(i))) -> Clear(); fNCPVImpacts[i] = 0 ; } } //_____________________________________________________________________________ void AliPHOSvImpacts::StepManager(void) { // Find impacts (tracks which enter the EMC, CPV) // and add them to to the impact lists AliPHOSv1::StepManager(); Float_t xyzm[3], xyzd[3], pm[3], pd[3]; TLorentzVector pmom ; // Lorentz momentum of the particle initiated hit TLorentzVector pos ; // Lorentz vector of the track current position Int_t copy ; Int_t tracknumber = gAlice->GetMCApp()->GetCurrentTrackNumber() ; Int_t primary = gAlice->GetMCApp()->GetPrimary( gAlice->GetMCApp()->GetCurrentTrackNumber() ); TString name = GetGeometry()->GetName() ; // Add impact to EMC static Int_t idPXTL = TVirtualMC::GetMC()->VolId("PXTL"); if( TVirtualMC::GetMC()->CurrentVolID(copy) == idPXTL && TVirtualMC::GetMC()->IsTrackEntering() ) { TVirtualMC::GetMC()->TrackMomentum(pmom); TVirtualMC::GetMC()->TrackPosition(pos) ; Int_t i; for (i=0; i<3; i++) xyzm[i] = pos[i]; for (i=0; i<3; i++) { xyzm[i] = pos[i] ; pm[i] = pmom[i]; } TVirtualMC::GetMC() -> Gmtod (xyzm, xyzd, 1); // transform coordinate from master to daughter system TVirtualMC::GetMC() -> 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)) { // Select tracks coming to the crystal from up or down sides Int_t pid = TVirtualMC::GetMC()->TrackPid(); Int_t module; TVirtualMC::GetMC()->CurrentVolOffID(10,module); module--; AddImpact("EMC ",fIshunt, primary,tracknumber, module, pid, pmom, xyzm); } } // Add impact to CPV static Int_t idPCPQ = TVirtualMC::GetMC()->VolId("PCPQ"); if( TVirtualMC::GetMC()->CurrentVolID(copy) == idPCPQ && TVirtualMC::GetMC()->IsTrackEntering() ) { TVirtualMC::GetMC()->TrackMomentum(pmom); TVirtualMC::GetMC()->TrackPosition(pos) ; Int_t i; for (i=0; i<3; i++) xyzm[i] = pos[i]; for (i=0; i<3; i++) { xyzm[i] = pos[i] ; pm[i] = pmom[i]; } Int_t pid = TVirtualMC::GetMC()->TrackPid(); Int_t module; TVirtualMC::GetMC()->CurrentVolOffID(3,module); module--; AddImpact("CPV ",fIshunt, primary,tracknumber, module, pid, pmom, xyzm); } }