// integrated ITS test beam of Nov. 04 //
// Author: mercedes.lopez.noriega@cern.ch //
////////////////////////////////////////////////////////
-#include "AliRun.h"
-#include "AliITSvBeamTestITS04.h"
-#include <TClonesArray.h>
-#include <TString.h>
+
+
#include "AliITS.h"
-#include "AliITSDetTypeSim.h"
#include "AliITSgeom.h"
#include "AliITShit.h"
-#include "AliITSresponseSDD.h"
-#include "AliITSCalibrationSDD.h"
-#include "AliITSCalibrationSPD.h"
-#include "AliITSCalibrationSSD.h"
-#include "AliITSsegmentationSDD.h"
-#include "AliITSsegmentationSPD.h"
-#include "AliITSsegmentationSSD.h"
-#include "AliMagF.h"
-#include "TVirtualMC.h"
+#include "AliITSvBeamTestITS04.h"
#include "AliMC.h"
+#include "AliMagF.h"
+#include "AliRun.h"
+#include "AliTrackReference.h"
+
+#include <TVirtualMC.h>
+#include <TClonesArray.h>
+#include <TGeoGlobalMagField.h>
+#include <TGeoManager.h>
+#include <TGeoMatrix.h>
+#include <TGeoPcon.h>
+#include <TGeoVolume.h>
+#include <TLorentzVector.h>
+#include <TString.h>
const Int_t AliITSvBeamTestITS04::fgkNumberOfSPD = 4;
const Int_t AliITSvBeamTestITS04::fgkNumberOfSDD = 2;
//===============================================================
-#include <TLorentzVector.h>
-#include "AliTrackReference.h"
-#include "AliITSDetTypeSim.h"
-#include "AliITSgeom.h"
-#include "AliITSgeomSDD.h"
-#include "AliITSgeomSPD.h"
-#include "AliITSgeomSSD.h"
-#include "AliITShit.h"
-#include "AliITSCalibrationSDD.h"
-#include "AliITSCalibrationSPD.h"
-#include "AliITSCalibrationSSD.h"
-#include "AliITSsegmentationSDD.h"
-#include "AliITSsegmentationSPD.h"
-#include "AliITSsegmentationSSD.h"
-
-#include <TGeoManager.h>
-#include <TGeoVolume.h>
-#include <TGeoPcon.h>
ClassImp(AliITSvBeamTestITS04)
//_____________________________________________________________
-AliITSvBeamTestITS04::AliITSvBeamTestITS04() : AliITS()
+ AliITSvBeamTestITS04::AliITSvBeamTestITS04() :
+AliITS(), // Base class
+fITSmotherVolume(0), // Pointer to ITS mother volume.
+fNspd(fgkNumberOfSPD), //Number of SPD modules
+fNsdd(fgkNumberOfSDD), //Number of SDD modules
+fNssd(fgkNumberOfSSD), //Number of SSD modules
+fGeomDetOut(kFALSE), // Flag to write .det file out
+fGeomDetIn(kFALSE), // Flag to read geometry file (JC)
+fWrite(), //! file name to write .det file
+fRead(), // file name to read .det file (JC)
+fMajorVersion(kvITS04),// Major Version
+fMinorVersion(1), // Minor Version
+fIgm(kvITS04) //! Init geometry object
{
//
// Constructor
//
-
- SetNumberOfSPD(fgkNumberOfSPD);
- SetNumberOfSDD(fgkNumberOfSDD);
- SetNumberOfSSD(fgkNumberOfSSD);
-
+
fIdN = 3;
fIdName = new TString[fIdN];
fIdName[0] = fgSPDsensitiveVolName;
fIdSens = new Int_t[fIdN];
for(Int_t i=0; i<fIdN; i++) fIdSens[i] = 0;
- //for writing out geometry
- fGeomDetOut = kFALSE;
-
- // for reading in geometry (JC)
- fGeomDetIn = kFALSE;
-
for(Int_t a=0;a<60;a++) fWrite[a] = '\0';
}
-
//_____________________________________________________________
-AliITSvBeamTestITS04::AliITSvBeamTestITS04(const char* name,const char *title)
- : AliITS(name,title)
+AliITSvBeamTestITS04::AliITSvBeamTestITS04(const char* name,const char *title):
+AliITS(name,title), // Base class
+fITSmotherVolume(0), // Pointer to ITS mother volume.
+fNspd(fgkNumberOfSPD), //Number of SPD modules
+fNsdd(fgkNumberOfSDD), //Number of SDD modules
+fNssd(fgkNumberOfSSD), //Number of SSD modules
+fGeomDetOut(kFALSE), // Flag to write .det file out
+fGeomDetIn(kFALSE), // Flag to read geometry file (JC)
+fWrite(), //! file name to write .det file
+fRead(), // file name to read .det file (JC)
+fMajorVersion(kvITS04),// Major Version
+fMinorVersion(1), // Minor Version
+fIgm(kvITS04) //! Init geometry object
{
//
// Constructor
//
-
- SetNumberOfSPD(fgkNumberOfSPD);
- SetNumberOfSDD(fgkNumberOfSDD);
- SetNumberOfSSD(fgkNumberOfSSD);
-
+
fIdN = 3;
fIdName = new TString[fIdN];
fIdName[0] = fgSPDsensitiveVolName;
fIdName[2] = fgSSDsensitiveVolName;
fIdSens = new Int_t[fIdN];
for(Int_t i=0; i<fIdN; i++) fIdSens[i] = 0;
-
- //for writing out geometry
- fGeomDetOut = kFALSE; // Don't write .det file
-
- // for reading in geometry (JC)
- fGeomDetIn = kFALSE;
-
for(Int_t a=0;a<60;a++) fWrite[a] = '\0';
}
-
-//______________________________________________________________________
-AliITSvBeamTestITS04::AliITSvBeamTestITS04(const AliITSvBeamTestITS04 &source) : AliITS(source){
- //Copy constructor (dummy)
- if(&source == this) return;
- Warning("Copy Constructor","Not allowed to copy AliITSvSDD03");
- return;
-}
-//______________________________________________________________________
-AliITSvBeamTestITS04& AliITSvBeamTestITS04::operator=(const AliITSvBeamTestITS04 &source){
-
-
- // This class is not to be copied. Function only dummy.
- if(&source == this) return *this;
- Warning("= operator","Not allowed to copy AliITSvSDD03");
- return *this;
-}
-
//__________________________________________________________________
-AliITSvBeamTestITS04::~AliITSvBeamTestITS04()
-{
+AliITSvBeamTestITS04::~AliITSvBeamTestITS04(){
//
// Destructor
//
}
-
//______________________________________________________________________
-void AliITSvBeamTestITS04::CreateMaterials()
-{
+void AliITSvBeamTestITS04::CreateMaterials(){
// Media defined here should correspond to the one defined in galice.cuts
// This file is read in (AliMC*) fMCApp::Init() { ReadTransPar(); }
-
// Create ITS materials
- Int_t ifield = gAlice->Field()->Integ();
- Float_t fieldm = gAlice->Field()->Max();
+ Int_t ifield = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
+ Float_t fieldm = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
Float_t tmaxfdSi = 0.1;
Float_t stemaxSi = 0.0075;
}
//______________________________________________________________________
-void AliITSvBeamTestITS04::CreateGeometry()
-{
- //Creates geometry
+void AliITSvBeamTestITS04::CreateGeometry(){
+ //Creates geometry
+ // These constant character strings are set by cvs during commit
+ // do not change them unless you know what you are doing!
+ const Char_t *cvsDate="$Date$";
+ const Char_t *cvsRevision="$Revision$";
TGeoManager *geoManager = gGeoManager;
TGeoVolume *vALIC = geoManager->GetTopVolume();
sITS->DefineSection(0,-100.0,0.01,100.0); // Units in cms
sITS->DefineSection(1,+100.0,0.01,100.0);
- TGeoMedium *air = gGeoManager->GetMedium("ITSair");
- fITSmotherVolume = new TGeoVolume("vITS",sITS,air);
+ TGeoMedium *air = geoManager->GetMedium("ITSair");
+ fITSmotherVolume = new TGeoVolume("ITSV",sITS,air);
+ const Int_t kLength=100;
+ Char_t vstrng[kLength];
+ if(fIgm.WriteVersionString(vstrng,kLength,(AliITSVersion_t)IsVersion(),
+ fMinorVersion,cvsDate,cvsRevision))
+ fITSmotherVolume->SetTitle(vstrng);
+ else Error("CreateGeometry","Error writing/setting version string");
+ //printf("Title set to %s\n",vstrng);
+ if(vALIC==0) {
+ Error("CreateGeometry","alic=0");
+ return;
+ } // end if
fITSmotherVolume->SetVisibility(kFALSE);
vALIC->AddNode(fITSmotherVolume,1,0);
void AliITSvBeamTestITS04::Init()
{
// Initialize the ITS after it has been created.
- Int_t i;
- for(i=0;i<20;i++) printf("*");
- printf( " ITSbeamtest_Init " );
- for(i=0;i<20;i++) printf("*"); printf("\n");
-
-// // Create geometry
-// if(!fGeomDetIn) this->InitAliITSgeom();
-
- // Initialize AliITS
+ // Inputs:
+ // none.
+ // Outputs:
+ // none.
+ // Return:
+ // none.
+
+ AliDebug(1,Form("Init: Major version %d Minor version %d",fMajorVersion,
+ fMinorVersion));
+ //
+ UpdateInternalGeometry();
AliITS::Init();
- for(i=0;i<40+16;i++) printf("*"); printf("\n");
+ if(fGeomDetOut) GetITSgeom()->WriteNewFile(fWrite);
+ //
}
-
+/*
//______________________________________________________________________
void AliITSvBeamTestITS04::InitAliITSgeom()
{
if (GetITSgeom()) SetITSgeom(0x0);
AliITSgeom* geom = new AliITSgeom(0,knlayers,nlad,ndet,nModTot);
SetITSgeom(geom);
- //*** Set default shapes
+ // *** Set default shapes
const Float_t kDxyzSPD[] = {fgkSPDwidthSens/2, fgkSPDthickSens/2,fgkSPDlengthSens/2};
if(!(GetITSgeom()->IsShapeDefined(kSPD)))
GetITSgeom()->ReSetShape(kSPD,new AliITSgeomSPD425Short(3,(Float_t *)kDxyzSPD));
return;
}
-
//______________________________________________________________________
void AliITSvBeamTestITS04::SetDefaults()
{
bz[160] = 425.0; // Set so that there is no zero pixel size for fNz.
seg0->SetBinSize(bx,bz); // Based on AliITSgeomSPD for now.
SetSegmentationModel(kSPD,seg0);
- // set digit and raw cluster classes to be used
- const char *kData0=(fDetTypeSim->GetCalibrationModel(kSPD))->DataType();
- if (strstr(kData0,"real"))
- fDetTypeSim->SetDigitClassName(kSPD,"AliITSdigit");
- else fDetTypeSim->SetDigitClassName(kSPD,"AliITSdigitSPD");
}
// SDD
seg1->SetDetSize(s1->GetDx()*kconv, // base this on AliITSgeomSDD
s1->GetDz()*4.*kconv, // for now.
s1->GetDy()*4.*kconv); // x,z,y full width in microns.
- seg1->SetDriftSpeed(resp1->GetDriftSpeed());
+ seg1->SetDriftSpeed(AliITSDriftSpeedSDD::DefaultDriftSpeed());
seg1->SetNPads(256,256);// Use AliITSgeomSDD for now
SetSegmentationModel(kSDD,seg1);
- const char *kData1=(fDetTypeSim->GetCalibrationModel(kSDD))->DataType();
- const char *kopt=resp1->GetZeroSuppOption();
- if((!strstr(kopt,"2D")) && (!strstr(kopt,"1D")) || strstr(kData1,"real") ){
- fDetTypeSim->SetDigitClassName(kSDD,"AliITSdigit");
- } else fDetTypeSim->SetDigitClassName(kSDD,"AliITSdigitSDD");
}
// SSD
seg2->SetPadSize(95.,0.); // strip x pitch in microns
seg2->SetNPads(768,0); // number of strips on each side.
seg2->SetAngles(0.0075,0.0275); // strip angels rad P and N side.
- seg2->SetAnglesLay5(0.0075,0.0275); // strip angels rad P and N side.
- seg2->SetAnglesLay6(0.0275,0.0075); // strip angels rad P and N side.
SetSegmentationModel(kSSD,seg2);
- const char *kData2=(fDetTypeSim->GetCalibrationModel(kSSD))->DataType();
- if(strstr(kData2,"real") ) fDetTypeSim->SetDigitClassName(kSSD,"AliITSdigit");
- else fDetTypeSim->SetDigitClassName(kSSD,"AliITSdigitSSD");
}
if(fgkNTYPES>3){Warning("SetDefaults","Only the four basic detector types are initialised!");}
return;
}
-
+*/
//______________________________________________________________________
void AliITSvBeamTestITS04::AddSPDGeometry(TGeoVolume *moth) const
{
// creator with the information to be recoreded about that hit.
// "Standard" StepManager. (Similar to AliITSv11) mln
- Int_t copy, id;
+ Int_t cpy0,mod,status,id,kk;
TLorentzVector position, momentum;
- static TLorentzVector position0;
- static Int_t stat0=0;
-
- if(!(this->IsActive())){
- return;
- } // end if !Active volume.
+ static AliITShit hit;// Saves on calls to construtors
+ if(!(this->IsActive())) return;
if(!(gMC->TrackCharge())) return;
+ //TClonesArray &lhits = *(GetDetTypeSim()->GetHits());
+ TClonesArray &lhits = *(Hits());
+ //
+ // Track status
+ // Track status
+ status = 0;
+ if(gMC->IsTrackInside()) status += 1;
+ if(gMC->IsTrackEntering()) status += 2;
+ if(gMC->IsTrackExiting()) status += 4;
+ if(gMC->IsTrackOut()) status += 8;
+ if(gMC->IsTrackDisappeared()) status += 16;
+ if(gMC->IsTrackStop()) status += 32;
+ if(gMC->IsTrackAlive()) status += 64;
- id=gMC->CurrentVolID(copy);
+ id=gMC->CurrentVolID(cpy0);
Bool_t sensvol = kFALSE;
- for(Int_t kk = 0; kk < fIdN; kk++)
- if(id == fIdSens[kk]) sensvol = kTRUE;
-
- if (sensvol && (gMC->IsTrackExiting())) {
- copy = fTrackReferences->GetEntriesFast();
- TClonesArray &lTR = *fTrackReferences;
- // Fill TrackReference structure with this new TrackReference.
- new(lTR[copy]) AliTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber());
- } // if Outer ITS mother Volume
-
- Int_t vol[5];
- TClonesArray &lhits = *fHits;
+ for(kk=0;kk<fIdN;kk++) if(id == fIdSens[kk]) sensvol = kTRUE;
+ if(!sensvol) return;
+
+ fIgm.DecodeDetector(mod,gGeoManager->GetMother(1)->GetNumber(),1,1,1);
//
- // Track status
- vol[3] = 0;
- vol[4] = 0;
// Fill hit structure.
- if(gMC->IsTrackInside()) vol[3] += 1;
- if(gMC->IsTrackEntering()) vol[3] += 2;
- if(gMC->IsTrackExiting()) vol[3] += 4;
- if(gMC->IsTrackOut()) vol[3] += 8;
- if(gMC->IsTrackDisappeared()) vol[3] += 16;
- if(gMC->IsTrackStop()) vol[3] += 32;
- if(gMC->IsTrackAlive()) vol[3] += 64;
-
- // Only entering charged tracks
- if(!(gMC->TrackCharge())) return;
-
- if( ((id = gMC->CurrentVolID(copy)) == fIdSens[0]) ||
- ((id = gMC->CurrentVolID(copy)) == fIdSens[1]) ||
- ((id = gMC->CurrentVolID(copy)) == fIdSens[2]) )
- {
- GetCurrentLayLaddDet(vol[0], vol[2], vol[1]);
- // vol[2], vol[1]) : in this order because the ladder
- // index and the det. index are exchanged in the constructor
- // of AliITShit...
- } else {
- return; // not an ITS volume?
- };
-
+ //
+ hit.SetModule(mod);
+ hit.SetTrack(gAlice->GetMCApp()->GetCurrentTrackNumber());
gMC->TrackPosition(position);
gMC->TrackMomentum(momentum);
- vol[4] = stat0;
+ hit.SetPosition(position);
+ hit.SetTime(gMC->TrackTime());
+ hit.SetMomentum(momentum);
+ hit.SetStatus(status);
+ hit.SetEdep(gMC->Edep());
+ hit.SetShunt(GetIshunt());
if(gMC->IsTrackEntering()){
- position0 = position;
- stat0 = vol[3];
- return;
+ hit.SetStartPosition(position);
+ hit.SetStartTime(gMC->TrackTime());
+ hit.SetStartStatus(status);
+ return; // don't save entering hit.
} // end if IsEntering
// Fill hit structure with this new hit.
- new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->GetMCApp()->GetCurrentTrackNumber(),
- vol, gMC->Edep(),gMC->TrackTime(),position,
- position0,momentum);
- //
- position0 = position;
- stat0 = vol[3];
+ //Info("StepManager","Calling Copy Constructor");
+ new(lhits[fNhits++]) AliITShit(hit); // Use Copy Construtor.
+ // Save old position... for next hit.
+ hit.SetStartPosition(position);
+ hit.SetStartTime(gMC->TrackTime());
+ hit.SetStartStatus(status);
return;
}
-
+/*
//______________________________________________________________________
Int_t AliITSvBeamTestITS04::GetCurrentLayLaddDet(Int_t &lay,Int_t &ladd, Int_t &det) const
{
return kTRUE;
}
-
+*/
//_____________________________________________________________
Int_t AliITSvBeamTestITS04::GetNumberOfSubDet(const TString& det) const{