KVBreakUp.cpp

//Created by KVClassFactory on Thu Mar 25 11:25:27 2010
//Author: bonnet
 
#include "KVBreakUp.h"
#include "TMath.h"
#include "TMethodCall.h"
#include "TRandom3.h"
#include "KVNucleusEvent.h"
 
using namespace std;
 
ClassImp(KVBreakUp)
 
 
 
 
 
void KVBreakUp::init(void)
{
   //Initialisation des variables
   //appele par le constructeur
 
   Ztotal = 0;
   bound = 0;
   SetConditions(0, 0, 0);
 
   //generateurs de nombre aleatoire
   alea = 0;
   RedefineTRandom("TRandom3");
 
   BreakUpMethod = "BreakUsingChain";
   SetBreakUpMethod(BreakUpMethod);
 
   lhisto = new KVHashList();
   lhisto->SetOwner(kTRUE);
   //Rempli dans les methodes type Break...()
   hzz = new TH1F("KVBreakUp_hzz", "distri z", 200, -0.5, 199.5);
   lhisto->Add(hzz);
   nraffine = 0;
 
   lobjects = new KVHashList();
   lobjects->SetOwner(kTRUE);
 
   DefineHistos();
 
   niter_tot = 0;
   tstart = tstop = tellapsed = 0;
 
   StorePartitions();
   current_event = KVEvent::Factory("SimEvent");
}
 
 
 
 
KVBreakUp::KVBreakUp(Int_t taille_max)
{
   //Constructeur, l'argument taille_max correspond
   //a la taille maximale du tableau ou sont enregistrees
   //les produits de la cassure (par defaut 5000)
   //Ce qui correspond a la cassure en 5000 morceaux max
   size_max = taille_max;
   size = new Int_t[size_max];
   init();
}
 
 
 
 
KVBreakUp::~KVBreakUp()
{
   //Destructeur
   delete alea;
   delete lhisto;
   delete lobjects;
 
   //delete current_event;
 
   if (bound) delete [] bound;
   delete [] size;
 
   TCanvas* c1 = 0;
   if ((c1 = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("BreakUp_Control"))) delete c1;
 
}
 
 
 
void KVBreakUp::Clear(Option_t* opt)
{
   //Remet a zero les compteurs
   //Efface les partitions enregistrees
   //Le contenu des histos
   KVPartitionList::Clear(opt);
   ResetTotalIterations();
   ResetHistos();
 
}
 
 
 
 
void KVBreakUp::DefineHistos()
{
   //Definition des histogrammes
   //A redefinir si besoin dans les classes filles
   //A Remplir dans la methode TreatePartition()
 
   hzt = new TH1F("KVBreakUp_hzt", "h ztotal", 200, -0.5, 199.5);
   lhisto->Add(hzt);
   hmt = new TH1F("KVBreakUp_hmt", "h mtotal", 200, -0.5, 199.5);
   lhisto->Add(hmt);
}
 
 
 
 
void KVBreakUp::SetZtot(Int_t zt)
{
   //Protected method
   if (zt != Ztotal) {
      Ztotal = zt;
      //zt correspond au nombre de possibles liens a casser
      if (bound) delete bound;
      bound = new Int_t[Ztotal];
      for (Int_t nn = 0; nn < Ztotal; nn += 1) bound[nn] = 1;
      nbre_nuc = Ztotal - Mtotal * Zmin;
   }
}
 
 
 
 
void KVBreakUp::SetMtot(Int_t mt)
{
   //Protected method
   if (mt > size_max) {
      Error("SetMtot", "%d -> La multiplicite max (%d) est depassee", mt, size_max);
      exit(EXIT_FAILURE);
   }
   Mtotal = mt;
   nbre_nuc = Ztotal - Mtotal * Zmin;
}
 
 
 
 
void KVBreakUp::SetZmin(Int_t zlim)
{
   //Protected method
 
   Zmin = zlim;
   nbre_nuc = Ztotal - Mtotal * Zmin;
 
}
 
 
 
 
void KVBreakUp::SetBreakUpMethod(KVString bup_method)
{
   //Protected method
   BreakUpMethod = bup_method;
 
}
 
 
 
 
void KVBreakUp::LinkEvent(KVEvent* evt)
{
 
   if (current_event) delete current_event;
   current_event = evt;
 
}
 
 
 
 
void KVBreakUp::RedefineTRandom(KVString TRandom_Method)
{
   //Permet de definir une classe pour le tirage aleatoire pour la
   //cassure en entiers
   //La class doit derivee de TRandom
 
   if (alea) delete alea;
   TClass* cl = new TClass(TRandom_Method.Data());
   alea = (TRandom*)cl->New();
   delete cl;
}
 
 
 
 
void KVBreakUp::SetConditions(Int_t zt, Int_t mt, Int_t zmin)
{
   //Definition du nombre a casser (zt), en plusieurs nombres entiers (mt)
   //avec une valeur minimale (zmin)
   SetZtot(zt);
   SetMtot(mt);
   SetZmin(zmin);
}
 
 
 
 
void KVBreakUp::DefineBreakUpMethod(KVString bup_method)
{
   //Permet de definir une methode
   //de cassure
 
   if (bup_method == "") {
      Info("DefineBreakUpMethod", "Available methods are");
      cout << "BreakUsingChain" << endl;
      cout << "BreakUsingPile" << endl;
      cout << "BreakUsingIndividual" << endl;
      cout << "BreakUsingLine" << endl;
      cout << "Make your choice" << endl;
   } else {
      SetBreakUpMethod(bup_method);
   }
 
}
 
 
 
 
Int_t KVBreakUp::GetZtot(void) const
{
   //Retourne le nombre entier a casser
   return Ztotal;
 
}
 
 
 
 
Int_t KVBreakUp::GetMtot(void) const
{
   //Retourne le nombre d'entiers apres cassure (la multiplicite)
   return Mtotal;
 
}
 
 
 
 
Int_t KVBreakUp::GetZmin(void) const
{
   //Retourne la taille minimale des entiers apres cassure
   return Zmin;
 
}
 
 
 
 
KVString KVBreakUp::GetBreakUpMethod(void) const
{
   //Retourne methode de cassure
   return BreakUpMethod;
 
}
 
 
 
 
void KVBreakUp::StorePartitions(Bool_t choix)
{
   //si choix=kTRUE, on enregistre les partitions
   SetBit(kStorePartitions, choix);
 
}
 
 
 
 
Int_t KVBreakUp::BreakUsingChain(void)
{
   //Methode de cassure
 
   //Conditions de depart
   //Mtotal clusters de taille minimale Zmin
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      size[mm] = Zmin - 1;
   }
   //On initilise la taille a Zmin-1
   //Le tirage aleatoire se fait
   //sur la taille restante Ztotal_corr
   Int_t Ztotal_corr = Ztotal - (Zmin - 1) * Mtotal;
 
   nl = "";
   Int_t bb = 0;
   Int_t mtire = 0;
   while (mtire < Mtotal) {
      //Tirage uniform d'entier entre 0 et Ztotal_corr - 1
      bb = TMath::Nint(alea->Uniform(0, (Ztotal_corr) - 0.5));
      //test si ce lien a deja ete casse
      if (bound[bb] == 1) {
         nl.Add(bb);
         mtire += 1;
         bound[bb] = 0;
      }
   }
 
   //La cassure a reussie
   //on recupere les tailles additionelles
   IntArray val = nl.GetArray();
   Int_t mtot_corr = val.size();
 
   Int_t zc = 0;
   //boucle sur les liens casses
   Int_t taille = 0;
   for (Int_t ii = 1; ii < mtot_corr; ii += 1) {
      taille = val[ii] - val[ii - 1];
      size[ii] += taille;  //mise a jour des tailles des clusters
      hzz->Fill(size[ii]);
      bound[val[ii]] = 1; //Reset du lien pour le prochain tirage
      zc += size[ii];   //incrementation de la taille totale de controle
   }
   taille = Ztotal_corr - val[mtot_corr - 1] + val[0];
 
   size[0] += taille;
   hzz->Fill(size[0]);
   bound[val[0]] = 1;
   zc += size[0];
 
   //cout << zc << " " << Ztotal << endl;
   if (zc == Ztotal && mtot_corr == Mtotal) return 1;
   else return 0;
 
}
 
 
 
 
Int_t KVBreakUp::BreakUsingLine(void)
{
   //Methode de cassure
 
   //Conditions de depart
   //Mtotal clusters de taille minimale Zmin
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      size[mm] = Zmin - 1;
   }
   //On initilise la taille a Zmin-1
   //Le tirage aleatoire se fait
   //sur la taille restante Ztotal_corr
   Int_t Ztotal_corr = Ztotal - (Zmin - 1) * Mtotal;
 
   nl = "";
   Int_t bb = 0;
   Int_t mtire = 0;
   while (mtire < Mtotal - 1) {
      //Tirage uniform d'entier entre 0 et Ztotal_corr - 1
      bb = TMath::Nint(alea->Uniform(1, (Ztotal_corr - 1) - 0.5));
      //test si ce lien a deja ete casse
      if (bound[bb] == 1) {
         //cout << "bb="<<bb << endl;
         nl.Add(bb);
         mtire += 1;
         bound[bb] = 0;
      }
   }
 
   //cout << "mtire="<<mtire << endl;
   //
   nl.Add(0);
   nl.Add(Ztotal_corr);
   //cout << "nl="<<nl.AsString() << endl;;
   //La cassure a reussie
   //on recupere les tailles additionelles
   IntArray val = nl.GetArray();
   Int_t mtot_corr = val.size();
   //cout << "mtot_corr="<<mtot_corr << endl;
   Int_t zc = 0;
   //boucle sur les liens casses
   Int_t taille = 0;
   for (Int_t ii = 1; ii < mtot_corr; ii += 1) {
      taille = val[ii] - val[ii - 1];
      //cout << "ii="<< ii <<" taille="<<taille<<endl;
      size[ii - 1] += taille; //mise a jour des tailles des clusters
      hzz->Fill(size[ii - 1]);
      bound[val[ii]] = 1; //Reset du lien pour le prochain tirage
      zc += size[ii - 1]; //incrementation de la taille totale de controle
   }
 
   //cout << "zc="<<zc << " " << Ztotal << endl;
   if (zc == Ztotal && mtot_corr - 1 == Mtotal) return 1;
   else return 0;
 
}
 
 
 
 
Int_t KVBreakUp::BreakUsingIndividual(void)
{
   //Methode de cassure
 
   //Conditions de depart
   // Mtotal clusters de taille minimale Zmin
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      size[mm] = Zmin;
   }
   //On initilise la taille a Zmin
   //Le tirage aleatoire se fait
   //sur la taille restante Ztotal_corr
   Int_t Ztotal_corr = Ztotal - (Zmin) * Mtotal;
 
   Double_t val[Mtotal], sum_val = 0;
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      //Tirage des Mtotal valeurs entre [0; 1[
      val[mm] = alea->Uniform(0, 1);
      sum_val += val[mm];
   }
   Int_t zc = 0;
   Int_t surplus[Mtotal];
   //normalisation a une charge totale de Ztotal_corr
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      surplus[mm] = TMath::Nint(val[mm] * Ztotal_corr / sum_val);
      zc += surplus[mm];
   }
   //Test effets spurieux pasage reel -> entier
   //Redistribution eventuelle de charge manquante ou en exces
   if (zc != Ztotal_corr) {
      Int_t diff = Ztotal_corr - zc;
      nraffine += 1;
      Int_t bb = 0;
      if (diff > 0) {
         for (Int_t mm = 0; mm < TMath::Abs(diff); mm += 1) {
            bb = TMath::Nint(alea->Uniform(0, (Mtotal) - 0.5));
            surplus[bb] += 1;
         }
      } else {
         for (Int_t mm = 0; mm < TMath::Abs(diff); mm += 1) {
            bb = TMath::Nint(alea->Uniform(0, (Mtotal) - 0.5));
            if (surplus[bb] > 0)
               surplus[bb] -= 1;
            else
               mm -= 1;
         }
      }
   }
   zc = 0;
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      size[mm] += surplus[mm];
      zc += size[mm];
      hzz->Fill(size[mm]);
   }
 
   if (zc == Ztotal) return 1;
   else return 0;
 
}
 
 
 
 
Int_t KVBreakUp::BreakUsingPile(void)
{
   //Methode de cassure
   Int_t bb = 0;
 
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      size[mm] = Zmin;
   }
   //on distribue aleatoirement les charges
   //disponibles une par une sur les differents fragments
   for (Int_t mm = 0; mm < nbre_nuc; mm += 1) {
      bb = TMath::Nint(alea->Uniform(0, (Mtotal) - 0.5));
      size[bb] += 1;
   }
 
   Int_t mc = 0, zc = 0;
   for (Int_t mm = 0; mm < Mtotal; mm += 1) {
      Int_t taille = size[mm];
      hzz->Fill(taille);
      mc += 1;
      zc += taille;
   }
   if (mc == Mtotal && zc == Ztotal) return 1;
   else return 0;
 
 
}
 
 
 
 
void KVBreakUp::TreatePartition()
{
   //Remplissage des histogrammes predefinis
   //A redefinir si besoin dans les classes filles
   //
   //Enregistrement de la partition (si demande, via StorePartitions(kTRUE) par defaut)
   //Si une partition identique est deja presente, on incremente sa population
   //sinon on enregistre celle-ci, voir KVPartitionList
   hmt->Fill(Mtotal);
   hzt->Fill(Ztotal);
 
 
   partition = new KVIntegerList();
   partition->Fill(size, Mtotal);
   Int_t* tab = partition->CreateTableOfValues();
 
   current_event->Clear();
   for (Int_t ii = 0; ii < Mtotal; ii += 1)
      current_event->AddNucleus()->SetZ(tab[ii]);
 
   delete [] tab;
 
   if (TestBit(kStorePartitions)) {
      if (Fill(partition))
         delete partition;
   } else {
      delete partition;
   }
}
 
 
 
 
 
KVEvent* KVBreakUp::GetCurrentEvent()
{
 
   return current_event;
 
}
 
 
 
 
void KVBreakUp::BreakNtimes(Int_t times)
{
   //On realise times fois la cassure
   //suivant les conditions definis vis la methode SetConditions et DefineBreakUpMethod
   //Plusieurs series de cassures peuvent etre ainsi realise
   //
   //Si l utilisateur veut changer de conditions et ne pas melanger les partitions
   //il faut appeler la methode Clear()
 
   Start();
 
   TMethodCall meth;
   meth.InitWithPrototype(this->IsA(), BreakUpMethod.Data(), "");
   Long_t ret;
   if (meth.IsValid() && meth.ReturnType() == TMethodCall::kLong) {
      for (Int_t nn = 0; nn < times; nn += 1) {
         if (nn % 10000 == 0) Info("BreakNtimes", "%d partitions generees sur %d", nn, times);
         meth.Execute(this, "", ret);
         if (ret == 1)
            TreatePartition();
         else {
            Info("BreakNtimes", "%s retourne %ld", BreakUpMethod.Data(), ret);
            nn -= 1;
         }
      }
   }
   Info("BreakNtimes", "Tirage termine");
   niter_tot += times;
 
   Stop();
   Info("BreakNtimes", "Temps ecoule en secondes : %d", GetDeltaTime());
 
}
 
 
 
 
void KVBreakUp::BreakNtimesOnGaussian(Int_t times, Double_t Ztot_moy, Double_t Ztot_rms, Double_t Mtot_moy, Double_t Mtot_rms, Int_t zmin)
{
   //On realise times fois la cassure sur une double gaussienne
   //la methode SetCondition n'a ici aucune influence
   //les valeurs ztot et mtot sont tirees aleatoirement sur une gaussienne
   //a chaque iteration
   //
   //Plusieurs series de cassures peuvent etre ainsi realise
   //
   //Si l utilisateur veut changer de conditions et ne pas melanger les partitions
   //il faut appeler la methode Clear()
 
   Start();
 
   TMethodCall meth;
   meth.InitWithPrototype(this->IsA(), BreakUpMethod.Data(), "");
   Long_t ret;
   TRandom3* gaus = new TRandom3();
   if (meth.IsValid() && meth.ReturnType() == TMethodCall::kLong) {
      for (Int_t nn = 0; nn < times; nn += 1) {
         if (nn % 1000 == 0) printf("%d partitions generees sur %d\n", nn, times);
 
         Int_t zt = TMath::Nint(gaus->Gaus(Ztot_moy, Ztot_rms));
         Int_t mt = TMath::Nint(gaus->Gaus(Mtot_moy, Mtot_rms));
 
         if (mt > 0 && zt > 0 && zt >= mt * zmin) {
            SetConditions(zt, mt, zmin);
            meth.Execute(this, "", ret);
            if (ret == 1)
               TreatePartition();
            else {
               //cout << BreakUpMethod << " retourne " << ret << endl;
            }
         } else {
            nn -= 1;
         }
      }
   }
   delete gaus;
   niter_tot += times;
 
   Stop();
   Info("BreakNtimesOnGaussian", "Temps ecoule en secondes : %d", GetDeltaTime());
 
}
 
 
 
 
KVEvent* KVBreakUp::BreakOnGaussian(Double_t Ztot_moy, Double_t Ztot_rms, Double_t Mtot_moy, Double_t Mtot_rms, Int_t zmin)
{
   //On realise la cassure sur une double gaussienne
   //
   //les valeurs ztot et mtot sont tirees aleatoirement sur une gaussienne
   //
   //Plusieurs series de cassures peuvent etre ainsi realise
   //
 
   TMethodCall meth;
   meth.InitWithPrototype(this->IsA(), BreakUpMethod.Data(), "");
   Long_t ret;
 
   if (meth.IsValid() && meth.ReturnType() == TMethodCall::kLong) {
      Int_t zt = -1;
      Int_t mt = -1;
      while (!(mt > 0 && zt > 0 && zt >= mt * zmin)) {
         zt = TMath::Nint(alea->Gaus(Ztot_moy, Ztot_rms));
         mt = TMath::Nint(alea->Gaus(Mtot_moy, Mtot_rms));
      }
 
      SetConditions(zt, mt, zmin);
      meth.Execute(this, "", ret);
      if (ret == 1) {
         TreatePartition();
         return current_event;
      } else {
         return 0;
      }
   }
   return 0;
}
 
 
 
 
void KVBreakUp::BreakFromHisto(TH2F* hh_zt_VS_mt, Int_t zmin)
{
   //On realise times fois la cassure suivant un histogramme a deux dimensions
   //avec en definition :
   // mtot -> axe X
   // ztot -> axe Y
   //la methode SetCondition n'a ici aucune influence
   //les valeurs ztot et mtot sont tirees aleatoirement sur cette histogramme
   //a chaque iteration
   //
   //Plusieurs series de cassures peuvent etre ainsi realise
   //
   //Si l utilisateur veut changer de conditions et ne pas melanger les partitions
   //il faut appeler la methode Clear()
 
   TH2F* h2 = hh_zt_VS_mt;
   if (!h2) return;
 
   Start();
 
   TMethodCall meth;
   meth.InitWithPrototype(this->IsA(), BreakUpMethod.Data(), "");
   Long_t ret;
 
   Int_t zt, mt;
   Int_t stat_tot = Int_t(h2->Integral());
   Int_t stat_par = 0;
   for (Int_t nx = 1; nx <= h2->GetNbinsX(); nx += 1)
      for (Int_t ny = 1; ny <= h2->GetNbinsY(); ny += 1) {
         Int_t stat = TMath::Nint(h2->GetBinContent(nx, ny));
         if (stat > 0) {
            mt = TMath::Nint(h2->GetXaxis()->GetBinCenter(nx));
            zt = TMath::Nint(h2->GetYaxis()->GetBinCenter(ny));
            if (mt > 0 && zt > 0 && zt >= mt * zmin) {
 
               SetConditions(zt, mt, zmin);
               for (Int_t nn = 0; nn < stat; nn += 1) {
                  meth.Execute(this, "", ret);
                  if (ret == 1)
                     TreatePartition();
                  stat_par += 1;
                  if (stat_par % 1000 == 0) printf("%d partitions generees sur %d\n", stat_par, stat_tot);
               }
            } else {
               cout << zt << " " << mt << endl;
            }
         }
      }
 
   niter_tot += stat_par;
 
   Stop();
   Info("BreakNtimesFromHisto", "Temps ecoule en secondes : %d", GetDeltaTime());
 
}
 
 
 
 
Int_t KVBreakUp::GetTotalIterations(void)
{
   //Retourne le nombre total d iterations
   //depuis le dernier clear
   return niter_tot;
 
}
 
 
 
 
KVHashList* KVBreakUp::GetHistos()
{
   //Retourne la liste des histogrammes
   //definis dans DefineHistos()
   //si l utilisateur a utilise lhisto->Add(TObject * )
   return lhisto;
 
}
 
 
 
 
KVHashList* KVBreakUp::GetObjects()
{
   //Retourne la liste d'objects autres ...
   return lobjects;
 
}
 
 
 
 
void KVBreakUp::ResetTotalIterations()
{
   //Comme c'est ىcrit
   niter_tot = 0;
 
}
 
 
 
 
void KVBreakUp::ResetHistos()
{
   //Met a zero le contenu des histogrammes
   lhisto->Execute("Reset", "");
 
}
 
 
 
 
void KVBreakUp::DrawPanel()
{
   //Trace les histos definis
   //A redefinir si besoin dans les classes filles
 
   TCanvas* c1 = 0;
   if (!(c1 = (TCanvas*)gROOT->GetListOfCanvases()->FindObject("BreakUp_Control"))) {
 
      c1 = new TCanvas("BreakUp_Control", "Control", 0, 0, 900, 900);
      c1->Divide(2, 2, 0.001, 0.001, 10);
      c1->cd(1);
      if (hzz->GetEntries() > 0)   gPad->SetLogy(1);
      hzz->Draw();
      c1->cd(2);
      if (hzt->GetEntries() > 0)   gPad->SetLogy(1);
      hzt->Draw();
      c1->cd(3);
      if (hmt->GetEntries() > 0)   gPad->SetLogy(1);
      hmt->Draw();
 
   }
}
 
 
 
 
void KVBreakUp::SaveHistos(KVString filename, KVString suff, Option_t* option)
{
   //Permet la sauvegarde des histogrammes
   //dans un fichier (option est l'option classique de TFile::TFile())
   //
   //Si filename=="" (defaut) -> nom du fichier = "KVBreakUp_Ouput.root";
   //Si suff=="" (defaut) -> definition d'un suffixe en utisant les arguments de SetConditions()
   //Ce suffixe est ensuite utilise dans la redefinition des noms des histogrammes
   //Exemple:
   // nom de l'histogramme : histo -> histo_suff
   //
   if (filename == "") filename = "KVBreakUp_Ouput.root";
   TFile* file = new TFile(filename.Data(), option);
 
   if (suff == "")
      suff.Form("Zt%d_Mt%d_Zm%d_%s", GetZtot(), GetMtot(), GetZmin(), GetBreakUpMethod().Data());
 
   KVString snom;
   for (Int_t nn = 0; nn < lhisto->GetEntries(); nn += 1) {
      snom.Form("%s_%s", lhisto->At(nn)->GetName(), suff.Data());
      lhisto->At(nn)->Write(snom.Data());
   }
   file->Close();
 
}
 
 
 
 
void KVBreakUp::PrintConfig() const
{
   // Comme c'est ىcrit
   // Why not "Print(Option_t*)" ?
   //   - because TCollection has several 'virtual' Print methods which
   //     are overloaded (i.e. have different arguments): BAD!
 
   Info("Print", "Configuration for the break up");
   printf(" Ztot=%d - Mtot=%d - Zmin=%d\n", GetZtot(), GetMtot(), GetZmin());
   printf(" Charge to be distributed %d - Biggest possible charge %d\n", nbre_nuc, Zmin + nbre_nuc);
   printf(" Methode de cassage aleatoire %s\n", GetBreakUpMethod().Data());
   alea->Print();
   printf(" Partition are stored via KVPartitionList : %d\n", Int_t(TestBit(kStorePartitions)));
   printf("------------------------------------------------------");
}
 
 
 
 
void KVBreakUp::Start()
{
   //protected method
   //Signal start
   TDatime time;
   tstart = time.GetHour() * 3600 + time.GetMinute() * 60 + time.GetSecond();
 
}
 
 
 
 
void KVBreakUp::Stop()
{
   //protected method
   //Signal stop
   TDatime time;
   tstop = time.GetHour() * 3600 + time.GetMinute() * 60 + time.GetSecond();
   tellapsed = tstop - tstart;
 
}
 
 
 
 
Int_t KVBreakUp::GetDeltaTime()
{
   //Retoune le temps ecoules (en seconde)
   //entre un appel Start() et un appel Stop()
   return tellapsed;
 
}