/**************************************************************************
 * 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.11  2000/06/27 07:31:07  morsch
fChambers = 0; deleted from constructor.

Revision 1.10  2000/06/26 14:02:38  morsch
Add class AliMUONConstants with MUON specific constants using static memeber data and access methods.

Revision 1.9  2000/06/15 07:58:49  morsch
Code from MUON-dev joined

Revision 1.8.4.9  2000/06/12 19:20:49  morsch
Constructor sets default geometry, segmentation and response parameters.

Revision 1.8.4.8  2000/06/09 21:55:28  morsch
Most coding rule violations corrected.

Revision 1.8.4.7  2000/05/02 13:15:18  morsch
Coding rule violations RS3, RN13 corected

Revision 1.8.4.6  2000/05/02 10:24:26  morsch
Public access to fdGas and fdAlu of AliMUONChamber replaced by getters.

Revision 1.8.4.5  2000/04/26 19:58:47  morsch
Obsolete reference to trig_ removed.

Revision 1.8.4.4  2000/04/19 19:42:47  morsch
change NCH to kNCH

Revision 1.8.4.3  2000/02/17 08:17:43  morsch
Gammas and neutrons are also scored in the stepmanager
*/

/////////////////////////////////////////////////////////
//  Manager and hits classes for set:MUON version 0    //
/////////////////////////////////////////////////////////

#include <TTUBE.h>
#include <TNode.h> 
#include <TRandom.h> 
#include <TLorentzVector.h> 
#include <iostream.h>

#include "AliMUONv0.h"
#include "AliMUONChamber.h"
#include "AliRun.h"
#include "AliMC.h"
#include "AliMagF.h"
#include "AliMUONHit.h"
#include "AliMUONPadHit.h"
#include "AliCallf77.h"
#include "AliConst.h" 
#include "AliMUONResponseV0.h"
#include "AliMUONResponseTrigger.h"
#include "AliMUONSegmentationV0.h"
#include "AliMUONSegmentationV01.h"
#include "AliMUONSegmentationV02.h"
#include "AliMUONSegmentationV04.h"
#include "AliMUONSegmentationV05.h"
#include "AliMUONSegmentationTrigger.h"
#include "AliMUONSegmentationTriggerX.h"
#include "AliMUONSegmentationTriggerY.h"
#include "AliMUONConstants.h"

ClassImp(AliMUONv0)
 
//___________________________________________
AliMUONv0::AliMUONv0() : AliMUON()
{
// Constructor
    fChambers = 0;
}
 
//___________________________________________
AliMUONv0::AliMUONv0(const char *name, const char *title)
       : AliMUON(name,title)
{
// Constructor
    SetIshunt(0);
    SetMaxStepGas(0.1);
    SetMaxStepAlu(0.1);
//
// Version 0
//
// First define the number of planes that are segmented (1 or 2) by a call
// to SetNsec. 
// Then chose for each chamber (chamber plane) the segmentation 
// and response model.
// They should be equal for the two chambers of each station. In a future
// version this will be enforced.
//
//  
    Int_t chamber;
    Int_t station;
// Default response
    AliMUONResponseV0* response0 = new AliMUONResponseV0;
    response0->SetSqrtKx3(0.7131);
    response0->SetKx2(1.0107);
    response0->SetKx4(0.4036);
    response0->SetSqrtKy3(0.7642);
    response0->SetKy2(0.9706);
    response0->SetKy4(0.3831);
    response0->SetPitch(0.25);
    response0->SetSigmaIntegration(10.);
    response0->SetChargeSlope(50);
    response0->SetChargeSpread(0.18, 0.18);
    response0->SetMaxAdc(4096);
    response0->SetZeroSuppression(6);
//--------------------------------------------------------
// Configuration for Chamber TC1/2  (Station 1) ----------           
//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    Float_t rseg1[4]={17.5, 55.2, 71.3, 95.5};
    Int_t   nseg1[4]={4, 4, 2, 1};
//
    chamber=1;
//^^^^^^^^^
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV01 *seg11=new AliMUONSegmentationV01;
    
    seg11->SetSegRadii(rseg1);
    seg11->SetPadSize(3, 0.5);
    seg11->SetDAnod(3.0/3./4);
    seg11->SetPadDivision(nseg1);
    
    SetSegmentationModel(chamber-1, 1, seg11);
//
    AliMUONSegmentationV02 *seg12=new AliMUONSegmentationV02;
    seg12->SetSegRadii(rseg1); 
    seg12->SetPadSize(0.75, 2.0);
    seg12->SetDAnod(3.0/3./4);
    seg12->SetPadDivision(nseg1);
    
    SetSegmentationModel(chamber-1, 2, seg12);
    
    SetResponseModel(chamber-1, response0);	    
    
    chamber=2;
//^^^^^^^^^
//
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV01 *seg21=new AliMUONSegmentationV01;
    seg21->SetSegRadii(rseg1);
    seg21->SetPadSize(3, 0.5);
    seg21->SetDAnod(3.0/3./4);
    seg21->SetPadDivision(nseg1);
    SetSegmentationModel(chamber-1, 1, seg21);
//
    AliMUONSegmentationV02 *seg22=new AliMUONSegmentationV02;
    seg22->SetSegRadii(rseg1); 
    seg22->SetPadSize(0.75, 2.);
    seg22->SetDAnod(3.0/3./4);
    seg22->SetPadDivision(nseg1);
    SetSegmentationModel(chamber-1, 2, seg22);
    
    SetResponseModel(chamber-1, response0);	    
//
//--------------------------------------------------------
// Configuration for Chamber TC3/4 -----------------------
//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    Float_t rseg2[4]={23.5, 47.1, 87.7, 122.5};
    Int_t   nseg2[4]={4, 4, 2, 1};
//
    chamber=3;
//^^^^^^^^^
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV01 *seg31=new AliMUONSegmentationV01;
    seg31->SetSegRadii(rseg2);
    seg31->SetPadSize(6, 0.5);
    seg31->SetDAnod(3.0/3./4);
    seg31->SetPadDivision(nseg2);
    SetSegmentationModel(chamber-1, 1, seg31);
//
    AliMUONSegmentationV02 *seg32=new AliMUONSegmentationV02;
    seg32->SetSegRadii(rseg2); 
    seg32->SetPadSize(0.75, 4.);
    seg32->SetPadDivision(nseg2);
    seg32->SetDAnod(3.0/3./4);
    
    SetSegmentationModel(chamber-1, 2, seg32);
    
    SetResponseModel(chamber-1, response0);	    
    
    chamber=4;
//^^^^^^^^^
//
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV01 *seg41=new AliMUONSegmentationV01;
    seg41->SetSegRadii(rseg2);
    seg41->SetPadSize(6, 0.5);
    seg41->SetDAnod(3.0/3./4);
    seg41->SetPadDivision(nseg2);
    SetSegmentationModel(chamber-1, 1, seg41);
//
    AliMUONSegmentationV02 *seg42=new AliMUONSegmentationV02;
    seg42->SetSegRadii(rseg2); 
    seg42->SetPadSize(0.75, 4.);
    seg42->SetPadDivision(nseg2);
    seg42->SetDAnod(3.0/3./4);
    
    SetSegmentationModel(chamber-1, 2, seg42);
    
    SetResponseModel(chamber-1, response0);	    


//--------------------------------------------------------
// Configuration for Chamber TC5/6 -----------------------
//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    chamber=5;
//^^^^^^^^^
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV01 *seg51=new AliMUONSegmentationV01;
    seg51->SetSegRadii(rseg2);
    seg51->SetPadSize(6, 0.5);
    seg51->SetDAnod(3.0/3./4);
    seg51->SetPadDivision(nseg2);
    SetSegmentationModel(chamber-1, 1, seg51);
//
    AliMUONSegmentationV02 *seg52=new AliMUONSegmentationV02;
    seg52->SetSegRadii(rseg2); 
    seg52->SetPadSize(0.75, 4.);
    seg52->SetPadDivision(nseg2);
    seg52->SetDAnod(3.0/3./4);
    
    SetSegmentationModel(chamber-1, 2, seg52);
    SetResponseModel(chamber-1, response0);	    
    
    chamber=6;
//^^^^^^^^^
//
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV01 *seg61=new AliMUONSegmentationV01;
    seg61->SetSegRadii(rseg2);
    seg61->SetPadSize(6, 0.5);
    seg61->SetDAnod(3.0/3./4);
    seg61->SetPadDivision(nseg2);
    SetSegmentationModel(chamber-1, 1, seg61);
//
    AliMUONSegmentationV02 *seg62=new AliMUONSegmentationV02;
    seg62->SetSegRadii(rseg2); 
    seg62->SetPadSize(0.75, 4.);
    seg62->SetPadDivision(nseg2);
    seg62->SetDAnod(3.0/3./4);
    
    SetSegmentationModel(chamber-1, 2, seg62);
    
    SetResponseModel(chamber-1, response0);	    


//
// Station 3
    station=3;
    SetPadSize(station, 1, 0.975, 0.55);

//--------------------------------------------------------
// Configuration for Chamber TC7/8  (Station 4) ----------           
//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

    Int_t   nseg4[4]={4, 4, 2, 1};

    chamber=7;
//^^^^^^^^^
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV04 *seg71=new AliMUONSegmentationV04;
    seg71->SetPadSize(10.,0.5);
    seg71->SetDAnod(0.25);
    seg71->SetPadDivision(nseg4);
    SetSegmentationModel(chamber-1, 1, seg71);
    AliMUONSegmentationV05 *seg72=new AliMUONSegmentationV05;
    seg72->SetPadSize(1,10);
    seg72->SetDAnod(0.25);
    seg72->SetPadDivision(nseg4);
    SetSegmentationModel(chamber-1, 2, seg72);
    
    SetResponseModel(chamber-1, response0);	    

    chamber=8;
//^^^^^^^^^
    SetNsec(chamber-1,2);
    AliMUONSegmentationV04 *seg81=new AliMUONSegmentationV04;
    seg81->SetPadSize(10., 0.5);
    seg81->SetPadDivision(nseg4);
    seg81->SetDAnod(0.25);
    SetSegmentationModel(chamber-1, 1, seg81);
    
    AliMUONSegmentationV05 *seg82=new AliMUONSegmentationV05;
    seg82->SetPadSize(1, 10);
    seg82->SetPadDivision(nseg4);
    seg82->SetDAnod(0.25);
    SetSegmentationModel(chamber-1, 2, seg82);
    
    SetResponseModel(chamber-1, response0);	    
//--------------------------------------------------------
// Configuration for Chamber TC9/10  (Station 5) ---------           
//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    chamber=9;
//^^^^^^^^^
    SetNsec(chamber-1,2);
//
    AliMUONSegmentationV04 *seg91=new AliMUONSegmentationV04;
    seg91->SetPadSize(10.,0.5);
    seg91->SetDAnod(0.25);
    seg91->SetPadDivision(nseg4);
    SetSegmentationModel(chamber-1, 1, seg91);
    
    AliMUONSegmentationV05 *seg92=new AliMUONSegmentationV05;
    seg92->SetPadSize(1,10);
    seg92->SetDAnod(0.25);
    seg92->SetPadDivision(nseg4);
    
    SetSegmentationModel(chamber-1, 2, seg92);
    
    SetResponseModel(chamber-1, response0);	    
    
    chamber=10;
//^^^^^^^^^
    SetNsec(chamber-1,2);
    AliMUONSegmentationV04 *seg101=new AliMUONSegmentationV04;
    seg101->SetPadSize(10., 0.5);
    seg101->SetPadDivision(nseg4);
    seg101->SetDAnod(0.25);
    SetSegmentationModel(chamber-1, 1, seg101);
    
    AliMUONSegmentationV05 *seg102=new AliMUONSegmentationV05;
    seg102->SetPadSize(1,10);
    seg102->SetPadDivision(nseg4);
    seg102->SetDAnod(0.25);
    SetSegmentationModel(chamber-1, 2, seg102);
    
    SetResponseModel(chamber-1, response0);	    

//--------------------------------------------------------
// Configuration for Trigger staions --------------------- 
//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
    AliMUONResponseTrigger* responseTrigger0 =  new AliMUONResponseTrigger;
    
    chamber=11;
    SetNsec(chamber-1,2);
    AliMUONSegmentationTriggerX *seg111=new AliMUONSegmentationTriggerX;
    SetSegmentationModel(chamber-1, 1, seg111);
    AliMUONSegmentationTriggerY *seg112=new AliMUONSegmentationTriggerY;
    SetSegmentationModel(chamber-1, 2, seg112);
    
    SetResponseModel(chamber-1, responseTrigger0);      
    
    chamber=12;
    SetNsec(chamber-1,2);
    AliMUONSegmentationTriggerX *seg121=new AliMUONSegmentationTriggerX;
    SetSegmentationModel(chamber-1, 1, seg121);
    AliMUONSegmentationTriggerY *seg122=new AliMUONSegmentationTriggerY;
    SetSegmentationModel(chamber-1, 2, seg122);
    
    SetResponseModel(chamber-1, responseTrigger0);      
    
    chamber=13;
    SetNsec(chamber-1,2);
    AliMUONSegmentationTriggerX *seg131=new AliMUONSegmentationTriggerX;
    SetSegmentationModel(chamber-1, 1, seg131);
    AliMUONSegmentationTriggerY *seg132=new AliMUONSegmentationTriggerY;
    SetSegmentationModel(chamber-1, 2, seg132);
    SetResponseModel(chamber-1, responseTrigger0);      
    
    chamber=14;
    SetNsec(chamber-1,2);
    AliMUONSegmentationTriggerX *seg141=new AliMUONSegmentationTriggerX;
    SetSegmentationModel(chamber-1, 1, seg141);
    AliMUONSegmentationTriggerY *seg142=new AliMUONSegmentationTriggerY;
    SetSegmentationModel(chamber-1, 2, seg142);
    
    SetResponseModel(chamber-1, responseTrigger0); 
}

void AliMUONv0::CreateGeometry()
{
// Creates coarse geometry for hit density simulations
    Int_t *idtmed = fIdtmed->GetArray()-1099;
//
     Float_t zpos, dAlu, tpar[3];
     Int_t idAir=idtmed[1100];
     Int_t idAlu=idtmed[1103];     

     AliMUONChamber *iChamber;
     // Loop over all chambers (tracking and trigger)
     for (Int_t ch = 0; ch < AliMUONConstants::NCh(); ch++) {
	 char alu[8];
	 char gas[8];
     
	 iChamber=(AliMUONChamber*) (*fChambers)[ch];
	 // Z of the chamber
	 zpos=iChamber->Z(); 
	 dAlu=iChamber->DAlu();
	 if (ch < AliMUONConstants::NTrackingCh()) {
	   // tracking chambers
	     sprintf(alu,"CA0%1d",ch);
	     sprintf(gas,"CG0%1d",ch);	 
	 } else {
	   // trigger chambers
	     sprintf(alu,"CA%2d",ch);
	     sprintf(gas,"CG%2d",ch);	 
	 }
//
	 printf("\n %d,  %s,  %s \n ", ch, alu, gas);
	 
	 tpar[0] = iChamber->RInner(); 
	 tpar[1] = iChamber->ROuter();
	 tpar[2] = (dAlu+0.2)/2.;
	 if (ch !=4 && ch !=5) {
	     gMC->Gsvolu(alu, "TUBE", idAlu, tpar, 3);
	     tpar[2] = 0.1;
	     gMC->Gsvolu(gas, "TUBE", idAir, tpar, 3);
	 } else {
	     gMC->Gsvolu(alu, "TUBE", idAlu, tpar, 3);
	     tpar[2] = 0.1;
	     gMC->Gsvolu(gas, "TUBE", idAir, tpar, 3);
	 }
	 gMC->Gspos(gas, 1, alu,  0., 0., 0., 0, "ONLY");
	 gMC->Gspos(alu, 1, "ALIC", 0., 0., zpos, 0, "ONLY");
	 iChamber->SetGid(gMC->VolId(gas));
     }
}

//___________________________________________
void AliMUONv0::CreateMaterials()
{
// Creates materials for coarse geometry
    AliMaterial(15, "AIR$      ", 14.61, 7.3, .001205, 30423.24, 67500);
    AliMaterial(9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);

    Float_t epsil  = .001; // Tracking precision, 
    Float_t stemax = -1.;  // Maximum displacement for multiple scat 
    Float_t tmaxfd = -20.; // Maximum angle due to field deflection 
    Float_t deemax = -.3;  // Maximum fractional energy loss, DLS 
    Float_t stmin  = -.8;
    Int_t isxfld   = gAlice->Field()->Integ();
    Float_t sxmgmx = gAlice->Field()->Max();

    //
    //    Air 
    AliMedium(1, "AIR_CH_US         ", 15, 1, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
    AliMedium(4, "ALU_CH_US          ", 9, 0, isxfld, sxmgmx, tmaxfd, fMaxStepAlu, 
	    fMaxDestepAlu, epsil, stmin);

}

void AliMUONv0::Init()
{
   // 
   // Initialize Tracking Chambers
   //

   printf("\n\n\n Start Init for version 0 - CPC chamber type\n\n\n");
   for (Int_t i=0; i<AliMUONConstants::NCh(); i++) {
       ( (AliMUONChamber*) (*fChambers)[i])->Init();
   }
}

void AliMUONv0::StepManager()
{
//
// Step manager for hit density simulations
  Int_t          copy, id;
  static Int_t   idvol;
  static Int_t   vol[2];
  Int_t          ipart;
  TLorentzVector pos;
  TLorentzVector mom;
  Float_t        theta,phi;
  
  //  modifs perso
  static Float_t hits[15];

  TClonesArray &lhits = *fHits;
  //
  // Only gas gap inside chamber
  // Tag chambers and record hits when track enters 
  idvol=-1;
  id=gMC->CurrentVolID(copy);
  
    for (Int_t i=1; i<=AliMUONConstants::NCh(); i++) {
      if(id==((AliMUONChamber*)(*fChambers)[i-1])->GetGid()){ 
	  vol[0]=i; 
	  idvol=i-1;
      }
    }
    if (idvol == -1) return;
  //
  // Get current particle id (ipart), track position (pos)  and momentum (mom) 
  gMC->TrackPosition(pos);
  gMC->TrackMomentum(mom);

  ipart  = gMC->TrackPid();
  //
  // record hits when track enters ...
//  if( !(gMC->TrackCharge()) ) return; 
  if( gMC->IsTrackEntering()) {
      Double_t tc = mom[0]*mom[0]+mom[1]*mom[1];
      Double_t rt = TMath::Sqrt(tc);
      theta   = Float_t(TMath::ATan2(rt,Double_t(mom[2])))*kRaddeg;
      phi     = Float_t(TMath::ATan2(Double_t(mom[1]),Double_t(mom[0])))*kRaddeg;
      hits[0] = Float_t(ipart);             // Geant3 particle type
      hits[1] = pos[0];                     // X-position for hit
      hits[2] = pos[1];                     // Y-position for hit
      hits[3] = pos[2];                     // Z-position for hit
      hits[4] = theta;                      // theta angle of incidence
      hits[5] = phi;                        // phi angle of incidence 
      hits[8] = -1;                         // first padhit
      hits[9] = -1;                         // last pad hit

      // modifs personel
      hits[10] = mom[3]; // hit Energy
      hits[11] = mom[0]; // Px
      hits[12] = mom[1]; // Py
      hits[13] = mom[2]; // Pz
      hits[14] = gMC->TrackTime();
      
      // fin modifs perso
      new(lhits[fNhits++]) 
	  AliMUONHit(fIshunt,gAlice->CurrentTrack(),vol,hits);

  }
}