doFastGlauberMcMaker.C

//****************************************************************************************************
//  Run Fast Monte Carlo Glauber simulation (StFastGlauberMcMaker)
//****************************************************************************************************
/****************************************************************************************************
 * $Id: doFastGlauberMcMaker.C,v 1.3 2012/04/25 05:29:34 hmasui Exp $
 * $Log: doFastGlauberMcMaker.C,v $
 * Revision 1.3  2012/04/25 05:29:34  hmasui
 * Added deformation flag
 *
 * Revision 1.2  2010/10/27 00:31:45  hmasui
 * Added loading relevant STAR libs
 *
 * Revision 1.1  2010/06/24 23:41:40  hmasui
 * Macro to run StFastGlauberMcMaker
 *
 ****************************************************************************************************/

//____________________________________________________________________________________________________
// Specify the system by TString, detailed parameters will be defined inside StFastGlauberMcMaker 
void doFastGlauberMcMaker
(
 const TString outputFileName = "fastmcglauber.root",
 const UInt_t nevents = 1000,
 const TString system = "AuAu",
 const Double_t energy = 200,
 const TString type = "default",
 const Bool_t isDeformed = kFALSE,
 const Double_t replusionDistance = 0.0
 )
{
        gBenchmark->Start("doFastGlauberMcMaker");

        gSystem->Load("St_base");
        gSystem->Load("StUtilities");
        gSystem->Load("StGlauberUtilities");
        gSystem->Load("StCentralityMaker");
        gSystem->Load("StGlauberTree");
        gSystem->Load("StFastGlauberMcMaker");

        StFastGlauberMcMaker* maker = new StFastGlauberMcMaker(outputFileName, system, energy, type, isDeformed);
        //  maker->DebugOn();

        // Print available types
        maker->Print("type");

        // Smearing option will be added in type later
        //  maker->DoHardCoreSmearing() ; // Hard-core smearing ON
        //  maker->DoGaussianSmearing() ; // Gaussian smearing ON

        // Repulsion option will be added in type later
        //  maker->SetRepulsionDistance(replusionDistance); // Repulsion distance between two nucleons

        // Collision profile were already added (default is hard-core)
        // specify "gauss" in type for Gaussian collision profile
        //  if( collisionProfile == 0 ) maker->DoHardCoreCollision() ;
        //  if( collisionProfile == 1 ) maker->DoGaussianCollision() ;

        // Run nevents
        maker->Run(nevents);

        // Close ROOT file etc
        maker->Finish();

        gBenchmark->Stop("doFastGlauberMcMaker");
        gBenchmark->Show("doFastGlauberMcMaker");
        gBenchmark->Reset();
}
//____________________________________________________________________________________________________
//  Specify all parameters by hand
void _doFastGlauberMcMaker
(
 const Char_t* outputFileName,
 const UInt_t nevents = 1000,
 const UInt_t A = 197,
 const Double_t R = 6.38,
 const Double_t d = 0.535,
 const Double_t beta2 = -0.131,
 const Double_t beta4 = -0.031,
 const Double_t sigma = 4.2,
 const Double_t energy = 200,
 const Double_t replusionDistance = 0.9,
 const UInt_t collisionProfile = 0 // 0=hard-core (default), 1=gaussian profile
 )
{
        gBenchmark->Start("doFastGlauberMcMaker");
        gSystem->Load("StGlauberUtilities");
        gSystem->Load("StCentralityMaker");
        gSystem->Load("StFastGlauberMcMaker");

        StFastGlauberMcMaker* maker = new StFastGlauberMcMaker(outputFileName, A, R, d, beta2, beta4, sigma, energy);

        // Hard-core smearing ON
        //  maker->DoHardCoreSmearing() ;

        // Gaussian smearing ON
        //  maker->DoGaussianSmearing() ;

        // Repulsion distance between two nucleons
        maker->SetRepulsionDistance(replusionDistance);

        // Collision profile
        if( collisionProfile == 0 ) maker->DoHardCoreCollision() ;
        if( collisionProfile == 1 ) maker->DoGaussianCollision() ;

        // Run nevents
        maker->Run(nevents);

        // Close ROOT file etc
        maker->Finish();

        gBenchmark->Stop("doFastGlauberMcMaker");
        gBenchmark->Show("doFastGlauberMcMaker");
        gBenchmark->Reset();
}