X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PMD%2FAliPMD.cxx;h=ba5f03387df89441704fc9d11f696d87aee1690d;hb=307b389fd7f04c274e3f7ca49cad4b62e8671b11;hp=4e93479d85dc4e2a030fa3591a545d2308e001a5;hpb=d3b63a89a43f2b37eb9ae649527d86304f47deb2;p=u%2Fmrichter%2FAliRoot.git diff --git a/PMD/AliPMD.cxx b/PMD/AliPMD.cxx index 4e93479d85d..ba5f03387df 100644 --- a/PMD/AliPMD.cxx +++ b/PMD/AliPMD.cxx @@ -13,27 +13,7 @@ * provided "as is" without express or implied warranty. * **************************************************************************/ -/* -$Log$ -Revision 1.10 2000/11/06 09:07:13 morsch -Set fRecPoints to zero in default constructor. - -Revision 1.9 2000/10/30 09:03:59 morsch -Prototype for PMD reconstructed hits (AliPMDRecPoint) added. - -Revision 1.8 2000/10/20 06:24:40 fca -Put the PMD at the right position in the event display - -Revision 1.7 2000/10/02 21:28:12 fca -Removal of useless dependecies via forward declarations - -Revision 1.6 2000/01/19 17:17:06 fca -Introducing a list of lists of hits -- more hits allowed for detector now - -Revision 1.5 1999/09/29 09:24:27 fca -Introduction of the Copyright and cvs Log - -*/ +/* $Id$ */ /////////////////////////////////////////////////////////////////////////////// // @@ -59,16 +39,22 @@ Introduction of the Copyright and cvs Log /////////////////////////////////////////////////////////////////////////////// #include +#include +#include #include #include -#include -#include +#include +#include "AliLog.h" +#include "AliLoader.h" +#include "AliPMDLoader.h" #include "AliPMD.h" #include "AliRun.h" -#include "AliMC.h" -#include "AliConst.h" -#include "AliPMDRecPoint.h" +#include "AliMC.h" +#include "AliPMDDigitizer.h" +#include "AliPMDhit.h" +#include "AliPMDDDLRawData.h" +#include "AliPMDRawToSDigits.h" ClassImp(AliPMD) @@ -79,7 +65,7 @@ AliPMD::AliPMD() // Default constructor // fIshunt = 0; - fRecPoints = NULL; + } //_____________________________________________________________________________ @@ -93,13 +79,10 @@ AliPMD::AliPMD(const char *name, const char *title) // // Allocate the array of hits fHits = new TClonesArray("AliPMDhit", 405); - gAlice->AddHitList(fHits); + gAlice->GetMCApp()->AddHitList(fHits); - fRecPoints = new TClonesArray("AliPMDRecPoint",10000); - fNRecPoints = 0; - - fIshunt = 1; + fIshunt = 0; fPar[0] = 1; fPar[1] = 1; @@ -119,13 +102,29 @@ AliPMD::AliPMD(const char *name, const char *title) fPadSize[3] = 1.5; } +AliLoader* AliPMD::MakeLoader(const char* topfoldername) +{ + // Makes PMD Loader + + fLoader = new AliPMDLoader(GetName(),topfoldername); + + if (fLoader) + { + AliDebug(100,"Success"); + } + else + { + AliError("Failure"); + } + + return fLoader; +} + AliPMD::~AliPMD() { // - // Default constructor + // Destructor // - delete fRecPoints; - fNRecPoints=0; } //_____________________________________________________________________________ @@ -136,11 +135,13 @@ void AliPMD::AddHit(Int_t track, Int_t *vol, Float_t *hits) // TClonesArray &lhits = *fHits; AliPMDhit *newcell, *curcell; -// printf("PMD++ Adding energy %f, prim %d, vol %d %d %d %d\n", -// hits[3],gAlice->GetPrimary(track-1),vol[0],vol[1],vol[2],vol[3]); + // printf("PMD++ Adding energy %f, prim %d, vol %d %d %d %d %d %d %d %d\n", + // hits[3],gAlice->GetPrimary(track-1),vol[0],vol[1],vol[2],vol[3], + // vol[4],vol[5],vol[6],vol[7]); + newcell = new AliPMDhit(fIshunt, track, vol, hits); - Int_t i; - for (i=0; iGetGeometry()->GetNode("alice"); - - // PMD - new TBRIK("S_PMD","PMD box","void",300,300,5); - Top->cd(); - Node = new TNode("PMD","PMD","S_PMD",0,0,-600,""); - Node->SetLineColor(kColorPMD); - fNodes->Add(Node); -} - -//_____________________________________________________________________________ -Int_t AliPMD::DistancetoPrimitive(Int_t , Int_t ) -{ - // - // Distance from mouse to detector on the screen - // dummy routine - // - return 9999; -} - //_____________________________________________________________________________ void AliPMD::SetPAR(Float_t p1, Float_t p2, Float_t p3,Float_t p4) { @@ -243,87 +213,117 @@ void AliPMD::StepManager() // } -void AliPMD::AddRecPoint(const AliPMDRecPoint &p) -{ - // - // Add a PMD reconstructed hit to the list - // - TClonesArray &lrecpoints = *fRecPoints; - new(lrecpoints[fNRecPoints++]) AliPMDRecPoint(p); -} - void AliPMD::MakeBranch(Option_t* option) { // Create Tree branches for the PMD - printf("Make Branch - TreeR address %p\n",gAlice->TreeR()); - const Int_t kBufferSize = 4000; - char branchname[30]; - + const char *cH = strstr(option,"H"); + if (cH && fLoader->TreeH() && (fHits == 0x0)) + fHits = new TClonesArray("AliPMDhit", 405); + AliDetector::MakeBranch(option); - - sprintf(branchname,"%sRecPoints",GetName()); - if (fRecPoints && gAlice->TreeR()) { - gAlice->TreeR()->Branch(branchname, &fRecPoints, kBufferSize); - printf("Making Branch %s for reconstructed hits\n",branchname); - } } void AliPMD::SetTreeAddress() { - // Set branch address for the TreeR - char branchname[30]; - AliDetector::SetTreeAddress(); + // Set branch address - TBranch *branch; - TTree *treeR = gAlice->TreeR(); + if (fLoader->TreeH() && fHits==0x0) + fHits = new TClonesArray("AliPMDhit", 405); + + AliDetector::SetTreeAddress(); +} - sprintf(branchname,"%s",GetName()); - if (treeR && fRecPoints) { - branch = treeR->GetBranch(branchname); - if (branch) branch->SetAddress(&fRecPoints); - } +void AliPMD::SetCpvOff() +{ + // Set the CPV plane off } +void AliPMD::SetPreOff() +{ + // Set the Preshower plane off -void AliPMD::ResetHits() +} +void AliPMD::SetModuleOff(Int_t /*imodule*/) { - // - // Reset number of hits and the hits array - // - AliDetector::ResetHits(); - fNRecPoints = 0; - if (fRecPoints) fRecPoints->Clear(); + // Set the desired module off + } -/////////////////////////////////////////////////////////////////////////////// -// // -// Photon Multiplicity Detector Version 1 // -// // -//Begin_Html -/* - -*/ -//End_Html -// // -/////////////////////////////////////////////////////////////////////////////// +//____________________________________________________________________________ +void AliPMD::Hits2SDigits() +{ +// create summable digits + AliRunLoader* runLoader = fLoader->GetRunLoader(); + AliPMDDigitizer* pmdDigitizer = new AliPMDDigitizer; + pmdDigitizer->OpengAliceFile(fLoader->GetRunLoader()->GetFileName().Data(), + "HS"); + pmdDigitizer->SetZPosition(361.5); + for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) { + pmdDigitizer->Hits2SDigits(iEvent); + } + fLoader->UnloadHits(); + fLoader->UnloadSDigits(); + delete pmdDigitizer; +} +//____________________________________________________________________________ +void AliPMD::SDigits2Digits() +{ + // creates sdigits to digits +} +//____________________________________________________________________________ +void AliPMD::Hits2Digits() +{ +// create digits + + AliRunLoader* runLoader = fLoader->GetRunLoader(); + AliPMDDigitizer* pmdDigitizer = new AliPMDDigitizer; + pmdDigitizer->OpengAliceFile(fLoader->GetRunLoader()->GetFileName().Data(), + "HD"); + pmdDigitizer->SetZPosition(361.5); + + for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) { + pmdDigitizer->Hits2Digits(iEvent); + } + fLoader->UnloadHits(); + fLoader->UnloadDigits(); + delete pmdDigitizer; -ClassImp(AliPMDhit) - -//_____________________________________________________________________________ -AliPMDhit::AliPMDhit(Int_t shunt,Int_t track, Int_t *vol, Float_t *hits): - AliHit(shunt, track) +} +// --------------------------------------------------------------------------- +AliDigitizer* AliPMD::CreateDigitizer(AliRunDigitizer* manager) const +{ + return new AliPMDDigitizer(manager); +} +// --------------------------------------------------------------------------- +void AliPMD::Digits2Raw() +{ +// convert digits of the current event to raw data + + fLoader->LoadDigits(); + TTree* digits = fLoader->TreeD(); + if (!digits) { + AliError("No digits tree"); + return; + } + + AliPMDDDLRawData rawWriter; + rawWriter.WritePMDRawData(digits); + + fLoader->UnloadDigits(); +} + +Bool_t AliPMD::Raw2SDigits(AliRawReader *rawReader) { - // - // Add a PMD hit - // - Int_t i; - for (i=0;i<5;i++) fVolume[i] = vol[i]; - fX=hits[0]; - fY=hits[1]; - fZ=hits[2]; - fEnergy=hits[3]; + // converts raw to sdigits + AliRunLoader* runLoader = fLoader->GetRunLoader(); + //runLoader->GetEvent(ievt); + + AliPMDRawToSDigits pmdr2sd; + pmdr2sd.Raw2SDigits(runLoader, rawReader); + fLoader->UnloadSDigits(); + return kTRUE; } - +