/************************************************************************** * 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. * **************************************************************************/ /* $Log$ Revision 1.26 2002/10/22 15:40:19 alibrary Introducing Riostream.h Revision 1.25 2002/10/14 14:57:32 hristov Merging the VirtualMC branch to the main development branch (HEAD) Revision 1.23.8.1 2002/07/24 09:50:10 alibrary Updating VirtualMC Revision 1.24 2002/07/23 11:48:05 alla new Digits structure Revision 1.23 2001/09/19 18:41:45 alla Asimmetric START geometry Revision 1.22 2001/07/27 13:03:12 hristov Default Branch split level set to 99 Revision 1.21 2001/06/27 16:06:59 hristov Rotation matrix in BuildGeometry has been changed to rotx999 Revision 1.20.2.1 2001/06/27 10:51:15 alla Rotation matrix in BuildGeometry has benn changed to rotx999 Revision 1.20 2001/05/16 14:57:21 alibrary New files for folders and Stack Revision 1.19 2001/04/04 12:10:18 alla changes according Coding Convension Revision 1.18 2001/03/12 17:46:43 hristov Changes needed on Sun with CC 5.0 Revision 1.17 2001/01/26 19:59:53 hristov Major upgrade of AliRoot code Revision 1.16 2001/01/17 10:56:08 hristov Corrections to destructors Revision 1.15 2001/01/01 13:10:42 hristov Local definition of digits removed Revision 1.14 2000/12/22 16:17:15 hristov Updated START code from Alla Revision 1.13 2000/12/18 11:39:41 alibrary Quick fix to avoid crash in display waiting for new version Revision 1.12 2000/12/04 08:48:19 alibrary Fixing problems in the HEAD Revision 1.11 2000/10/13 13:14:08 hristov Bug fixes and code cleaning Revision 1.10 2000/10/02 21:28:13 fca Removal of useless dependecies via forward declarations Revision 1.9 2000/07/13 16:41:29 fca New START corrected for coding conventions Revision 1.8 2000/03/27 17:24:25 alla Modifing geometry Revision 1.6 2000/01/21 15:45:23 fca New Version from Alla Revision 1.5 2000/01/19 17:17:15 fca Introducing a list of lists of hits -- more hits allowed for detector now Revision 1.4 1999/11/12 15:04:00 fca Modifications from A.Maevskaya Revision 1.3 1999/09/29 09:24:29 fca Introduction of the Copyright and cvs Log */ /////////////////////////////////////////////////////////////////////////////// // // // START (T-Zero) Detector // // This class contains the base procedures for the START // // detector // // // //Begin_Html /*

The responsible person for this module is Alla Maevskaia.

*/
//End_Html
//                                                                           //
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include 
#include 

#include "TMath.h"
#include "TTUBE.h"
#include "TNode.h"
#include "TRandom.h"
#include "TGeometry.h"
#include "TFile.h"
#include "TParticle.h"

#include "AliRun.h"
#include "AliSTART.h"
#include "AliSTARTdigit.h"
#include "AliSTARThit.h"
#include "AliSTARThitPhoton.h"
#include "AliSTARTvertex.h"

ClassImp(AliSTART)

static  AliSTARTdigit *digits; 

//_____________________________________________________________________________
AliSTART::AliSTART()
{
  //
  // Default constructor for class AliSTART
  //
  fIshunt   = 1;
  fHits     = 0;
  fDigits   = 0;
  fPhotons  = 0;
}
 
//_____________________________________________________________________________
AliSTART::AliSTART(const char *name, const char *title)
       : AliDetector(name,title)
{
  //
  // Standard constructor for START Detector
  //

  
  //
  // Initialise Hit array
  fHits       = new TClonesArray("AliSTARThit",  405);
  gAlice->AddHitList(fHits);

  fPhotons  = new TClonesArray("AliSTARThitPhoton", 10000);
  gAlice->AddHitList (fPhotons);
  
  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)
{
  //
  // Add a START hit
  //
  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;
  const int kColorSTART  = 19;

  top=gAlice->GetGeometry()->GetNode("alice");

  // START define the different volumes
  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)
{
  //
  // Calculate the distance from the mouse to the START on the screen
  // Dummy routine
  //
  return 9999;
}
 
//-------------------------------------------------------------------------
void AliSTART::Init()
{
  //
  // Initialis the START after it has been built
  Int_t i;
  //
  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, const char *file)
{
  //
  // Specific START branches
  //
  // Create Tree branches for the START.
  Int_t buffersize = 4000;
  char branchname[20];
  sprintf(branchname,"%s",GetName());

  AliDetector::MakeBranch(option,file);

  const char *cD = strstr(option,"D");
  const char *cH = strstr(option,"H");
  
  if (cH)
  {
     sprintf (branchname, "%shitPhoton", GetName());
     MakeBranchInTree (gAlice->TreeH(), branchname, &fPhotons, 50000, file);
  } 

  
  if (cD) {
    digits = new AliSTARTdigit();
    MakeBranchInTree(gAlice->TreeD(), 
                     branchname, "AliSTARTdigit", digits, buffersize, 1, file);
  } 
/*
  char *cR = strstr(option,"R");
  
  if (cR)   {  
    MakeBranchInTree(gAlice->TreeR(), 
                     branchname, "Int_t", &fZposit, buffersize, 1, file);
  }
  */
}    

//_____________________________________________________________________________
void AliSTART::ResetHits()
{
  AliDetector::ResetHits();
  
  fNPhotons = 0;
  if (fPhotons)  fPhotons->Clear();
}

//_____________________________________________________________________________
void AliSTART::SetTreeAddress()
{
  TBranch	*branch;
  TTree		*treeH;

  AliDetector::SetTreeAddress();
  treeH = gAlice->TreeH();
  
  if (treeH)
    if (fPhotons)
    {
       branch = treeH->GetBranch ("STARThitPhoton");
       if (branch)  branch->SetAddress (&fPhotons);
    }
}


//_____________________________________________________________________________

void AliSTART::Hit2digit(Int_t evnum) 
{
  //
  // From hits to digits
  //
  /*
  Float_t x,y,e;
  Int_t nbytes = 0;
  Int_t hit;
  Int_t nhits;
  Int_t volume,pmt;
  char nameTH[8],nameTD[8];
  Float_t timediff,timeright,timeleft,timeav;
  Float_t besttimeright,besttimeleft,meanTime;
  Int_t channelWidth=10;

  TParticle *particle;
  AliSTARThit  *startHit;

  Int_t buffersize=256;

  digits= new AliSTARTdigit();
  TBranch *bDig=0;

  
 // Event ------------------------- LOOP  
 
    sprintf(nameTD,"TreeD%d",evnum);
    TTree *td = new TTree(nameTD,"START");
    bDig = td->Branch("START","AliSTARTdigit",&digits,buffersize);

    besttimeright=9999.;
    besttimeleft=9999.;
    Int_t timeDiff=0;
    Int_t timeAv=0;

    Int_t nparticles = gAlice->GetEvent(evnum);
    if (nparticles <= 0) return;
    printf("\n nparticles %d\n",nparticles);
    
   
    sprintf(nameTH,"TreeH%d",evnum);
    printf("%s\n",nameTH);
    TTree *th = gAlice->TreeH();
    Int_t ntracks    = (Int_t) th->GetEntries();
    if (ntracks<=0) return;
    // Start loop on tracks in the hits containers
    for (Int_t track=0; trackResetHits();
      nbytes += th->GetEvent(track);
      particle=gAlice->Particle(track);
      nhits =fHits->GetEntriesFast();
      
      for (hit=0;hitUncheckedAt(hit);
	pmt=startHit->fPmt;
	e=startHit->fEtot;
	x=startHit->X();
	y=startHit->Y();
	volume = startHit->fVolume;
	if(volume==1){
	  timeright = startHit->fTime;
	  if(timerightfTime;
	  //                printf("timeleft %f\n",timeleft);
	  if(timeleftGaus(besttimeright,0.05);
    Float_t besttimeleftGaus=gRandom->Gaus(besttimeleft,0.05);
    timediff=besttimerightGaus-besttimeleftGaus;
    meanTime=(besttimerightGaus+besttimeleftGaus)/2.;
    if ( TMath::Abs(timediff)<2. && meanTime<3.) 
     {
       //we assume centre of bunch is 5ns after TTS signal
       //TOF values are relative of the end of bunch
       Float_t ppBunch=25;
    
       ppBunch=ppBunch-10/2;
       Float_t t1=1000.*besttimeleftGaus;
       Float_t t2=1000.*besttimerightGaus;
       t1=t1/channelWidth+ppBunch; //time in ps to channelWidth
       t2=t2/channelWidth+ppBunch; //time in ps to channelWidth
     
       timeav=(t1+t2)/2.;
     
       // Time to TDC signal
       // 256 channels for timediff, range 1ns
       
       timediff=128+1000*timediff/channelWidth; // time in ps 

       timeAv = (Int_t)(timeav);   // time (ps) channel numbres
       timeDiff = (Int_t)(timediff); // time ( ps) channel numbres
       digits->Set(timeAv,timeDiff);
     }
    else
      {timeAv=999999; timeDiff=99999;}
      
    td->Fill();
    printf("digits-> %d \n",digits->GetTime());
    td->Write();
  */
}