[u/mrichter/AliRoot.git] / RICH / AliRICHv0.cxx

/**************************************************************************
 * 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$ */


/////////////////////////////////////////////////////////////
//  Manager and hits classes for set: RICH default version //
/////////////////////////////////////////////////////////////

#include "Riostream.h"

#include <TNode.h> 
#include <TRandom.h> 
#include <TTUBE.h>
#include <TVirtualMC.h>

#include "AliConst.h" 
#include "AliPDG.h" 
#include "AliRICHGeometry.h"
#include "AliRICHResponse.h"
#include "AliRICHResponseV0.h"
#include "AliRICHSegmentationV0.h"
#include "AliRICHv0.h"
#include "AliRun.h"

ClassImp(AliRICHv0)
    
//___________________________________________
AliRICHv0::AliRICHv0() : AliRICH()
{

// Default constructor

    //fChambers = 0;
}

//___________________________________________
AliRICHv0::AliRICHv0(const char *name, const char *title)
    : AliRICH(name,title)
{
    //
// Version 0
// Default Segmentation, no hits
    AliRICHSegmentationV0* segmentation = new AliRICHSegmentationV0;
//
//  Segmentation parameters
    segmentation->SetPadSize(0.84,0.80);
    segmentation->SetDAnod(0.84/2);
//
//  Geometry parameters
    AliRICHGeometry* geometry = new AliRICHGeometry;
    geometry->SetGapThickness(8);
    geometry->SetProximityGapThickness(.4);
    geometry->SetQuartzLength(131);
    geometry->SetQuartzWidth(126.2);
    geometry->SetQuartzThickness(.5);
    geometry->SetOuterFreonLength(131);
    geometry->SetOuterFreonWidth(40.3);
    geometry->SetInnerFreonLength(131);
    geometry->SetInnerFreonWidth(40.3);
    geometry->SetFreonThickness(1.5);
//
//  Response parameters
    AliRICHResponseV0*  response   = new AliRICHResponseV0;
    response->SetSigmaIntegration(5.);
    response->SetChargeSlope(27.);
    response->SetChargeSpread(0.18, 0.18);
    response->SetMaxAdc(4096);
    response->SetAlphaFeedback(0.036);
    response->SetEIonisation(26.e-9);
    response->SetSqrtKx3(0.77459667);
    response->SetKx2(0.962);
    response->SetKx4(0.379);
    response->SetSqrtKy3(0.77459667);
    response->SetKy2(0.962);
    response->SetKy4(0.379);
    response->SetPitch(0.25);
    response->SetWireSag(0);                     // 1->On, 0->Off
    response->SetVoltage(2150);                  // Should only be 2000, 2050, 2100 or 2150
//
//    AliRICH *RICH = (AliRICH *) gAlice->GetDetector("RICH"); 
    
    fCkovNumber=0;
    fFreonProd=0;
    Int_t i=0;
    
    fChambers = new TObjArray(kNCH);
    for (i=0; i<kNCH; i++) {
      
      //PH      (*fChambers)[i] = new AliRICHChamber();  
      fChambers->AddAt(new AliRICHChamber(), i);  
      
    }
  
    for (i=0; i<kNCH; i++) {
      SetGeometryModel(i,geometry);
      SetSegmentationModel(i, segmentation);
      SetResponseModel(i, response);
      SetDebugLevel(0);
    }
}


//___________________________________________

void AliRICHv0::Init()
{

  if(fDebug) {
    printf("%s: *********************************** RICH_INIT ***********************************\n",ClassName());
    printf("%s: *                                                                               *\n",ClassName());
    printf("%s: *                       AliRICHv0 Default version started                       *\n",ClassName());
    printf("%s: *                                                                               *\n",ClassName());
  }

  
  AliSegmentation*  segmentation;
  AliRICHGeometry*  geometry;
  AliRICHResponse*  response;


    // 
    // Initialize Tracking Chambers
    //
    for (Int_t i=0; i<kNCH; i++) {
	//printf ("i:%d",i);
      //PH	( (AliRICHChamber*) (*fChambers)[i])->Init(i);  
	( (AliRICHChamber*) fChambers->At(i))->Init(i);  
    }  
    
    //
    // Set the chamber (sensitive region) GEANT identifier
    
    //PH    ((AliRICHChamber*)(*fChambers)[0])->SetGid(1);  
    //PH    ((AliRICHChamber*)(*fChambers)[1])->SetGid(2);  
    //PH    ((AliRICHChamber*)(*fChambers)[2])->SetGid(3);  
    //PH    ((AliRICHChamber*)(*fChambers)[3])->SetGid(4);  
    //PH    ((AliRICHChamber*)(*fChambers)[4])->SetGid(5);  
    //PH    ((AliRICHChamber*)(*fChambers)[5])->SetGid(6);  
    //PH    ((AliRICHChamber*)(*fChambers)[6])->SetGid(7); 

    ((AliRICHChamber*)fChambers->At(0))->SetGid(1);  
    ((AliRICHChamber*)fChambers->At(1))->SetGid(2);  
    ((AliRICHChamber*)fChambers->At(2))->SetGid(3);  
    ((AliRICHChamber*)fChambers->At(3))->SetGid(4);  
    ((AliRICHChamber*)fChambers->At(4))->SetGid(5);  
    ((AliRICHChamber*)fChambers->At(5))->SetGid(6);  
    ((AliRICHChamber*)fChambers->At(6))->SetGid(7);  

    segmentation=Chamber(0).GetSegmentationModel(0);
    geometry=Chamber(0).GetGeometryModel();
    response=Chamber(0).GetResponseModel();
    
    Float_t offset       = 490 + 1.276 - geometry->GetGapThickness()/2;        //distance from center of mother volume to methane
    Float_t deltaphi     = 19.5;                                               //phi angle between center of chambers - z direction
    Float_t deltatheta   = 20;                                                 //theta angle between center of chambers - x direction
    Float_t cosphi       = TMath::Cos(deltaphi*TMath::Pi()/180);
    Float_t sinphi       = TMath::Sin(deltaphi*TMath::Pi()/180);
    Float_t costheta     = TMath::Cos(deltatheta*TMath::Pi()/180);
    Float_t sintheta     = TMath::Sin(deltatheta*TMath::Pi()/180);

    Float_t pos1[3]={0.                , offset*cosphi         , offset*sinphi};
    Float_t pos2[3]={offset*sintheta   , offset*costheta       , 0. };
    Float_t pos3[3]={0.                , offset                , 0.};
    Float_t pos4[3]={-offset*sintheta  , offset*costheta       , 0.};
    Float_t pos5[3]={offset*sinphi     , offset*costheta*cosphi, -offset*sinphi};
    Float_t pos6[3]={0.                , offset*cosphi         , -offset*sinphi};
    Float_t pos7[3]={ -offset*sinphi   , offset*costheta*cosphi, -offset*sinphi};

    Chamber(0).SetChamberTransform(pos1[0],pos1[1],pos1[2],new TRotMatrix("rot993","rot993",90., 0.               , 90. - deltaphi, 90.             , deltaphi, -90.           ));
    Chamber(1).SetChamberTransform(pos2[0],pos2[1],pos2[2],new TRotMatrix("rot994","rot994",90., -deltatheta      , 90.           , 90.- deltatheta , 0.      , 0.             ));
    Chamber(2).SetChamberTransform(pos3[0],pos3[1],pos3[2],new TRotMatrix("rot995","rot995",90., 0.               , 90.           , 90.             , 0.      , 0.             ));
    Chamber(3).SetChamberTransform(pos4[0],pos4[1],pos4[2],new TRotMatrix("rot996","rot996",90.,  deltatheta      , 90.           , 90 + deltatheta , 0.      , 0.             ));
    Chamber(4).SetChamberTransform(pos5[0],pos5[1],pos5[2],new TRotMatrix("rot997","rot997",90., 360. - deltatheta, 108.2         , 90.- deltatheta ,18.2     , 90 - deltatheta));
    Chamber(5).SetChamberTransform(pos6[0],pos6[1],pos6[2],new TRotMatrix("rot998","rot998",90., 0.               , 90 + deltaphi , 90.             , deltaphi, 90.            ));
    Chamber(6).SetChamberTransform(pos7[0],pos7[1],pos7[2],new TRotMatrix("rot999","rot999",90., deltatheta       , 108.2         , 90.+ deltatheta ,18.2     , 90 + deltatheta));   
     
    if(fDebug) {
      printf("%s: *                            Pads            : %3dx%3d                          *\n",
	     ClassName(),segmentation->Npx(),segmentation->Npy());
      printf("%s: *                            Pad size        : %5.2f x%5.2f mm2                 *\n",
	     ClassName(),segmentation->Dpx(),segmentation->Dpy()); 
      printf("%s: *                            Gap Thickness   : %5.1f cm                         *\n",
	     ClassName(),geometry->GetGapThickness());
      printf("%s: *                            Radiator Width  : %5.1f cm                         *\n",
	     ClassName(),geometry->GetQuartzWidth());
      printf("%s: *                            Radiator Length : %5.1f cm                         *\n",
	     ClassName(),geometry->GetQuartzLength());
      printf("%s: *                            Freon Thickness : %5.1f cm                         *\n",
	     ClassName(),geometry->GetFreonThickness());
      printf("%s: *                            Charge Slope    : %5.1f ADC                        *\n",
	     ClassName(),response->ChargeSlope());
      printf("%s: *                            Feedback Prob.  : %5.2f %%                          *\n",
	     ClassName(),response->AlphaFeedback()*100);
      printf("%s: *                                                                               *\n",
	     ClassName());
      printf("%s: *********************************************************************************\n",
	     ClassName());
    }
}
Commit	Line	Data
4c039060	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
88cb7938	16	/* $Id$ */
6e36c0f2	17
6e36c0f2	18
c71fc6da	19	/////////////////////////////////////////////////////////////
	20	// Manager and hits classes for set: RICH default version //
	21	/////////////////////////////////////////////////////////////
ddae0931	22
88cb7938	23	#include "Riostream.h"
88cb7938	24
ddae0931	25	#include <TNode.h>
ddae0931	26	#include <TRandom.h>
88cb7938	27	#include <TTUBE.h>
88cb7938	28	#include <TVirtualMC.h>
ddae0931	29
88cb7938	30	#include "AliConst.h"
	31	#include "AliPDG.h"
	32	#include "AliRICHGeometry.h"
237c933d	33	#include "AliRICHResponse.h"
237c933d	34	#include "AliRICHResponseV0.h"
88cb7938	35	#include "AliRICHSegmentationV0.h"
88cb7938	36	#include "AliRICHv0.h"
ddae0931	37	#include "AliRun.h"
ddae0931	38
	39	ClassImp(AliRICHv0)
	40
	41	//___________________________________________
	42	AliRICHv0::AliRICHv0() : AliRICH()
	43	{
237c933d	44
	45	// Default constructor
	46
6e36c0f2	47	//fChambers = 0;
ddae0931	48	}
	49
	50	//___________________________________________
	51	AliRICHv0::AliRICHv0(const char name, const char title)
	52	: AliRICH(name,title)
	53	{
c71fc6da	54	//
c71fc6da	55	// Version 0
237c933d	56	// Default Segmentation, no hits
c4384528	57	AliRICHSegmentationV0* segmentation = new AliRICHSegmentationV0;
c71fc6da	58	//
c71fc6da	59	// Segmentation parameters
c4384528	60	segmentation->SetPadSize(0.84,0.80);
c4384528	61	segmentation->SetDAnod(0.84/2);
c71fc6da	62	//
c71fc6da	63	// Geometry parameters
237c933d	64	AliRICHGeometry* geometry = new AliRICHGeometry;
	65	geometry->SetGapThickness(8);
	66	geometry->SetProximityGapThickness(.4);
	67	geometry->SetQuartzLength(131);
	68	geometry->SetQuartzWidth(126.2);
	69	geometry->SetQuartzThickness(.5);
	70	geometry->SetOuterFreonLength(131);
	71	geometry->SetOuterFreonWidth(40.3);
	72	geometry->SetInnerFreonLength(131);
	73	geometry->SetInnerFreonWidth(40.3);
a2db5cde	74	geometry->SetFreonThickness(1.5);
c71fc6da	75	//
c71fc6da	76	// Response parameters
c4384528	77	AliRICHResponseV0* response = new AliRICHResponseV0;
	78	response->SetSigmaIntegration(5.);
	79	response->SetChargeSlope(27.);
	80	response->SetChargeSpread(0.18, 0.18);
	81	response->SetMaxAdc(4096);
	82	response->SetAlphaFeedback(0.036);
	83	response->SetEIonisation(26.e-9);
	84	response->SetSqrtKx3(0.77459667);
	85	response->SetKx2(0.962);
	86	response->SetKx4(0.379);
	87	response->SetSqrtKy3(0.77459667);
	88	response->SetKy2(0.962);
	89	response->SetKy4(0.379);
	90	response->SetPitch(0.25);
	91	response->SetWireSag(0); // 1->On, 0->Off
	92	response->SetVoltage(2150); // Should only be 2000, 2050, 2100 or 2150
c71fc6da	93	//
	94	// AliRICH RICH = (AliRICH ) gAlice->GetDetector("RICH");
	95
237c933d	96	fCkovNumber=0;
237c933d	97	fFreonProd=0;
7f36f6bc	98	Int_t i=0;
c71fc6da	99
237c933d	100	fChambers = new TObjArray(kNCH);
7f36f6bc	101	for (i=0; i<kNCH; i++) {
c71fc6da	102
2682e810	103	//PH (*fChambers)[i] = new AliRICHChamber();
2682e810	104	fChambers->AddAt(new AliRICHChamber(), i);
c71fc6da	105
	106	}
	107
7f36f6bc	108	for (i=0; i<kNCH; i++) {
237c933d	109	SetGeometryModel(i,geometry);
c4384528	110	SetSegmentationModel(i, segmentation);
c4384528	111	SetResponseModel(i, response);
8fa8fc03	112	SetDebugLevel(0);
6e36c0f2	113	}
ddae0931	114	}
	115
	116
6e36c0f2	117
ddae0931	118	//___________________________________________
	119
	120	void AliRICHv0::Init()
	121	{
c71fc6da	122
9e1a0ddb	123	if(fDebug) {
	124	printf("%s: ********************************* RICH_INIT *********************************\n",ClassName());
	125	printf("%s: * *\n",ClassName());
	126	printf("%s: * AliRICHv0 Default version started *\n",ClassName());
	127	printf("%s: * *\n",ClassName());
	128	}
c71fc6da	129
c71fc6da	130
a2f7eaf6	131	AliSegmentation* segmentation;
c71fc6da	132	AliRICHGeometry* geometry;
	133	AliRICHResponse* response;
	134
	135
ddae0931	136	//
	137	// Initialize Tracking Chambers
	138	//
2fc6393f	139	for (Int_t i=0; i<kNCH; i++) {
6e36c0f2	140	//printf ("i:%d",i);
2682e810	141	//PH ( (AliRICHChamber) (fChambers)[i])->Init(i);
2682e810	142	( (AliRICHChamber*) fChambers->At(i))->Init(i);
6e36c0f2	143	}
ddae0931	144
	145	//
	146	// Set the chamber (sensitive region) GEANT identifier
	147
2682e810	148	//PH ((AliRICHChamber)(fChambers)[0])->SetGid(1);
	149	//PH ((AliRICHChamber)(fChambers)[1])->SetGid(2);
	150	//PH ((AliRICHChamber)(fChambers)[2])->SetGid(3);
	151	//PH ((AliRICHChamber)(fChambers)[3])->SetGid(4);
	152	//PH ((AliRICHChamber)(fChambers)[4])->SetGid(5);
	153	//PH ((AliRICHChamber)(fChambers)[5])->SetGid(6);
	154	//PH ((AliRICHChamber)(fChambers)[6])->SetGid(7);
	155
	156	((AliRICHChamber*)fChambers->At(0))->SetGid(1);
	157	((AliRICHChamber*)fChambers->At(1))->SetGid(2);
	158	((AliRICHChamber*)fChambers->At(2))->SetGid(3);
	159	((AliRICHChamber*)fChambers->At(3))->SetGid(4);
	160	((AliRICHChamber*)fChambers->At(4))->SetGid(5);
	161	((AliRICHChamber*)fChambers->At(5))->SetGid(6);
	162	((AliRICHChamber*)fChambers->At(6))->SetGid(7);
6e36c0f2	163
c71fc6da	164	segmentation=Chamber(0).GetSegmentationModel(0);
	165	geometry=Chamber(0).GetGeometryModel();
	166	response=Chamber(0).GetResponseModel();
	167
449e6743	168	Float_t offset = 490 + 1.276 - geometry->GetGapThickness()/2; //distance from center of mother volume to methane
	169	Float_t deltaphi = 19.5; //phi angle between center of chambers - z direction
	170	Float_t deltatheta = 20; //theta angle between center of chambers - x direction
	171	Float_t cosphi = TMath::Cos(deltaphi*TMath::Pi()/180);
	172	Float_t sinphi = TMath::Sin(deltaphi*TMath::Pi()/180);
	173	Float_t costheta = TMath::Cos(deltatheta*TMath::Pi()/180);
	174	Float_t sintheta = TMath::Sin(deltatheta*TMath::Pi()/180);
	175
	176	Float_t pos1[3]={0. , offsetcosphi , offsetsinphi};
	177	Float_t pos2[3]={offsetsintheta , offsetcostheta , 0. };
	178	Float_t pos3[3]={0. , offset , 0.};
	179	Float_t pos4[3]={-offsetsintheta , offsetcostheta , 0.};
	180	Float_t pos5[3]={offsetsinphi , offsetcosthetacosphi, -offsetsinphi};
	181	Float_t pos6[3]={0. , offsetcosphi , -offsetsinphi};
	182	Float_t pos7[3]={ -offsetsinphi , offsetcosthetacosphi, -offsetsinphi};
	183
	184	Chamber(0).SetChamberTransform(pos1[0],pos1[1],pos1[2],new TRotMatrix("rot993","rot993",90., 0. , 90. - deltaphi, 90. , deltaphi, -90. ));
	185	Chamber(1).SetChamberTransform(pos2[0],pos2[1],pos2[2],new TRotMatrix("rot994","rot994",90., -deltatheta , 90. , 90.- deltatheta , 0. , 0. ));
	186	Chamber(2).SetChamberTransform(pos3[0],pos3[1],pos3[2],new TRotMatrix("rot995","rot995",90., 0. , 90. , 90. , 0. , 0. ));
	187	Chamber(3).SetChamberTransform(pos4[0],pos4[1],pos4[2],new TRotMatrix("rot996","rot996",90., deltatheta , 90. , 90 + deltatheta , 0. , 0. ));
	188	Chamber(4).SetChamberTransform(pos5[0],pos5[1],pos5[2],new TRotMatrix("rot997","rot997",90., 360. - deltatheta, 108.2 , 90.- deltatheta ,18.2 , 90 - deltatheta));
	189	Chamber(5).SetChamberTransform(pos6[0],pos6[1],pos6[2],new TRotMatrix("rot998","rot998",90., 0. , 90 + deltaphi , 90. , deltaphi, 90. ));
	190	Chamber(6).SetChamberTransform(pos7[0],pos7[1],pos7[2],new TRotMatrix("rot999","rot999",90., deltatheta , 108.2 , 90.+ deltatheta ,18.2 , 90 + deltatheta));
c71fc6da	191
9e1a0ddb	192	if(fDebug) {
	193	printf("%s: * Pads : %3dx%3d *\n",
	194	ClassName(),segmentation->Npx(),segmentation->Npy());
	195	printf("%s: * Pad size : %5.2f x%5.2f mm2 *\n",
	196	ClassName(),segmentation->Dpx(),segmentation->Dpy());
	197	printf("%s: * Gap Thickness : %5.1f cm *\n",
	198	ClassName(),geometry->GetGapThickness());
	199	printf("%s: * Radiator Width : %5.1f cm *\n",
	200	ClassName(),geometry->GetQuartzWidth());
	201	printf("%s: * Radiator Length : %5.1f cm *\n",
	202	ClassName(),geometry->GetQuartzLength());
	203	printf("%s: * Freon Thickness : %5.1f cm *\n",
	204	ClassName(),geometry->GetFreonThickness());
	205	printf("%s: * Charge Slope : %5.1f ADC *\n",
	206	ClassName(),response->ChargeSlope());
	207	printf("%s: * Feedback Prob. : %5.2f %% *\n",
	208	ClassName(),response->AlphaFeedback()*100);
	209	printf("%s: * *\n",
	210	ClassName());
	211	printf("%s: *********************************************************************************\n",
	212	ClassName());
	213	}
ddae0931	214	}
ddae0931	215