* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-*/
+/* $Id$ */
///////////////////////////////////////////////////////////////////////////////
// //
// //
// //
///////////////////////////////////////////////////////////////////////////////
-#include <stdlib.h>
-#include <TTUBE.h>
-#include <TNode.h>
+#include <Riostream.h>
+
+#include <TFile.h>
#include <TGeometry.h>
+#include <TMath.h>
+#include <TNode.h>
+#include <TParticle.h>
+#include <TRandom.h>
+#include <TTUBE.h>
+#include <TVirtualMC.h>
+#include <AliESD.h>
+
+#include "AliLoader.h"
#include "AliRun.h"
#include "AliSTART.h"
-#include <iostream.h>
-#include <fstream.h>
+#include "AliSTARTLoader.h"
+#include "AliSTARTdigit.h"
+#include "AliSTARThit.h"
+#include "AliSTARThitPhoton.h"
#include "AliMC.h"
-//#include "TGeant3.h"
-
+#include "AliSTARTDigitizer.h"
+
ClassImp(AliSTART)
-
+
+static AliSTARTdigit *digits;
+
//_____________________________________________________________________________
AliSTART::AliSTART()
{
//
// Default constructor for class AliSTART
//
- fIshunt = 0;
+ fIshunt = 1;
+ fHits = 0;
+ fDigits = 0;
+ fPhotons = 0;
}
//_____________________________________________________________________________
//
// Standard constructor for START Detector
//
-
- AliModule *fmd = gAlice->GetModule("FMD");
- if(fmd) {
- Int_t fmdversion = fmd->IsVersion();
- if(fmdversion==0 || fmdversion==1) {
- Error("ctor","Versions 0 and 1 of FMD incompatible with START\n");
- exit(1);
- }
- }
+
+
//
// Initialise Hit array
- fHits = new TClonesArray("AliSTARThit", 405);
-
- fIshunt = 0;
- fIdSens1 = 0;
+ fHits = new TClonesArray("AliSTARThit", 405);
+ gAlice->GetMCApp()->AddHitList(fHits);
+ fPhotons = new TClonesArray("AliSTARThitPhoton", 10000);
+ gAlice->GetMCApp()->AddHitList (fPhotons);
+ if (GetDebug()>2) cout<<" Debug "<<endl;
+ fIshunt = 1;
+ fIdSens = 0;
+ fNPhotons = 0;
SetMarkerColor(kRed);
}
+
+//_____________________________________________________________________________
+AliSTART::~AliSTART() {
+ if (fHits) {
+ fHits->Delete();
+ delete fHits;
+ }
+ if (fPhotons) {
+ fPhotons->Delete();
+ delete fPhotons;
+ }
+}
//_____________________________________________________________________________
void AliSTART::AddHit(Int_t track, Int_t *vol, Float_t *hits)
TClonesArray &lhits = *fHits;
new(lhits[fNhits++]) AliSTARThit(fIshunt,track,vol,hits);
}
-
+
+//_____________________________________________________________________________
+void AliSTART::AddHitPhoton(Int_t track, Int_t *vol, Float_t *hits)
+{
+ // Add a START hit of photons
+
+ TClonesArray &lhits = *fPhotons;
+ new(lhits[fNPhotons++]) AliSTARThitPhoton(fIshunt,track,vol,hits);
+}
+
+//_____________________________________________________________________________
+
+void AliSTART::AddDigit(Int_t * /*tracks*/, Int_t * /*digits*/)
+{
+
+ // Add a START digit to the list. Dummy function.
+
+}
+
//_____________________________________________________________________________
void AliSTART::BuildGeometry()
{
//
// Build simple ROOT TNode geometry for event display
//
- TNode *Node, *Top;
+ TNode *node, *top;
const int kColorSTART = 19;
- //
- Top=gAlice->GetGeometry()->GetNode("alice");
+
+ top=gAlice->GetGeometry()->GetNode("alice");
// START define the different volumes
- new TRotMatrix("rot999","rot999", 90,0,90,90,180,0);
-
- new TTUBE("S_STR1","START volume 1","void",5.,10.7,5.3);
- Top->cd();
- Node = new TNode("STR1","STR1","S_STR1",0,0,75.,"");
- Node->SetLineColor(kColorSTART);
- fNodes->Add(Node);
-
- new TTUBE("S_STR2","START volume 2","void",5.,10.7,5.3);
- Top->cd();
- Node = new TNode("STR2","STR2","S_STR2",0,0,-75,"rot999");
- Node->SetLineColor(kColorSTART);
- fNodes->Add(Node);
+ new TRotMatrix("rotx999","rot999", 90,0,90,90,180,0);
+
+ new TTUBE("S_0ST1","START volume 1","void",5.,10.7,5.3);
+ top->cd();
+ node = new TNode("0ST1","0ST01","S_0ST1",0,0,-69.7,"");
+ node->SetLineColor(kColorSTART);
+ fNodes->Add(node);
+
+ new TTUBE("S_0ST2","START volume 2","void",5.,10.7,5.3);
+ top->cd();
+ node = new TNode("0ST2","0ST2","S_0ST2",0,0,350,"rotx999");
+ node->SetLineColor(kColorSTART);
+ fNodes->Add(node);
}
//_____________________________________________________________________________
-Int_t AliSTART::DistanceToPrimitive(Int_t px, Int_t py)
+Int_t AliSTART::DistanceToPrimitive(Int_t /*px*/, Int_t /*py*/)
{
//
// Calculate the distance from the mouse to the START on the screen
return 9999;
}
-//_____________________________________________________________________________
-
//-------------------------------------------------------------------------
void AliSTART::Init()
{
//
// Initialis the START after it has been built
Int_t i;
- AliMC* pMC = AliMC::GetMC();
- //
- printf("\n");
- for(i=0;i<35;i++) printf("*");
- printf(" START_INIT ");
- for(i=0;i<35;i++) printf("*");
- printf("\n");
//
- // Here the START initialisation code (if any!)
- for(i=0;i<80;i++) printf("*");
- printf("\n");
- //
- //
- fIdSens1=pMC->VolId("PTOP");
-
+ if(fDebug) {
+ printf("\n%s: ",ClassName());
+ for(i=0;i<35;i++) printf("*");
+ printf(" START_INIT ");
+ for(i=0;i<35;i++) printf("*");
+ printf("\n%s: ",ClassName());
+ //
+ // Here the START initialisation code (if any!)
+ for(i=0;i<80;i++) printf("*");
+ printf("\n");
+ }
}
//---------------------------------------------------------------------------
void AliSTART::MakeBranch(Option_t* option)
{
-
+ //
+ // Specific START branches
+ //
// Create Tree branches for the START.
Int_t buffersize = 4000;
- char branchname[10];
+ char branchname[20];
sprintf(branchname,"%s",GetName());
+
+ const char *cD = strstr(option,"D");
+ const char *cH = strstr(option,"H");
+
+ if (cH && fLoader->TreeH())
+ {
+ if (fPhotons == 0x0) fPhotons = new TClonesArray("AliSTARThitPhoton", 10000);
+ sprintf (branchname, "%shitPhoton", GetName());
+ MakeBranchInTree (fLoader->TreeH(), branchname, &fPhotons, 50000, 0);
+ if (fHits == 0x0) fHits = new TClonesArray("AliSTARThit", 405);
+ }
+
AliDetector::MakeBranch(option);
- char *D = strstr(option,"D");
+ if (cD) {
+ digits = new AliSTARTdigit();
+ MakeBranchInTree(fLoader->TreeD(), branchname, "AliSTARTdigit", digits, buffersize, 1, 0);
+ }
+}
- if (fDigits && gAlice->TreeD() && D) {
- gAlice->TreeD()->Branch(branchname,&fDigits, buffersize);
- printf("Making Branch %s for digits\n",branchname);
- }
+//_____________________________________________________________________________
+void AliSTART::ResetHits()
+{
+ AliDetector::ResetHits();
-}
-
-ClassImp(AliSTARThit)
-
+ fNPhotons = 0;
+ if (fPhotons) fPhotons->Clear();
+}
+
//_____________________________________________________________________________
-AliSTARThit::AliSTARThit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits):
- AliHit(shunt, track)
+void AliSTART::SetTreeAddress()
{
- //
- // Add a START hit
- //
+ TBranch *branch;
+ TTree *treeH;
+
+
+ treeH = TreeH();
+
+ if (treeH)
+ {
+ if (fPhotons == 0x0) fPhotons = new TClonesArray("AliSTARThitPhoton", 10000);
+ branch = treeH->GetBranch("STARThitPhoton");
+ if (branch) branch->SetAddress(&fPhotons);
+ if (fHits == 0x0) fHits = new TClonesArray("AliSTARThit", 405);
+ }
+
+ AliDetector::SetTreeAddress();
-// Int_t i;
- fVolume = vol[0];
- fPmt=vol[1];
-//printf("fvolume %d\n",fVolume);
-//printf("fpmt %d\n",fPmt);
-
- fX=hits[0];
- fY=hits[1];
- fZ=hits[2];
- fEdep=hits[3];
- fEtot=hits[4];
- fParticle=Int_t (hits[5]);
- fTime=hits[6];
-
-//for (i=0; i<=6; i++) {printf("Hits up %f\n",hits[i]);}
}
-// ClassImp(AliSTARTdigit)
+//______________________________________________________________________
+AliLoader* AliSTART::MakeLoader(const char* topfoldername)
+{
+ Info("MakeLoader", "Creating AliSTARTLoader. Top folder is %s.", topfoldername);
+ fLoader = new AliSTARTLoader(GetName(), topfoldername);
+ return fLoader;
+}
+
+//_____________________________________________________________________________
+AliDigitizer* AliSTART::CreateDigitizer(AliRunDigitizer* manager) const
+{
+ return new AliSTARTDigitizer(manager);
+}